ea520a.m

来自「用MATLAB程序解决高阶问题；通过20多个实例来说明线性代数在后续课和工程中的」· M 代码 · 共 24 行

% 《工程线性代数（MATLAB版）》第五章例题5.20程序ea520a
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 二次型方程拟合问题
%
clear
x=[-1;0;2;2;0];y=[0;2;3;-2;-3]; 
A=[x.^2,x.*y,y.^2,x,y],
b=-ones(5,1)
K=A\b
ezplot('2*X.^2-X.*Y+Y.^2-4*X+Y-6*ones(size(X))=0'),hold on
plot(x,y,'sk')
axis([-4 4 -4 4]),grid on
set(gcf,'color','w'),pause
disp('增加一个点的超定情况'),pause
x1=[x;1];y1=[y;2]; 
A1=[x1.^2,x1.*y1,y1.^2,x1,y1],
b=-ones(6,1)
K1=pinv(A1)*b,pause
disp('将所得的K1各项代入椭圆方程，再用ezplot画图'),pause
echo on
ezplot('32/117*X.^2-89/468*X.*Y+43/234*Y.^2-265/468*X+107/468*Y-ones(size(X))=0')
pause, echo off
plot(x1,y1,'*r')