m0725.m

来自「MATLAB仿真技术与实例应用教程」· M 代码 · 共 13 行

x  = 0.1*randn(1,500); % 生成机信号.
b  = fir1(32,0.55);     % 待鉴别的滤波器
d  = filter(b,1,x);    % 自适应滤波器的输入
w0 = zeros(1,33);      % 初始化滤波器的系数
p0 = 5*eye(33);        % 初始化相关矩阵
lambda = 1.0;         % 权重因子
s = initrls(w0,p0,lambda);
[y,e,s] = adaptrls(x,d,s);
stem([b.',s.coeffs.']);
legend('实际值','估计值');
title('RLS自适应滤波'); grid on;