采用高斯拟合时显示精度溢出或得到结果显然存在明显的较大偏差

Question

数据如下

data1 := {{6.2, 575.8}, {6.3, 1025.0}, {6.4, 1168.3}, {6.5,840.3}, {6.6, 485.5}}

采用FindFit，FindFit[data1, a*Exp[-(x - b)^2/c^2], {a, b, c}, x]，然后显示

LinearAlgebra`BLAS`TRSV::oflow: 在计算过程中，碰到机器精度溢出.

求问这是为什么

上面的问题已被解答，然而又有一个相似的问题。。。

data2 = {{7.0, 647.0}, {7.1, 1909.3}, {7.2, 2352.8}, {7.3, 
    1682.0}, {7.4, 633.0}};
FindFit[data2, a*Exp[-(x - b)^2/c^2], {a, b, c}, x, 
 MaxIterations -> 10000]

 

所得的结果明显有问题

NonlinearModelFit[data2, a*Exp[-(x - b)^2/c^2], {a, b, c}, x, 
 MaxIterations -> Infinity]

分析之后得到

{
 {"", "Estimate", "Standard Error", "t\[Hyphen]Statistic", 
  "P\[Hyphen]Value"},
 {a, 4182.48, 1.67714*10^24, 2.49381*10^-21, 1.},
 {b, 352.755, 4.2033*10^21, 8.39233*10^-20, 1.},
 {c, 60.1843, 3.66042*10^20, 1.64419*10^-19, 1.}
}

显然存在问题

Answer 1

data1 = {{6.2, 575.8}, {6.3, 1025.0}, {6.4, 1168.3}, {6.5, 840.3}, {6.6, 485.5}};
FindFit[data1, a*Exp[-(x - b)^2/c^2], {a, b, c}, x, MaxIterations -> 1000]
Plot[a*Exp[-(x - b)^2/c^2] /. %, {x, 6.2, 6.6}, Epilog -> Point[data1]]

V10.4  Win7

没出问题。

Comments

谢谢，原来提高迭代次数就可以了。。

Answer 2

已解决。，大致估算出b的值时即可得到较好的结果