gnuplotで・・・

質問

0.0000 5.501251e+02
0.5000 2.937885e+02
1.0000 3.253154e+02
1.5000 3.603518e+02
2.0000 1.771118e+02
2.5000 2.438578e+02
3.0000 1.212656e+02
3.5000 -7.970373e+01
4.0000 -3.334714e+01
4.5000 -1.631910e+02
5.0000 -6.219916e+01
5.5000 -2.850575e+02
・・・・
・・・・
を s(t) の測定データ（ t の単位は秒、s の単位は Pa ）とみなし、グラフを書き、
σ(t)=σ0cos(ωｔ+δ)で fit し、σ0、δを求めるのですが、ωのパラメーターの初期値をどのくらいにすればいいのかわかりません。こういう場合、viaの後に3つパラメータ初期値の数字をおくのでしょうか。
わかる方詳しく教えてください。fitting自体がありえるのかよくわからないところがあるので、のコツみたいのがあったらお願いします。

回答No.1

0123456789A · ベストアンサー

gnuplotのfittingは最小二乗法を使っています。
fitting自体はありえると思いますが、非線形なものではあまり信用できないかもしれません。
最小二乗fitでの正当性は考察しなければなりません。
さて、ご質問の件ですが、初期値を入れる必要はありません。
gnuplotが適当に値を入れてiterationすることで最適な
パラメータの値σ0,ω,δと更に最終的に得られたもののエラーを返してくれるはずです。

念のためやり方を書いておきます。お持ちのデータファイル名をhoge.datとします。

>f(x)=a*cos(b*x+c)
>fit f(x) "hoge.dat" via a,b,c

これでfit完了.同じディレクトリにfit.logというファイルが出来ているので
それを見ればどれだけの精度でfittingできたのかもわかります。

お礼

ありがとうございました。