1.算法运行效果图预览
2.算法运行软件版本
matlab2022a
3.算法理论概述
基于OFDM通信系统的PAPR抑制算法是降低OFDM信号峰均比（Peak-to-Average Power Ratio，PAPR）的技术，以提高通信系统的性能和稳定性。其中，IPTS（选择性映射迭代削峰）、OPTS（优化的PTS）和CEPTS（压缩扩展变换选择性映射）是三种常见的PAPR抑制算法。下面将详细介绍这三种算法的原理和数学公式。
3.1、IPTS算法
IPTS算法是一种基于选择性映射迭代削峰的方法，通过迭代削峰和选择性映射，降低OFDM信号的PAPR。具体步骤如下：
对OFDM信号进行IFFT变换，得到时域信号。
对时域信号进行削峰处理，将峰值超过一定阈值的信号进行削减。
对削减后的信号进行IFFT变换，得到新的频域信号。
通过选择性映射，选择PAPR最低的频域信号作为输出信号。
3.2、OPTS算法
OPTS算法是一种基于优化的PTS方法，通过将OFDM信号分成多个子块，对每个子块进行相位旋转和幅度调整，以降低PAPR。具体步骤如下：
将OFDM信号分成多个子块。
对每个子块进行相位旋转和幅度调整，使得子块的PAPR最低。
将调整后的子块重新组合成完整的OFDM信号。
3.3、CEPTS算法
CEPTS算法是一种基于压缩扩展变换选择性映射的方法，通过对OFDM信号进行压缩扩展变换和选择性映射，降低PAPR。具体步骤如下：
对OFDM信号进行压缩扩展变换，得到多个变换后的信号。
对每个变换后的信号进行PAPR计算，选择PAPR最低的信号作为输出信号。
将选择的信号进行逆变换，得到最终的OFDM信号。
综上所述，IPTS、OPTS和CEPTS算法都是通过不同的方式对OFDM信号进行处理，以降低PAPR，提高通信系统的性能和稳定性。具体选择哪种算法需要根据实际应用场景和性能需求来决定。
4.部分核心程序
```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) = 10log10(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)) + ijimag(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
```
