基于OFDM通信系统的PAPR抑制算法matlab仿真,对比IPTS,OPTS,CEPTS三种算法

news2024/11/21 1:47:41

1.算法运行效果图预览

2.算法运行软件版本

3.部分核心程序

4.算法理论概述

4.1、IPTS算法

4.2、OPTS算法

4.3、CEPTS算法

5.算法完整程序工程

1.算法运行效果图预览

2.算法运行软件版本

matlab2022a

3.部分核心程序

for k=1:Nframes
    if mod(k,10) == 0
       k/10
    end
    %产生数据源
    QPSK_Ind     = floor(length(Map_qpsk)*rand(1,Nfft)) + 1;
    %调制，这里为了研究PAPR性能，所以不加入编码模块和交织模块
    Qpsk_mod     = Map_qpsk(QPSK_Ind(1,:));   
    %进行IFFT变换
    Dat_Ifft     = ifft(Qpsk_mod,[],2); 
    %计算功率和PAPR
    Signal_Power = abs(Dat_Ifft.^2);
    Peak_Power   = max(Signal_Power,[],2);
    Mean_Power   = mean(Signal_Power,2);
    PAPRo(k)     = 10*log10(Peak_Power./Mean_Power);
    %随机分块
    QPSK_Ind     = randperm(Nfft);
    A            = zeros(Npts,Nfft);
    for v=1:Npts
        A(v,QPSK_Ind(v:Npts:Nfft)) = Qpsk_mod(QPSK_Ind(v:Npts:Nfft));
    end
    a            = ifft(A,[],2); 
    %限幅
    Tho          = mean2(abs(a));
    [rr,cc] = size(a);
    for i = 1:rr
        for j = 1:cc
            if abs(a(i,j)) > Tho
               a(i,j) = Tho*(real(a(i,j)) + ij*imag(a(i,j)))/abs(a(i,j));
            end
        end
    end
    %PCME算法
    P0   = 0.5*ones(1,Npts);%初始概率为0.5
    Ps   = zeros(Iter,Npts);
    P    = zeros(Iter,Npts);
    for iter = 1:Iter
        %根据随机分布，产生一组序列c
        ....................................................................
        
        for j = 1:J
            Phase_Factor = repmat(1-2*c(j,:)',1,Nfft);  
            aa           = sum(a.*Phase_Factor);
            Signal_Power = abs(aa.^2);
            Peak_Power   = max(Signal_Power,[],2);
            Mean_Power   = mean(Signal_Power,2);
            F(j)         = 10*log10(Peak_Power./Mean_Power);
        end
        
        %对当前迭代产生的J个F进行增序排序
        [F2,IND]      = sort(F);
        %计算rj
        r(iter)   = sum(F2(1:ceil(Po*J)))/ceil(Po*J);
        IND2      = find(F <= r(iter));
        
        if isempty(IND2) == 1
           IND2 = 1;
        else
           IND2 = IND2;  
        end
        
        for pp = 1:Npts
            
            for s1 = 1:J
                I(pp,s1) = c(IND(s1),pp);
                tmp11s(s1) = I(pp,s1)*exp(-1*F2(s1)); 
                tmp12s(s1) = exp(-1*F2(s1)); 
            end
            P(iter,pp) =  sum(tmp11s)/sum(tmp12s);
        end        
        %更新概率P
        ......................................................
    end
    %根据PMCE计算得到的相位因子来计算PAPR值
    aa                = sum(a.*repmat(sign((1-2*Ps(iter,:)))',1,Nfft));
    Signal_Power      = abs(aa.^2);
    Peak_Power_temp   = max(Signal_Power,[],2);
    Mean_Power_temp   = mean(Signal_Power,2);
    PAPR_temp(k)      = 10*log10(Peak_Power_temp./Mean_Power_temp);
end
01_062m

4.算法理论概述

基于OFDM通信系统的PAPR抑制算法是降低OFDM信号峰均比（Peak-to-Average Power Ratio，PAPR）的技术，以提高通信系统的性能和稳定性。其中，IPTS（选择性映射迭代削峰）、OPTS（优化的PTS）和CEPTS（压缩扩展变换选择性映射）是三种常见的PAPR抑制算法。下面将详细介绍这三种算法的原理和数学公式。