我需要使用polyfit命令来确定以下二阶两项方程的最佳指数系数:
y= a1 * exp(b1*x1^2+b2*x2^2+b3*x1*x2+b4*x1+b5*x2)
对于最简单的情况(y= a* exp(kx)
),我会这样做:
yeqn = @(coefs, x) coefs(1).*exp(coefs(2).*x);
xi = x;
eta = log(y);
P = polyfit(xi, eta, 1);
MyCoefs(1) = exp(P(2));
MyCoefs(2) = P(1)
有没有机会为我的情况做同样的事情,或者我应该使用另一种方法?