統計学 では 、一般化された マーカムQ関数 は次 のように定義される。
ν
{\displaystyle \nu}
質問
ν
(
1つの
、
b
)
=
1
1つの
ν
−
1
∫
b
∞
x
ν
経験
(
−
x
2
+
1つの
2
2
)
私
ν
−
1
(
1つの
x
)
d
x
{\displaystyle Q_{\nu }(a,b)={\frac {1}{a^{\nu -1}}}\int _{b}^{\infty }x^{\nu }\exp \left(-{\frac {x^{2}+a^{2}}{2}}\right)I_{\nu -1}(ax)\,dx}
ここで 、 およびは の 次数 の第一種修正ベッセル関数 です 。 のとき 、積分は任意の に対して収束します 。 マーカムの Q 関数は、 非心カイ 分布 、 非心カイ二乗分布 、および ライス分布 に対する 相補累積分布関数 として現れます。工学では、この関数はレーダーシステム、通信システム、待ち行列システム、および信号処理の研究に登場します。この関数は 、パルスレーダーの Jess Marcum によって に対して初めて研究され、その名が付けられました。 [1]
b
≥
0
{\displaystyle b\geq 0}
1つの
、
ν
>
0
{\displaystyle a,\nu >0}
私
ν
−
1
{\displaystyle I_{\nu -1}}
ν
−
1
{\displaystyle \nu -1}
b
>
0
{\displaystyle b>0}
ν
{\displaystyle \nu}
ν
=
1
{\displaystyle \nu =1}
プロパティ
有限積分表現
という事実を用いて 、一般化されたマーカムQ関数は、次のように有限積分として定義することもできる。
質問
ν
(
1つの
、
0
)
=
1
{\displaystyle Q_{\nu }(a,0)=1}
質問
ν
(
1つの
、
b
)
=
1
−
1
1つの
ν
−
1
∫
0
b
x
ν
経験
(
−
x
2
+
1つの
2
2
)
私
ν
−
1
(
1つの
x
)
d
x
。
{\displaystyle Q_{\nu }(a,b)=1-{\frac {1}{a^{\nu -1}}}\int _{0}^{b}x^{\nu }\exp \left(-{\frac {x^{2}+a^{2}}{2}}\right)I_{\nu -1}(ax)\,dx.}
しかし、マルカムQ関数の積分表現は、(i)積分の極限が関数の引数に依存しないこと、(ii)極限が有限であること、(iii)被積分関数がこれらの引数のガウス関数であることなどが好ましい。の正の整数値に対して、そのような表現は三角積分 [2] [3] で与えられる。
ν
=
ん
{\displaystyle \nu =n}
質問
ん
(
1つの
、
b
)
=
{
H
ん
(
1つの
、
b
)
1つの
<
b
、
1
2
+
H
ん
(
1つの
、
1つの
)
1つの
=
b
、
1
+
H
ん
(
1つの
、
b
)
1つの
>
b
、
{\displaystyle Q_{n}(a,b)=\left\{{\begin{array}{lr}H_{n}(a,b)&a<b,\\{\frac {1}{2}}+H_{n}(a,a)&a=b,\\1+H_{n}(a,b)&a>b,\end{array}}\right.}
どこ
H
ん
(
1つの
、
b
)
=
ζ
1
−
ん
2
π
経験
(
−
1つの
2
+
b
2
2
)
∫
0
2
π
コス
(
ん
−
1
)
θ
−
ζ
コス
ん
θ
1
−
2
ζ
コス
θ
+
ζ
2
経験
(
1つの
b
コス
θ
)
d
θ
、
{\displaystyle H_{n}(a,b)={\frac {\zeta ^{1-n}}{2\pi }}\exp \left(-{\frac {a^{2}+b^{2}}{2}}\right)\int _{0}^{2\pi }{\frac {\cos(n-1)\theta -\zeta \cos n\theta }{1-2\zeta \cos \theta +\zeta ^{2}}}\exp(ab\cos \theta )\mathrm {d} \theta ,}
そしてその比率は 定数です。
ζ
=
1つの
/
b
{\displaystyle \zeta =a/b}
任意の実数に対して 、そのような有限三角積分は [4]で与えられる。
ν
>
0
{\displaystyle \nu >0}
質問
ν
(
1つの
、
b
)
=
{
H
ν
(
1つの
、
b
)
+
C
ν
(
1つの
、
b
)
1つの
<
b
、
1
2
+
H
ν
(
1つの
、
1つの
)
+
C
ν
(
1つの
、
b
)
1つの
=
b
、
1
+
H
ν
(
1つの
、
b
)
+
C
ν
(
1つの
、
b
)
1つの
>
b
、
{\displaystyle Q_{\nu }(a,b)=\left\{{\begin{array}{lr}H_{\nu }(a,b)+C_{\nu }(a,b)&a<b,\\{\frac {1}{2}}+H_{\nu }(a,a)+C_{\nu }(a,b)&a=b,\\1+H_{\nu }(a,b)+C_{\nu }(a,b)&a>b,\end{array}}\right.}
ここで、 は前に定義したとおりであり 、追加の補正項は次のように与えられる。
H
n
(
a
,
b
)
{\displaystyle H_{n}(a,b)}
ζ
=
a
/
b
{\displaystyle \zeta =a/b}
C
ν
(
a
,
b
)
=
sin
(
ν
π
)
π
exp
(
−
a
2
+
b
2
2
)
∫
0
1
(
x
/
ζ
)
ν
−
1
ζ
+
x
exp
[
−
a
b
2
(
x
+
1
x
)
]
d
x
.
{\displaystyle C_{\nu }(a,b)={\frac {\sin(\nu \pi )}{\pi }}\exp \left(-{\frac {a^{2}+b^{2}}{2}}\right)\int _{0}^{1}{\frac {(x/\zeta )^{\nu -1}}{\zeta +x}}\exp \left[-{\frac {ab}{2}}\left(x+{\frac {1}{x}}\right)\right]\mathrm {d} x.}
の整数値の場合 、補正項は 消える傾向があります。
ν
{\displaystyle \nu }
C
ν
(
a
,
b
)
{\displaystyle C_{\nu }(a,b)}
単調性と対数凹性
一般化されたマーカムQ関数は、 すべての およびに対して および で厳密に増加し 、 すべての およびに対して で厳密に減少します [5]
Q
ν
(
a
,
b
)
{\displaystyle Q_{\nu }(a,b)}
ν
{\displaystyle \nu }
a
{\displaystyle a}
a
≥
0
{\displaystyle a\geq 0}
b
,
ν
>
0
{\displaystyle b,\nu >0}
b
{\displaystyle b}
a
,
b
≥
0
{\displaystyle a,b\geq 0}
ν
>
0.
{\displaystyle \nu >0.}
この関数は、すべての [5] に対して 対数凹関数 で ある。
ν
↦
Q
ν
(
a
,
b
)
{\displaystyle \nu \mapsto Q_{\nu }(a,b)}
[
1
,
∞
)
{\displaystyle [1,\infty )}
a
,
b
≥
0.
{\displaystyle a,b\geq 0.}
関数は、 すべての および に対して厳密に対数凹であり 、これは一般化されたMarcum Q関数が新しいものが中古品よりも優れているという性質を満たすことを意味します。 [6]
b
↦
Q
ν
(
a
,
b
)
{\displaystyle b\mapsto Q_{\nu }(a,b)}
(
0
,
∞
)
{\displaystyle (0,\infty )}
a
≥
0
{\displaystyle a\geq 0}
ν
>
1
{\displaystyle \nu >1}
この関数は、すべての [5] に対して 対数凹関数である。
a
↦
1
−
Q
ν
(
a
,
b
)
{\displaystyle a\mapsto 1-Q_{\nu }(a,b)}
[
0
,
∞
)
{\displaystyle [0,\infty )}
b
,
ν
>
0.
{\displaystyle b,\nu >0.}
シリーズ表現
一般化されたマルカムQ関数は 不完全ガンマ関数を使って次のように表される [7] [8] [9]
ν
>
0
{\displaystyle \nu >0}
Q
ν
(
a
,
b
)
=
1
−
e
−
a
2
/
2
∑
k
=
0
∞
1
k
!
γ
(
ν
+
k
,
b
2
2
)
Γ
(
ν
+
k
)
(
a
2
2
)
k
,
{\displaystyle Q_{\nu }(a,b)=1-e^{-a^{2}/2}\sum _{k=0}^{\infty }{\frac {1}{k!}}{\frac {\gamma (\nu +k,{\frac {b^{2}}{2}})}{\Gamma (\nu +k)}}\left({\frac {a^{2}}{2}}\right)^{k},}
ここで、 は 下側不完全ガンマ関数 です。これは通常、 次の一般化 Marcum Q 関数 の標準表現と呼ばれます。
γ
(
s
,
x
)
{\displaystyle \gamma (s,x)}
ν
{\displaystyle \nu }
一般化されたマルカムQ関数は、 一般化されたラゲール多項式を 使用して次のように 表すこともできる [9]
ν
>
0
{\displaystyle \nu >0}
Q
ν
(
a
,
b
)
=
1
−
e
−
a
2
/
2
∑
k
=
0
∞
(
−
1
)
k
L
k
(
ν
−
1
)
(
a
2
2
)
Γ
(
ν
+
k
+
1
)
(
b
2
2
)
k
+
ν
,
{\displaystyle Q_{\nu }(a,b)=1-e^{-a^{2}/2}\sum _{k=0}^{\infty }(-1)^{k}{\frac {L_{k}^{(\nu -1)}({\frac {a^{2}}{2}})}{\Gamma (\nu +k+1)}}\left({\frac {b^{2}}{2}}\right)^{k+\nu },}
ここで、 次 、 位 の一般化ラゲール多項式です 。
L
k
(
α
)
(
⋅
)
{\displaystyle L_{k}^{(\alpha )}(\cdot )}
k
{\displaystyle k}
α
{\displaystyle \alpha }
一般化されたマルカムQ関数は ノイマン級数展開としても表すことができる [4] [8]
ν
>
0
{\displaystyle \nu >0}
Q
ν
(
a
,
b
)
=
e
−
(
a
2
+
b
2
)
/
2
∑
α
=
1
−
ν
∞
(
a
b
)
α
I
−
α
(
a
b
)
,
{\displaystyle Q_{\nu }(a,b)=e^{-(a^{2}+b^{2})/2}\sum _{\alpha =1-\nu }^{\infty }\left({\frac {a}{b}}\right)^{\alpha }I_{-\alpha }(ab),}
1
−
Q
ν
(
a
,
b
)
=
e
−
(
a
2
+
b
2
)
/
2
∑
α
=
ν
∞
(
b
a
)
α
I
α
(
a
b
)
,
{\displaystyle 1-Q_{\nu }(a,b)=e^{-(a^{2}+b^{2})/2}\sum _{\alpha =\nu }^{\infty }\left({\frac {b}{a}}\right)^{\alpha }I_{\alpha }(ab),}
ここで、合計は 1 ずつ増加します。 が 整数値を想定する場合、 となることに注意してください 。
α
{\displaystyle \alpha }
I
α
(
a
b
)
=
I
−
α
(
a
b
)
{\displaystyle I_{\alpha }(ab)=I_{-\alpha }(ab)}
非負の半整数値に対して 、一般化されたマーカムQ関数の閉じた形式の表現は [8] [10]である。
ν
=
n
+
1
/
2
{\displaystyle \nu =n+1/2}
Q
n
+
1
/
2
(
a
,
b
)
=
1
2
[
e
r
f
c
(
b
−
a
2
)
+
e
r
f
c
(
b
+
a
2
)
]
+
e
−
(
a
2
+
b
2
)
/
2
∑
k
=
1
n
(
b
a
)
k
−
1
/
2
I
k
−
1
/
2
(
a
b
)
,
{\displaystyle Q_{n+1/2}(a,b)={\frac {1}{2}}\left[\mathrm {erfc} \left({\frac {b-a}{\sqrt {2}}}\right)+\mathrm {erfc} \left({\frac {b+a}{\sqrt {2}}}\right)\right]+e^{-(a^{2}+b^{2})/2}\sum _{k=1}^{n}\left({\frac {b}{a}}\right)^{k-1/2}I_{k-1/2}(ab),}
ここで は 相補誤差関数 である 。半整数パラメータを持つベッセル関数は有限和展開を持つため、 [4]
e
r
f
c
(
⋅
)
{\displaystyle \mathrm {erfc} (\cdot )}
I
±
(
n
+
0.5
)
(
z
)
=
1
π
∑
k
=
0
n
(
n
+
k
)
!
k
!
(
n
−
k
)
!
[
(
−
1
)
k
e
z
∓
(
−
1
)
n
e
−
z
(
2
z
)
k
+
0.5
]
,
{\displaystyle I_{\pm (n+0.5)}(z)={\frac {1}{\sqrt {\pi }}}\sum _{k=0}^{n}{\frac {(n+k)!}{k!(n-k)!}}\left[{\frac {(-1)^{k}e^{z}\mp (-1)^{n}e^{-z}}{(2z)^{k+0.5}}}\right],}
ここで、は 非負の整数であり、半整数パラメータを持つ一般化マーカムQ関数を正確に表すことができる。より正確には、 [4]
n
{\displaystyle n}
Q
n
+
1
/
2
(
a
,
b
)
=
Q
(
b
−
a
)
+
Q
(
b
+
a
)
+
1
b
2
π
∑
i
=
1
n
(
b
a
)
i
∑
k
=
0
i
−
1
(
i
+
k
−
1
)
!
k
!
(
i
−
k
−
1
)
!
[
(
−
1
)
k
e
−
(
a
−
b
)
2
/
2
+
(
−
1
)
i
e
−
(
a
+
b
)
2
/
2
(
2
a
b
)
k
]
,
{\displaystyle Q_{n+1/2}(a,b)=Q(b-a)+Q(b+a)+{\frac {1}{b{\sqrt {2\pi }}}}\sum _{i=1}^{n}\left({\frac {b}{a}}\right)^{i}\sum _{k=0}^{i-1}{\frac {(i+k-1)!}{k!(i-k-1)!}}\left[{\frac {(-1)^{k}e^{-(a-b)^{2}/2}+(-1)^{i}e^{-(a+b)^{2}/2}}{(2ab)^{k}}}\right],}
非負整数 に対して 、 は ガウスQ関数 です 。あるいは、半整数のベッセル関数を双曲正弦関数と双曲余弦関数の和としてより簡潔に表現することもできます。 [11]
n
{\displaystyle n}
Q
(
⋅
)
{\displaystyle Q(\cdot )}
I
n
+
1
2
(
z
)
=
2
z
π
[
g
n
(
z
)
sinh
(
z
)
+
g
−
n
−
1
(
z
)
cosh
(
z
)
]
,
{\displaystyle I_{n+{\frac {1}{2}}}(z)={\sqrt {\frac {2z}{\pi }}}\left[g_{n}(z)\sinh(z)+g_{-n-1}(z)\cosh(z)\right],}
ここで 、 の任意の整数値に対して 、 、 となります 。
g
0
(
z
)
=
z
−
1
{\displaystyle g_{0}(z)=z^{-1}}
g
1
(
z
)
=
−
z
−
2
{\displaystyle g_{1}(z)=-z^{-2}}
g
n
−
1
(
z
)
−
g
n
+
1
(
z
)
=
(
2
n
+
1
)
z
−
1
g
n
(
z
)
{\displaystyle g_{n-1}(z)-g_{n+1}(z)=(2n+1)z^{-1}g_{n}(z)}
n
{\displaystyle n}
再帰関係と生成関数
部分積分すると、一般化マーカムQ関数は次の再帰関係を満たすことが分かる [8] [10]
Q
ν
+
1
(
a
,
b
)
−
Q
ν
(
a
,
b
)
=
(
b
a
)
ν
e
−
(
a
2
+
b
2
)
/
2
I
ν
(
a
b
)
.
{\displaystyle Q_{\nu +1}(a,b)-Q_{\nu }(a,b)=\left({\frac {b}{a}}\right)^{\nu }e^{-(a^{2}+b^{2})/2}I_{\nu }(ab).}
Q
ν
−
n
(
a
,
b
)
=
Q
ν
(
a
,
b
)
−
(
b
a
)
ν
e
−
(
a
2
+
b
2
)
/
2
∑
k
=
1
n
(
a
b
)
k
I
ν
−
k
(
a
b
)
,
{\displaystyle Q_{\nu -n}(a,b)=Q_{\nu }(a,b)-\left({\frac {b}{a}}\right)^{\nu }e^{-(a^{2}+b^{2})/2}\sum _{k=1}^{n}\left({\frac {a}{b}}\right)^{k}I_{\nu -k}(ab),}
Q
ν
+
n
(
a
,
b
)
=
Q
ν
(
a
,
b
)
+
(
b
a
)
ν
e
−
(
a
2
+
b
2
)
/
2
∑
k
=
0
n
−
1
(
b
a
)
k
I
ν
+
k
(
a
b
)
,
{\displaystyle Q_{\nu +n}(a,b)=Q_{\nu }(a,b)+\left({\frac {b}{a}}\right)^{\nu }e^{-(a^{2}+b^{2})/2}\sum _{k=0}^{n-1}\left({\frac {b}{a}}\right)^{k}I_{\nu +k}(ab),}
正の整数 に対して、 と なります。前者の再帰性は、負の に対して、一般化 Marcum Q 関数を正式に定義するために使用できます 。 に対して およびを取ると 、一般化 Marcum Q 関数のノイマン級数表現が得られます。
n
{\displaystyle n}
ν
{\displaystyle \nu }
Q
∞
(
a
,
b
)
=
1
{\displaystyle Q_{\infty }(a,b)=1}
Q
−
∞
(
a
,
b
)
=
0
{\displaystyle Q_{-\infty }(a,b)=0}
n
=
∞
{\displaystyle n=\infty }
Q
ν
+
1
(
a
,
b
)
−
(
1
+
c
ν
(
a
,
b
)
)
Q
ν
(
a
,
b
)
+
c
ν
(
a
,
b
)
Q
ν
−
1
(
a
,
b
)
=
0
,
{\displaystyle Q_{\nu +1}(a,b)-(1+c_{\nu }(a,b))Q_{\nu }(a,b)+c_{\nu }(a,b)Q_{\nu -1}(a,b)=0,}
どこ
c
ν
(
a
,
b
)
=
(
b
a
)
I
ν
(
a
b
)
I
ν
+
1
(
a
b
)
.
{\displaystyle c_{\nu }(a,b)=\left({\frac {b}{a}}\right){\frac {I_{\nu }(ab)}{I_{\nu +1}(ab)}}.}
ベッセル関数の出現を消去すると、3次の再帰関係が得られる [7]。
a
2
2
Q
ν
+
2
(
a
,
b
)
=
(
a
2
2
−
ν
)
Q
ν
+
1
(
a
,
b
)
+
(
b
2
2
+
ν
)
Q
ν
(
a
,
b
)
−
b
2
2
Q
ν
−
1
(
a
,
b
)
.
{\displaystyle {\frac {a^{2}}{2}}Q_{\nu +2}(a,b)=\left({\frac {a^{2}}{2}}-\nu \right)Q_{\nu +1}(a,b)+\left({\frac {b^{2}}{2}}+\nu \right)Q_{\nu }(a,b)-{\frac {b^{2}}{2}}Q_{\nu -1}(a,b).}
導関数と関連したもう一つの再帰関係は次のように表される。
Q
ν
+
1
(
a
,
b
)
=
Q
ν
(
a
,
b
)
+
1
a
∂
∂
a
Q
ν
(
a
,
b
)
,
{\displaystyle Q_{\nu +1}(a,b)=Q_{\nu }(a,b)+{\frac {1}{a}}{\frac {\partial }{\partial a}}Q_{\nu }(a,b),}
Q
ν
−
1
(
a
,
b
)
=
Q
ν
(
a
,
b
)
+
1
b
∂
∂
b
Q
ν
(
a
,
b
)
.
{\displaystyle Q_{\nu -1}(a,b)=Q_{\nu }(a,b)+{\frac {1}{b}}{\frac {\partial }{\partial b}}Q_{\nu }(a,b).}
積分に対する 通常の生成関数は [10] である。
Q
ν
(
a
,
b
)
{\displaystyle Q_{\nu }(a,b)}
ν
{\displaystyle \nu }
∑
n
=
−
∞
∞
t
n
Q
n
(
a
,
b
)
=
e
−
(
a
2
+
b
2
)
/
2
t
1
−
t
e
(
b
2
t
+
a
2
/
t
)
/
2
,
{\displaystyle \sum _{n=-\infty }^{\infty }t^{n}Q_{n}(a,b)=e^{-(a^{2}+b^{2})/2}{\frac {t}{1-t}}e^{(b^{2}t+a^{2}/t)/2},}
どこ
|
t
|
<
1.
{\displaystyle |t|<1.}
対称関係
2つのノイマン級数表現を用いると、正の積分に対して次の対称関係が得られる。
ν
=
n
{\displaystyle \nu =n}
Q
n
(
a
,
b
)
+
Q
n
(
b
,
a
)
=
1
+
e
−
(
a
2
+
b
2
)
/
2
[
I
0
(
a
b
)
+
∑
k
=
1
n
−
1
a
2
k
+
b
2
k
(
a
b
)
k
I
k
(
a
b
)
]
.
{\displaystyle Q_{n}(a,b)+Q_{n}(b,a)=1+e^{-(a^{2}+b^{2})/2}\left[I_{0}(ab)+\sum _{k=1}^{n-1}{\frac {a^{2k}+b^{2k}}{(ab)^{k}}}I_{k}(ab)\right].}
特に、 私たちは
n
=
1
{\displaystyle n=1}
Q
1
(
a
,
b
)
+
Q
1
(
b
,
a
)
=
1
+
e
−
(
a
2
+
b
2
)
/
2
I
0
(
a
b
)
.
{\displaystyle Q_{1}(a,b)+Q_{1}(b,a)=1+e^{-(a^{2}+b^{2})/2}I_{0}(ab).}
特別な値
マルカムQ関数の具体的な値は [6]
Q
ν
(
0
,
0
)
=
1
,
{\displaystyle Q_{\nu }(0,0)=1,}
Q
ν
(
a
,
0
)
=
1
,
{\displaystyle Q_{\nu }(a,0)=1,}
Q
ν
(
a
,
+
∞
)
=
0
,
{\displaystyle Q_{\nu }(a,+\infty )=0,}
Q
ν
(
0
,
b
)
=
Γ
(
ν
,
b
2
/
2
)
Γ
(
ν
)
,
{\displaystyle Q_{\nu }(0,b)={\frac {\Gamma (\nu ,b^{2}/2)}{\Gamma (\nu )}},}
Q
ν
(
+
∞
,
b
)
=
1
,
{\displaystyle Q_{\nu }(+\infty ,b)=1,}
Q
∞
(
a
,
b
)
=
1
,
{\displaystyle Q_{\infty }(a,b)=1,}
については 、ノイマン級数表現の2つの形式を減算することにより、 [10]
a
=
b
{\displaystyle a=b}
Q
1
(
a
,
a
)
=
1
2
[
1
+
e
−
a
2
I
0
(
a
2
)
]
,
{\displaystyle Q_{1}(a,a)={\frac {1}{2}}[1+e^{-a^{2}}I_{0}(a^{2})],}
これを再帰式と組み合わせると、
Q
n
(
a
,
a
)
=
1
2
[
1
+
e
−
a
2
I
0
(
a
2
)
]
+
e
−
a
2
∑
k
=
1
n
−
1
I
k
(
a
2
)
,
{\displaystyle Q_{n}(a,a)={\frac {1}{2}}[1+e^{-a^{2}}I_{0}(a^{2})]+e^{-a^{2}}\sum _{k=1}^{n-1}I_{k}(a^{2}),}
Q
−
n
(
a
,
a
)
=
1
2
[
1
+
e
−
a
2
I
0
(
a
2
)
]
−
e
−
a
2
∑
k
=
1
n
I
k
(
a
2
)
,
{\displaystyle Q_{-n}(a,a)={\frac {1}{2}}[1+e^{-a^{2}}I_{0}(a^{2})]-e^{-a^{2}}\sum _{k=1}^{n}I_{k}(a^{2}),}
任意の非負整数に対して 。
n
{\displaystyle n}
については 、一般化マーカムQ関数の基本的な積分定義を用いると、 [8] [10]
ν
=
1
/
2
{\displaystyle \nu =1/2}
Q
1
/
2
(
a
,
b
)
=
1
2
[
e
r
f
c
(
b
−
a
2
)
+
e
r
f
c
(
b
+
a
2
)
]
.
{\displaystyle Q_{1/2}(a,b)={\frac {1}{2}}\left[\mathrm {erfc} \left({\frac {b-a}{\sqrt {2}}}\right)+\mathrm {erfc} \left({\frac {b+a}{\sqrt {2}}}\right)\right].}
については 、
ν
=
3
/
2
{\displaystyle \nu =3/2}
Q
3
/
2
(
a
,
b
)
=
Q
1
/
2
(
a
,
b
)
+
2
π
sinh
(
a
b
)
a
e
−
(
a
2
+
b
2
)
/
2
.
{\displaystyle Q_{3/2}(a,b)=Q_{1/2}(a,b)+{\sqrt {\frac {2}{\pi }}}\,{\frac {\sinh(ab)}{a}}e^{-(a^{2}+b^{2})/2}.}
なぜなら 私たちは
ν
=
5
/
2
{\displaystyle \nu =5/2}
Q
5
/
2
(
a
,
b
)
=
Q
3
/
2
(
a
,
b
)
+
2
π
a
b
cosh
(
a
b
)
−
sinh
(
a
b
)
a
3
e
−
(
a
2
+
b
2
)
/
2
.
{\displaystyle Q_{5/2}(a,b)=Q_{3/2}(a,b)+{\sqrt {\frac {2}{\pi }}}\,{\frac {ab\cosh(ab)-\sinh(ab)}{a^{3}}}e^{-(a^{2}+b^{2})/2}.}
が固定され 大きいと 仮定すると 、一般化されたマルカムQ関数は次の漸近形を持つ [7]
ν
{\displaystyle \nu }
a
b
{\displaystyle ab}
ζ
=
a
/
b
>
0
{\displaystyle \zeta =a/b>0}
Q
ν
(
a
,
b
)
∼
∑
n
=
0
∞
ψ
n
,
{\displaystyle Q_{\nu }(a,b)\sim \sum _{n=0}^{\infty }\psi _{n},}
ここで 、
ψ
n
{\displaystyle \psi _{n}}
ψ
n
=
1
2
ζ
ν
2
π
(
−
1
)
n
[
A
n
(
ν
−
1
)
−
ζ
A
n
(
ν
)
]
ϕ
n
.
{\displaystyle \psi _{n}={\frac {1}{2\zeta ^{\nu }{\sqrt {2\pi }}}}(-1)^{n}\left[A_{n}(\nu -1)-\zeta A_{n}(\nu )\right]\phi _{n}.}
関数 と関数 は次のように与えられる。
ϕ
n
{\displaystyle \phi _{n}}
A
n
{\displaystyle A_{n}}
ϕ
n
=
[
(
b
−
a
)
2
2
a
b
]
n
−
1
2
Γ
(
1
2
−
n
,
(
b
−
a
)
2
2
)
,
{\displaystyle \phi _{n}=\left[{\frac {(b-a)^{2}}{2ab}}\right]^{n-{\frac {1}{2}}}\Gamma \left({\frac {1}{2}}-n,{\frac {(b-a)^{2}}{2}}\right),}
A
n
(
ν
)
=
2
−
n
Γ
(
1
2
+
ν
+
n
)
n
!
Γ
(
1
2
+
ν
−
n
)
.
{\displaystyle A_{n}(\nu )={\frac {2^{-n}\Gamma ({\frac {1}{2}}+\nu +n)}{n!\Gamma ({\frac {1}{2}}+\nu -n)}}.}
この関数は 再帰を満たす
A
n
(
ν
)
{\displaystyle A_{n}(\nu )}
A
n
+
1
(
ν
)
=
−
(
2
n
+
1
)
2
−
4
ν
2
8
(
n
+
1
)
A
n
(
ν
)
,
{\displaystyle A_{n+1}(\nu )=-{\frac {(2n+1)^{2}-4\nu ^{2}}{8(n+1)}}A_{n}(\nu ),}
および
n
≥
0
{\displaystyle n\geq 0}
A
0
(
ν
)
=
1.
{\displaystyle A_{0}(\nu )=1.}
ϕ
0
=
2
π
a
b
b
−
a
e
r
f
c
(
b
−
a
2
)
.
{\displaystyle \phi _{0}={\frac {\sqrt {2\pi ab}}{b-a}}\mathrm {erfc} \left({\frac {b-a}{\sqrt {2}}}\right).}
したがって、と仮定すると 、一般化マルカムQ関数の第一項の漸近近似は [7]である。
b
>
a
{\displaystyle b>a}
Q
ν
(
a
,
b
)
∼
ψ
0
=
(
b
a
)
ν
−
1
2
Q
(
b
−
a
)
,
{\displaystyle Q_{\nu }(a,b)\sim \psi _{0}=\left({\frac {b}{a}}\right)^{\nu -{\frac {1}{2}}}Q(b-a),}
ここで は ガウスQ関数 です 。ここ で
Q
(
⋅
)
{\displaystyle Q(\cdot )}
Q
ν
(
a
,
b
)
∼
0.5
{\displaystyle Q_{\nu }(a,b)\sim 0.5}
a
↑
b
.
{\displaystyle a\uparrow b.}
の場合には 、 [7]
a
>
b
{\displaystyle a>b}
Q
ν
(
a
,
b
)
∼
1
−
ψ
0
=
1
−
(
b
a
)
ν
−
1
2
Q
(
a
−
b
)
.
{\displaystyle Q_{\nu }(a,b)\sim 1-\psi _{0}=1-\left({\frac {b}{a}}\right)^{\nu -{\frac {1}{2}}}Q(a-b).}
ここ でも
Q
ν
(
a
,
b
)
∼
0.5
{\displaystyle Q_{\nu }(a,b)\sim 0.5}
a
↓
b
.
{\displaystyle a\downarrow b.}
差別化
の偏微分は 次 のよう に与えられる [12] [13]
Q
ν
(
a
,
b
)
{\displaystyle Q_{\nu }(a,b)}
a
{\displaystyle a}
b
{\displaystyle b}
∂
∂
a
Q
ν
(
a
,
b
)
=
a
[
Q
ν
+
1
(
a
,
b
)
−
Q
ν
(
a
,
b
)
]
=
a
(
b
a
)
ν
e
−
(
a
2
+
b
2
)
/
2
I
ν
(
a
b
)
,
{\displaystyle {\frac {\partial }{\partial a}}Q_{\nu }(a,b)=a\left[Q_{\nu +1}(a,b)-Q_{\nu }(a,b)\right]=a\left({\frac {b}{a}}\right)^{\nu }e^{-(a^{2}+b^{2})/2}I_{\nu }(ab),}
∂
∂
b
Q
ν
(
a
,
b
)
=
b
[
Q
ν
−
1
(
a
,
b
)
−
Q
ν
(
a
,
b
)
]
=
−
b
(
b
a
)
ν
−
1
e
−
(
a
2
+
b
2
)
/
2
I
ν
−
1
(
a
b
)
.
{\displaystyle {\frac {\partial }{\partial b}}Q_{\nu }(a,b)=b\left[Q_{\nu -1}(a,b)-Q_{\nu }(a,b)\right]=-b\left({\frac {b}{a}}\right)^{\nu -1}e^{-(a^{2}+b^{2})/2}I_{\nu -1}(ab).}
2つの偏導関数は次のように関係づけられる。
1
a
∂
∂
a
Q
ν
(
a
,
b
)
+
1
b
∂
∂
b
Q
ν
+
1
(
a
,
b
)
=
0.
{\displaystyle {\frac {1}{a}}{\frac {\partial }{\partial a}}Q_{\nu }(a,b)+{\frac {1}{b}}{\frac {\partial }{\partial b}}Q_{\nu +1}(a,b)=0.}
のn 次偏 微分は 次のように与えられる [10]
Q
ν
(
a
,
b
)
{\displaystyle Q_{\nu }(a,b)}
∂
n
∂
a
n
Q
ν
(
a
,
b
)
=
n
!
(
−
a
)
n
∑
k
=
0
[
n
/
2
]
(
−
2
a
2
)
−
k
k
!
(
n
−
2
k
)
!
∑
p
=
0
n
−
k
(
−
1
)
p
(
n
−
k
p
)
Q
ν
+
p
(
a
,
b
)
,
{\displaystyle {\frac {\partial ^{n}}{\partial a^{n}}}Q_{\nu }(a,b)=n!(-a)^{n}\sum _{k=0}^{[n/2]}{\frac {(-2a^{2})^{-k}}{k!(n-2k)!}}\sum _{p=0}^{n-k}(-1)^{p}{\binom {n-k}{p}}Q_{\nu +p}(a,b),}
∂
n
∂
b
n
Q
ν
(
a
,
b
)
=
n
!
a
1
−
ν
2
n
b
n
−
ν
+
1
e
−
(
a
2
+
b
2
)
/
2
∑
k
=
[
n
/
2
]
n
(
−
2
b
2
)
k
(
n
−
k
)
!
(
2
k
−
n
)
!
∑
p
=
0
k
−
1
(
k
−
1
p
)
(
−
a
b
)
p
I
ν
−
p
−
1
(
a
b
)
.
{\displaystyle {\frac {\partial ^{n}}{\partial b^{n}}}Q_{\nu }(a,b)={\frac {n!a^{1-\nu }}{2^{n}b^{n-\nu +1}}}e^{-(a^{2}+b^{2})/2}\sum _{k=[n/2]}^{n}{\frac {(-2b^{2})^{k}}{(n-k)!(2k-n)!}}\sum _{p=0}^{k-1}{\binom {k-1}{p}}\left(-{\frac {a}{b}}\right)^{p}I_{\nu -p-1}(ab).}
不平等
Q
ν
2
(
a
,
b
)
>
Q
ν
−
1
(
a
,
b
)
+
Q
ν
+
1
(
a
,
b
)
2
>
Q
ν
−
1
(
a
,
b
)
Q
ν
+
1
(
a
,
b
)
{\displaystyle Q_{\nu }^{2}(a,b)>{\frac {Q_{\nu -1}(a,b)+Q_{\nu +1}(a,b)}{2}}>Q_{\nu -1}(a,b)Q_{\nu +1}(a,b)}
すべて およびについて 。
a
≥
b
>
0
{\displaystyle a\geq b>0}
ν
>
1
{\displaystyle \nu >1}
境界
単調性と対数凹性に基づく
一般化 Marcum-Q 関数のさまざまな上限と下限は、関数の単調性と対数凹性、および が半整数値の
場合 に の閉じた形式表現があることを利用して取得できます。
ν
↦
Q
ν
(
a
,
b
)
{\displaystyle \nu \mapsto Q_{\nu }(a,b)}
Q
ν
(
a
,
b
)
{\displaystyle Q_{\nu }(a,b)}
ν
{\displaystyle \nu }
とを、 それぞれ実数を最も近い左半奇数と右半奇数に
マッピングする半整数丸め演算子のペアとします 。
⌊
x
⌋
0.5
{\displaystyle \lfloor x\rfloor _{0.5}}
⌈
x
⌉
0.5
{\displaystyle \lceil x\rceil _{0.5}}
x
{\displaystyle x}
⌊
x
⌋
0.5
=
⌊
x
−
0.5
⌋
+
0.5
{\displaystyle \lfloor x\rfloor _{0.5}=\lfloor x-0.5\rfloor +0.5}
⌈
x
⌉
0.5
=
⌈
x
+
0.5
⌉
−
0.5
{\displaystyle \lceil x\rceil _{0.5}=\lceil x+0.5\rceil -0.5}
ここで 、およびは 整数の床関数と天井関数を表します。
⌊
x
⌋
{\displaystyle \lfloor x\rfloor }
⌈
x
⌉
{\displaystyle \lceil x\rceil }
全ての 関数の単調性は 、 次のような単純な境界を与える [14] [8] [15]
ν
↦
Q
ν
(
a
,
b
)
{\displaystyle \nu \mapsto Q_{\nu }(a,b)}
a
≥
0
{\displaystyle a\geq 0}
b
>
0
{\displaystyle b>0}
Q
⌊
ν
⌋
0.5
(
a
,
b
)
<
Q
ν
(
a
,
b
)
<
Q
⌈
ν
⌉
0.5
(
a
,
b
)
.
{\displaystyle Q_{\lfloor \nu \rfloor _{0.5}}(a,b)<Q_{\nu }(a,b)<Q_{\lceil \nu \rceil _{0.5}}(a,b).}
しかし、この境界の相対誤差は の場合にはゼロに近づかない 。 [5] の整数値の場合 、この境界は次のように減少する。
b
→
∞
{\displaystyle b\to \infty }
ν
=
n
{\displaystyle \nu =n}
Q
n
−
0.5
(
a
,
b
)
<
Q
n
(
a
,
b
)
<
Q
n
+
0.5
(
a
,
b
)
.
{\displaystyle Q_{n-0.5}(a,b)<Q_{n}(a,b)<Q_{n+0.5}(a,b).}
整数値に対する一般化マルカムQ関数の非常に良い近似は、 上限と下限の算術平均を取ることによって得られる [15]
ν
=
n
{\displaystyle \nu =n}
Q
n
(
a
,
b
)
≈
Q
n
−
0.5
(
a
,
b
)
+
Q
n
+
0.5
(
a
,
b
)
2
.
{\displaystyle Q_{n}(a,b)\approx {\frac {Q_{n-0.5}(a,b)+Q_{n+0.5}(a,b)}{2}}.}
の対数凹性を利用することで、より厳密な境界を得ることができる 。 [ 5]
ν
↦
Q
ν
(
a
,
b
)
{\displaystyle \nu \mapsto Q_{\nu }(a,b)}
[
1
,
∞
)
{\displaystyle [1,\infty )}
Q
ν
1
(
a
,
b
)
ν
2
−
v
Q
ν
2
(
a
,
b
)
v
−
ν
1
<
Q
ν
(
a
,
b
)
<
Q
ν
2
(
a
,
b
)
ν
2
−
ν
+
1
Q
ν
2
+
1
(
a
,
b
)
ν
2
−
ν
,
{\displaystyle Q_{\nu _{1}}(a,b)^{\nu _{2}-v}Q_{\nu _{2}}(a,b)^{v-\nu _{1}}<Q_{\nu }(a,b)<{\frac {Q_{\nu _{2}}(a,b)^{\nu _{2}-\nu +1}}{Q_{\nu _{2}+1}(a,b)^{\nu _{2}-\nu }}},}
ここで、 については となる 。この境界の厳密さは 、 または が増加するにつれて向上する。この境界の相対誤差は のときに 0 に収束する 。 [5] の整数値については 、この境界は次のように減少する。
ν
1
=
⌊
ν
⌋
0.5
{\displaystyle \nu _{1}=\lfloor \nu \rfloor _{0.5}}
ν
2
=
⌈
ν
⌉
0.5
{\displaystyle \nu _{2}=\lceil \nu \rceil _{0.5}}
ν
≥
1.5
{\displaystyle \nu \geq 1.5}
a
{\displaystyle a}
ν
{\displaystyle \nu }
b
→
∞
{\displaystyle b\to \infty }
ν
=
n
{\displaystyle \nu =n}
Q
n
−
0.5
(
a
,
b
)
Q
n
+
0.5
(
a
,
b
)
<
Q
n
(
a
,
b
)
<
Q
n
+
0.5
(
a
,
b
)
Q
n
+
0.5
(
a
,
b
)
Q
n
+
1.5
(
a
,
b
)
.
{\displaystyle {\sqrt {Q_{n-0.5}(a,b)Q_{n+0.5}(a,b)}}<Q_{n}(a,b)<Q_{n+0.5}(a,b){\sqrt {\frac {Q_{n+0.5}(a,b)}{Q_{n+1.5}(a,b)}}}.}
コーシー・シュワルツ境界
整数値に対する三角関数の積分表現を用いると 、次のコーシー・シュワルツ境界が得られる [3]。
ν
=
n
{\displaystyle \nu =n}
e
−
b
2
/
2
≤
Q
n
(
a
,
b
)
≤
exp
[
−
1
2
(
b
2
+
a
2
)
]
I
0
(
2
a
b
)
2
n
−
1
2
+
ζ
2
(
1
−
n
)
2
(
1
−
ζ
2
)
,
ζ
<
1
,
{\displaystyle e^{-b^{2}/2}\leq Q_{n}(a,b)\leq \exp \left[-{\frac {1}{2}}(b^{2}+a^{2})\right]{\sqrt {I_{0}(2ab)}}{\sqrt {{\frac {2n-1}{2}}+{\frac {\zeta ^{2(1-n)}}{2(1-\zeta ^{2})}}}},\qquad \zeta <1,}
1
−
Q
n
(
a
,
b
)
≤
exp
[
−
1
2
(
b
2
+
a
2
)
]
I
0
(
2
a
b
)
ζ
2
(
1
−
n
)
2
(
ζ
2
−
1
)
,
ζ
>
1
,
{\displaystyle 1-Q_{n}(a,b)\leq \exp \left[-{\frac {1}{2}}(b^{2}+a^{2})\right]{\sqrt {I_{0}(2ab)}}{\sqrt {\frac {\zeta ^{2(1-n)}}{2(\zeta ^{2}-1)}}},\qquad \zeta >1,}
どこ 。
ζ
=
a
/
b
>
0
{\displaystyle \zeta =a/b>0}
指数型境界
解析的な目的のためには、単純な指数形式で境界値を持つことがしばしば有用であるが、それが達成可能な最も厳しい境界値ではないかもしれない。 とすると 、整数値 に対するそのような境界値の1つは 次のように与えられる [16] [3]
ζ
=
a
/
b
>
0
{\displaystyle \zeta =a/b>0}
ν
=
n
{\displaystyle \nu =n}
e
−
(
b
+
a
)
2
/
2
≤
Q
n
(
a
,
b
)
≤
e
−
(
b
−
a
)
2
/
2
+
ζ
1
−
n
−
1
π
(
1
−
ζ
)
[
e
−
(
b
−
a
)
2
/
2
−
e
−
(
b
+
a
)
2
/
2
]
,
ζ
<
1
,
{\displaystyle e^{-(b+a)^{2}/2}\leq Q_{n}(a,b)\leq e^{-(b-a)^{2}/2}+{\frac {\zeta ^{1-n}-1}{\pi (1-\zeta )}}\left[e^{-(b-a)^{2}/2}-e^{-(b+a)^{2}/2}\right],\qquad \zeta <1,}
Q
n
(
a
,
b
)
≥
1
−
1
2
[
e
−
(
a
−
b
)
2
/
2
−
e
−
(
a
+
b
)
2
/
2
]
,
ζ
>
1.
{\displaystyle Q_{n}(a,b)\geq 1-{\frac {1}{2}}\left[e^{-(a-b)^{2}/2}-e^{-(a+b)^{2}/2}\right],\qquad \zeta >1.}
のとき 、境界は次のように簡略化される。
n
=
1
{\displaystyle n=1}
e
−
(
b
+
a
)
2
/
2
≤
Q
1
(
a
,
b
)
≤
e
−
(
b
−
a
)
2
/
2
,
ζ
<
1
,
{\displaystyle e^{-(b+a)^{2}/2}\leq Q_{1}(a,b)\leq e^{-(b-a)^{2}/2},\qquad \zeta <1,}
1
−
1
2
[
e
−
(
a
−
b
)
2
/
2
−
e
−
(
a
+
b
)
2
/
2
]
≤
Q
1
(
a
,
b
)
,
ζ
>
1.
{\displaystyle 1-{\frac {1}{2}}\left[e^{-(a-b)^{2}/2}-e^{-(a+b)^{2}/2}\right]\leq Q_{1}(a,b),\qquad \zeta >1.}
コーシー・シュワルツ不等式によって得られる別の境界は次のように与えられる [3]。
e
−
b
2
/
2
≤
Q
n
(
a
,
b
)
≤
1
2
2
n
−
1
2
+
ζ
2
(
1
−
n
)
2
(
1
−
ζ
2
)
[
e
−
(
b
−
a
)
2
/
2
+
e
−
(
b
+
a
)
2
/
2
]
,
ζ
<
1
{\displaystyle e^{-b^{2}/2}\leq Q_{n}(a,b)\leq {\frac {1}{2}}{\sqrt {{\frac {2n-1}{2}}+{\frac {\zeta ^{2(1-n)}}{2(1-\zeta ^{2})}}}}\left[e^{-(b-a)^{2}/2}+e^{-(b+a)^{2}/2}\right],\qquad \zeta <1}
Q
n
(
a
,
b
)
≥
1
−
1
2
ζ
2
(
1
−
n
)
2
(
ζ
2
−
1
)
[
e
−
(
b
−
a
)
2
/
2
+
e
−
(
b
+
a
)
2
/
2
]
,
ζ
>
1.
{\displaystyle Q_{n}(a,b)\geq 1-{\frac {1}{2}}{\sqrt {\frac {\zeta ^{2(1-n)}}{2(\zeta ^{2}-1)}}}\left[e^{-(b-a)^{2}/2}+e^{-(b+a)^{2}/2}\right],\qquad \zeta >1.}
チェルノフ型境界
一般化マルカムQ関数のチェルノフ型境界は 整数で、次のように与えられる [16] [3]
ν
=
n
{\displaystyle \nu =n}
(
1
−
2
λ
)
−
n
exp
(
−
λ
b
2
+
λ
n
a
2
1
−
2
λ
)
≥
{
Q
n
(
a
,
b
)
,
b
2
>
n
(
a
2
+
2
)
1
−
Q
n
(
a
,
b
)
,
b
2
<
n
(
a
2
+
2
)
{\displaystyle (1-2\lambda )^{-n}\exp \left(-\lambda b^{2}+{\frac {\lambda na^{2}}{1-2\lambda }}\right)\geq \left\{{\begin{array}{lr}Q_{n}(a,b),&b^{2}>n(a^{2}+2)\\1-Q_{n}(a,b),&b^{2}<n(a^{2}+2)\end{array}}\right.}
ここで、チェルノフパラメータ の
最適値は
(
0
<
λ
<
1
/
2
)
{\displaystyle (0<\lambda <1/2)}
λ
0
{\displaystyle \lambda _{0}}
λ
0
=
1
2
(
1
−
n
b
2
−
n
b
2
1
+
(
a
b
)
2
n
)
.
{\displaystyle \lambda _{0}={\frac {1}{2}}\left(1-{\frac {n}{b^{2}}}-{\frac {n}{b^{2}}}{\sqrt {1+{\frac {(ab)^{2}}{n}}}}\right).}
半線形近似
1次のマルカムQ関数は[17] で半線形近似できる。
Q
1
(
a
,
b
)
=
{
1
,
i
f
b
<
c
1
−
β
0
e
−
1
2
(
a
2
+
(
β
0
)
2
)
I
0
(
a
β
0
)
(
b
−
β
0
)
+
Q
1
(
a
,
β
0
)
,
i
f
c
1
≤
b
≤
c
2
0
,
i
f
b
>
c
2
{\displaystyle {\begin{aligned}Q_{1}(a,b)={\begin{cases}1,~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\mathrm {if} ~b<c_{1}\\-\beta _{0}e^{-{\frac {1}{2}}\left(a^{2}+\left(\beta _{0}\right)^{2}\right)}I_{0}\left(a\beta _{0}\right)\left(b-\beta _{0}\right)+Q_{1}\left(a,\beta _{0}\right),~~~~~\mathrm {if} ~c_{1}\leq b\leq c_{2}\\0,~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\mathrm {if} ~b>c_{2}\end{cases}}\end{aligned}}}
どこ
β
0
=
a
+
a
2
+
2
2
,
{\displaystyle {\begin{aligned}\beta _{0}={\frac {a+{\sqrt {a^{2}+2}}}{2}},\end{aligned}}}
c
1
(
a
)
=
max
(
0
,
β
0
+
Q
1
(
a
,
β
0
)
−
1
β
0
e
−
1
2
(
a
2
+
(
β
0
)
2
)
I
0
(
a
β
0
)
)
,
{\displaystyle {\begin{aligned}c_{1}(a)=\max {\Bigg (}0,\beta _{0}+{\frac {Q_{1}\left(a,\beta _{0}\right)-1}{\beta _{0}e^{-{\frac {1}{2}}\left(a^{2}+\left(\beta _{0}\right)^{2}\right)}I_{0}\left(a\beta _{0}\right)}}{\Bigg )},\end{aligned}}}
そして
c
2
(
a
)
=
β
0
+
Q
1
(
a
,
β
0
)
β
0
e
−
1
2
(
a
2
+
(
β
0
)
2
)
I
0
(
a
β
0
)
.
{\displaystyle {\begin{aligned}c_{2}(a)=\beta _{0}+{\frac {Q_{1}\left(a,\beta _{0}\right)}{\beta _{0}e^{-{\frac {1}{2}}\left(a^{2}+\left(\beta _{0}\right)^{2}\right)}I_{0}\left(a\beta _{0}\right)}}.\end{aligned}}}
マーカムQ関数は次のように再表現すると便利である [18]
P
N
(
X
,
Y
)
=
Q
N
(
2
N
X
,
2
Y
)
.
{\displaystyle P_{N}(X,Y)=Q_{N}({\sqrt {2NX}},{\sqrt {2Y}}).}
は、正規化された検出閾値 を持つ、 一定の受信信号対雑音比 の非コヒーレントに積分された受信信号サンプル の検出確率として解釈できます 。 この Marcum Q 関数の同等の形式では、 および が与えられ た場合、 および が成り立ちます 。 を表す式は多数存在します 。しかし、数値計算に最も信頼性が高く、正確で、効率的な 5 つの式を以下に示します。それらは形式 1 です。 [18]
P
N
(
X
,
Y
)
{\displaystyle P_{N}(X,Y)}
N
{\displaystyle N}
X
{\displaystyle X}
Y
{\displaystyle Y}
a
{\displaystyle a}
b
{\displaystyle b}
X
=
a
2
/
2
N
{\displaystyle X=a^{2}/2N}
Y
=
b
2
/
2
{\displaystyle Y=b^{2}/2}
P
N
(
X
,
Y
)
{\displaystyle P_{N}(X,Y)}
P
N
(
X
,
Y
)
=
∑
k
=
0
∞
e
−
N
X
(
N
X
)
k
k
!
∑
m
=
0
N
−
1
+
k
e
−
Y
Y
m
m
!
,
{\displaystyle P_{N}(X,Y)=\sum _{k=0}^{\infty }e^{-NX}{\frac {(NX)^{k}}{k!}}\sum _{m=0}^{N-1+k}e^{-Y}{\frac {Y^{m}}{m!}},}
フォーム2: [18]
P
N
(
X
,
Y
)
=
∑
m
=
0
N
−
1
e
−
Y
Y
m
m
!
+
∑
m
=
N
∞
e
−
Y
Y
m
m
!
(
1
−
∑
k
=
0
m
−
N
e
−
N
X
(
N
X
)
k
k
!
)
,
{\displaystyle P_{N}(X,Y)=\sum _{m=0}^{N-1}e^{-Y}{\frac {Y^{m}}{m!}}+\sum _{m=N}^{\infty }e^{-Y}{\frac {Y^{m}}{m!}}\left(1-\sum _{k=0}^{m-N}e^{-NX}{\frac {(NX)^{k}}{k!}}\right),}
フォーム3: [18]
1
−
P
N
(
X
,
Y
)
=
∑
m
=
N
∞
e
−
Y
Y
m
m
!
∑
k
=
0
m
−
N
e
−
N
X
(
N
X
)
k
k
!
,
{\displaystyle 1-P_{N}(X,Y)=\sum _{m=N}^{\infty }e^{-Y}{\frac {Y^{m}}{m!}}\sum _{k=0}^{m-N}e^{-NX}{\frac {(NX)^{k}}{k!}},}
フォーム4: [18]
1
−
P
N
(
X
,
Y
)
=
∑
k
=
0
∞
e
−
N
X
(
N
X
)
k
k
!
(
1
−
∑
m
=
0
N
−
1
+
k
e
−
Y
Y
m
m
!
)
,
{\displaystyle 1-P_{N}(X,Y)=\sum _{k=0}^{\infty }e^{-NX}{\frac {(NX)^{k}}{k!}}\left(1-\sum _{m=0}^{N-1+k}e^{-Y}{\frac {Y^{m}}{m!}}\right),}
そして5つ目は次のようになる: [18]
1
−
P
N
(
X
,
Y
)
=
e
−
(
N
X
+
Y
)
∑
r
=
N
∞
(
Y
N
X
)
r
/
2
I
r
(
2
N
X
Y
)
.
{\displaystyle 1-P_{N}(X,Y)=e^{-(NX+Y)}\sum _{r=N}^{\infty }\left({\frac {Y}{NX}}\right)^{r/2}I_{r}(2{\sqrt {NXY}}).}
これら5つの形態のうち、2番目の形態が最も堅牢である。 [18]
アプリケーション
一般化された Marcum Q 関数は、多くのランダム変数の累積分布関数 (cdf) を表すために使用できます。
が速度パラメータ を持つ 指数分布 である 場合 、その累積分布関数は次のように表される。
X
∼
E
x
p
(
λ
)
{\displaystyle X\sim \mathrm {Exp} (\lambda )}
λ
{\displaystyle \lambda }
F
X
(
x
)
=
1
−
Q
1
(
0
,
2
λ
x
)
{\displaystyle F_{X}(x)=1-Q_{1}\left(0,{\sqrt {2\lambda x}}\right)}
が形状パラメータ と速度パラメータを持つ アーラン分布 である 場合 、そのcdfは次のように与えられる。
X
∼
E
r
l
a
n
g
(
k
,
λ
)
{\displaystyle X\sim \mathrm {Erlang} (k,\lambda )}
k
{\displaystyle k}
λ
{\displaystyle \lambda }
F
X
(
x
)
=
1
−
Q
k
(
0
,
2
λ
x
)
{\displaystyle F_{X}(x)=1-Q_{k}\left(0,{\sqrt {2\lambda x}}\right)}
が自由度を持つ カイ2乗分布 である 場合 、その累積分布関数は次のように表される。
X
∼
χ
k
2
{\displaystyle X\sim \chi _{k}^{2}}
k
{\displaystyle k}
F
X
(
x
)
=
1
−
Q
k
/
2
(
0
,
x
)
{\displaystyle F_{X}(x)=1-Q_{k/2}(0,{\sqrt {x}})}
が形状パラメータ と速度パラメータを持つ ガンマ分布 である 場合 、その累積分布関数は次のように表される。
X
∼
G
a
m
m
a
(
α
,
β
)
{\displaystyle X\sim \mathrm {Gamma} (\alpha ,\beta )}
α
{\displaystyle \alpha }
β
{\displaystyle \beta }
F
X
(
x
)
=
1
−
Q
α
(
0
,
2
β
x
)
{\displaystyle F_{X}(x)=1-Q_{\alpha }(0,{\sqrt {2\beta x}})}
が形状パラメータ と尺度パラメータを持つ ワイブル分布 である 場合 、そのcdfは次のように与えられる。
X
∼
W
e
i
b
u
l
l
(
k
,
λ
)
{\displaystyle X\sim \mathrm {Weibull} (k,\lambda )}
k
{\displaystyle k}
λ
{\displaystyle \lambda }
F
X
(
x
)
=
1
−
Q
1
(
0
,
2
(
x
λ
)
k
2
)
{\displaystyle F_{X}(x)=1-Q_{1}\left(0,{\sqrt {2}}\left({\frac {x}{\lambda }}\right)^{\frac {k}{2}}\right)}
がパラメータ を持つ 一般化ガンマ分布 である 場合 、その累積分布関数は次のように表される。
X
∼
G
G
(
a
,
d
,
p
)
{\displaystyle X\sim \mathrm {GG} (a,d,p)}
a
,
d
,
p
{\displaystyle a,d,p}
F
X
(
x
)
=
1
−
Q
d
p
(
0
,
2
(
x
a
)
p
2
)
{\displaystyle F_{X}(x)=1-Q_{\frac {d}{p}}\left(0,{\sqrt {2}}\left({\frac {x}{a}}\right)^{\frac {p}{2}}\right)}
が非心パラメータ と自由度を持つ 非心カイ二乗分布 である 場合 、そのcdfは次のように与えられる。
X
∼
χ
k
2
(
λ
)
{\displaystyle X\sim \chi _{k}^{2}(\lambda )}
λ
{\displaystyle \lambda }
k
{\displaystyle k}
F
X
(
x
)
=
1
−
Q
k
/
2
(
λ
,
x
)
{\displaystyle F_{X}(x)=1-Q_{k/2}({\sqrt {\lambda }},{\sqrt {x}})}
がパラメータ を持つ レイリー分布 である 場合 、そのcdfは次のように与えられる。
X
∼
R
a
y
l
e
i
g
h
(
σ
)
{\displaystyle X\sim \mathrm {Rayleigh} (\sigma )}
σ
{\displaystyle \sigma }
F
X
(
x
)
=
1
−
Q
1
(
0
,
x
σ
)
{\displaystyle F_{X}(x)=1-Q_{1}\left(0,{\frac {x}{\sigma }}\right)}
がパラメータ を持つ マクスウェル・ボルツマン分布 である 場合 、その累積分布関数は次のように与えられる。
X
∼
M
a
x
w
e
l
l
(
σ
)
{\displaystyle X\sim \mathrm {Maxwell} (\sigma )}
σ
{\displaystyle \sigma }
F
X
(
x
)
=
1
−
Q
3
/
2
(
0
,
x
σ
)
{\displaystyle F_{X}(x)=1-Q_{3/2}\left(0,{\frac {x}{\sigma }}\right)}
が自由度を持つ カイ分布 である 場合 、その累積分布関数は次のように与えられる。
X
∼
χ
k
{\displaystyle X\sim \chi _{k}}
k
{\displaystyle k}
F
X
(
x
)
=
1
−
Q
k
/
2
(
0
,
x
)
{\displaystyle F_{X}(x)=1-Q_{k/2}(0,x)}
が形状パラメータとして、広がりパラメータとして 仲上分布 で ある 場合 、その累積分布関数は次のように表される。
X
∼
N
a
k
a
g
a
m
i
(
m
,
Ω
)
{\displaystyle X\sim \mathrm {Nakagami} (m,\Omega )}
m
{\displaystyle m}
Ω
{\displaystyle \Omega }
F
X
(
x
)
=
1
−
Q
m
(
0
,
2
m
Ω
x
)
{\displaystyle F_{X}(x)=1-Q_{m}\left(0,{\sqrt {\frac {2m}{\Omega }}}x\right)}
がパラメータ と を持つ ライス 分布 である 場合 、その累積分布関数は次のように表される。
X
∼
R
i
c
e
(
ν
,
σ
)
{\displaystyle X\sim \mathrm {Rice} (\nu ,\sigma )}
ν
{\displaystyle \nu }
σ
{\displaystyle \sigma }
F
X
(
x
)
=
1
−
Q
1
(
ν
σ
,
x
σ
)
{\displaystyle F_{X}(x)=1-Q_{1}\left({\frac {\nu }{\sigma }},{\frac {x}{\sigma }}\right)}
が非心パラメータ と自由度を持つ 非心カイ分布 である 場合 、その累積分布関数は次のように表される。
X
∼
χ
k
(
λ
)
{\displaystyle X\sim \chi _{k}(\lambda )}
λ
{\displaystyle \lambda }
k
{\displaystyle k}
F
X
(
x
)
=
1
−
Q
k
/
2
(
λ
,
x
)
{\displaystyle F_{X}(x)=1-Q_{k/2}(\lambda ,x)}
^ JI Marcum (1960). パルスレーダーによるターゲット検出の統計理論: 数学的付録、 IRE Trans. Inform. Theory、 vol. 6、59-267。
^ MK SimonとM.-S. Alouini (1998)。一般化フェージングチャネルを介したデジタル通信のパフォーマンスに対する統一的なアプローチ、 IEEEの議事録 、86(9)、1860-1877。
^ abcde A. Annamalai および C. Tellambura (2001). 一般化 Marcum-Q 関数の Cauchy-Schwarz 境界とその応用、 Wireless Communications and Mobile Computing 、1(2)、243-253。
^ abcd A. Annamalai および C. Tellambura (2008)。 実数次 M に対する一般化 Marcum Q 関数 Q M ( a , b ) の単純な指数積分表現とその応用。2008 IEEE 軍事通信会議 、サンディエゴ、カリフォルニア州、米国
^ abcdefg Y. Sun、A. Baricz、S. Zhou (2010)。一般化された Marcum および Nuttall Q 関数の単調性、対数凹性、および厳密な境界について。IEEE Transactions on Information Theory 、56(3)、1166–1186、 ISSN 0018-9448
^ ab Y. Sun および A. Baricz (2008) 一般化 Marcum Q 関数の不等式 応用数学および計算 203(2008) 134-141。
^ abcdef NM Temme (1993). 非心カイ二乗分布の漸近的および数値的側面。Computers Math. Applic. 、25(5)、55-63。
^ abcdef A. Annamalai、C. Tellambura、John Matyjas (2009)。「 分数次数 M を持つ一般化 Marcum Q 関数 Q M ( a 、 b ) の新しい展開とその応用」 2009 第 6 回 IEEE コンシューマー通信およびネットワーキング会議 、1–5、 ISBN 978-1-4244-2308-8
^ ab S. Andras、A. Baricz、Y. Sun (2011)「一般化Marcum Q関数:直交多項式アプローチ」 Acta Univ. Sapientiae Mathematica 、3(1)、60-76。
^ abcdefg YA Brychkov (2012). Marcum Q関数のいくつかの特性について。 積分変換と特殊関数 23(3), 177-182。
^ M. Abramowitz と IA Stegun (1972)。公式 10.2.12、修正球面ベッセル関数、 数学関数ハンドブック 、p. 443
^ WK Pratt (1968). MarcumのQ関数の偏微分. IEEE紀要 , 56(7), 1220-1221.
^ R. Esposito (1968). MarcumのQ関数の偏微分に関するコメント。IEEE 紀要 、56(12)、2195-2195。
^ VM Kapinas、SK Mihos、GK Karagiannidis (2009)。一般化されたMarcumおよびNuttal Q関数の単調性について。IEEE Transactions on Information Theory 、55(8)、3701-3710。
^ ab R. Li, PY Kam, H. Fu (2010). 幾何学的アプローチによる一般化Marcum Q関数の新しい表現と境界、およびその応用。IEEE Trans. Commun. 、58(1)、157-169。
^ ab MK Simon および M.-S. Alouini (2000)。一般化 Marcum Q 関数の指数型境界とフェージング チャネル上の誤り確率解析への応用。IEEE Trans. Commun. 48(3), 359-366。
^ H. Guo、B. Makki、M.-S. Alouini、T. Svensson、「予測アンテナシステムへの応用を伴う1次Marcum Q関数の半線形近似」、 IEEE Open Journal of the Communications Society 、vol.2、pp.273-286、2021年、doi:10.1109/OJCOMS.2021.3056393。
^ abcdefg DA Shnidman (1989). 検出確率の計算と一般化Marcum Q関数。IEEE Transactions on Information Theory、 35(2)、389-400。
参考文献
Marcum, JI (1950)「Q 関数表」。 米国空軍 RAND 研究メモ M-339 。カリフォルニア州サンタモニカ: Rand Corporation、1950 年 1 月 1 日。
ナットル、アルバート H. (1975): QM 関数を含むいくつかの積分 、 IEEE Transactions on Information Theory 、21(1)、95–96、 ISSN 0018-9448
Shnidman, David A. (1989): 検出確率の計算と一般化Marcum Q関数、 IEEE Transactions on Information Theory、 35(2)、389-400。
ワイスタイン、エリック・W・ マーカムQ関数。MathWorld —Wolfram Webリソースより。[1]