matlab解幂函数或对数方程组的参数,方程如下：y=a*x^(b)(x1,y1)=(1400,25),(x2,y2)=(5000,10)

Question

matlab解幂函数或对数方程组的参数,方程如下：
y=a*x^(b)
(x1,y1)=(1400,25),(x2,y2)=(5000,10)

Answer 1

>> x=[1400 5000];
>> y=[25 10];
>> fun=inline('a(1).*x.^a(2)',‘a','x');
>> a=lsqcurvefit(fun,[1;1],x,y);
>> vpa(a)
ans =
4597.7756792206319005344994366169 a值
-0.71980788585361032971832173643634 b值
a=4597.77568
b= -0.71981