すみませんがお助けください。
任意の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ランキング
-
いちおう元彼に解いてもらった...
-
等分除と包含除について 割り算...
-
数学の反例について。 P⇒Qの反...
-
数学のユークリッドの原論はな...
-
箱ひげ図についての質問です。 ...
-
「帰納法とは、本来前提となる...
-
ぜみ
-
スカラー場とベクトル場
-
数学のやり直しをしている社会...
-
②の後、「よって、」の直後がわ...
-
√0.25=±0.5である。 これはなぜ...
-
仕事をクビになり会社の門で憔...
-
今は、割合が入った数学の問題...
-
フルランク
-
(1)長さlの棒の中点を通り、棒...
-
ガンマ関数Γ(x)は、階乗からど...
-
微分方程式の計算でわからなく...
-
京都大学で出題された次の問題...
-
学生時代、偏差値60以上の高校...
-
社会人で数3Cをやり直したいの...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
今更で申し訳ないのですが、疑...
-
数学の法則を発見しました
-
画像の質問①〜③に答えてわかり...
-
小数点の引き算
-
わかりませんでした。
-
√2の値の計算方法を教えてくだ...
-
f(z)=tan(z)のマクローリン展開...
-
数学がわからない。
-
複素数平面
-
2乗を平方。3乗を立方。それじ...
-
a, bがa>0, b>0,1/a+2/b=3を満...
-
なにしてるかわからない。。
-
写真の様な解き方はおかしいで...
-
数学Aの合同式について質問です...
-
級数の係数を求める
-
不完全定理により、「ある命題...
-
複素数平面について質問です。 ...
-
新高一です。 数1の因数分解が...
-
ならなくない??
-
三角関数
おすすめ情報