matlab绘制二元二次曲线图,急求拟合二元二次方程，和作MATLAB三维图 - 数学 - 小木虫 - 学术 科研 互动社区...

function eryuanihe

A=[87.5                3.5                2389.433                2600.933                2921.633

162.5        3.5                2856.667                3351.3                3240.867

162.5        6.5                4617.7                4261.767                4013.833

87.5                6.5                2679.3                2707.533                2926.867

125                2                1617                1417.533                1369.6

125                8                4293.867                4343.6                4837.967

50                5                1971.467                2241.433                2196.967

200                5                2661.133                2578.9                2457.667

125                5                3782.667                4267.667                3247.633

125                5                3818.2                3022.1                3221.533

125                5                3366.533                4261.733                3618.533

125                5                4570.233                3907.933                3291.5

];

subplot(1,3,1)

nihe(A(:,1),A(:,2),A(:,3),20,20)

xlabel('x1');ylabel('x2');zlabel('y1');

hold on

subplot(1,3,2)

xlabel('x1');ylabel('x2');zlabel('y2');

nihe(A(:,1),A(:,2),A(:,4),20,20)

subplot(1,3,3)

nihe(A(:,1),A(:,2),A(:,5),20,20)

xlabel('x1');ylabel('x2');zlabel('y3');

function nihe(x1,x2,y,M,N)

B(:,1)=x1.^2;

B(:,2)=x2.^2;

B(:,3)=x1.*x2;

B(:,4)=x1;

B(:,5)=x2;

B(:,6)=1;

B(:,7)=y;

C=B(:,[1:6]);

xs=inv(C'*C)*(C'*y);

x1t=linspace(min(x1),max(x1),M);

x2t=linspace(min(x2),max(x2),N);

[x1t,x2t]=meshgrid(x1t,x2t);

yt=xs(1)*x1t.^2+xs(2)*x2t.^2+xs(3)*x1t.*x2t+xs(4)*x1t+xs(5)*x2t+xs(6);

mesh(x1t,x2t,yt);，