m0726.m

来自「matlab的经典讲义」· M 代码 · 共 11 行

x  = 0.1*randn(1,500); % 生成信号
b  = fir1(31,0.5);     % 待鉴别的滤波器
d  = filter(b,1,x);    % 预测信号
w0 = zeros(1,32);      % 初始化滤波器的系数
mu = 0.8;              % LMS步长
s = initlms(w0,mu);
[y,e,s] = adaptlms(x,d,s);
stem([b.',s.coeffs.']);
legend('实际值','估计值');
title('LMS自适应滤波'); grid on;