exm041_2.m

来自「北京航空航天大学出版社张志涌、徐彦琴主编的matlab教程附带的mfile例程压」· M 代码 · 共 19 行

%exm041_2.m

%（1）生成矩阵 及 ，并用三种方法求解
A=gallery(5)				
A(:,1)=[];
y=[1.7 7.5 6.3 0.83 -0.082]';
x=inv(A'*A)*A'*y,
xx=pinv(A)*y,
xxx=A\y

%（2）计算三个解的范数
nx=norm(x),
nxx=norm(xx),
nxxx=norm(xxx)

%（3）比较三种解法的方程误差
e=norm(y-A*x),
ee=norm(y-A*xx),
eee=norm(y-A*xxx)