91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 2KB
    文件類(lèi)型: .rar
    金幣: 2
    下載: 3 次
    發(fā)布日期: 2021-04-21
  • 語(yǔ)言: Matlab
  • 標(biāo)簽: PAPR??TR??OFDM??

資源簡(jiǎn)介

預(yù)留子載波法(TR)降低PAPR的matlab算法

資源截圖

代碼片段和文件信息

clear;
clc;
carr=1024;???%子載波個(gè)數(shù)
sym_c=10;???%每個(gè)子載波含有的OFDM符號(hào)數(shù)
bit_sym=2;??%每個(gè)符號(hào)含有的比特?cái)?shù),QPSK調(diào)制
IFFT_n=256;?%IFFT點(diǎn)數(shù)
r=1/10;?????%保護(hù)間隔和OFDM數(shù)據(jù)的比例;
%產(chǎn)生信號(hào)
sum=carr*sym_c*bit_sym;
colume=sum/(2*carr);
signal=rand(1sum)>0.5;???%初始信號(hào)
%QPSK調(diào)制QPSK_sig里面存放的是調(diào)制后的信號(hào),數(shù)目sumQ
sumQ=sum/2;%OFDM符號(hào)總數(shù)
imag=sqrt(-1);??????????????????????????%?虛部?j
QPSK=[-1+imag-1-imag1+imag1-imag];???%創(chuàng)建QPSK?映射表
SIGNAL=zeros(1sumQ);??????%計(jì)算并存放調(diào)制前的十進(jìn)制數(shù)據(jù)(將二進(jìn)制轉(zhuǎn)為十進(jìn)制數(shù)據(jù),第20行)
QPSK_sig=zeros(1sumQ);????%存放調(diào)制后的QPSK信號(hào)
for?n=1:sumQ
?SIGNAL(n)=signal(2*n-1)*2+signal(2*n);?%將二進(jìn)制換算成十進(jìn)制
end
for?i=1:sumQ
????if?SIGNAL(i)==0;
????????QPSK_sig(i)=QPSK(1);
????elseif?SIGNAL(i)==1;
????????????QPSK_sig(i)=QPSK(2);
????elseif?SIGNAL(i)==2;
?????????QPSK_sig(i)=QPSK(3);
????elseif?SIGNAL(i)==3;
???????????????QPSK_sig(i)=QPSK(4);
????end
end?????????????????????
%串/并轉(zhuǎn)換??????計(jì)算第i個(gè)載波上面的信號(hào)to_par(i,:)
colume=sumQ/carr;%每個(gè)子載波的符號(hào)數(shù)
for?i=1:carr??%?carr載波個(gè)數(shù)
????for?j=1:colume;??
????????to_par(ij)=QPSK_sig(carr*(j-1)+i);
????end
end
colume=sumQ/carr;
%?%?to_par=reshape(QPSK_sigcarrcolume);
%每個(gè)子載波上進(jìn)行?IFFT變換??(調(diào)制后的QPSK信號(hào)進(jìn)行IFFT)
for?j=1:colume
??????y(:j)=ifft(to_par(:j));?%每一列進(jìn)行一次IFFT?
end
%?%?y=ifft(to_par);
yr=real(y);???%實(shí)部
yi=(y-yr)*sqrt(-1)*(-1);%虛部??
y_=yr+yi*sqrt(-1);
y_series=reshape(y_1(colume)*carr);
y1=y_(:1);
y2=abs(y1);%采樣點(diǎn)幅值
y2_power=y2.^2;%采樣點(diǎn)功率
y2_power_mean=mean(y2_power);
sigma=sqrt(y2_power_mean);
y_PAPR=10*log10(y2_power./sigma^2);
PAPR_max=max(y_PAPR);
PAPR_mean=mean(y_PAPR);
%*******************************************進(jìn)行TR消峰處理***********************
L=32;%預(yù)留子載波數(shù)目
u=0.3;%設(shè)置步進(jìn)參數(shù)
for?i=1:carr;?????????????????%生成IFFT變化矩陣
????for?j=1:carr;
????????Q(ij)=exp(sqrt(-1)*2*pi*(i-1)*(j-1)/carr)/sqrt(carr);
????end
end
Ql=Q(:carr-L+1:carr);???????????%矩陣Q尖
x(:1)=y1;%步驟一:算法初始化
Imax=16;????????????????????%迭代次數(shù)
i=0;
A=y2_power_mean;%門(mén)限值
U=zeros(carr1);
flag=1;
while((i????flag=0;
????????for?n=1:carr
????????????if?abs(x(n))>A
????????????????q_nrow=Ql(n:);
????????????????U=U+(x(n)-A*exp(sqrt(-1)*phase(x(n))))*Ql*q_nrow‘;
????????????????flag=1;
????????????else
????????????????U=U+zeros(carr1);
????????????end
????????end
????????x=x-u*U;
????i=i+1;
end
x_power=abs(x).^2;
x_power_mean=mean(x_power);
x_sigma=sqrt(x_power_mean);
x_PAPR=10*log10(x_power./x_sigma^2);
subplot(211);
plot(y2_power‘b-‘);
xlim([01024]);
ylim([00.02]);
grid?on;
subplot(212);
plot(x_power‘r-‘);
xlim([01024]);
ylim([00.02]);
grid?on;
figure;
%************************************統(tǒng)計(jì)CCDF特性*************************************
PAPR0=2:0.5:9;
for?k=1:length(PAPR0)
????count(k)=0;
????x_count(k)=0;
end
for?l=1:carr;
????for?k=1:length(PAPR0)
????????if?y_PAPR(l)>PAPR0(k)
????????????count(k)=count(k)+1;
????????end
????????if?x_PAPR(l)>PAPR0(k)
????????????x_count(k)=x_count(k)+1;
????????end
????end
end
CCDF=count/carr;
x_CCDF=x_count/carr;
subplot;
semilogy(PAPR0CCDF‘.

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????3401??2012-12-11?16:45??ofdm_tr.m

-----------?---------??----------?-----??----

?????????????????3401????????????????????1


評(píng)論

共有 條評(píng)論