• 締切済み

gnuplotでの最小二乗法について。

他のカテで質問したのですが、 こちらの方が関係性が深いと思い移動しました。 y,x1,x2 を測定データとして、線形関数 y=a*x1-b*x2 におけるパラメータ a,b をgnuplotを使って最小二乗法で求めたいのですが、参列のデータ(x1,x2,y)を用意して、 f(x)=a*x1-b*x2 fit f(x) "data.dat" via a,b とうった時点で、undefined variable: x1 とエラーが出てしまいます。 どうすればよいでしょうか。よろしくお願いします。

みんなの回答

  • Sbacteria
  • ベストアンサー率42% (55/129)
回答No.1

ぱっと、見ただけですが、2変数関数なのですね? それじゃ、f(x)の定義では、だめなのでは無いですか?

karin19
質問者

お礼

お礼が遅れてしまって申し訳ございません。 ありがとうございました。

関連するQ&A

  • gnuplotで原点を通る最小二乗法

    gnuplotで原点を通る最小二乗法を行う方法が知りたいです. 現在は以下のようなコマンドで直線を引いています. fit a*x+b "データが入ったファイル.txt" every ::2::5 via a,b replot a*x+b title"least-squares method" よろしくお願い致します.

  • 最小二乗法の問題

    (x、y)についてのん組のデータ(x1、y1).......(xn、yn)が直線y=ax+bの近似できるとき、aとbを最小二乗法により求めよ。

  • 円の最小二乗法の公式

    いくつかのデータから最小二乗法で近似曲線を求めたいのですが、よくわかりません。そのデータ集の近似曲線は円になります。 最小二乗法を調べ、1次、2次関数についてはわかりました。ある点の座標を(x1,y1), (x2,y2)…、近似曲線上の座標を(x1,y’1),(x2,y’2)… とした時、 (y’1-y1)^2 + (y’2-y2)^2 … が最小となるような係数a,b などを偏微分 → 連立方程式で求めるという方法でした。 円についても、同様の方法で r^2 = (x-a)^2 + (y-b)^2 のような近似曲線の式が求められるのでしょうか?1次関数などのように、y’1-y1を求めようとすると、±√ が出てきてしまい、ややこしくなる気がしますが、これを解くしかないのでしょうか?もしくは別の解法があるのでしょうか?詳しく教えていただけたらと思います。よろしくお願いします。

  • 最小二乗法

    n組のデータ (xi, yi) を,特定点(X0, Y0) を通る直線 y = ax+b でフィッティングしたい。最小二乗法で係数a,bを求めるため の式を導きなさい。 という問題で 各データの残差を二乗した和が最小になるときのa,bを求めるのですが 特定点(X0,Y0)を通るにはどうすればよいでしょうか? ただ単に、特定点を通らずフィッティングするやりかたはわかるのですが・・・。 よろしくお願いします。

  • 3次元の最小二乗法

    A(x1,y1,z1)とB(x2,y2,z2)とC(x3,y3,z3)とD(x4,y4,z4)の点で 最小2乗法ほうを使い直線近似したいのですが、どのような式になりますかご存知の方教えていただけますか?

  • 重みつき最小二乗法

    Yiに重みWiがあるとき、Y=aXに対して線形回帰を行うと、正規方程式より a = [WXY]/[WX^2] となるとおもうのですが、このときaの分散はどうなりますか? また、Xが正規分布に従うとき、c/Xはどのような分布に従うのでしょうか?やはり正規分布のような気はしているのですが・・・ よろしくおねがいします。

  • 最小二乗法について

    いくつか教えてほしいです。 min f(x)=||Ax-b||2・・・・(I) (1)||Ax-b||2ならわかるんですが||Ax-b||の上と下に2がつくのはどういったことでしょうか?? (2)(I)式を行列表現で表すと (A^T)Ax=(A^T)bとなりますが過程は <Ax,b><Ax,b> =<x,(A^T)Ax>-2<(A^T)b,x>+||b|となり最小化するので偏微分をし=0とおくのですが <A,x>の微分はAというのを用いれば-2<(A^T)b,x>の微分が2(A^T)bとなるのはわかるんですが <x,(A^T)Ax>の微分が2(A^T)Axとなるのがわかりません。 内積の微分は講義ではやってないので独学になります。 教授に聞いたらx=(x1,x2) (A^T)Aを2×2{{a,b},{c,d}}とおいて実際に計算し微分すれば理解できるかもと言われ計算してみたのですがよくわかりませんでした。 どうかご教授お願いします。

  • gnuplot 4次元データ(3次元座標における値)の表示

    まず、次のようなデータファイルがあります。 ----- x1, y1, z1, f1 x2, y1, z1, f2 x3, y1, z1, f3 x1, y2, z1, f4 x2, y2, z1, f5 x3, y2, z1, f6 以下、延々と続く。 ----- gnuplotで、zをz1に固定して、xy平面上に、fを立体的に表示する 方法があれば教えてください。 なお、pm3d で平面にカラーマップを表示する方法は分っています。

  • Gnuplotについて エラーメッセージ

    Gnuplotの使い方を覚えたくて勉強しています。 オプションを指定した上でプロットしようと、以下のように入力したところ、 plot [-5:5] f(x) axes x1y1 title "fx" with lines, a=1.0, b=0.16, f(x)=a*sin(x)+b*sin(3*x) 以下のようなエラーメッセージが表示されました。 function to plot expected プロットの文法が間違っていたのでしょうか。 このメッセージの意味するところがわかりません。 また、このようなエラーメッセージに遭遇した際に どのように対処すべきか、書かれたサイトなどが ありましたらお教えください。 (わたしが探した限りではみつかりませんでした) よろしくおねがいいたします。

  • GNUPLOTで点にラベルをつけて描きたい

    UNIXのGNUPLOTについて質問します。 例えば以下のようなデータがあったとします。 1.2 3.4 A 3.2 6.7 B 5.3 9.8 C この場合、1列目と2列目をそれぞれX,Yとし、3列目をその点または点につけるラベルとして2次元の散布図をGNUPLOTで書かせることは可能ですか?