main.m

anoldi算法,有运行结果,效果还行,可自己改进,能对大型矩阵操作

tic
clear;clc;
n=2000;
a=hilb(n);
b1=a*ones(n,1);
aph=1e-1;
A=a'*a+aph*eye(n);
b=a'*b1;
m=25;
x=zeros(n,1);
[x, error] = gm( A, x, b, m)
t1=toc