clear;
clc;
carr=1024;???%子載波個數
sym_c=10;???%每個子載波含有的OFDM符號數
bit_sym=2;??%每個符號含有的比特數，QPSK調制
IFFT_n=256;?%IFFT點數
r=1/10;?????%保護間隔和OFDM數據的比例；
%產生信號
sum=carr*sym_c*bit_sym;
colume=sum/（2*carr）;
signal=rand（1sum）>0.5;???%初始信號
%QPSK調制QPSK_sig里面存放的是調制后的信號，數目sumQ
sumQ=sum/2;%OFDM符號總數
imag=sqrt（-1）;??????????????????????????%?虛部?j
QPSK=[-1+imag-1-imag1+imag1-imag];???%創建QPSK?映射表
SIGNAL=zeros（1sumQ）;??????%計算并存放調制前的十進制數據（將二進制轉為十進制數據，第20行）
QPSK_sig=zeros（1sumQ）;????%存放調制后的QPSK信號
for?n=1:sumQ
?SIGNAL（n）=signal（2*n-1）*2+signal（2*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?????????????????????
%串/并轉換??????計算第i個載波上面的信號to_par（i，：）
colume=sumQ/carr;%每個子載波的符號數
for?i=1:carr??%?carr載波個數
????for?j=1:colume;??
????????to_par（ij）=QPSK_sig（carr*（j-1）+i）;
????end
end
colume=sumQ/carr;
%?%?to_par=reshape（QPSK_sigcarrcolume）;
%每個子載波上進行?IFFT變換??（調制后的QPSK信號進行IFFT）
for?j=1:colume
??????y（:j）=ifft（to_par（:j））;?%每一列進行一次IFFT?
end
%?%?y=ifft（to_par）;
yr=real（y）;???%實部
yi=（y-yr）*sqrt（-1）*（-1）;%虛部??
y_=yr+yi*sqrt（-1）;
y_series=reshape（y_1（colume）*carr）;
y1=y_（:1）;
y2=abs（y1）;%采樣點幅值
y2_power=y2.^2;%采樣點功率
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）;
%*******************************************進行TR消峰處理***********************
L=32;%預留子載波數目
u=0.3;%設置步進參數
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;????????????????????%迭代次數
i=0;
A=y2_power_mean;%門限值
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;
%************************************統計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‘.

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

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

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

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