a=27.971;
b=-12.167;
c=48.01;
d=-48.904;
p=25.986;
x=[0:0.1:1.0];
y=[0:0.1:0.6];
[X,Y]=meshgrid(x,y);
Z=a.*X.^2+b.*X+c.*Y.^2+d.*Y+p;
mesh(X,Y,Z)
surf(X,Y,Z)
hold on
a=29.126;
b=-10.517;
c=49.002;
d=-50.541;
p=25.867;
x=[0:0.1:1.0];
y=[0:0.1:0.6];
[X,Y]=meshgrid(x,y);
Z=a.*X.^2+b.*X+c.*Y.^2+d.*Y+p;
mesh(X,Y,Z)
surf(X,Y,Z)
hold on
a=28.506;
b=-9.44;
c=52.392;
d=-53.711;
p=27.044;
x=[0:0.1:1.0];
y=[0:0.1:0.6];
[X,Y]=meshgrid(x,y);
Z=a.*X.^2+b.*X+c.*Y.^2+d.*Y+p;
mesh(X,Y,Z)
surf(X,Y,Z)
hold on
a=28.741;
b=-10.01;
c=50.792;
d=-51.924;
p=26.204;
x=[0:0.1:1.0];
y=[0:0.1:0.6];
[X,Y]=meshgrid(x,y);
Z=a.*X.^2+b.*X+c.*Y.^2+d.*Y+p;
mesh(X,Y,Z)
surf(X,Y,Z)
hold on
a=31.441;
b=-10.304;
c=53.486;
d=-54.264;
p=27.847;
x=[0:0.1:1.0];
y=[0:0.1:0.6];
[X,Y]=meshgrid(x,y);
Z=a.*X.^2+b.*X+c.*Y.^2+d.*Y+p;
mesh(X,Y,Z)
surf(X,Y,Z)
title('由下到上：主-主、主-次、主-支、次-次、次-支'）
xlabel('V/C')
ylabel('g/c')
zlabel('FC1')
axis([0 1.0 0 0.6 10 42])