本科毕设，做基于802.11ad的信道估计，想实现基于LS的信道估计算法，信道估计序列是有64位格蕾a序列做实部，b序列做虚部的UW，重复9次，在接受段移除cp，对剩余8段做时域平均。后利用ls估计通用公式得到频域的信道估计CFRE=Ufcn*yfc，ufcn为64位UWfft后作为对角线矩阵的逆矩阵，yfc为平均后的接收信号的fft，具体代码如下，代码的信道估计结果完全不对，但自己实在找不出是哪里有问题，求大神指教

clear all
clc
seed = 0;

j = sqrt(-1);
randn('state',seed);
rand('state',seed);

[imp_res] = cr_ch_model;
[imp_reszz]=[imp_res]';

scpreamble=testPHY_Preamble_Gen('SCPHY');%sc的发送序列
scchipout=conv(scpreamble,imp_reszz)';%经过信道，1024*1
sccesr=scchipout(65:576);%提取ces
UW1=sccesr(1:64);
UW2=sccesr(65:128);
UW3=sccesr(129:192);
UW4=sccesr(193:256);
UW5=sccesr(257:320);
UW6=sccesr(321:384);
UW7=sccesr(385:448);
UW8=sccesr(449:512);%8个重复的UW
ycaverage=zeros(64,1);
for n=1:64
ycaverage(n)=(UW1(n)+UW2(n)+UW3(n)+UW4(n)+UW5(n)+UW6(n)+UW7(n)+UW8(n))/8;
end%做时域平均
yfc=fft(ycaverage,64);%做FFT
UW64=UW(64);
SCPHY_CEF_Seq=[UW64];
UW64s=SCPHY_CEF_Seq.*exp(j*pi/2.*(0:length(SCPHY_CEF_Seq)-1));
ufcfft=fft(UW64s,64);
Ufc=diag(ufcfft);
Ufcn=inv(Ufc);
CFRE=Ufcn*yfc;%频域
CIRE=fft(CFRE);