積分の項を単純化する
T=e−12((zy−μxσx)2−y)yk/2−2
多項式検索ようにp(y)
[p(y)e−12((zy−μxσx)2−y)]′=p′(y)e−12((zy−μxσx)2−y)+p(y)[−12((zy−μxσx)2−y)]′e−12((zy−μxσx)2−y)=T
次のようなを見つけることになりますp(y)
p′(y)+p(y)[−12((zy−μxσx)2−y)]′=yk/2−2
または
p′(y)−12p(y)(zμxσ2xy−2z2σ2xy−3−1)=yk/2−2
すべてのべき乗を別々に評価して行うことができますy
コメント後に編集する
上記のソリューションは、発散するため機能しません。
しかし、他の何人かはこのタイプの製品に取り組んできました。
Fourrier変換の使用:
Schoenecker、Steven、Tod Luginbuhl。「2つのガウス確率変数の積とガウスおよびガンマ確率変数の積の特性関数。」IEEE信号処理レター23.5(2016):644-647。
http://ieeexplore.ieee.org/document/7425177/#full-text-section
およびの積に対して、特性関数を取得しました。Z=XYX∼N(0,1)Y∼Γ(α,β)
φZ=1βα|t|−αexp(14β2t2)D−α(1β|t|)
ウィテカーの機能(http://people.math.sfu.ca/~cbm/aands/page_686.htm)Dα
Mellin変換の使用:
スプリンガーとトムソンは、ベータ、ガンマ、およびガウスの分布確率変数の積の評価をより一般的に説明しています。
スプリンガー、MD、WEトンプソン。「ベータ、ガンマ、ガウス確率変数の積の分布。」SIAM Journal on Applied Mathematics 18.4(1970):721-737。
http://epubs.siam.org/doi/10.1137/0118065
Mellin積分変換を使用します。メリンは、変換のメリン変換の積である及び(参照http://epubs.siam.org/doi/10.1137/0118065又はhttps://projecteuclid.org/euclid.aoms/1177730201を)。研究された製品の場合、この製品の逆変換は、計算方法を提供および証明するメイジャーG関数として表現できます。ZXY
同じ手法を使用できる可能性はありますが、ガウス分布変数とガンマ分布変数の積は分析しませんでした。すぐにこれを行おうとすると、H関数(https://en.wikipedia.org/wiki/Fox_H-function)を取得することが可能になるはずです。ただし、G-機能するか、他の単純化を行います。
M{fY(x)|s}=2s−1Γ(12k+s−1)/Γ(12k)
そして
M{fX(x)|s}=1π2(s−1)/2σs−1Γ(s/2)
あなたが得る
M{fZ(x)|s}=1π232(s−1)σs−1Γ(s/2)Γ(12k+s−1)/Γ(12k)
の分布は次のとおりです。Z
fZ(y)=12πi∫c+i∞c−i∞y−sM{fZ(x)|s}ds
少なくともH関数として項を削除するために変数を変更した後)232(s−1)
残っているのは、この逆メリン変換をG関数として表現するパズルです。と両方出現はこれを複雑にします。ガウス分布変数のみの積の別のケースでは、変数代入することにより、をに変換できます。しかし、カイ二乗分布の条件のために、これはもう機能しません。たぶんこれが、このケースの解決策を誰も提供していない理由です。ss/2s/2sx=w2