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

    %=====================SCIFDMA  PAPR =========================%
%============initialization ========================%
datablock=256;
N = 512;        %- Subcarrier number
nd=5;          % Number of information OFDM symbol for one loop
M=16;
Q=N/datablock;              %bandwide spreading factor
%==============transmitter=====================%
    seridata=randsrc(1,datablock*nd,0:M-1);
    paradata=reshape(seridata,datablock,nd);      %size(51*nd*ml)    
%===========ml modulation ================%
x= qammod(paradata,M);
Y=fft(x);
%------------IFFT-------------------------%
Y1=IFDMA(Y,Q,datablock,nd);
id=ifft(Y1);
%==============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)');