%**********File：系統(tǒng)辨識******************
clear?all
SNR_DB=[0:1:12];
sum=1000;
%********************產(chǎn)生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）;

%*******************系統(tǒng)辨識（LMS算法）***********************

N=13;%濾波器長度
umax=2/（N*mean（（abs（x））.^2））;%步長應(yīng)小于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