exm04022_1.m

来自「这本书是matlab经典书籍」· M 代码 · 共 15 行

randn('state',0);
A=gallery('randsvd',100,2e13,2);	%产生条件数为2e13的100阶随机矩阵
x=ones(100,1);				%指定真解
b=A*x;							%为使Ax=b方程一致，用A和 x 生成 b 向量。
cond(A)

tic								%启动计时器Stopwatch Timer
xi=inv(A)*b;					% xi 是用“求逆”法解恰定方程所得的解。
ti=toc							%关闭计时器，并显示解方程所用的时间。
eri=norm(x-xi)				%解向量 xi 与真解向量 x 的范-2误差。
rei=norm(A*xi-b)/norm(b)	%方程的范-2相对残差  

tic;xd=A\b;			%启动计时器，“左除”求方程解。
td=toc,erd=norm(x-xd),red=norm(A*xd-b)/norm(b)