exm05422_1.m

来自「书名应为《掌握和精通MATLAB》附有完整的源程序」· M 代码 · 共 21 行

%exm05422_1.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)