現在、ガウシアン関数y=a+b*exp(-(x-c)^2/d^2)に下記のようなデータを使用しフィッティングを行いたいのですが、
手法やパラメータa,b,c,dの求め方がわかりません。
どなたか教えていただけませんか。
よろしくお願いいたします。
(x,y)={
48.8006092
48.8056105
48.8105942
48.8156000
48.8206021
48.8256127
48.8306131
48.8356169
48.8406146
48.8456077
48.8506141
48.8556236
48.8606115
48.8656179
48.8706296
48.8756176
48.8806272
48.8856294
.....}
No.1ベストアンサー
- 回答日時:
Excelしか使えないのであれば、ソルバーを使って、以下の手順で「残差2乗和」を最小とするパラメータ a, b, c, d を探すと良いでしょう。
ただし、パラメータの初期値があまりかけ離れていると変な値に収束するか解が見つかりません。a, b, c は元のデータのグラフから見当をつけられると思います( a はベースライン高さ、b はピーク高さ、 c はピーク位置x)。d は直感では見当をつけられないので、(ピークの半値全幅)/1.67 で計算してそれを初期値としてください。【ソルバーを使った最小自乗法】
(1) Excelのメニューの [ツール] → [アドイン] で [ソルバーアドイン] の左側の□をチェックして OK
(2) A列に x データ、B列に y データを書き込む( x は A1 から、y は B1 から下方向に書き込む)
(3) D1からD4にフィッティングパラメータの初期値を書き込む(a → D1、b → D2、c → D3、d → D4)
(4) セル C1 に以下の式を貼り付ける(これをCopy&Paste)
=($D$1+$D$2*EXP(-1*(A1-$D$3)^2/($D$4)^2)-B1)^2
(5) セルC1をコピーして、C2以下の全データ分のC列にペースト(これでC列=「残差2乗」となる)
(6) セルD5 に =sum(C1:C??)と書く(??はC列最後の行番号)
(7) Excelのメニューの [ツール] → [ソルバー] で [目的セル」を $D$5、[目標値] を 最小、[変化させるセル]を $D$1:$D$4 とする。この意味は、「セルD1~D4に書かれている数値を変化させて、 D5セルを最小となるようにする」ということです。
(8) ソルバーウィンドウのオプションボタンをクリック → 制限時間を 1000、反復回数を 1000、精度・公差・収束をすべて 1e-10 とし、OKをクリック
(9) ソルバーウィンドウの実行ボタンをクリック
(10) フィッティングパラメータが見つかったら、「最適解が見つかりました・・」と出るので、OKをクリック
(11) ExcelのセルD1~D4にフィッティングパラメータが書き込まれている
なお、精度・公差・収束の値をあまり小さくすると収束しないので、もしデータのばらつきが大きくて収束しないときは、これらの値を適宜、大きくしていってみてください。ご質問のデータはGaussianのほんの1部でしたのでこちらで実験することはできませんでした。ちなみに私は通常、カレイダグラフというグラフソフトで任意関数のフィッティングをやっています。
この回答への補足
inara様 御回答ありがとうございます。私の質問内容が少し抜けておりました。Excel、gnuplot等のソフトを使用せず、計算結果が出力されるまでの途中のアルゴリズムを教えていただけませんかと記載すればよかったのですが質問が不十分でした。大変御迷惑をお掛けし、申し訳ありませんでした。
補足日時:2007/11/15 09:17お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
質問したい事が2つあります。 ①...
-
これなぜせんぶんAB上だったり...
-
√2が無理数であることの証明で...
-
以下の問題が示している領域が...
-
数学の式変形を教えてください!
-
Aさんがいたとします。Aさんの...
-
数学 算数の通分について 分数...
-
数学の問題を始めると何故か涙...
-
締め上げ暗号解析に耐性のある...
-
数学が大好きな国の国旗のデザイン
-
水溶液の希釈。なぜこれで解け...
-
中3 連立方程式 代金について...
-
これの(2)でan>0を確認してい...
-
√2の証明によく出てくる言葉で...
-
隣り合う平方数の大きい数から...
-
これなぜ最後の不定形が0に収束...
-
なぜこのように置換しようと思...
-
これの(3)はどういった発想で解...
-
図の問題の|z|=1の写像を計算し...
-
このシグマのしきはなにを表し...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
小学生の時(40年前)に、18÷...
-
この問題の解答が(写真) A. a <...
-
私がばかなのか?
-
素数(合成数の並びの最大数)...
-
楕円と回転行列について
-
背理法について
-
数学が本当にできないので相談...
-
数学の問題ですが、わかりません
-
ぜっったいちがくないですか?...
-
一般相対論でのローレンツ条件...
-
cos^2θ/tanθ=1でθを出すことは...
-
a^3+b^3=(a+b)(a^2-ab+b^2)とな...
-
当選確率5%(5/100)のくじを32...
-
つかれた
-
a^2+b^2+c^2-ab-bc-ca を因数分...
-
a≥1とする。曲線y=x(x-a)(x-3a)...
-
数学わからないので、ムカつい...
-
河野玄斗の数学力をどう思いま...
-
ライプニッツ記法
-
有界はどうつかいますか?
おすすめ情報