clear?all;
close?all;
fprintf（‘OFDM信道估計仿真\n\n‘）;
carrier_cnt=64;%-----------載波數(shù)目
num_symbol=50;%--------------OFDM符號個數(shù)
Guard=8;%--------------------循環(huán)前綴
pilot_Inter=8;%--------------導(dǎo)頻間隔
modulation_mode=16;%---------調(diào)制方式
SNR=[10:10:200];%-------------信噪比取值
NumLoop=15;%-----------------循環(huán)次數(shù)
num_bit_err=zeros（length（SNR）NumLoop）;
num_bit_err_dft=zeros（length（SNR）NumLoop）;
num_bit_err_ls=zeros（length（SNR）NumLoop）;
MSE=zeros（length（SNR）NumLoop）;
MSE1=zeros（length（SNR）NumLoop）;
MSE2=zeros（length（SNR）NumLoop）;
%%%%%%%%%%%%%%%%%%%%%%%主程序循環(huán)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for?c1=1:length（SNR）
????fprintf（‘\n\n\n仿真信噪比為%f\n\n‘SNR（c1））;
????for?num1=1:NumLoop
????????%---------------產(chǎn)生發(fā)送的隨機序列——————————————-
????????BitsLen=carrier_cnt*num_symbol;
????????BitsTx=randint（1BitsLen）;
????????%---------------符號調(diào)制---------------------------------------
????????Modulated_Sequence=qammod（BitsTxmodulation_mode）;
????????%---------------導(dǎo)頻格式---------------------------------------
????????pilot_len=carrier_cnt;???????????????????%導(dǎo)頻長度等于載波數(shù)目
????????pilot_symbols=round（rand（1pilot_len））;??%產(chǎn)生導(dǎo)頻符號長度為64的01序列
????????for?i=1:pilot_len
????????????if?pilot_symbols（1i）==0
????????????????pilot_symbols（1i）=pilot_symbols（1i）-1;
????????????else
????????????????pilot_symbols（1i）=pilot_symbols（1i）;
????????????end
????????end
????????pilot_symbols=pilot_symbols‘;%產(chǎn)生一列導(dǎo)頻數(shù)據(jù)，共有64個
????????%----------------計算導(dǎo)頻和數(shù)據(jù)數(shù)目----------------------------
????????num_pilot=ceil（num_symbol/pilot_Inter）;%導(dǎo)頻序列的數(shù)目
????????if?rem（num_symbolpilot_Inter）==0
????????????num_pilot=num_pilot+1;
????????end%如果符號數(shù)目可以被導(dǎo)頻間隔整除，還要在末尾加一個導(dǎo)頻
????????num_data=num_symbol+num_pilot;%總數(shù)據(jù)量為符號數(shù)目+導(dǎo)頻符號數(shù)目
????????%----------------導(dǎo)頻位置計算----------------------------------
????????pilot_Indx=zeros（1num_pilot）;%7段導(dǎo)頻?1?10?19?28?37?46?57?
????????Data_Indx=zeros（1num_pilot*（pilot_Inter+1））;
????????for?i=1:num_pilot-1
????????????pilot_Indx（1i）=（i-1）*（pilot_Inter+1）+1;
????????end
????????pilot_Indx（1num_pilot）=num_data;
????????%pilot_Indx完成了導(dǎo)頻位置的選定1?10?19?28?37?46?57
????????for?j=0:num_pilot
????????????Data_Indx（1（1+j*pilot_Inter）:（j+1）*pilot_Inter）=（2+j*（pilot_Inter+1））:（（j+1）*（pilot_Inter+1））;
????????end
????????Data_Indx=Data_Indx（11:num_symbol）;%指定了數(shù)據(jù)的位置%完成了數(shù)據(jù)位置的設(shè)定
????????%----------------導(dǎo)頻插入-------------------------------------
????????piloted_ofdm_syms=zeros（carrier_cntnum_data）;%一個二維數(shù)組，行數(shù)為載波數(shù)目，列數(shù)為符號數(shù)+導(dǎo)頻數(shù)
????????piloted_ofdm_syms（:Data_Indx）=reshape（Modulated_Sequencecarrier_cntnum_symbol）;
????????p