%要解决的问题：拟合Z=a*X^b*Y^c，其中a,b,c为要确定的参数。

%思路：将上式转化为ln(Z)=ln(a)+b*ln(X)+c*ln(Y)，这样就变成了简单的多元回归。

Data=xlsread('原始数据.xlsx');

x=Data(:,1);

y=Data(:,2);

z=Data(:,3);

X=log(x);Y=log(y);Z=log(z);

m=size(x);n=max(m);

XX=[ones(b,1) X Y];

[p,bint,r,rint,s]=regress(Z,XX);

a=p(1);b=p(2);c=p(3);

[x1,x2]=meshgrid(x,y);

z1=a*(x1.^b).*x2.^c;

mesh(x1,x2,z1);