3

我正在研究由两个数组组成的曲线拟合数据:

t: 1, 3, 4, 7, 8, 10

P: 2.1, 4.6, 5.4, 6.1, 6.4, 6.6

两个变量之间的关系由 给出P = mt/(b+t)。我被告知通过将方程曲线拟合到数据点来确定常数 m 和 b。这应该通过写出方程的倒数并使用一阶多项式来完成。这是我的代码:

t = [1 3 4 7 8 10];
P = [2.1 4.6 5.4 6.1 6.4 6.6];

p = polyfit(t, t./P, 1);


m = 1/p(1)
b = p(2)*m

tm = 1:0.01:10;
Pm = (m*tm)./(b+tm);

plot(t,P, 'o', tm, Pm)

书中的答案是m = 9.4157b = 3.4418。上面的代码产生m = 8.4807b = 2.6723。我的错误是什么?任何建议将不胜感激。感谢您的时间。

4

1 回答 1

1

要跟进@David_G 的评论,您似乎有更好的答案。事实上,如果你通过 MATLAB 中的曲线拟合工具箱运行数据,你会得到:

General model:
  f(t) = m*t/(b+t)
Coefficients (with 95% confidence bounds):
   b =       2.587  (1.645, 3.528)
   m =       8.448  (7.453, 9.443)

Goodness of fit:
  SSE: 0.1594
  R-square: 0.9888
  Adjusted R-square: 0.986
  RMSE: 0.1996

您的解决方案几乎一样好:

Goodness of fit:
  SSE: 0.1685
  R-square: 0.9881
  Adjusted R-square: 0.9852
  RMSE: 0.2053

而且他们两个都比书中的那个好:

Goodness of fit:
  SSE: 0.404
  R-square: 0.9716
  Adjusted R-square: 0.9645
  RMSE: 0.3178
于 2013-05-24T10:40:56.080 回答