0

メスバウアー分光データのカスタムフィットにcftoolを使用しています。ガンマとN0の2つの係数があります。

N = f(v)
  = (299792458000^2*Gamma^2*N0)/(299792458000^2*Gamma^2+4*v^2*(4.29383292e-15)^2)

デフォルト設定（信頼領域、ロバストオフなど）を使用すると、次のようになります。

Fit computation did not converge:
Fitting stopped because the number of iterations or function evaluations exceeded the specified maximum.

Fit found when optimization terminated:

General model:
     f(v) = (299792458000^2*Gamma^2*N0)/(299792458000^2*Gamma^2+4*v^2*(4.29383292e-
                    15)^2)
Coefficients (with 95% confidence bounds):
       Gamma =      0.9137  (-Inf, Inf)
       N0 =   2.454e+04  (2.059e+04, 2.849e+04)

Goodness of fit:
  SSE: 6.41e+11
  R-square: -2068
  Adjusted R-square: -2073
  RMSE: 4.013e+04

Warning: A negative R-square is possible if the model does not contain a constant term and the fit is poor (worse than just fitting the mean). Try changing the model or using a different StartPoint.

Levenberg-Marquardtに切り替えると、データを直線で表すことができます。

General model:
     f(v) = (299792458000^2*Gamma^2*N0)/(299792458000^2*Gamma^2+4*v^2*(4.29383292e-
                    15)^2)
Coefficients (with 95% confidence bounds):
       Gamma =       0.793  (-Inf, Inf)
       N0 =   6.456e+04  (6.447e+04, 6.465e+04)

Goodness of fit:
  SSE: 3.098e+08
  R-square: 2.22e-16
  Adjusted R-square: -0.002513
  RMSE: 882.3

ここに画像の説明を入力してください

どちらの場合も、なぜこれがそれほどひどく失敗するのですか？

4

1 に答える 1

1

f（v）f(v)=N0/(1+(2.8645e-26*(v/Gamma))^2)は、（v / Gamma）が10 ^ 25まで大きくなり始めるまで、分母の1が支配的になるように単純化されます。ガンマが0.793、vが|15|の場合 matlabはN0以外に収束するのに苦労するかもしれないと思います

于 2013-01-20T19:01:53.273 に答える