%***************************************************
% SLP-OFDM系統PARP理論分布曲線 
% N 子載波數  
% M 隨機相位序列個數 
%---------------------------------------------------
clc; clear
N = 128; 
x = 0:0.1:16;  % x為PAPR值，單位dB
M = 1;  p1 = (1-(1- exp(-10.^(x./10))).^N).^M;
M = 2;  p2 = (1-(1- exp(-10.^(x./10))).^N).^M;
M = 4;  p3 = (1-(1- exp(-10.^(x./10))).^N).^M;
M = 8;  p4 = (1-(1- exp(-10.^(x./10))).^N).^M;
M = 16;  p5 = (1-(1- exp(-10.^(x./10))).^N).^M;
M = 32;  p6 = (1-(1- exp(-10.^(x./10))).^N).^M;
figure(1);
semilogy(x,p1);hold on;
semilogy(x,p2);hold on;
semilogy(x,p3);hold on;
semilogy(x,p4);hold on;
semilogy(x,p5);hold on;
semilogy(x,p6);hold off;
ylim([10E-6 1]);xlabel('峰平功率比PAPR0(dB)');ylabel('CCDF');