粗大メモ置き場

個人用,たまーに来訪者を意識する雑記メモ

お手軽関数フィッティングチャレンジ

過去記事でフィッティングチャレンジしたんですけど,
今一度復習してみようという内容です。

過去記事との変更点

ossyaritoori.hatenablog.com
ossyaritoori.hatenablog.com

  • Levenberg-Marquardt法を用いた二変数関数,6変数のフィッティング。
  • jacob関数を用いてヤコビ行列計算を簡略化。
  • 関数の値を計算し,ノイズを載せることで観測データをシミュレートした。
  • ezsurfを用いた図示(fsurf)の方が良いらしい。

プログラム

解説(略)


結果

付与したガウスノイズの標準偏差が1のとき。
f:id:ossyaritoori:20170408005014p:plain

付与したガウスノイズの標準偏差が5のとき。
f:id:ossyaritoori:20170408005108p:plain

初期条件が良かったのもあるが高々~5%くらいの誤差なら十分フィッテイング出来そう。