频偏纠正matlab实现,Matlab关于OFDM同步定时估计和频偏估计的算法

各位高手你们好，在下有理了。我最近在写有关ofdm的同步的问题，主要是SC算法，

我已经把M(d)的图像仿真出来了，也找到了M(d)取最大时的那个点，还有小数频偏也算算出来了，但是我不会算'定时估计均方差'和‘频率估计均方差’，请给位高手指点一下 在下感激不尽!!!

我写的代码 ，请高手们看一下那儿出错了：

前半部分仿真可能是对的 ，但是后面求频偏估计均方误差与定时估计均方误差的图像时，就不对了

希望告诉们能给我解答一下

clear

clc

gilen=32;

para=512;   % Number of parallel channel to transmit (points)

fftlen=512; % FFT length

noc=512;    % Number of carrier

nd=6;       % Number of information OFDM symbol for one loop

ml=2;       % Modulation level : QPSK

sr=250000;  % Symbol rate

br=sr.*ml;  % Bit rate per carrier

%gilen=32;   %Length of guard interval (points)

fullen=gilen+para;  %ofdm符号总长

deltad=250;  %时延

deltaf=0.25;  %频偏

deltaderr=0;%时延均方误差

deltaferr=0;%频偏均方误差

derr=zeros(1,11);%时延均方误差

ferr=zeros(1,11);%频偏均方误差

%f=-1;

error=zeros(1,11);

o=0

%************************** main loop part **************************

for ebn0=0:2:20

o=o+1;

%************************** main loop part **************************

nloop=100;  % Number of simulation loops

noe = 0;    % Number of error data

nod = 0;    % Number of transmitted data

eop=0;      % Number of error packet

nop=0;