用MATLAB实现ofdm的ls估计算法

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