LMMSE算法实现

len=500;
R=randint（1 len）;
I=randint（1 len）;
Real=2*R-1;
Imag=2*I-1;
X=（Real+1j*Imag）/sqrt（2）;
scatterplot（X）;
title（‘发送端QPSK信号星座图‘）;

h=[0.04 -0.05 0.07 -0.21 -0.5 0.72 0.36 0 0.21 0.03 0.07];
X1=conv（Xh）;
X2=zeros（1len）;
for i=1:len
    X2（i）=X1（i）;
end
scatterplot（X2）;
title（‘经过线性系统QPSK信号星座图‘）;

SNR=20;
Y_R=awgn（real（X2）SNR）;
Y_I=awgn（imag（X2）SNR）;
Y=Y_R+1j*Y_I;
scatterplot（Y）;
title（‘叠加噪声后QPSK信号星座图‘）;

N=length（h）;
H1 = zero