• 最小二乗法によるフィッティング

    結構ニーズがあるようでして、どうやって使うのか?といったお便りも少なくない。ネットで検索すれば山ほど情報が得られるだろうがその多くは直線の当てはめや特定の関数の当てはめに特化している。しかしいつでもそういうパタンにはならないのが悲しいところである。全くである。自分もそう感じていたりする。そういった悲しさからCCTLには最小二乗法で任意の(非線形の)関数の当てはめができる関数がある。nlsqfit がそれである。計算もそこそこ早い。しかしもっと強力な関数がfindminNdfindminGANdであるがこれらは最適化関数なのだが目的関数が異なるだけで最小二乗法にも使えるのです。ただ扱いがちょっと面倒な点があるのでうまく使えない、といった感じのお便りもありもう少し便利に使えるようなパッケージを作っている。
    例えば次のような使い方が可能になる(予定である)。
    #%
    # 非線形式に対応した最小二乗フィッティングパッケージ
    %data "Optimiz.pkg";

    # 非線形式に対応した最小二乗フィッティングを使う
    call Optimiz:Fitting(n, P, m, param, expr)->(xmin, xmax, param);


    実行例(まだパッケージは未完なのでちまちましたスクリプトを書いて実行した例です)

    グラフィックが貧弱なのでgnuplotで表示しています。

2006年05月04日 20時30分54秒

inserted by FC2 system