が均一かつ通常選択されている場合、球面仰角はどのように分布しますか?


7

フォローアップとして 極座標方法、、分散されたときにとIF?θ(x,y)U(1,1)×U(1,1)(x,y)N(0,1)×N(0,1)

仮定どのようにしている及び分散しますか?(x,y,z)U(10,10)×U(10,10)×U(10,10)θϕ

が次のようになるのは、前の質問のすばらしい回答から明らかです。 θここに画像の説明を入力してください

しかし、なぜがで最大尤度を取得しないのですか?ϕϕ=π/4

ここに画像の説明を入力してください

正規分布でを選択すると次の2つのpdfが得られます。x,y,z

ここに画像の説明を入力してください

および分布の名前はどちらの場合にもありますか?私にとっては、区間分布のように見えます。θϕβ[90,90]


1
最後の1つはベータ版のディストリビューションではありません。トリガー関数で考えます。角度がデカルト座標にどのように関連しているかに注意してください。
Glen_b-2017

2
角度のみを考慮しているため、球に投影できます。ガウスの場合、球上で均一な分布になりますが、その状況では緯度が均一ではないことに注意してください。サイト上のいくつかの質問はいくつかの問題を議論しますが、前の質問と大きく異ならない基本的な議論はあなたをそこに導くことができます。たとえば、ここに
Glen_b -Reinstate Monica

先に質問するつもりでしたが、と定義する必要があります。私は彼らがあなたの議論から何であるかを仮定しましたが、私たちが同じページにいることを確認するのは素晴らしいことです。ベータ版で配布されない理由を明確にする回答を掲載しました。θϕ
Glen_b-2017

私はあなたの2番目のグラフィックを信じていません。ポイントが立方体に均一に分布している場合、緯度()の密度はゼロになりません。最後の2つのヒストグラムは確かにベータ分布です。その結果は、より高い次元に一般化されます。その説明はここにあります。「sphere」と「Beta」を一緒に検索してください。ϕ±π/2
whuber

2
申し訳ありません:あなたは正しいです。ヒストグラムの詳細を確認できました。私は直感を再調整する必要があります-したがって、あなたが指摘したこと(+1)に感謝しています。
whuber

回答:


7

ここでの議論では、は事実上経度であり、は事実上緯度であると想定しています。おそらく、より一般的な球面座標は、赤道からではなく北極からの角度を使用し、そこから2つのシンボルの役割を入れ替えますが、どちらにしても問題はないので、あなたの表記はそうです。θϕ

質問で使用されているように見える角度を示すグラフィック

半径の分布はここでは重要ではなく、角度のみに注意してください。したがって、角度を変更せずにすべてを単位球に投影できます。これは、通常の場合に非常に役立ちます。

3次元標準法線のような球対称の分布では、傾斜の分布の外観は、極の近くよりも赤道の近くの球の表面の方がはるかに広いという事実に関係しています。

赤道付近の緯度でより多くの領域を示すグラフィック

数学をたどる(または、前の2D質問と同様の確率の要素で幾何学的な引数を書く)と、傾斜の密度は比例するはずです。これが「確率の要素」という用語で動機付けられるべき幾何学的な議論です:cos(ϕ)

緯度での半径を示す画像phi is cos(phi)

赤道半径が1と緯度における半径であるので IS、で円周緯度 IS比例する、および密度そうで IS比例する。ϕcos(ϕ)ϕcos(ϕ)ϕcos(ϕ)


均一なケース:3D均一で一定の半径に正規化されている場合、2Dの場合と同じ理由で、球上で密度の均一性はありません。球に投影すると、さらに多くのものが得られます "角が角よりも角の近くにある球の(密度はエッジの中央付近の部分が間にあります)-角よりも角に近い角の方が角よりも立方体の体積が多いため顔の中央近く。

これは、立方体で多数のランダムな値を均一に生成し、それらを球に投影することで確認できます。キューブの面の近くよりもコーナーの方がボリュームが多いので、面よりもコーナーから「内側」に見える密度が高くなります。経度に対して赤道の上に高さをプロットすると(これは投影されたz値であり、ここで)、以下の一番上のプロットを取得します。z=z/rr=x2+y2+z2

単位球に投影された立方体$ [-1,1] ^ 3 $の多くのランダムな均一値のプロット。(i)高さ/経度と(ii)緯度/経度に変換

その高さは、前の図の直角三角形の垂直辺に対応しています。その高さはの()。これを緯度()に変換するには、投影された垂直方向の高さのアークシンを取得します。これは、下のプロットに表示されているものです。これにより、極に近づくほど「伸縮」し、緯度の関数としての密度が北極と南極で0に低下します(均一と通常の両方の場合)。sinϕz=sin(ϕ)ϕ

以下のための密度、その後オーバーその二変量密度の不可欠になります。ϕθ

ファイの限界を計算するための統合を示す2変量シータファイ密度の画像

の周辺(つまり、固定値で実行されるストリップ)を見ると、の密度に4つのピークが見られます。赤道から離れたの密度のペアのピークを作成します。これは、立方体のコーナーと上/下のエッジが投影する単位球の表面上の領域に対応します。θθθϕ


好奇心から、絵を描くのにどんなプログラムを使いますか?
0x90、2017年

4
@ 0x90これらの3D図面をペイントで作成しました。手で。(括弧の誤った配置に基づいているため、質問の最後のプロットはおそらく削除できます。)
Glen_b -Reinstate Monica

一様な場合は、質問に示されているのヒストグラムを提供できない可能性があります。したがって、質問のその部分に対する定量的な応答が最も役立ちます。ϕ
whuber

@whuberは実際に私は0x90のプロットに同意します ϕ。球を介したスライスの描画を見て、(球形座標で)単位球上のポイントに変換した後、赤道からの角度はZ座標のアークサインになるでしょうか?これにより、2つのヒストグラムの表示が生成されます。ϕ質問
Glen_b-モニカの復活2017

6

球面緯度の補完的な累積分布 ϕ 立方体のランダムな点が [1,1]3 関数をグラフ化する円錐の上にあります z=cot(ϕ)x2+y2。これらの点は、立方体全体に均一に分布しているためです(立方体の体積は8)、このチャンスは、円錐と立方体の上部の間の体積の1/8です。緯度が超えたときπ/4、このボリュームは高さのある右円錐のボリュームです 1 そしてベース cot(ϕ)、 に等しい

F+(ϕ)=18π3cot2(ϕ).

図の2つの右側のプロットを参照してください。

緯度が以下の場合 arctan(1/2)、これは半無限の円錐と立方体の交差の体積です。極座標での積分は次の式を与えます

F(ϕ)=18(443tan(ϕ)(2+2tanh1(tan(π8)))).

図の左端の2つのプロットを参照してください。

図

これらの式の負の導関数が密度を与えます。の間にarctan(1/2)π/5 そして π/4円錐と立方体の交差が複雑な遷移領域です。正確な表現を作成することはできますが、面倒です。私たちが知っていることは、密度はF の派生物に F+ なので ϕそれらのポイント間で異なります。これは、100万個のシミュレーション値のヒストグラムで示されます(立方体の上半分のみから、下半分は鏡像になります)。ゴールドカーブは、ddϕF 右側の赤い曲線は ddϕF+.

ヒストグラム

これにより、モードが ϕ=±π/4、ただしこれらの値と ±arctan(1/2)

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.