function PAPR_16QAM_Clipping2 Num_Loop=10000;%总共模拟10000个点,数字越大越精确 Num_Subcarrier=64;%子载波个数 Mapper=[3+3i 1+3i -1+3i -3+3i 3+1i 1+1i -1+1i -3+1i 3-1i 1-1i -1-1i -3-1i 3-3i 1-3i -1-3i -3-3i]/sqrt(4);%QPSK J=5;0%Oversampling factor过采样因子 CR=1;%Clipping Ratio裁剪比 PAPR_Original=zeros( 1,Num_Loop); PAPR_Clipping=zeros(1,Num_Loop); for n=1:Num_Loop InputSymbolIndex=randi([1 16],1,Num_Subcarrier); OFDM_Freq=Mapper(InputSymbolIndex); OFDM_Time=sqrt(J*Num_Subcarrier)*ifft([OFDM_Freq(1:Num_Subcarrier/2) zeros(1,(J-1)*Num_Subcarrier) OFDM_Freq(Num_Subcarrier/2+1:Num_Subcarrier)]); PAPR_Original(n)=10*log10(max(abs(OFDM_Time))^2./mean(abs(OFDM_Time).^2)); for p=1:J*Num_Subcarrier if abs(OFDM_Time(p))>CR OFDM_Time(p)=CR*cos(angle(OFDM_Time(p)))+1i*CR*sin(angle(OFDM_Time(p))); end end PAPR_Clipping(n)=10*log10(max(abs(OFDM_Time)).^2./mean(abs(OFDM_Time).^2)); end Pr_Ori=[]; Pr_Clip=[]; PAPRO_Start=2; PAPRO_End=12; step=0.2; count=1; for m=PAPRO_Start:step:PAPRO_End temp_Ori=sum(PAPR_Original(:)>m)/Num_Loop; temp_Clip=sum(PAPR_Clipping(:)>m)/Num_Loop; Pr_Ori(count)=temp_Ori; Pr_Clip(count)=temp_Clip; count=count+1; end semilogy(PAPRO_Start:step:PAPRO_End,Pr_Ori, '-b*',PAPRO_Start:step:PAPRO_End,Pr_Clip, '-r*'); axis([PAPRO_Start PAPRO_End 10^-4 1]); grid on; hold on:
时间: 2024-01-30 10:02:36 浏览: 21
这段代码是一个用MATLAB编写的仿真程序,用于比较原始OFDM信号的峰均功率比(PAPR)和经过裁剪后的OFDM信号的PAPR。程序分为两个主要部分:第一部分是生成随机的OFDM符号,计算原始OFDM信号的PAPR;第二部分是对每个OFDM符号进行裁剪,并计算裁剪后的OFDM信号的PAPR。最后,程序将原始OFDM信号和裁剪后的OFDM信号的PAPR与一系列阈值进行比较,并绘制两个PAPR曲线的图像,以便比较它们在不同的PAPR阈值下的性能。
相关问题
function PAPR_QPSK_Clipping Num_Loop=10000;%总共模拟10000个点,数字越大越精确 Num_Subcarrier=64;%子载波个数 Mapper=[1+1i -1+1i 1-1i -1-1i]/sqrt(2);%QPSK J=5;0%Oversampling factor过采样因子 CR=1;%Clipping Ratio裁剪比 PAPR_Original=zeros( 1,Num_Loop); PAPR_Clipping=zeros(1,Num_Loop); for n=1:Num_Loop InputSymbolIndex=randi([1 4],1,Num_Subcarrier); OFDM_Freq=Mapper(InputSymbolIndex); OFDM_Time=sqrt(J*Num_Subcarrier)*ifft([OFDM_Freq(1:Num_Subcarrier/2) zeros(1,(J-1)*Num_Subcarrier) OFDM_Freq(Num_Subcarrier/2+1:Num_Subcarrier)]); PAPR_Original(n)=10*log10(max(abs(OFDM_Time))^2./mean(abs(OFDM_Time).^2)); for p=1:J*Num_Subcarrier if abs(OFDM_Time(p))>CR OFDM_Time(p)=CR*cos(angle(OFDM_Time(p)))+1i*CR*sin(angle(OFDM_Time(p))); end end PAPR_Clipping(n)=10*log10(max(abs(OFDM_Time)).^2./mean(abs(OFDM_Time).^2)); end Pr_Ori=[]; Pr_Clip=[]; PAPRO_Start=2; PAPRO_End=12; step=0.2; count=1; for m=PAPRO_Start:step:PAPRO_End temp_Ori=sum(PAPR_Original(:)>m)/Num_Loop; temp_Clip=sum(PAPR_Clipping(:)>m)/Num_Loop; Pr_Ori(count)=temp_Ori; Pr_Clip(count)=temp_Clip; count=count+1; end semilogy(PAPRO_Start:step:PAPRO_End,Pr_Ori, '-y*',PAPRO_Start:step:PAPRO_End,Pr_Clip, '-g*'); axis([PAPRO_Start PAPRO_End 10^-4 1]); grid on; hold on:
这段代码是一个用MATLAB实现的PAPR(峰均比)模拟程序,用于比较原始OFDM信号和进行裁剪处理后的信号的PAPR性能。其中QPSK调制信号通过FFT变换转换成OFDM信号,然后进行裁剪处理,最终统计PAPR性能。在程序的最后,使用semilogy函数将原始信号和裁剪后信号的PAPR性能进行比较并绘制成图像。
papr_eight3dinsert(m_2d,m_3d,num_subcarrier)
`papr_eight3dinsert(m_2d, m_3d, num_subcarrier)` 是一个函数,用于将2D的信号转换为3D的信号。函数的输入参数为`m_2d`(2D的信号矩阵)、`m_3d`(输出的3D信号矩阵)和`num_subcarrier`(子信道的数量)。
函数的功能是通过将原始的2D信号分割成多个子信道,然后在垂直方向上复制这些子信道,以形成3D信号。具体的步骤如下:
1. 首先,获取2D信号矩阵的行数和列数,并计算出每个子信道的行数。
2. 在`m_3d`中创建一个新的矩阵,行数为`num_subcarrier`,列数与2D信号的列数相同。
3. 使用循环遍历子信道的数量,并进行以下操作:
- 在`m_3d`中选择第`i`个子信道,并将其与`m_2d`中的对应的子信道复制。
- 复制完成后,就可以在垂直方向上形成3D信号的主体部分。
4. 最后,返回经过转换的3D信号矩阵`m_3d`。
这个函数的作用在于将2D信号转换为3D信号,使得信号可以在不同的维度上进行处理和分析。这在一些通信和图像处理应用中,能够更好地利用信号的空间相关性,提高信号处理的效果。