main.m

来自「anoldi算法,有运行结果,效果还行,可自己改进,能对大型矩阵操作」· M 代码 · 共 12 行

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