ccdf.m

改变子载波的数目n

% Program 2-1-2
% Pr1.m
% Function to calculate Pr(papr_in_dB>papr0_in_dB)
% programmed by phui
function Pprr=ccdf(pr)

%****************** variables *************************
% pr     : Input data (para-by-nd matrix)
% para   : Number of paralell channels
% fftlen : FFT length
% nloop  : Number of simulation loops
% *****************************************************
nloop=1000;
PPP=12000;
for ppp=1:1:PPP;
   papr0_in_dB=2+0.001*ppp;                  %set a papr0_in_dB
Pr(ppp)=(sum(pr>papr0_in_dB))/nloop;         %calculate Pr(papr_in_dB>papr0_in_dB)
end 
Pprr=Pr;