x=[70, 70, 70, 80, 80, 80, 90, 90, 90, 100, 100, 100, 200, 200, 200, 300, 300, 300]';
y=[25, 50, 100, 25, 50, 100, 25, 50, 100, 25, 50, 100, 25, 50, 100, 25, 50, 100]';
z=[3.9, 8.5, 25, 5.5, 12, 34, 7.5, 18, 46, 10, 22, 60, 65, 150, 400, 200, 450, 1200]';
X=[x y];n=length(z);
fX=@(b,X)b(1)*X(:,1).^b(2).*X(:,2).^b(3);
b=[0.000000429798, 2.7107579, 1.3650115]
for l=1:5
    b=lsqcurvefit(fX,b,X,z);
    b=nlinfit(X,z,fX,b);
end
vpa(b,10)
SSz=var(z)*(n-1)
zhat=fX(b,X);
RSS=(z-yhat)'*(z-zhat)
Rsquare=(SSz-RSS)/SSz

我是一个MATLAB初学者，MATLAB很多功能都不了解，对于这个问题，我说一下我的想法。
可以方程两边先取对数，将其化为线性方程，在用解线性方程组的方法，将a，b，c看做未知数，用MATLAB解出最小二乘解，也许可以解出答案。