clc;
clear?all;
%生成訓(xùn)練序列
%用BPSK調(diào)制
X=zeros（6464）;
d=rand（641）;
??????for?i=1:64
???????if（d（i）>=0.5）
???????????d（i）=+1;
???????else
???????????d（i）=-1;
???????end
????end
?for?i=1:64
?????X（ii）=d（i）;
?end
%計(jì)算出信道向量G
?%信道特性
??tau=[0.5?3.5];
for?k=1:64
??????s=0;
??????for?m=1:2
?????????s=s+（exp（-j*pi*（1/64）*（k+63*tau（m）））?*?（（?sin（pi*tau（m））?/?sin（pi*（1/64）*（tau（m）-k）））））;
??????end
g（k）=s/sqrt（64）;
end
G=g‘;%信道向量
H=fft（G）;%?頻域
XFG=X*H;
n1=ones（641）;
n1=n1*0.000000000000000001i;%加入復(fù)高斯白噪聲
noise=awgn（n18）;%設(shè)置SNR為8db
variance=var（noise）;
N=fft（noise）;
Y=XFG+N;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%求出G的自相關(guān)矩陣Rgg
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
gg=zeros（6464）;
for?i=1:64
????gg（