lms.m

最小均方误差自适应算法

clc;
clear all;
close all;
hold off
hh=5;
M=300;                  
w=randn(1,M);
x=zeros(1,M);           %期望响应
xx=zeros(1,M+4);
xx(1,1)=0;xx(1,2)=0;xx(1,3)=0;xx(1,4)=0;
for i=5:M+4
    xx(1,i)=w(1,i-4)+1.352*xx(1,i-1)-1.338*xx(1,i-2)+0.662*xx(1,i-3)-0.240*xx(1,i-4);
end
x(1,1:M)=xx(1,5:M+4);
v=randn(1,M);            %白噪声 
y=x+v;                   %滤波器输入
w=zeros( hh, 1 ) ;%初始化
for n=hh : M 
u=x(n:-1:n-hh+1)';
s(n)=w'*u;%系统输出
e(n) =x(n)-s(n) ;%误差
mu=0.0020
w = w + mu * u * e(n) ;%迭代方程
end
%画仿真图
plot(s);
hold on
plot(x,'r');
hold on
title('系统输出和期望输出对比') ;
figure
plot((abs(e)));
title('估计误差');
hold off