[slgu]

原则上a,b应一起估计，方法也很简单，即使用 lsqcurvefit or nlinfit 等命令。
b=lsqcurvefit(fun,b,x,y)
b=nlinfit(x,y,fun,b)
fun为拟合函数，如例其函数可表示为：
fun=@(b,x)b(1)*x.^b(2);
即可。
圆括号中的b为初值，一般情况下，象这种简单的幂函数方程，b=rand(1,2)随机初值即可。