exm040422_1.m

来自「MATLAB 6.5 从基础讲解其功能」· M 代码 · 共 6 行

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  

nx=norm(x),nxx=norm(xx),nxxx=norm(xxx)  

e=norm(y-A*x),ee=norm(y-A*xx),eee=norm(y-A*xxx)