>> 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