exm0522_1.m

matlab程序

%exm0522_1.m

clear
rand('seed',12);	%选定随机因子, 使得可重复产生随机矩阵A
A=rand(10)     	%生成100X100均匀分布的随机矩阵

A=rand(10)+1.e8;	%生成100X100均匀分布的随机矩阵
x=ones(10,1);		%令解向量x为全1的100元列向量
b=A*x;			    %为使Ax=b方程一致, 用A和x生成b向量

%（2）"求逆"法解恰定方程的误差、残差、运算次数和所用时间
tic				            %浮点运算计数器置0 ；启动计时器Stopwatch Timer
xi=inv(A)*b;				% xi 是用"求逆"法解恰定方程所得的解。
ti=toc						%关闭计时器，并显示解方程所用的时间。

eri=norm(x-xi)				%解向量 xi 与真解向量 x 的范-2误差。
rei=norm(A*xi-b)/norm(b)	%方程的范-2相对残差  

%（3）"左除"法解恰定方程的误差、残差、运算次数和所用时间
tic;
xd=A\b;		%是用"左除"法解恰定方程所得的解。
td=toc,
erd=norm(x-xd),
red=norm(A*xd-b)/norm(b)