递推最小二乘法

主要用于求解超定方程的未知解，递推最小二乘法将用于实现机器人动力学参数辨识等方面，相关程序如下：

相关代码说明

使用递推最小二乘法，求解Ax=b的解

A为m*x维的矩阵，元素服从独立同分布的正态分布

b为m维的已知向量，元素也是服从独立同分布的正态分布

相关代码引用（第五清风创作）

clear;
format long;
A = randn([10000 10]);
b = randn([10000 1]);
x = rand(10,1);
I = eye(10, 10);
P = (10^6) * I;

for k = 1:10000
    Ak = A(k,:);
    Q1 = P*(Ak');
    Q2 = 1 + Ak * P * (Ak');
    Q = Q1/Q2;
    x = x + Q * (b(k) - Ak*x);
    P = (I - Q*Ak)*P;
    result2(:,k) = x;
    result1(k) = k;
end
result1 = result1';
%result = [result1; result2];
plot(result1, result2);