%IEEE?802.11a系統的仿真
clear?all?
%%%參數設置部分%%%

Nsp?=?52;?%系統子載波數，（不包括直流載波）
Nfft?=?64;?%FFT長度
Ncp?=?16;?%循環前綴長度
Ns?=?Nfft+Ncp;?%一個完整OFDM符號長度
noc?=?53;?%包含直流載波的總的子載波數
Nd?=?6;?%每幀包含的OFDM符號數（不包括訓練符號）
M1?=?4;?%QPSK調制
M2?=?16;?%16-QAM調制
sr?=?250000;?%OFDM符號速率
EbN0?=?0:2:30;?%歸一化信噪比
Nfrm?=?10000;?%每種信噪比下的仿真幀數
ts?=?1/sr/Ns;?%OFDM符號抽樣時間間隔
t?=?0:ts:（Ns*（Nd+1）*Nfrm-1）*ts;?%抽樣時刻
fd?=?100;?%最大多普勒頻移
h?=?rayleigh（fdt）;?%生成單徑Rayleigh衰落信道

%訓練符號頻域數據，采用802.11a中的長訓練符號數據
Preamble=[1?1?-1?-1?1?1?-1?1?-1?1?1?1?1?1?1?-1?-1?1?1?-1?1?-1?1?1?1?1?...?
????1?-1?-1?1?1?-1?1?-1?1?-1?-1?-1?-1?-1?1?1?-1?-1?1?-1?1?-1?1?1?1?1];?
Preamble1=zeros（1Nfft）;?
Preamble1（2:27）=Preamble（27:end）;?%訓練符號重排后的數據????
Preamble1（39:end）=Preamble（1:26）;?
preamble1=ifft（Preamble1）;?%訓練符合時域數據????
preamble1=[preamble1（Nfft-Ncp+1:end）?preamble1];?%加入循環前綴
?
%%%仿真循環%%%
for?ii=1:length（EbN0）?
?????
????%%%發射機部分%%%
????msg1=randsrc（NspNd*Nfrm[0:M1-1]）;?%QPSK信息數據???
????msg2=randsrc（NspNd*Nfrm[0:M2-1]）;?%16-QAM信息數據???
????data1=pskmod（msg1M1pi/4）;?%QPSK調制???
????data2=qammod（msg2M2）/sqrt（10）;?%16-QAM調制并歸一化
?
????data3=zeros（NfftNd*Nfrm）;?%根據FFT要求，對數據重排???????????
????data4=zeros（NfftNd*Nfrm）;?
?
????data3（2:27:）=data1（27:end:）;?%重排后的數據
????data3（39:end:）=data1（1:26:）;?
?
????data4（2:27:）=data2（27:end:）;?%重排后的數據
????data4（39:end:）=data2（1:26:）;?
?????????
????clear?data1?data2;?%清除不需要的臨時變量??
?
????data3=ifft（data3）;?%IFFT變換
????data4=ifft（data4）;?
?
????data3=[data3（Nfft-Ncp+1:end:）;data3];?%加入循環前綴
????data4=[data4（Nfft-Ncp+1:end:）;data4];?
?
????spow1=norm（data3‘fro‘）.^2/（Nsp*Nd*Nfrm）;?%計算數據符號能量????????????
????spow2=norm（data4‘fro‘）.^2/（Nsp*Nd*Nfrm）;?
?????????
????data5=zeros（Ns（Nd+1）*Nfrm）;?%加入訓練符號????
????data6=data5;?
????for?indx=1:Nfrm?
????????data5（:（indx-1）*（Nd+1）+1）=preamble1.‘;?
????????data5（:（indx-1）*（Nd+1）+2:indx*（Nd+1））=data3（:（indx-1）*Nd+1:indx*Nd）;?
?????????????
????????data6（:（indx-1）*（Nd+1）+1）=preamble1.‘;?
????????data6（:（indx-1）*（Nd+1）+2:indx*（Nd+1））=data4（:（indx-1）*Nd+1:indx*Nd）;?
????end?
?????????
????clear?data3?data4?
?????????
????data5=reshape（data51Ns*（Nd+1）*Nfrm）;?%并串變換??
????data6=reshape（data61Ns*（Nd+1）*Nfrm）;?
???????
????sigma1=sqrt（1/2*spow1/log2（M1）*10.^（-EbN0（ii）/10））;?%根據EbN0計算噪聲標準差????
????sigma2=sqrt（1/2*spow2/log2（M2）*10.^（-EbN0（ii）/10））;?
?????????
?????????
????for?indx=1:Nfrm?
????????dd1=data5（（indx-1）*Ns*（Nd+1）+1:indx*Ns*（Nd+1））;?%當前幀的發射數據??????
????????dd2=data6（（indx-1）*Ns*（Nd+1）+1:indx*Ns*（Nd+1））;?
???????????
????????hh=h（（indx-1）*Ns*（Nd+1）+1:indx*Ns*（Nd+1））;?%當前幀對應的信道參數????
?????????
????????%信號通過單徑Rayleigh衰落信道，并加入高斯白噪聲?????????????????
????????r1=hh.*dd1+sigma1*（randn（1length（dd1））+j*randn（1length（dd1）））;????
????????r2=hh.*dd2+sigma2*（randn（1length（dd2））+j*randn（1length（dd2）））;?
?????????????
????????r1=reshape（r1NsNd+1）;?%串并變換
????????r2=reshape（r2NsNd+1）;?
?
????????r1=r1（Ncp+1:end:）;?%移除循環前綴
????????r2=r2（Ncp+1:end:）;??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2020-04-09?13:13??Matlab仿真IEEE?802.11a系統\
?????文件????????5352??2020-04-09?12:50??Matlab仿真IEEE?802.11a系統\OFDM_802.11a.m
?????文件?????????479??2020-04-09?13:06??Matlab仿真IEEE?802.11a系統\rayleigh.m
?????文件???????32891??2020-04-09?13:05??Matlab仿真IEEE?802.11a系統\Simulation?results.png
?????文件??????228126??2020-04-09?13:11??Matlab仿真IEEE?802.11a系統\System?parameters.png