最小二乗フィッティングは他のグラフツールよりも使いやすいと思います。フィッティングする関数を直接数式で書くことができるので、とても直感的です。下のようなデータをフィッティングする例を考えます。フィッティングするgnuplotのコマンドはfitです。a,bのように書きます。f(x)はフィッティングする関数です。2のように指定します。指定しないと、デフォルト値の1からフィッティングを開始します。非線形の最小二乗フィッティングは一般に初期値に依存しますから、初期値を指定するほうが良いでしょう。fitの後ろの[][]はそれぞれxとyの範囲を指定します。データのどの列を使うか指定します。a,bはフィッティングするパラメタを指定します。'hoge.log'はフィッティング結果を書き込むログファイルを指定します。d$$でフィッティングしています。outputフィッティング結果フィッティングした結果を下に示します。フィッティング後にフィッティング関数をプロットするとフィッティングした結果(フィッティング・パラメタの値)を参照してプロットしてくれます。フィッティングする範囲を指定するには(xでもyでも)xの範囲を[-1:1]に限定するにはfitの後ろに[-1:1]をつけます。b$$の形でフィッティングしています。outputデータがある条件を満足する時だけフィッティングしたいデータがある条件を満足する時だけフィッティングしたい時はusingの部分に条件(三項演算子)をいれることで実現できます。上の例だとデータ1列目($1)が正の時に2列目($2)をy軸とします。1/0とすると無視されます。$$でフィッティングしています。フィッティングを制御する変数がいくつかあります。普通はあまり気にしなくて良いでしょう。FIT_LIMIT収束条件を指定します。デフォルト値は1e-05です。一回前と今の二乗残差の比がこの値より小さいと収束と見なされます。FIT_MAXITER収束計算の繰り返し回数最大値を指定します。0だと制限なしで繰り返します。法で使うラムダ値の開始値を指定します。この値を指定しない場合のデフォルト値は自動計算されます。FIT_LAMBDA_FACTORMarquardt-Levenberg法で使うラムダ値の倍率を指定します。デフォルト値は10.0になっています。フィッティング・ログの例フィッティングした際のログの例です。収束の様子や、フィッティング・パラメタの値が記録されています。matrix”というものが表示されます。対角成分はいつも1です。非対角要素は、そのフィッティング・パラメタ同士が独立であれば、ゼロに近くなります。残差ながら、Excelなどで見るピアソンの相関係数とは同じものではないです。コメントをキャンセルメールアドレスが公開されることはありません。date*日本語が含まれない投稿は無視されますのでご注意ください。コメントデータの処理方法の詳細はこちらをご覧ください。使うRECOMMENDこちらの記事も人気です。