A 回答 (5件)
- 最新から表示
- 回答順に表示
No.1
- 回答日時:
最小二乗法とは残差二乗和を各係数で偏微分した式を=0とした連立方程式を立てて解くことですので、
残差二乗和をS=Σ(ax^2+by^2+2cz-1)^2 として、
各係数a,b,c毎に偏微分した式で連立方程式を立てると
(Σx^4)a+(Σx^2 y^2)b+2(Σx^2 z)c = Σx^2
(Σx^2 y^2)a+(Σy^4)b+2(Σy^2 z)c = Σy^2
(Σx^2 z)a+(Σy^2 z)b+2(Σz^2)c = Σz
となります。
この連立方程式を解けば係数a,b,cが求まります。
参考URL:http://www.eli.hokkai-s-u.ac.jp/~kikuchi/ma2/cha …
この回答への補足
回答ありがとうございます。
例えば(x1,y1,z1)における誤差dは
d=z1-(aXk^2+bYk^2+1)になるのではないのですか?
ようするに
残和二乗和Σd=Σ{Zk-(aXk^2+bYk^2+1)}になると思うのですが、
どうでしょうか?
ご意見お願いします。
No.2
- 回答日時:
こんにちは.
データの分布と誤差の乗り方によってはかなり難しくなりますよ.
というのも,ax^2+by^2+2cz=1は楕円の「標準型」で,
適切な直交変換をかけてxy平面の断面が楕円に,長軸と短軸がx軸とy軸になるようにしたものだからです.
入力する観測データは上記のことが期待できるのでしょうか?
また,誤差の分布はどのようになっていますか?
観測データの真値(x^,y^,z^)が楕円放物面に拘束されているとして,
観測データ(x,y,z)が(x^+u,y^+v,z^+w)という誤差を含んで観測されるとなると,
パラメターの現在値(a,b,c)で書ける楕円放物面への直交射影を計算せねばなりません.
さらにパラメターはa>0,b>0なる非負制約のもとで更新する必要があります
(単純な最小二乗法ではa>0,b>0となる保証がなく,「楕円放物面」以外の二次曲面が当てはまる恐れがあります).
これらの情報が不足しています.問題設定をもう少しきちんとした方が回答もつきやすいと思います.
No.3
- 回答日時:
普通なら(線形だろうが非線形だろうが)エクセルソルバで最小二乗法は解くのですが、数千となるとだめなのでしょうか。
となるとパラメータa, b, cについて線形なのでlp_solveはだめでしょうか?http://lpsolve.sourceforge.net/5.5/
あるいはoctaveを使うとか。
http://ja.wikipedia.org/wiki/GNU_Octave
No.4
- 回答日時:
おはようございます.
まだ閉じていなかったのですね.
さて,k番目のデータの真値(xk,yk,zk)が(xk,yk,zk+dk)と観測されれば,
残差は次のように書けます.
dk=zk-(axk^2+byk^2+1)/2c ...(1)
これはzのみに誤差がのっている場合に対応します.
とりあえず,これの平方和を最小化する方法を,と言うことであれば,
cの逆数を定義するか何かすれば,
非負制約つき二次計画(QP)辺りを利用できるかもしれません.
線型計画(LP)でないのは「二乗誤差」の最小化が目的だからです.
裸のoctaveでは結構難しいと思います.
この意味ではExcelソルバのGRG2は最も良い選択肢の一つでしょう.
しかしながら,z以外にも誤差がある場合は事情が異なります.
これは(xk+dxk, yk+dyk, zk+dzk)を元の式に代入すればすぐに分かりますが,
式(1)のように誤差関数を簡単には書けません.
したがって,パラメターの現在値で書ける曲面上にあって,
観測データに最も近い点を(場合によっては数値的に)
計算する必要が生じるわけです(直交射影).
誤差の分布によっては「評価関数の計算」からして骨が折れますが,
これが計算できなければどんな方法でもパラメターは更新できません.
Excelソルバも同様でしょう.
だから必要以上に問題が難しくならないように「補足要求」をしたのです.
参考までに述べると「最も一般的に」考えると,
修士以上の学位論文が執筆可能な内容になると思われます.
No.5
- 回答日時:
No.3です。
No.1で例示されたURLには線形最小二乗法とかかれており、モデル関数が線形の場合は連立方程式を解けば求まりますが、本問は非線形なので数値的に解くしかないと思います。私がNo.3でa, b, cについて線形と申し上げたのは誤りなので訂正させてください。lp_solveは使えません。ちなみに経験者も一般人の誤りでした(寝ぼけてたのかな)。
さて、非線形のカーブフィッティングはoctaveでもできるかとは思いますが、私もここで質問して(答えは得られませんでしたが)調べたところwnlibというAPIで解くことができました。私の場合はlogistic curveへのフィッティングでしたが、ご参考まで。
参考URL:http://www.willnaylor.com/mantext/wnnlp.txt
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 どうして放物線ですか? 15 2023/06/11 09:53
- 数学 初歩的な質問で申し訳ないのですが、 平面における直線→y=ax+b 空間における直線→ax+b=y= 2 2022/04/01 13:22
- 数学 ラグランジュの未定乗数法を用いる問題 3 2023/05/15 14:48
- 物理学 歌口と楕円形の太鼓 1 2023/05/15 23:21
- 日用品・生活雑貨 (昔の?)レンタル用DVDケースはなぜあの形状? 1 2023/03/29 19:41
- DIY・エクステリア VVF2.0 2Cを通す為の楕円の穴を天井(2cm程度の板)に開けるには? 4 2022/10/06 18:38
- 数学 画像の問題について質問です。問題式を楕円の式に変形して、積分範囲を0<=x<=a √(z^2-1) 3 2022/08/29 13:44
- JavaScript 最小二乗法 2 2023/01/01 20:57
- 数学 平面と円柱の交点の求め方を教えていただきたいです. 6 2023/08/12 18:38
- Word(ワード) ワード。フリーフォームの使い方が分かりません。 1 2022/10/06 16:18
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
近似曲線の数式を手計算で出し...
-
Σの添え字について
-
Π←これは一体?
-
シグマの記号の読み方
-
2重ΣΣのΣ記号は交換可能でしょ...
-
Σの下にくるk=1のkってなに...
-
平面の計算方法
-
Σk(k+1) k=1 式を教えて下さい ...
-
Σの上が2n
-
Σx^2と(Σx)^2の違いは?
-
数列の問題です。次の数列の和...
-
二重和(ΣΣ)の計算方法について
-
19 Σk k=6 の和を求めろという...
-
理系数学プラチカの45(2)のまた...
-
最小二乗法における有効数字に...
-
a1=1,an+1=an+3n-1 この条...
-
分散について
-
Σのk=2
-
数学で答えを教えて欲しいので...
-
調和数列の和なんですが。。。。。
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報