linear_function.m

来自「最小二乘法(Least-square method)在matlab平台下的实现程」· M 代码 · 共 21 行

clear all, clc
% A=[1,2,3;2,14,4;3,4,5];  % A对称正定
% x=[2,3,4.5]'
% b=A*x;
A=[4,-1,1;-1,4.25,2.75;1,2.75,3.5];
x=[2,-1,-1]'
b=[6,-0.5,1.25]';
% 最速下降法------------------------%
x0=[1,-1,-1]';
r0=b-A*x0;
k=0;
while sum(abs(r0))>1e-8
    k=k+1;
    alfa=r0'*r0/(r0'*A*r0);
    x0=x0+alfa*r0;
    r0=b-A*x0;
    if k> 1000
        break;
    end
end
x0