OFDM同步算法中的最大似然函数法程序

%%%%最大似然ML  最大似然方法联合实现符号定时同步和载波同步仿真
clear all;
%********************** preparation part ***************************
para=128;   % Number of parallel channel to transmit （points） %并行信道传输个数
fftlen=128; % FFT length %fft 长度
noc=1024;    % Number of carrier %载波的个数
nd=6;       % Number of information OFDM symbol for one loop%OFDM symbol 循环个数的说明？
ml=2;       % Modulation level : QPSK %调制方式：QPSK
sr=250000;  % Symbol rate  %符号的比特率
br=sr.*ml;  % Bit rate per carrier %每个载波的比特率?
gilen=128;   % Length of guard interval （points）%%保护间隔的长度
ebn0=3;     % Eb/N0  %Eb为单位比特的平均信号的能量；n0为噪声的单边功率普密度
SNR=15;     % ENR  %信噪比
T=10^（-6）;   % Sampling time : 1us%%采样时间1us 抽样频率1024kHz
%************************** transmitter发射机 *********************************
%************************** Data generation 数据生成程序**************************** 
seldata=rand（1para*nd*ml）>0.5;  %  rand : built in function  %产生1x（128*6*2）=1X1536的数据
M=length（seldata）/noc;
%****************** Serial to parallel conversion串并转换 ***********************
paradata=reshape（seldataparand*ml）; %  reshape : built in function
%************************** QPSK modulation %QPSK 调制***************************** 
[ichqch]=qpskmod（paradataparandml）;
kmod=1/sqrt（2）; %  sqrt : built in function
ich1=ich.*kmod;%实部的数据
qch1=qch.*kmod;%虚部的数据
%******************* IFFT ************************
x=ich1+qch1.*i;
y=ifft（x）;      %  ifft : built in function
ich2=real（y）;   %  real : buil