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