![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?8acaa2e)
すみませんがお助けください。
任意の3次元座標(Xi,Yi,Zi)がn個あるとします。
Xi>0,Yi>0,Ziは正負混合です。
上記n個の座標から近似平面式Z=aX+bY+cを構成する乗数a,b,cを求めたいのですが下記理由により最小二乗法は使えません。
最小二乗法の場合は
D=Σ(Zi-aXi-bYi-c)^2が最小になるa,b,cを求めることになりますが、今回求めたいa,b,cは、
各座標のdi=(Zi-aXi-bYi-c)、すなわちn個のdiについて『diの最大値-diの最小値』が最も小さくなる場合のa,b,cです。
なお『diの最大値-diの最小値』についての補足ですが
diの最大値=5,diの最小値=2の場合
『diの最大値-diの最小値』は5-2=3となり、
diの最大値=5,diの最小値=-2の場合
『diの最大値-diの最小値』は5-(-2)=7となります
お手数かけて恐縮ですが宜しくお願い申し上げます。
A 回答 (2件)
- 最新から表示
- 回答順に表示
No.2
- 回答日時:
>ここの再度「ツール→ソルバー~」の操作を不要にして、Xi、Yi、Zi の一部数値を変更すると、即a,b,cの数値が変更される方法はありますでしょうか?
これはわかりません。マクロとかVBAとかで出来るのかどうか??ですが、エクセルに詳しい方の回答に期待するか、あるいは再度ここで質問されてみたらいかがでしょうか?
プログラムを組めばエクセルソルバーのような手作業は省けますが、この問題は最小自乗法ではないのでニュートン法などが使えず結構難しそうです。多少時間がかかっても問題なければ乱数を使ったモンテカルロ法か遺伝的アルゴリズムのような手法が使えるかもしれませんが。。。
最初のご回答が期待以上の内容で、大いに業務効率UPにつながりました。ここから先は自分で色々とやってみることにします。
本当にありがとうございました。
No.1
- 回答日時:
エクセルのソルバーを使えば下記の手順で a、b、c が求まります。
以下は EXCEL2002の場合で、データ数は n = 20 とします。
1)B~D 列にデータ Xi、Yi、Zi を入力する。
(第1行は xi 等のラベルを入れ、データは第2行から第23行まで)
2)K2~K4 に可変パラメータ a、b、c の初期値(例えばゼロ)を入れる。
3)E2 に =D2-$K$2*B2-$K$3*C2-$K$4 を入れて、E3~E23 にそれをコピーする。
4)G2 に =MAX(E2:E23) と入力する。
5)G3 に =MIN(E2:E23) と入力する。
6)G4 に =G2-G3 と入力する。
7)G4 を選択した状態でツール→ソルバーをクリックする。
(ここでソルバーのウインドウが開き、目的セルが $G$4 となっている。)
8)ソルバーウインドウの目標値を◎最小値に選ぶ。
9)ソルバーウインドウの変化させるセルを選択した後にワークシートのK2~K4を選択する。
10)ソルバーウインドウの実行ボタンをクリックする。
11)最適解が見つかったら◎解を記入するを選択してOKボタンをクリックする。
12)上記の結果 K2~k4 に解 a、b、c が得られる。
蛇足ながらツールの中にソルバーが見当たらない時はアドインする必要があります。
またこの問題では問題ないようですが最適解に収束しない場合は K2~K4 の初期値を変えます。
ありがとうございました。おかげさまで解決しました。
ところで一旦「ツール→ソルバー~」の操作をしてa、b、cを求めた後に、Xi、Yi、Zi の一部数値を変更した場合は再度「ツール→ソルバー~」の操作をしてa,b,cを求めることになりますが、ここの再度「ツール→ソルバー~」の操作を不要にして、Xi、Yi、Zi の一部数値を変更すると、即a,b,cの数値が変更される方法はありますでしょうか?もしあればご教示いただきたく存じます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 工学 過渡現象 RL回路で 回路方程式 Ri+L(di/dt)=E 定常解 is=E/R 過渡解 Ri+L 1 2022/06/08 21:47
- Word(ワード) ネットの文章をワードに貼り付けるとき、幅を広げる方法 6 2022/05/30 11:04
- 数学 『iの微分』 7 2023/02/06 19:45
- JavaScript 最小二乗法 2 2023/01/01 20:57
- 数学 整数問題5 類難題 6 2023/04/08 00:05
- 数学 条件付き極値問題といわれる問題です。ラグランジュの乗数法 について、質問したいことがあります。 条件 3 2023/05/15 21:38
- デジタルカメラ 動画撮影に使うレンズ。どちらのレンズがいいか悩んでいます 1 2023/01/02 18:21
- 楽器・演奏 RolandのJuno di で、デュアルモードで作った音をスプリットモードで使用することはできない 1 2023/08/10 17:57
- モデル 外国人のモデルさんについて 2 2022/09/25 13:29
- オープンソース イラストのDLサイト 1 2023/05/11 17:33
このQ&Aを見た人はこんなQ&Aも見ています
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
小学生の時(40年前)に、18÷...
-
見づらいですけど、同時型微分...
-
同一地点からAが時速 15km、Bが...
-
ドモアブルの定理はサインまた...
-
写真のように
-
数Bの階差数列についててす。 ...
-
微分方程式 について d²y/dx² ...
-
数Ⅲです 写真の問題がわかりま...
-
フーリエサイン変換の
-
高校一年生です。 私はどうして...
-
1+2+3+…=?
-
数学Aの問題です。 A高校の男子...
-
整数問題です。
-
複素関数論のローラン展開について
-
√-1 は、何になるのでしょうか
-
数Ⅲ極限 写真の問題の途中計算...
-
fft
-
どうしても欲しい本が()理工書
-
数Ⅲ極限です。 limx→-♾️ sinx/x...
-
偏微分方程式って
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報