すみませんがお助けください。
任意の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ランキング
-
確率の問題 数学と実生活と
-
過去に 「ii) f(z)=1/(z^2-1) r...
-
整数問題 兎に角 難問です 千葉...
-
「an=(n-1)/(n+1)のときlim[n→∞...
-
こちらの式はtan(z)のローラン...
-
これの極限値を求める問題で、 ...
-
指数関数と階乗。グラフで表し...
-
△ABCの辺BCの延長上の点Dを通る...
-
数学Aについて、4でも6でも割り...
-
返信の続きはありますか
-
iphone の自分の声が高すぎるん...
-
30分の動画を2倍速で見たら、3...
-
2024.4.7 03:42の質問に対する2...
-
逆三角関数の方程式の問題です...
-
連続的ポストごめんなさい
-
2022年 東京理科大 難易度判定
-
方程式の計算について
-
背景は何
-
【数学・標準偏差σ】標準偏差の...
-
ベクトル3重積
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
整数問題 兎に角 難問です 千葉...
-
えこれわかるひといますか?
-
長方形の分類(幾何学)
-
京都大学理系 過去問 整数問題
-
天孫降臨の神武天皇のY染色体...
-
中学数学の図形の問題です。
-
数学I アホらしい質問なのでそ...
-
正方行列Aについて
-
このルートを外す計算どうすれ...
-
素数についての一考察
-
これて最後どうやりますか??
-
この数学の問題はどうやって解...
-
高校数学の整数問題です。
-
暗闇で2人が出会うには両方動...
-
なんでですか?
-
大学入試の数学で、解答を進め...
-
n^2+n-4032はどうやって解くん...
-
(-1) ^2πってなんで1じゃないん...
-
標準偏差
-
一橋大学過去問 整数 素数 かな...
おすすめ情報