標準正規分布に従う2個の確率変数の積が従う確率密度関数を計算してみた。

はじめに

かなり前の記事で、ちょいと野暮用で正規分布に従う2個の確率変数の和と差が従う確率密度関数を計算してみました。

この記事はその続編のようなものです。

調子に乗って~~あまり使わないかもしれない~~標準正規分布関数に従う2個の確率変数の積を計算しようと試みたわけですが、本Webサイトの管理人たるpandaの計算の能力が未熟なせいもあってか、2018年の年末にものすごい勢いでハマってしまったので、メモっておくことにしました。

問題の定義

まず、解くべき問題を以下のように定義します。

標準正規分布に従う2個の確率密度変数がある。このとき、によって定義される確率変数が従う確率密度関数を求めよ。ただし、とする。

なお、なぜとしたかというと、

「場合分けが面倒だからです。」

と潔く(本Webサイト比)言い切っておくことにします。

[2019/08/17 追記] の場合についても別途書きました。こちらをご参照いただけると幸いです。

サクサクと計算。

の符号による場合分け。

問題を定義したところで、サクサクと計算します。

まず、積分の領域を定めます。

積分を実行するの範囲を2次元平面上にプロットすると、以下の図の水色で示した部分になります。

この記事の最初の方でといきなり変数の範囲を限定したにもかかわらず、の符号によって場合分けをする必要がありそうです。(´・ω・｀)

の符号を考慮しつつ、の積で定義される確率変数のとる値が以下となる確率は( $1$ )式で表されます。

[2019/01/12追記]不等号の向きが間違っていたので、修正しました。

(1)

( $1$ )式の右辺の第1項は( $2$ )式のように計算できます。上図の(a)の向きに最初に積分するイメージになります。なお、( $2$ )式の右辺ではをと書き換えています。

(2)

ここで、ちょっと気になるのはの場合のところですが、この場合は( $2$ )式のうち変数による積分区間がになるので、

(3)

と計算できて積分可能ですので、(

2

)式は

の場合にも成り立つと考えることができます。

さらに、( $2$ )式でと置きます。すると、( $2$ )式の右辺は、

(4)

と計算できます。

次に、( $1$ )式の右辺の第2項を計算してみます。

( $1$ )式の右辺の第2項は、

(5)

と計算できます。上図の(b)の向きに最初に積分するイメージになります。これを(

4

)式と似たような式に変形できないか考えてみます。

まず、は偶関数ですので、( $5$ )式の右辺は( $6$ )式のように書き換えることができます。

(6)

次に、

と置きます。すると、(

6

)式の右辺は、

(7)

となりますので、さらに

と置きます。すると、(

7

)式の右辺は、

(8)

となります。

( $8$ )式のをに置き換えると( $4$ )式の右辺と同じ式が得られます。

したがって、は( $9$ )式で表されます(をに置き換えています)。

(9)

( $9$ )式の及びについての積分の順序は入れ替えることができます。そこで、積分の順序を入れ替えてからで微分すると、における確率密度関数が以下のように計算できます。

(10)

さらに変数変換。

ここからは結論がある程度見えていないと計算をするモチベーションが起こりにくいところではありますが、( $10$ )式をすでによく知られた形に変形できないか考えてみます。

( $10$ )式でとおきます。すると、となりますので( $11$ )式のように変形できます。

(11)

最後の仕上げに、とおきます。( $11$ )式の積分区間はですが、の積分区間はとなることと、であることから、

(12)

となります。

ここで、( $12$ )式の右辺ののの部分をじーっと見ると、はに依存しない定数であり、は偶関数であることから、( $12$ )式の右辺は、

12 式 の 右 辺 (13)

となります。

( $13$ )式は0次の第2種変形Bessel関数(Modified Bessel Function of the Second Kind, )の積分表示をで割ったものに等しくなりますので、確率密度関数は( $14$ )式のように表すことができます。

(14)

おまけ

この記事を最初に書いた時点(2019年1月)では実際に計算してみたわけではありませんが、の場合は前節の同様の手順で( $15$ )式となることを示すことができます。 [2019/01/12追記] こちらのページに書きました。

(15)

そこで、(

14

)式及び(

15

)式をまとめて(

16

)式と書くこともできます。

(16)

まとめ

2019年は新年最初の記事からいきなり最初から数式満載の記事でお送りしてみました。

[2019/01/12,2019/02/06追記]ちょこちょこと間違っているところがあったので、修正しました。

を0次の第2種変形Bessel関数で表すために変数の変換を何回か繰り返していますが、最後のとおいたところで置き換え後の積分の範囲が誤っていたことに気づかず、どうしても分母の2が取れずにハマってしまいました。

今後気を付けます。

Bessel関数は日常生活においてはそう簡単に登場することはないと思いますが、何かの参考にしていただけると幸いです。

この記事は以上です。

References / 参考文献

Bessel function (英語版Wikipedia)

はじめに

問題の定義

サクサクと計算。

の符号による場合分け。

さらに変数変換。

おまけ

まとめ

References / 参考文献

Related Posts / 関連ページ

Ads Blocker Detected!!! / 広告ブロッカー的なものを検知しました!!