plot.m

多项式二次曲线拟合程序！用于实验设计！另有画图函数并且可以进行方差分析1

a=-1:0.05:1;
b=-1:0.05:1;
[X,Y]=meshgrid(a,b);
[m,n]=size(X)

for i=1:m
    for j=1:m
        XY=[X(i,j),Y(i,j),0];
        Z(i,j)=A+XY*B+XY*C*XY';
    end
end
 
size(Z)
figure(1)
meshc(X,Y,Z)
figure(2)
[CC,h]=contour(Z);
set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2)
colormap cool
clear Z X Y