%OFDM系統(tǒng)參數(shù)設(shè)置
clc;
clear?all;
NFFT=64;?????????????????????????????%FFT長(zhǎng)度
G=0;?????????????????????????????????%保護(hù)間隔長(zhǎng)度
M_ary=4;?????????????????????????????%進(jìn)制數(shù)
P_A=sqrt（2）;?????????????????????????%導(dǎo)頻符號(hào)的幅度
D_t=4;???????????????????????????????%時(shí)域內(nèi)導(dǎo)頻序列的間隔，即一個(gè)訓(xùn)練序列與另一個(gè)訓(xùn)練序列的間隔
t_a=50*10^（-9）;??????????????????????%HiperLAN/2中的采樣間隔
%蒙特卡洛信道參數(shù)設(shè)置
symbol_duration=NFFT*t_a;????????????%一個(gè)OFDM信號(hào)的時(shí)間
number_of_summations=40;?????????????%蒙特卡洛方法中的調(diào)和系數(shù)
f_dmax=50.0;?????????????????????????%最大多普勒頻移
load?h_decimation.am?-ascii;?????????%數(shù)據(jù)文件，存儲(chǔ)多徑信道p[k]與套k參數(shù)
h11_initial=h_decimation;????????????%h11:發(fā)送天線1到接收天線1信道參數(shù)
h12_initial=h_decimation;????????????%h12:發(fā)送天線1到接收天線2信道參數(shù)
h21_initial=h_decimation;????????????%h21:發(fā)送天線2到接收天線1信道參數(shù)
h22_initial=h_decimation;????????????%h22:發(fā)送天線2到接收天線2信道參數(shù)
N_P=length（h_decimation）;
NofOFDMSymbol=1000;??????????????????%OFDM符號(hào)個(gè)數(shù)，包括導(dǎo)頻OFDM符號(hào)與數(shù)據(jù)OFDM符號(hào)
No_Of_OFDM_Data_Symbol=NofOfdmsYMBOLceil（NofOFDMSymbol/D_t）;%數(shù)據(jù)OFDM符號(hào)個(gè)數(shù)
length_data=（No_Of_OFDM_Data_Symbol）*NFFT;%總的數(shù)據(jù)長(zhǎng)度，例如采用4QAM則是四進(jìn)制數(shù)據(jù)的長(zhǎng)度
Number_Relz=100;?????????????????????%循環(huán)數(shù)
ser_relz=[];?????????????????????????%存儲(chǔ)每個(gè)循環(huán)的誤碼率
for?number_of_relialization=1:Number_Relz;%循環(huán)開始
%產(chǎn)生信道中的隨機(jī)系數(shù)
u11=rand（N_Pnumber_of_summations）;
u12=rand（N_Pnumber_of_summations）;
u21=rand（N_Pnumber_of_summations）;
u22=rand（N_Pnumber_of_summations）;
%發(fā)送的二進(jìn)制比特流
source_data1=randint（length_data2）;
source_data2=randint（length_data2）;
%轉(zhuǎn)換為四進(jìn)制信號(hào)
symbols1=bi2de（source_data1）;
symbols2=bi2de（source_data2）;
%進(jìn)行4QAM映射
QASK_Symbol1=dmodce（symbols111‘qask‘M_ary）;
%存放第一個(gè)天線發(fā)送的數(shù)據(jù)
QASK_Symbol2=dmodce（symbols211‘qask‘M_ary）;
%存放第二個(gè)天線發(fā)送的串行數(shù)據(jù)
%天線IIFFT前的信號(hào)，為IFFT點(diǎn)的并行數(shù)據(jù)，即將串行數(shù)據(jù)轉(zhuǎn)換為并行數(shù)據(jù)
Data_Patterm1=[];%為IFFT轉(zhuǎn)唄輸入信號(hào)
m=0;
for?i=0:No_Of_OFDM_Data_Symbol-1;
????QASK_tem=[];
????for?n=1:NFFT;
????????QASK_tem=[QASK_temQASK_Symbol1（i*NFFT+n）];
????end;
????Data_Pattern1=[Data_Pattern1;QASK_tem];
????clear?QASK_tem;
end;
%天線2IFFT前的信號(hào)，為IFFT點(diǎn)的并行數(shù)據(jù)，即將串行數(shù)據(jù)轉(zhuǎn)換為并行數(shù)據(jù)
Data_Pattern2=[];
m=0;
for?i=0:No_Of_OFDM_Data_Symbol-1;
????QASK_tem=[];
????for?n=1:NFFT;
????????QASK_tem=[QASK_temQASK_Symbol2（i*NFFT+n）];
????end;
????Data_Pattern2=[Data_Pattern2;QASK_tem];
????clear?QASK_tem;
end;
%天線1的導(dǎo)頻信號(hào)
PP_A1=[];
for?m=0:NFFT-1;
????PP_A1=[PP_A1P_A*exp（j*D_f*pi*（m）^2/NFFT）];
end;
%FFT矩陣
F=[];
for?k=0:NFFT-1;
????W_tem=[];
????for?n=0:NFFT-1;
????????W_tem=[W_temexp（-j*2*pi*n*k/NFFT）];
????end;
????F=[F;W_tem];
end;
%LS估計(jì)系數(shù)
PP=[diag（PP_A1）*F（:1:n_p）diag（PP_A2）*F（:1:N_P）];
%參看Q
Q=inv（PP*PP）;%%%%%%%%%%%%%%%%%%%%%%%%參看式5—80中的QQ的逆
%天線1發(fā)送信號(hào)。將導(dǎo)頻插入并行OFDM數(shù)據(jù)，一個(gè)塊狀導(dǎo)頻OFDM符號(hào)后3個(gè)數(shù)據(jù)
%OFDM符號(hào)
TS1_BeforeIFFT=Insert_PilotSymbol（PP_A2Data_Pattern2D_tNofOFDMSymbolNFFT）;
ser_without_isic=[];
snr_min=0;%最小性噪比
snr_max=50;%最大性噪比
step=5;%信噪比步長(zhǎng)
for?snr=snr_min:step:snr_max;%信噪比循環(huán)
????%對(duì)第一個(gè)發(fā)送天線
????rsl1_frame=[];
????%存放并行的經(jīng)過多徑與加性高斯白噪聲信道的第一個(gè)發(fā)送天線到第一個(gè)接收天線的每個(gè)OFDM符號(hào)
????rsl2_frame=[];
????%存放并行的經(jīng)過多徑與加性高斯白噪聲信道的第一個(gè)發(fā)送天線到第二個(gè)接收天線的每個(gè)OFDM符號(hào)
????initial_time=0;%初始時(shí)間
????for?i=0:NofOFD