自适应滤波算法LMS matlab

matlab 自适应滤波LMS算法 绘制等值线，绘制学习曲线，包含实验报告

clear allclose allclc;
%产生r信号%
pr=0.01;
r=randn（150000）;
r=r*sqrt（pr）;
%产生信号xk与期望信号dk%
N=1000;%训练长度
k=0:1:N-1;
xk=sin（2*pi*k/16）;
dk=2*cos（2*pi*k/16）;
r=r（1:N）;

%计算自相关矩阵%
R0=（1/N）*xk*xk‘+pr;
R1=（1/N）*xk*[0xk（1:N-1）]‘;
R=[R0 R1;R1 R0];
%计算互相关矩阵%
Rdx0=（1/N）*dk*xk‘;
Rdx1=（1/N）*dk*[0xk（1:N-1）]‘;
Rdx=[Rdx0Rdx1]‘;
%计算dk的均方%
Edj2=（1/N）*dk*dk‘;
%产生输入信号%
xk=xk+r;%输入信号

%计算最佳权矢量%
Wopt=inv（R）*Rdx;
%计算最小均方误差%
ksimin=Edj2+Wopt‘*R*Wopt-2*Rdx‘*Wopt;
%初始化均方误差%
ksi=zeros（1N）;

%LMS迭代算法%
M=2;%FIR滤波器的长度
u1=0.1;
w0=[0 0];
[hyw]=lms（xkdku1Mw0）;
figure（1）plot（w（1:）w（2:））axis（[-28-100]）hold ongrid on;
text（w0（1）w0（2）‘u=0.1k=1000‘）;
for n=1:1:N
    ksi（n）=ksimin+（w（:n）-Wopt）‘*R*（w（:n）-Wopt）;
end;
figure（2）plot（kksi‘r‘）hold ongrid on;%绘制学

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-11-08 22:27  adaptive filteryf\
     文件        2316  2012-11-05 16:55  adaptive filteryf\lab_2.asv
     文件        2334  2012-11-01 21:00  adaptive filteryf\lab_2.m
     文件        1425  2012-11-01 15:14  adaptive filteryf\lms.asv
     文件        1470  2012-11-01 15:15  adaptive filteryf\lms.m
     文件      229888  2012-11-08 22:27  adaptive filteryf\自适应滤波实验报告.doc