3次元での点群に対する最小二乗法での平面の算出について(点と平面の距離。残差ではない。)
--
点と平面のZ軸方向の距離(残差)の二乗和を最小とする場合には、
平面をax+by+c=zとして、Σ(ax+by+c-z)^2をa,b,cのそれぞれで偏微分して
それを=0とした連立方程式を解くことで解を得ることが出来ました。
また、式の形も、ある点のxとyを平面の式へ代入した際の値と、点のz値の差分を見ており、
簡単に納得のできるものとなりました。
これに対して、点と平面の距離(空間的な最小距離)の二乗和を最小とする場合には、
どのような流れで計算すれば良いのでしょうか?
点と平面の距離は|Ax+By+Cz+D| (A,B,Cは単位ベクトル)として求まりますが、
これをどう使うのかが分かりません。
Σ(Ax+By+Cz+D)^2をA,B,C,Dのそれぞれで偏微分して=0としても、
定数項が無いため、連立方程式の解がすべてゼロとなってしまいます。
強引に、Σ(A'x+B'y+C'z+1)^2として変形させて解いてみましたが、
得られたA',B',C'からA,B,C,Dに戻すと、Dがきちんと出ませんでした。(他についても怪しい。)
こういった状況に迷い込んでしまい、どう考えるのが良いのか分からなくなってしまいました。
指南いただけませんでしょうか?
No.2ベストアンサー
- 回答日時:
平面の式は、単に Ax+By+Cz+D=0 としたのでは、一意に決まりません。
同じ平面が、 2Ax+2By+2Cz+2D=0 とでも 3Ax+3By+3Cz+3D=0 とでも
書けるからです。
そのために、「(A,B,C) は単位ベクトル」としたのではありませんか?
だから、Σ(Ax+By+Cz+D)^2 を最小化するときに、単なる最小値でなく、
A^2+B^2+C^2=1 という制約下での最小値を探せばよいのです。
ラグランジュの未定乗数法が使えます。
あるいは、制約なしで、Σ(Ax+By+Cz+D)^2/√(A^2+B^2+C^2) を最小化
してもよいのだけれど。
この回答への補足
とても納得できました!
ラグランジュの未定乗数法を手計算でやってみたのですが、λの導出が3次方程式になってしまったのと、係数が多くなり過ぎて力尽きました。。。間違えたのでしょうか?
Σ(Ax+By+Cz+D)^2/√(A^2+B^2+C^2) を最小化する場合には、これも非線形であるために、
偏微分=0として解くには難解となりそうだな。。。と躊躇しております。
非線形最小化であるため、ニュートン・ラフソン法であればすんなりと出来たのですが、
この問題に対して妥当な方法と言えるでしょうか?
No.4
- 回答日時:
No.3です。
3次元で考えれば各点から平面に下ろした垂線の足の長さになりますから、Z方向の差分ではないと思いますよ。
この回答への補足
追記ありがとうございます。
Σ(ax+by+c-z)^2として解くと、平面の方程式に点のxとyを代入して算出される本来のz'値ax+by+cと、点のz値との差分を二乗しているため、各点から平面のZ方向の差分であると考えています。
どこに誤りがあるのでしょうか?いただいた説明では分かりませんでした。。。
No.1
- 回答日時:
平面の方程式をax+by+cz+d=0とすると、
点(xi,yi,zi)と平面の距離=|axi+byi+czi+d|/√(a^2+b^2+c^2)
したがって、
二乗和=Σ[i=1~n](axi+byi+czi+d)^2/(a^2+b^2+c^2)
これをa,b,c,dで偏微分した式を連立させてa,b,c,dを求めることに
なりますが、分母にも変数があることから、非線形連立方程式を
解くことになろうと思います。
最小二乗直交距離フィッティング等のキーワードで検索して
みてもよいと思います。
この回答への補足
上の補足にも記載したのですが、非線形連立方程式という形で一発で解けるのでしょうか?
非線形最小化であるため、ニュートン・ラフソン法であればすんなりと出来たのですが、
この問題に対して妥当な方法と言えるでしょうか?
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 初歩的な質問で申し訳ないのですが、 平面における直線→y=ax+b 空間における直線→ax+b=y= 2 2022/04/01 13:22
- 数学 ベクトル方程式(ヘッセの標準形)についての質問 2 2022/04/23 18:00
- 数学 3次元実ベクトル空間において, 平面 P:x-y+z+1=0 と直線 L:2(x-1)=-y=-z 3 2022/10/29 14:39
- その他(教育・科学・学問) 関数、写像について 1 2022/04/10 23:45
- 数学 (1)の平面の式を求める問題で ABベクトルとACベクトルの外積が平面の法線になるから ax+by+ 2 2023/04/13 13:50
- 数学 球面と接する直線の軌跡が表す領域 4 2023/07/30 12:37
- 数学 都市経済学の問題です。わかる方教えてください。 地主が得る地代を増やすために、都市の端までの距離 x 2 2023/07/18 17:41
- 数学 条件付き極値問題といわれる問題です。ラグランジュの乗数法 について、質問したいことがあります。 条件 3 2023/05/15 21:38
- 数学 数学得意な方教えてください 平面上の任意の点から同一平面上にある多角形の各頂点までの距離の平均は、多 1 2023/08/28 16:36
- 数学 数2Bの数列の問題です。 自分は、 まず数列 an=ar^(n-1)と置き こちらの問題の、y= の 1 2022/07/07 16:26
このQ&Aを見た人はこんなQ&Aも見ています
-
プロが教えるわが家の防犯対策術!
ホームセキュリティのプロが、家庭の防犯対策を真剣に考える 2組のご夫婦へ実際の防犯対策術をご紹介!どうすれば家と家族を守れるのかを教えます!
-
平面の計算方法
数学
-
近似式(z=ax+by+c)を取得したい
その他(コンピューター・テクノロジー)
-
XYZ平面の近似式についてご教示ください
数学
-
-
4
最小二乗平面
C言語・C++・C#
-
5
3次元の近似直線
数学
-
6
仮想平面の計算
数学
-
7
任意の面内にある点の座標から面の傾きを求める方法を教えて下さい。
数学
-
8
EXCELにてローパスフィルタを作成する
その他(教育・科学・学問)
-
9
Excelで平坦度
その他(Microsoft Office)
-
10
平面の方程式、dが持つ意味?
数学
-
11
3次元における近似平面と点Pの交点を求めたい
数学
-
12
三次元形状曲面の導出法
その他(プログラミング・Web制作)
-
13
平面度
数学
-
14
N点間の中心と重心の求め方
数学
-
15
空間上の測定された点群から最小二乗法により円(中心、直径)を算出する方法
数学
-
16
3次元空間で3点を通る平面を2次元座標で表すには
数学
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
全員と同じグループを経験でき...
-
5406を13で割ったときの絶対値...
-
3次元での点群に対する最小二...
-
数学 3次関数の最小値・最大値...
-
高校数学1の問題集に、2次関数...
-
距離の和を最小にする点を求め...
-
数学の公式の実践的な使い方を...
-
高校数学で分からない問題があ...
-
数A 整数の性質 x.yを整数とす...
-
数学の平面ベクトルの問題なの...
-
アルキメデスの公理について 任...
-
楕円のxの最大値と最小値
-
2進数のバイアス表現について
-
2つの放物線間の最短距離
-
急いでいます 数学の問題
-
至急!!二次関数について aは...
-
以下の説明文の変数係数とは何...
-
「無限の一つ前の数字は何?」...
-
年代と年台・・・どちらが正し...
-
フーリエ変換後の負の周波数成...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
5406を13で割ったときの絶対値...
-
全員と同じグループを経験でき...
-
3次元での点群に対する最小二...
-
斜線D中を通る直線の傾きkの最...
-
1/x+1/y≦1/2 , 2<x,2<yのとき、...
-
3で割ると2余り、7で割ると4余...
-
中学受験用の小5算数の問題です
-
至急!!二次関数について aは...
-
高校数学1の問題集に、2次関数...
-
Gnuplot 最小二乗フィッティン...
-
2進数のバイアス表現について
-
数学2です x>0のとき、x + 16/(...
-
EXCEL ドラッグしたセル...
-
正の約数の個数が20個である最...
-
間違いの理由を教えてください...
-
1/x+1/y+2/z=1を満たす自然数解
-
数学 3次関数の最小値・最大値...
-
y=x^xの最小値
-
数学Aの確率
-
2つの放物線間の最短距離
おすすめ情報