LMMSE

线性最小均方误差估计LMMSE的matlab算法

%**********File：系统辨识******************
clear all
SNR_DB=[0:1:12];
sum=1000;
%********************产生QPSK信号***********************
data=randsrc（sum2[0 1]）;
[a1b1]=find（data（:1）==0&data（:2）==0）;
message（a1）=（-1-j）*sqrt（1/2）;
[a2b2]=find（data（:1）==0&data（:2）==1）;
message（a2）=（-1+j）*sqrt（1/2）;
[a3b3]=find（data（:1）==1&data（:2）==0）;
message（a3）=（1-j）*sqrt（1/2）;
[a4b4]=find（data（:1）==1&data（:2）==1）;
message（a4）=（1+j）*sqrt（1/2）;
x=message;
%scatterplot（x）

h=[0.04-0.050.07-0.21-0.50.720.3600.210.030.07]; 
d=conv（xh）;
%scatterplot（d）;
d=awgn（d20）;%叠加高斯白噪声 信噪比为50dB
%scatterplot（d）;

%*******************系统辨识（LMS算法）***********************

N=13;%滤波器长度
umax=2/（N*mean（（abs（x））.^2））;%步长应小于umax，计算出umax=0.13333
u=0.06;%取u=0.06
W=zeros（1N）;
e=zeros（1su

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1759  2010-12-16 22:26  LMMSE.m

----------- ---------  ---------- -----  ----

                 1759                    1