MATLAB：功率函数曲线拟合

Question

我想知道幂多项式 y = x ^m 的MATLAB代码形式。<登记/> 通过将该曲线拟合到一组数据，我需要 m 的数值。我的初步计划如下：

load A   
N = A(:, 1);
t = A(:, 2);    
p = polyfit(log(N), log(t), 1);    
f = (p, t);    
plot(t, N, '-', t, f, '-');

我唯一怀疑的是用于曲线拟合的上述polyfit函数用于等式 y = bx ^m 这不符合我的数据集，我需要 b = 1。

Answer 1

polyfit的第三个参数是多项式次数， m 在你的情况下。

您可以将大于 m 的参数传递给 polyfit 。然后，如果你对模型的假设是正确的 p [m + 1] 将是 p （并接近1）的最大组成部分，而其他的应该接近于零