OFDM峰均比限幅clipping抑制matlab仿真程序

限幅(clipping)，是降低OFDM系统PAPR最直接的方法。根据峰均比的统计特性可知，高峰平比出现的概率极小，削去过高的瞬时高幅值,降低整个系统的误比特率性能，改善CCDF曲线

clear all;
clc;


fl=512;  % 设置FFT长度
Num=100;   
ca=256; %设置并行传输的子载波个数
NN=1:.1:14; % CCDF的门限值
ccdf0=zeros（1131）;
ccdf1=ccdf0;
ccdf2=ccdf0;
ccdf3=ccdf0;
ccdf4=ccdf0;

 for i=1:Num   
Signal（（i*（i-1）/2）*ca*4+1:（i*（i-1）/2+i）*ca*4）=rand（1ca*i*4）>0.5;%产生0，1 随机序列，符号数为ca*（1+Num）*Num*2
 end

 Ns=（1+Num）*Num;  %设置一个帧结构中OFDM信号的个数
 
for i=1:ca
    for j=1:Ns*2
        SigPara（ij）=Signal（i*j）;%串并变换
    end
end
%QPSK调制，将数据分为两个通道
for j=1:Ns
    ich（:j）=SigPara（:2*j-1）;
    qch（:j）=SigPara（:2*j）;
end
kmod=1./sqrt（2）;
ich1=ich.*kmod;
qch1=qch.*kmod;
x=ich1+qch1.*sqrt（-1）; 

x1=zeros（4*caNs）;   %   过采样 
   for i=1:ca/2
   x1（i:）=x（i:）;
   x1（i+7/2*ca:）=x（i+ca/2:）;
   end
   
 y=ifft（x[]2）; %频域数据变时域
y1=（abs（y））.^2;  
power_average=mean（y1）;
max_power=max（y1）;
papr0=10*log10（max_power./power_average）;  % 原来的信号

% 限幅
le=0.5*10^-2;  %限幅电平模值  
y2=y;
ccdf1=ccdf0;
for i=1:ca
    for j=1:Ns
        if abs（y（ij））>le;
            y2（ij）=le.*exp（angle（y（ij））*sqrt（-1））; 
        end      
    end
end 
y3=（abs（y2））.^2;
power_average_le=mean（y3）;
max_power_le=max（y3）;
papr1=10*log10（max_power_le./power_average_le）;




% 限幅
le=0.8*10^-2;  %限幅电平模值 
y3=y;
ccdf1=ccdf0;
for i=1:ca
    for j=1:Ns
        if abs（y（ij））>le;
            y3（ij）=le.*exp（angle（y（ij））*sqrt（-1））; 
        end