>> z = [9.73875 20.75 36.59875

13.5725 29.6325 50.93875

18.97875 36.59875 80.13875

20.75125 38.22125 90.925

22.055 44.58 104.7725];

>> [p,t] = meshgrid([0.8 1 1.2],[60:60:300]);

>> stats = reglm(z(:),[p(:), t(:)],'quadratic',{'p','t'});

>> [pnew,tnew] = meshgrid(linspace(0.8,1.2,20),linspace(60,300,20));

>> pp = pnew(:);

>> tt = tnew(:);

>> zhat = [ones(400,1) pp tt pp.*tt pp.^2 tt.^2]*stats.beta;

>> mesh(pnew,tnew,reshape(zhat,[20,20]));

>> hold on

>> plot3(p(:),t(:),z(:),'k*')

拟合结果：

------------------------------------方差分析表------------------------------------

方差来源    自由度            平方和             均方             F值          p值

回归       5.0000       11548.9147        2309.7829         93.4739      0.0000

残差       9.0000         222.3942          24.7105

总计      14.0000       11771.3089

均方根误差(Root MSE)         4.9710           判定系数(R-Square)    0.9811

因变量均值(Dependent Mean)        41.2168        调整的判定系数(Adj R-Sq)    0.9706

-----------------------------------参数估计-----------------------------------

变量               估计值            标准误             t值          p值

常数项            242.6188          69.0439           3.5140      0.0066

p           -513.7781         137.3777          -3.7399      0.0046

t             -0.3637           0.1212          -3.0002      0.0150

p*t              0.6022           0.0926           6.5010      0.0001

p*p            272.2625          68.0677           3.9999      0.0031

t*t             -0.0003           0.0002          -1.1946      0.2628