compare.m

MATLAB科学计算与工程分析源代码源程序4

% 该M文件用来演示求逆法与除法求解线性方程组在时间与精度上的区别
% 编写日期：2007-5-14

A=1000*rand(1000,1000);   %随机生成一个1000维的系数矩阵
x=ones(1000,1);
b=A*x;
disp('利用矩阵的逆求解所用时间及误差为：');
tic
y=inv(A)*b;
t1=toc
error1=norm(y-x)     %利用2-范数来刻画结果与精确解的误差

disp('利用除法求解所用时间及误差为：')
tic
y=A\b;
t2=toc
error2=norm(y-x)