www.gusucode.com > 对fdma和scfdma仿真性能的比较，可以清晰具体的比较峰均比的性能 > code/ofdm_papr.m

    %=====================OFDM  PAPR =========================%
%============initialization ========================%
para=1024;       %Number of parallel channel to transmit
N = 2048;        %- Subcarrier number
nd=50;          % Number of information OFDM symbol for one loop
M=4;
%==============transmitter=====================%
    seridata=randsrc(1,para*nd,0:M-1);
    paradata=reshape(seridata,para,nd);      %size(51*nd*ml)    
%===========ml modulation ================%
x=pskmod(paradata, M);
x1=IFDMA(x,N/para,para,nd);
%------------IFFT-------------------------%
Y=ifft(x1);
id=Y;
%==============CCDF calculation====================%
peak = abs(id).^2 ;                     %Peak Power
avg = norm(id)^2./length(id);           %Average Power  largest singular value
papr = peak / avg;                      %PAPR
papr_dB = 10*log10(papr);               %PAPR in dB
PAPR_new1=papr_dB;
[y, PAPRo,PAPR_thr, new_ccdf,MAXIMUM]=ccdf_lee(PAPR_new1,N);
figure(3),  semilogy(PAPR_thr, new_ccdf, 'b*-');hold on;grid on;
axis([0 15   10^-3   10^0]);
xlabel('PAPRo(dB)');
ylabel('Pr(PAPR>PAPRo)');