任意の3点、点A(x1,y1,z1)、点B(x2,y2,z2)、点C(x3,y3,z3)を通る球体半径rの中心点D(a,b,c)を求めるために、エクセルで常に計算出来るように書式を作りたいのですが、ソルバーを使ってもうまくいきません。
どのように制約条件を定めたら良いのでしょうか。
分かる方ご教授のほどよろしくお願いします。
計算の参考にしたのはこちらの方法です。
https://detail.chiebukuro.yahoo.co.jp/qa/questio …
A 回答 (1件)
- 最新から表示
- 回答順に表示
No.1
- 回答日時:
こんにちは
>ソルバーを使ってもうまくいきません。
ソルバーの場合、初期値や変化量によっては収束しない可能性も考えられるので、例えば、初期値を何点か変えてみるなどしてみてもダメでしょうか?
一方で、ご提示の題意であれば計算でも求められるはずなので、(かなり面倒そうではありますが)そちらに関してを以下に。
現在、スマホからなので、文章のみで考え方だけの回答になりますが、ご容赦ください。
まず、考えを簡単にするために、3点で決まる平面をXY平面となるように座標変換します。
求めた座標を逆変換すれば元の座標系での位置にもどすことができますので。
(座標変換はご理解なさっているものと仮定。エクセルを利用するのであれば、行列計算もサポートしていますので、計算自体は比較的簡単にできると思います。)
XY平面上で3点の外心を求めると、この半径(Rとします)が、与えられる半径rの取り得る最小値であることは理解できるものと思います。
求める点(中心点)は外心から3点のなす平面に対する法線(=z軸に並行な線)上にありますので、その長さをdとすれば、
d^2 = r^2 - R^2
となるので、dを求めることができます。
dは正負の2点存在しますが、XY平面に対して上下対象な位置に題意を満たす中心が存在することを考えれば理解できると思います。
また、r<Rの場合は、実数解が存在しませんが、これは「rが外心半径よりも小さな値になり得ない」ことと一致します。
一番面倒そうなのが、外心の座標を算出する計算ですが、ゴリゴリと計算してくださっているサイトがありましたのでそちらにお任せすることにしておきます。
https://qiita.com/tatesuke/items/59133758ec25146 …
外心の座標x,yと上記のdから、求める中心座標は(x, y, ±d)として算出できるはずです。
エクセルをご利用になるのであれば、計算を分解して、それぞれの(途中)結果をセルに表示させながら、順次計算をして行くようにすれば、計算式が長く複雑になることを防げることができると思います。
若干の計算誤差が蓄積される可能性はありますが、それは、まとめて計算した場合でも同様であろうと推測します。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# ある線が円の範囲に入っているかの計算 1 2022/12/07 16:14
- 物理学 中心を同じに点に持つ半径aの導体球(導体1)、内半径b、外半径Cの導体球殻(導体2)があるとして、導 1 2023/08/12 23:36
- 数学 写真の図は中心(a,b)半径rの円とその円周上の(x1,y1)における接線lと円の中心とlを結ぶ任意 4 2023/08/08 16:20
- 物理学 導体球殻 電場・電位 2 2023/01/28 11:51
- 数学 数学ベクトルに関しての質問 3 2022/05/25 23:21
- 数学 球面と接する直線の軌跡が表す領域 4 2023/07/30 12:37
- 統計学 統計学、エクセルがわかりません!解答と詳しい解説をお願いします! (1)それぞれの地域別に記述統計量 9 2022/08/21 16:30
- 数学 数Bです。 定点O、Aと動点Pがある。ベクトルOA=ベクトルa、ベクトルop=ベクトルPとするとき、 3 2022/07/04 23:12
- 物理学 物理の問題です 2 2022/12/17 22:43
- 物理学 原点中心とする半径10cmの演習上、質点が1分間に600回の割合で反時計回りに運動している。 (1) 4 2023/05/29 12:46
このQ&Aを見た人はこんなQ&Aも見ています
関連するカテゴリからQ&Aを探す
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
YUVって?
-
Excelである地点に一番近い座標...
-
彼女とハグをする時胸は当たる...
-
[乗換]
-
緯度、経度 1分は何km?
-
昔の住所表記を現在の住所表記...
-
Excelでの緯度経度入力方法
-
東西南北の方角がどうしてもわ...
-
Excelについての質問です。 2点...
-
google map 地図のスケールが分...
-
50・60・70・80年代の世界地図
-
北海道から沖縄までの距離は何...
-
北緯・緯度の読み方について
-
地図の縮尺について
-
町域とは
-
郵便番号の下四桁が「0000」と...
-
「大字」は、もう必要ないので...
-
この地図記号は何の意味ですか...
-
絶対値についてです |a+b|^2 ...
-
記号の意味 GH、FH
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
Excelである地点に一番近い座標...
-
YUVって?
-
ノーズR補正計算
-
エクセルを使って円周率を作る方法
-
ある三点を通る球体の半径rの...
-
Excel リストボックスについて...
-
彼女とハグをする時胸は当たる...
-
メルカトル図法での大圏コース...
-
Excelについての質問です。 2点...
-
メルカトル図法
-
緯度、経度 1分は何km?
-
東京(北緯34度41分、東経139度...
-
この地図記号は何の意味ですか...
-
昔の住所表記を現在の住所表記...
-
google map 地図のスケールが分...
-
[乗換]
-
東西南北の方角がどうしてもわ...
-
北海道から沖縄までの距離は何...
-
町域とは
-
縮尺5万分の1の地図上で1センチ...
おすすめ情報