lmsfilter.m

来自「各种滤波器」· M 代码 · 共 19 行

function [y,e,w]=LMSfilter(x,d,mu,N)
% 用LMS算法滤波
% 输入参数分别是:输入信号,期望信号,搜索步长,滤波器长度
% 输出参数分别是:输出信号,误差信号,滤波器系数
L=length(x);
y=zeros(1,L);
% LMS滤波器 
w=zeros(1,N);   % 滤波器系数
X=zeros(1,N);   % 延时区
e=0;
% 滤波，迭代
for i=1: L
    X=[x(i),X(1:N-1)];
    y(i)=w*X';
    e=d(i)-y(i);
    w=w+mu*e*X;
end
% 误差信号
e=d-y;