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

    function [y, PAPRo,PAPR_thr, new_ccdf,MAXIMUM]=ccdf_lee(PAPR_new1,N)
%---------------------------------------------------------------%
PAPR_new=vec2mat(PAPR_new1, N);
MAXIMUM = 10*log10(N);
PAPR_thr = [0:1:MAXIMUM];

% method 1
for n = 1 : length(PAPR_thr)
   PAPR = PAPR_thr(n);
   count1 = 0;
   for m = 1 : size(PAPR_new,1)
       PAPR_neww(m)=max(PAPR_new(m,:));
      if PAPR < PAPR_neww(m)
         count1 = count1 + 1;
      end      
   end
   new_ccdf(n) = count1/length(PAPR_neww);  
end

%CCDF theory=======================================================
PAPRo = [0:0.1:MAXIMUM];
alpha = 2.8;
K = 1;
PAPRx = 10.^(PAPRo/10);
y = (1 - (1 - exp( -PAPRx )).^(alpha*(N))).^K;  % CCDF