递归最小二乘matlab

以下是递归最小二乘的 MATLAB 代码：

function [theta, P] = rls(x, y, lambda, P)
% x: 输入数据
% y: 输出数据
% lambda: 遗忘因子
% P: 初始P矩阵

N = length(x);
theta = zeros(size(x, 2), 1);

for i = 1:N
    k = P * x(i, :)' / (lambda + x(i, :) * P * x(i, :)');
    theta = theta + k * (y(i) - x(i, :) * theta);
    P = (P - k * x(i, :) * P) / lambda;
end
end