自适应滤波器Matlab实现
直接上代码吧,原理一大堆,自行了解即可!这是LMS即自适应滤波器的Matlab代码。
function [yn,W,en]=LMS(xt,dn,M,mu)
itr = length(xt);
en = zeros(itr,1); % 误差序列,en(k)表示第k次迭代时预期输出与实际输入的误差
W = zeros(M,itr);% 每一行代表一个加权参量,每一列代表-次迭代,初始为0
% 迭代计算
for k = M:itr % 第k次迭代
x = xt(k:-1:k-M+1) ; % 滤波器M个抽头的输入
y = W(:,k-1).' * x ; % 滤波器的输出
en(k