ea520.m

matlab 在线性代数中的应用 利用matlab实现可以提高对所学知识的理解 更有助于你考研

% 《工程线性代数（MATLAB版）》第五章例题5.20程序ea520
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 二次型方程插值问题
%
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
pause
[X,Y]=meshgrid(-4:0.1:4);
Z=2*X.^2-X.*Y+Y.^2-4*X+Y-6*ones(size(X));
contour(X,Y,Z,[0 0]),hold on
plot(x,y,'s')