m0727.m

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

x  = 0.1*randn(1,500); % 生成信号
b  = fir1(31,0.5);     % 待检验的滤波器
d  = filter(b,1,x);    % 预测信号
w0 = zeros(1,32);      % 初始化滤波器的系数
k0 = 0.5*eye(32);      % 初始化状态误差矩阵
qm = 2;                % 测量噪声协方差
qp = 0.1*eye(32);      % 处理噪声协方差  
s = initkalman(w0,k0,qm,qp);
[y,e,s] = adaptkalman(x,d,s);
stem([b.',s.coeffs.']);
legend('实际值','估计值');
title('Kalman自适应滤波'); 
grid on;