Articulo de referencia

Función Q de Marcum

En estadística , la función Q de Marcum generalizada de orden ν {\displaystyle \nu } se define como Q ν ( a , b ) = 1 a ν − 1 ∫ b ∞ incógnita ν exp ⁡ ( − incógnita 2 + a 2 2 ) I...

En estadística , la función Q de Marcum generalizada de ordenν{\displaystyle \nu }se define como

Qν(a,b)=1aν1bincógnitaνexp(incógnita2+a22)Iν1(aincógnita)dincógnita{\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}

dóndeb0{\displaystyle b\geq 0}ya,ν>0{\displaystyle a,\nu >0}yIν1{\displaystyle I_{\nu -1}}es la función de Bessel modificada de primer tipo de ordenν1{\displaystyle \nu -1}. Sib>0{\displaystyle b>0}, la integral converge para cualquierν{\displaystyle \nu }La función Q de Marcum aparece como una función de distribución acumulativa complementaria para las distribuciones chi no central , chi cuadrado no central y Rice . En ingeniería, esta función aparece en el estudio de sistemas de radar, sistemas de comunicación, sistemas de colas y procesamiento de señales. Esta función se estudió por primera vez paraν=1{\displaystyle \nu =1}por Jess Marcum, y por lo tanto nombrado en su honor, para radares pulsados. [ 1 ]

Propiedades

Representación integral finita

Utilizando el hecho de queQν(a,0)=1{\displaystyle Q_{\nu }(a,0)=1}, la función Q de Marcum generalizada puede definirse alternativamente como una integral finita como

Qν(a,b)=11aν10bincógnitaνexp(incógnita2+a22)Iν1(aincógnita)dincógnita.{\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.}

Sin embargo, es preferible tener una representación integral de la función Q de Marcum tal que (i) los límites de la integral sean independientes de los argumentos de la función, (ii) y que los límites sean finitos, (iii) y que el integrando sea una función gaussiana de estos argumentos. Para valores enteros positivos deν=norte{\displaystyle \nu =n}, dicha representación viene dada por la integral trigonométrica [ 2 ] [ 3 ]

Qnorte(a,b)={Hnorte(a,b)a<b,12+Hnorte(a,a)a=b,1+Hnorte(a,b)a>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.}

dónde

Hnorte(a,b)=ζ1norte2πexp(a2+b22)02πporque(norte1)θζporquenorteθ12ζporqueθ+ζ2exp(abporqueθ)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 ,}

y la proporciónζ=a/b{\displaystyle \zeta =a/b}es una constante.

Para cualquier realν>0{\displaystyle \nu >0}, dicha integral trigonométrica finita viene dada por [ 4 ]

Qν(a,b)={Hν(a,b)+doν(a,b)a<b,12+Hν(a,a)+doν(a,b)a=b,1+Hν(a,b)+doν(a,b)a>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.}

dóndeHnorte(a,b){\displaystyle H_{n}(a,b)}es como se definió anteriormente,ζ=a/b{\displaystyle \zeta =a/b}y el término de corrección adicional viene dado por

doν(a,b)=pecado(νπ)πexp(a2+b22)01(incógnita/ζ)ν1ζ+incógnitaexp[ab2(incógnita+1incógnita)]dincógnita.{\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.}

Para valores enteros deν{\displaystyle \nu }, el término de correccióndoν(a,b){\displaystyle C_{\nu }(a,b)}tienden a desaparecer.

Monotonicidad y concavidad logarítmica

  • La función Q de Marcum generalizadaQν(a,b){\displaystyle Q_{\nu }(a,b)}está aumentando estrictamente enν{\displaystyle \nu }ya{\displaystyle a}a pesar dea0{\displaystyle a\geq 0}yb,ν>0{\displaystyle b,\nu >0}y está disminuyendo estrictamente enb{\displaystyle b}a pesar dea,b0{\displaystyle a,b\geq 0}yν>0.{\displaystyle \nu >0.}[ 5 ]
  • La funciónνQν(a,b){\displaystyle \nu \mapsto Q_{\nu }(a,b)}es logarítmicamente cóncava en[1,){\displaystyle [1,\infty )}a pesar dea,b0.{\displaystyle a,b\geq 0.}[ 5 ]
  • La funciónbQν(a,b){\displaystyle b\mapsto Q_{\nu }(a,b)}es estrictamente logarítmicamente cóncava en(0,){\displaystyle (0,\infty )}a pesar dea0{\displaystyle a\geq 0}yν>1{\displaystyle \nu >1}, lo que implica que la función Q de Marcum generalizada satisface la propiedad de que lo nuevo es mejor que lo usado. [ 6 ]
  • La funcióna1Qν(a,b){\displaystyle a\mapsto 1-Q_{\nu }(a,b)}es logarítmicamente cóncava en[0,){\displaystyle [0,\infty )}a pesar deb,ν>0.{\displaystyle b,\nu >0.}[ 5 ]

Representación de la serie

  • La función Q generalizada de Marcum de ordenν>0{\displaystyle \nu >0}puede representarse utilizando la función Gamma incompleta como [ 7 ] [ 8 ] [ 9 ]
Qν(a,b)=1mia2/2k=01k¡γ(ν+k,b22)Γ(ν+k)(a22)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},}
dóndeγ(s,incógnita){\displaystyle \gamma (s,x)}es la función Gamma incompleta inferior . Esto se suele llamar la representación canónica de laν{\displaystyle \nu }Función Q generalizada de Marcum de orden -ésimo.
Qν(a,b)=1mia2/2k=0(1)kLk(ν1)(a22)Γ(ν+k+1)(b22)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 },}
dóndeLk(α)(){\displaystyle L_{k}^{(\alpha )}(\cdot )}es el polinomio de Laguerre generalizado de gradok{\displaystyle k}y de ordenα{\displaystyle \alpha }.
Qν(a,b)=mi(a2+b2)/2α=1ν(ab)αIα(ab),{\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),}
1Qν(a,b)=mi(a2+b2)/2α=ν(ba)αIα(ab),{\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),}
donde las sumas están en incrementos de uno. Tenga en cuenta que cuandoα{\displaystyle \alpha }asume un valor entero, tenemosIα(ab)=Iα(ab){\displaystyle I_{\alpha }(ab)=I_{-\alpha }(ab)}.
  • Para valores semi-enteros no negativosν=norte+1/2{\displaystyle \nu =n+1/2}, tenemos una expresión de forma cerrada para la función Q de Marcum generalizada como [ 8 ] [ 10 ]
Qnorte+1/2(a,b)=12[mirFdo(ba2)+mirFdo(b+a2)]+mi(a2+b2)/2k=1norte(ba)k1/2Ik1/2(ab),{\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),}
dóndemirFdo(){\displaystyle \mathrm {erfc} (\cdot )}es la función de error complementaria . Dado que las funciones de Bessel con parámetro semi-entero tienen expansiones de suma finita como [ 4 ]
I±(norte+0,5)(z)=1πk=0norte(norte+k)¡k¡(nortek)¡[(1)kmiz(1)nortemiz(2z)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],}
dóndenorte{\displaystyle n}es un entero no negativo, podemos representar exactamente la función Q de Marcum generalizada con un parámetro semi-entero. Más precisamente, tenemos [ 4 ]
Qnorte+1/2(a,b)=Q(ba)+Q(b+a)+1b2πi=1norte(ba)ik=0i1(i+k1)¡k¡(ik1)¡[(1)kmi(ab)2/2+(1)imi(a+b)2/2(2ab)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],}
para enteros no negativosnorte{\displaystyle n}, dóndeQ(){\displaystyle Q(\cdot )}es la función Q gaussiana . Alternativamente, también podemos expresar de forma más compacta las funciones de Bessel con semi-entero como suma de funciones seno y coseno hiperbólicas: [ 11 ]
Inorte+12(z)=2zπ[gramonorte(z)sinh(z)+gramonorte1(z)aporrear(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],}
dóndegramo0(z)=z1{\displaystyle g_{0}(z)=z^{-1}},gramo1(z)=z2{\displaystyle g_{1}(z)=-z^{-2}}, ygramonorte1(z)gramonorte+1(z)=(2norte+1)z1gramonorte(z){\displaystyle g_{n-1}(z)-g_{n+1}(z)=(2n+1)z^{-1}g_{n}(z)}para cualquier valor entero denorte{\displaystyle n}.

Relación de recurrencia y función generadora

  • Integrando por partes, podemos demostrar que la función Q de Marcum generalizada satisface la siguiente relación de recurrencia [ 8 ] [ 10 ].
Qν+1(a,b)Qν(a,b)=(ba)νmi(a2+b2)/2Iν(ab).{\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).}
  • La fórmula anterior se generaliza fácilmente como [ 10 ].
Qνnorte(a,b)=Qν(a,b)(ba)νmi(a2+b2)/2k=1norte(ab)kIνk(ab),{\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ν+norte(a,b)=Qν(a,b)+(ba)νmi(a2+b2)/2k=0norte1(ba)kIν+k(ab),{\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),}
para entero positivonorte{\displaystyle n}. La recurrencia anterior puede utilizarse para definir formalmente la función Q de Marcum generalizada para valores negativosν{\displaystyle \nu }. TomandoQ(a,b)=1{\displaystyle Q_{\infty }(a,b)=1}yQ(a,b)=0{\displaystyle Q_{-\infty }(a,b)=0}paranorte={\displaystyle n=\infty }, obtenemos la representación en serie de Neumann de la función Q generalizada de Marcum.
Qν+1(a,b)(1+doν(a,b))Qν(a,b)+doν(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,}
dónde
doν(a,b)=(ba)Iν(ab)Iν+1(ab).{\displaystyle c_{\nu }(a,b)=\left({\frac {b}{a}}\right){\frac {I_{\nu }(ab)}{I_{\nu +1}(ab)}}.}
Podemos eliminar la aparición de la función de Bessel para obtener la relación de recurrencia de tercer orden [ 7 ].
a22Qν+2(a,b)=(a22ν)Qν+1(a,b)+(b22+ν)Qν(a,b)b22Qν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).}
  • Otra relación de recurrencia, que la relaciona con sus derivados, viene dada por
Qν+1(a,b)=Qν(a,b)+1aaQν(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)+1bbQν(a,b).{\displaystyle Q_{\nu -1}(a,b)=Q_{\nu }(a,b)+{\frac {1}{b}}{\frac {\partial }{\partial b}}Q_{\nu }(a,b).}
  • La función generadora ordinaria deQν(a,b){\displaystyle Q_{\nu }(a,b)}para integralν{\displaystyle \nu }es [ 10 ]
norte=tnorteQnorte(a,b)=mi(a2+b2)/2t1tmi(b2t+a2/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},}
dónde|t|<1.{\displaystyle |t|<1.}

Relación de simetría

  • Utilizando las dos representaciones de la serie de Neumann, podemos obtener la siguiente relación de simetría para integrales positivos.ν=norte{\displaystyle \nu =n}
Qnorte(a,b)+Qnorte(b,a)=1+mi(a2+b2)/2[I0(ab)+k=1norte1a2k+b2k(ab)kIk(ab)].{\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].}
En particular, paranorte=1{\displaystyle n=1}tenemos
Q1(a,b)+Q1(b,a)=1+mi(a2+b2)/2I0(ab).{\displaystyle Q_{1}(a,b)+Q_{1}(b,a)=1+e^{-(a^{2}+b^{2})/2}I_{0}(ab).}

Valores especiales

Algunos valores específicos de la función Marcum-Q son [ 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)=Γ(ν,b2/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,}
  • Paraa=b{\displaystyle a=b}, al restar las dos formas de representaciones de series de Neumann, tenemos [ 10 ]
Q1(a,a)=12[1+mia2I0(a2)],{\displaystyle Q_{1}(a,a)={\frac {1}{2}}[1+e^{-a^{2}}I_{0}(a^{2})],}
que cuando se combina con la fórmula recursiva da
Qnorte(a,a)=12[1+mia2I0(a2)]+mia2k=1norte1Ik(a2),{\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}),}
Qnorte(a,a)=12[1+mia2I0(a2)]mia2k=1norteIk(a2),{\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}),}
para cualquier entero no negativonorte{\displaystyle n}.
  • Paraν=1/2{\displaystyle \nu =1/2}, utilizando la definición integral básica de la función Q generalizada de Marcum, tenemos [ 8 ] [ 10 ]
Q1/2(a,b)=12[mirFdo(ba2)+mirFdo(b+a2)].{\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].}
  • Paraν=3/2{\displaystyle \nu =3/2}, tenemos
Q3/2(a,b)=Q1/2(a,b)+2πsinh(ab)ami(a2+b2)/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}.}
  • Paraν=5/2{\displaystyle \nu =5/2}tenemos
Q5/2(a,b)=Q3/2(a,b)+2πabaporrear(ab)sinh(ab)a3mi(a2+b2)/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}.}

Formas asintóticas

  • Arroganteν{\displaystyle \nu }ser arreglado yab{\displaystyle ab}grande, dejaζ=a/b>0{\displaystyle \zeta =a/b>0}, entonces la función generalizada de Marcum-Q tiene la siguiente forma asintótica [ 7 ]
Qν(a,b)norte=0ψnorte,{\displaystyle Q_{\nu }(a,b)\sim \sum _{n=0}^{\infty }\psi _{n},}
dóndeψnorte{\displaystyle \psi _{n}}es dado por
ψnorte=12ζν2π(1)norte[Anorte(ν1)ζAnorte(ν)]ϕnorte.{\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}.}
Las funcionesϕnorte{\displaystyle \phi _{n}}yAnorte{\displaystyle A_{n}}son dados por
ϕnorte=[(ba)22ab]norte12Γ(12norte,(ba)22),{\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),}
Anorte(ν)=2norteΓ(12+ν+norte)norte¡Γ(12+νnorte).{\displaystyle A_{n}(\nu )={\frac {2^{-n}\Gamma ({\frac {1}{2}}+\nu +n)}{n!\Gamma ({\frac {1}{2}}+\nu -n)}}.}
La funciónAnorte(ν){\displaystyle A_{n}(\nu )}satisface la recursión
Anorte+1(ν)=(2norte+1)24ν28(norte+1)Anorte(ν),{\displaystyle A_{n+1}(\nu )=-{\frac {(2n+1)^{2}-4\nu ^{2}}{8(n+1)}}A_{n}(\nu ),}
paranorte0{\displaystyle n\geq 0}yA0(ν)=1.{\displaystyle A_{0}(\nu )=1.}
  • En el primer término de la aproximación asintótica anterior, tenemos
ϕ0=2πabbamirFdo(ba2).{\displaystyle \phi _{0}={\frac {\sqrt {2\pi ab}}{b-a}}\mathrm {erfc} \left({\frac {b-a}{\sqrt {2}}}\right).}
Por lo tanto, suponiendob>a{\displaystyle b>a}, la aproximación asintótica del primer término de la función generalizada de Marcum-Q es [ 7 ]
Qν(a,b)ψ0=(ba)ν12Q(ba),{\displaystyle Q_{\nu }(a,b)\sim \psi _{0}=\left({\frac {b}{a}}\right)^{\nu -{\frac {1}{2}}}Q(b-a),}
dóndeQ(){\displaystyle Q(\cdot )}es la función Q gaussiana . AquíQν(a,b)0,5{\displaystyle Q_{\nu }(a,b)\sim 0.5}comoab.{\displaystyle a\uparrow b.}
Para el caso en quea>b{\displaystyle a>b}, tenemos [ 7 ]
Qν(a,b)1ψ0=1(ba)ν12Q(ab).{\displaystyle Q_{\nu }(a,b)\sim 1-\psi _{0}=1-\left({\frac {b}{a}}\right)^{\nu -{\frac {1}{2}}}Q(a-b).}
Aquí tambiénQν(a,b)0,5{\displaystyle Q_{\nu }(a,b)\sim 0.5}comoab.{\displaystyle a\downarrow b.}

Diferenciación

aQν(a,b)=a[Qν+1(a,b)Qν(a,b)]=a(ba)νmi(a2+b2)/2Iν(ab),{\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),}
bQν(a,b)=b[Qν1(a,b)Qν(a,b)]=b(ba)ν1mi(a2+b2)/2Iν1(ab).{\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).}
Podemos relacionar las dos derivadas parciales como
1aaQν(a,b)+1bbQν+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.}
  • La n -ésima derivada parcial deQν(a,b){\displaystyle Q_{\nu }(a,b)}con respecto a sus argumentos viene dado por [ 10 ]
norteanorteQν(a,b)=norte¡(a)nortek=0[norte/2](2a2)kk¡(norte2k)¡pag=0nortek(1)pag(nortekpag)Qν+pag(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),}
nortebnorteQν(a,b)=norte¡a1ν2nortebnorteν+1mi(a2+b2)/2k=[norte/2]norte(2b2)k(nortek)¡(2knorte)¡pag=0k1(k1pag)(ab)pagIνpag1(ab).{\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).}

Desigualdades

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 pesar deab>0{\displaystyle a\geq b>0}yν>1{\displaystyle \nu >1}.

Límites

Basado en la monotonicidad y la concavidad logarítmica

Se pueden obtener varios límites superiores e inferiores de la función generalizada de Marcum-Q utilizando la monotonicidad y la concavidad logarítmica de la función.νQν(a,b){\displaystyle \nu \mapsto Q_{\nu }(a,b)}y el hecho de que tenemos una expresión de forma cerrada paraQν(a,b){\displaystyle Q_{\nu }(a,b)}cuandoν{\displaystyle \nu }tiene un valor semi-entero.

Dejarincógnita0,5{\displaystyle \lfloor x\rfloor _{0.5}}yincógnita0,5{\displaystyle \lceil x\rceil _{0.5}}denotan el par de operadores de redondeo semi-entero que mapean un número realincógnita{\displaystyle x}a su semi-entero impar más cercano a la izquierda y a la derecha, respectivamente, según las relaciones

incógnita0,5=incógnita0,5+0,5{\displaystyle \lfloor x\rfloor _{0.5}=\lfloor x-0.5\rfloor +0.5}
incógnita0,5=incógnita+0,50,5{\displaystyle \lceil x\rceil _{0.5}=\lceil x+0.5\rceil -0.5}

dóndeincógnita{\displaystyle \lfloor x\rfloor }yincógnita{\displaystyle \lceil x\rceil }denotan las funciones de piso y techo enteras.

  • La monotonicidad de la funciónνQν(a,b){\displaystyle \nu \mapsto Q_{\nu }(a,b)}a pesar dea0{\displaystyle a\geq 0}yb>0{\displaystyle b>0}nos da la siguiente cota simple [ 14 ] [ 8 ] [ 15 ]
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).}
Sin embargo, el error relativo de esta cota no tiende a cero cuandob{\displaystyle b\to \infty }. [ 5 ] Para valores enteros deν=norte{\displaystyle \nu =n}, este límite se reduce a
Qnorte0,5(a,b)<Qnorte(a,b)<Qnorte+0,5(a,b).{\displaystyle Q_{n-0.5}(a,b)<Q_{n}(a,b)<Q_{n+0.5}(a,b).}
Una muy buena aproximación de la función Q generalizada de Marcum para valores enteros.ν=norte{\displaystyle \nu =n}se obtiene tomando la media aritmética del límite superior e inferior [ 15 ]
Qnorte(a,b)Qnorte0,5(a,b)+Qnorte+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}}.}
  • Se puede obtener una cota más ajustada explotando la log-concavidad deνQν(a,b){\displaystyle \nu \mapsto Q_{\nu }(a,b)}en[1,){\displaystyle [1,\infty )}como [ 5 ]
Qν1(a,b)ν2vQν2(a,b)vν1<Qν(a,b)<Qν2(a,b)ν2ν+1Qν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 }}},}
dóndeν1=ν0,5{\displaystyle \nu _{1}=\lfloor \nu \rfloor _{0.5}}yν2=ν0,5{\displaystyle \nu _{2}=\lceil \nu \rceil _{0.5}}paraν1.5{\displaystyle \nu \geq 1.5}. La firmeza de este límite mejora a medida quea{\displaystyle a}oν{\displaystyle \nu }aumenta. El error relativo de esta cota converge a 0 a medida queb{\displaystyle b\to \infty }. [ 5 ] Para valores enteros deν=norte{\displaystyle \nu =n}, este límite se reduce a
Qnorte0,5(a,b)Qnorte+0,5(a,b)<Qnorte(a,b)<Qnorte+0,5(a,b)Qnorte+0,5(a,b)Qnorte+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)}}}.}

Cauchy-Schwarz vinculado

Utilizando la representación integral trigonométrica para valores enterosν=norte{\displaystyle \nu =n}, se puede obtener la siguiente cota de Cauchy-Schwarz [ 3 ]

mib2/2Qnorte(a,b)exp[12(b2+a2)]I0(2ab)2norte12+ζ2(1norte)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,}
1Qnorte(a,b)exp[12(b2+a2)]I0(2ab)ζ2(1norte)2(ζ21),ζ>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,}

dóndeζ=a/b>0{\displaystyle \zeta =a/b>0}.

Límites de tipo exponencial

Para fines analíticos, a menudo es útil tener límites en forma exponencial simple, aunque no sean los límites más ajustados que se puedan lograr.ζ=a/b>0{\displaystyle \zeta =a/b>0}, uno de esos límites para valores enteros ν=norte{\displaystyle \nu =n}se da como [ 16 ] [ 3 ]

mi(b+a)2/2Qnorte(a,b)mi(ba)2/2+ζ1norte1π(1ζ)[mi(ba)2/2mi(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,}
Qnorte(a,b)112[mi(ab)2/2mi(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.}

Cuandonorte=1{\displaystyle n=1}, el límite se simplifica a dar

mi(b+a)2/2Q1(a,b)mi(ba)2/2,ζ<1,{\displaystyle e^{-(b+a)^{2}/2}\leq Q_{1}(a,b)\leq e^{-(b-a)^{2}/2},\qquad \zeta <1,}
112[mi(ab)2/2mi(a+b)2/2]Q1(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.}

Otro límite de este tipo obtenido mediante la desigualdad de Cauchy-Schwarz se da como [ 3 ].

mib2/2Qnorte(a,b)122norte12+ζ2(1norte)2(1ζ2)[mi(ba)2/2+mi(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}
Qnorte(a,b)112ζ2(1norte)2(ζ21)[mi(ba)2/2+mi(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.}

Límite de tipo Chernoff

Límites de tipo Chernoff para la función Q de Marcum generalizada, dondeν=norte{\displaystyle \nu =n}es un número entero, está dado por [ 16 ] [ 3 ]

(12λ)norteexp(λb2+λnortea212λ){Qnorte(a,b),b2>norte(a2+2)1Qnorte(a,b),b2<norte(a2+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.}

donde el parámetro de Chernoff(0<λ<1/2){\displaystyle (0<\lambda <1/2)}tiene valor óptimoλ0{\displaystyle \lambda _{0}}de

λ0=12(1norteb2norteb21+(ab)2norte).{\displaystyle \lambda _{0}={\frac {1}{2}}\left(1-{\frac {n}{b^{2}}}-{\frac {n}{b^{2}}}{\sqrt {1+{\frac {(ab)^{2}}{n}}}}\right).}

Aproximación semilineal

La función Marcum-Q de primer orden puede aproximarse semilinealmente mediante [ 17 ].

Q1(a,b)={1,                                                                                   iF b<do1β0mi12(a2+(β0)2)I0(aβ0)(bβ0)+Q1(a,β0),     iF do1bdo20,                                                                                   iF b>do2{\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}}}

dónde

β0=a+a2+22,{\displaystyle {\begin{aligned}\beta _{0}={\frac {a+{\sqrt {a^{2}+2}}}{2}},\end{aligned}}}
do1(a)=máximo(0,β0+Q1(a,β0)1β0mi12(a2+(β0)2)I0(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}}}

y

do2(a)=β0+Q1(a,β0)β0mi12(a2+(β0)2)I0(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}}}

Formas equivalentes para un cálculo eficiente

Es conveniente volver a expresar la función Q de Marcum como [ 18 ].

PAGnorte(incógnita,Y)=Qnorte(2norteincógnita,2Y).{\displaystyle P_{N}(X,Y)=Q_{N}({\sqrt {2NX}},{\sqrt {2Y}}).}

ElPAGnorte(incógnita,Y){\displaystyle P_{N}(X,Y)}puede interpretarse como la probabilidad de detección denorte{\displaystyle N}muestras de señal recibidas integradas de forma incoherente con una relación señal-ruido constante,incógnita{\displaystyle X}, con un umbral de detección normalizadoY{\displaystyle Y}. En esta forma equivalente de la función Q de Marcum, para un valor dadoa{\displaystyle a}yb{\displaystyle b}, tenemosincógnita=a2/2norte{\displaystyle X=a^{2}/2N}yY=b2/2{\displaystyle Y=b^{2}/2}Existen muchas expresiones que pueden representarPAGnorte(incógnita,Y){\displaystyle P_{N}(X,Y)}Sin embargo, a continuación se presentan los cinco métodos más fiables, precisos y eficientes para el cálculo numérico. Son de la forma uno: [ 18 ]

PAGnorte(incógnita,Y)=k=0minorteincógnita(norteincógnita)kk¡metro=0norte1+kmiYYmetrometro¡,{\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!}},}

segundo formulario: [ 18 ]

PAGnorte(incógnita,Y)=metro=0norte1miYYmetrometro¡+metro=nortemiYYmetrometro¡(1k=0metronorteminorteincógnita(norteincógnita)kk¡),{\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),}

tercer formulario: [ 18 ]

1PAGnorte(incógnita,Y)=metro=nortemiYYmetrometro¡k=0metronorteminorteincógnita(norteincógnita)kk¡,{\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!}},}

formulario cuatro: [ 18 ]

1PAGnorte(incógnita,Y)=k=0minorteincógnita(norteincógnita)kk¡(1metro=0norte1+kmiYYmetrometro¡),{\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),}

y forma cinco: [ 18 ]

1PAGnorte(incógnita,Y)=mi(norteincógnita+Y)r=norte(Ynorteincógnita)r/2Ir(2norteincógnitaY).{\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}}).}

Entre estas cinco formas, la segunda es la más robusta. [ 18 ]

Aplicaciones

La función Q generalizada de Marcum se puede utilizar para representar la función de distribución acumulativa (FDA) de muchas variables aleatorias:

  • Siincógnitamiincógnitapag(λ){\displaystyle X\sim \mathrm {Exp} (\lambda )}es una distribución exponencial con parámetro de tasaλ{\displaystyle \lambda }, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Q1(0,2λincógnita){\displaystyle F_{X}(x)=1-Q_{1}\left(0,{\sqrt {2\lambda x}}\right)}
  • Siincógnitamirlanortegramo(k,λ){\displaystyle X\sim \mathrm {Erlang} (k,\lambda )}es una distribución de Erlang con parámetro de formak{\displaystyle k}y parámetro de velocidadλ{\displaystyle \lambda }, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qk(0,2λincógnita){\displaystyle F_{X}(x)=1-Q_{k}\left(0,{\sqrt {2\lambda x}}\right)}
  • Siincógnitaχk2{\displaystyle X\sim \chi _{k}^{2}}es una distribución chi-cuadrado conk{\displaystyle k}grados de libertad, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qk/2(0,incógnita){\displaystyle F_{X}(x)=1-Q_{k/2}(0,{\sqrt {x}})}
  • SiincógnitaGRAMOametrometroa(α,β){\displaystyle X\sim \mathrm {Gamma} (\alpha ,\beta )}es una distribución gamma con parámetro de formaα{\displaystyle \alpha }y parámetro de velocidadβ{\displaystyle \beta }, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qα(0,2βincógnita){\displaystyle F_{X}(x)=1-Q_{\alpha }(0,{\sqrt {2\beta x}})}
  • SiincógnitaWmiibll(k,λ){\displaystyle X\sim \mathrm {Weibull} (k,\lambda )}es una distribución de Weibull con parámetros de formak{\displaystyle k}y parámetro de escalaλ{\displaystyle \lambda }, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Q1(0,2(incógnitaλ)k2){\displaystyle F_{X}(x)=1-Q_{1}\left(0,{\sqrt {2}}\left({\frac {x}{\lambda }}\right)^{\frac {k}{2}}\right)}
  • SiincógnitaGRAMOGRAMO(a,d,pag){\displaystyle X\sim \mathrm {GG} (a,d,p)}es una distribución gamma generalizada con parámetrosa,d,pag{\displaystyle a,d,p}, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qdpag(0,2(incógnitaa)pag2){\displaystyle F_{X}(x)=1-Q_{\frac {d}{p}}\left(0,{\sqrt {2}}\left({\frac {x}{a}}\right)^{\frac {p}{2}}\right)}
  • Siincógnitaχk2(λ){\displaystyle X\sim \chi _{k}^{2}(\lambda )}es una distribución chi-cuadrado no central con parámetro de no centralidadλ{\displaystyle \lambda }yk{\displaystyle k}grados de libertad, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qk/2(λ,incógnita){\displaystyle F_{X}(x)=1-Q_{k/2}({\sqrt {\lambda }},{\sqrt {x}})}
  • SiincógnitaRaylmiigramoh(σ){\displaystyle X\sim \mathrm {Rayleigh} (\sigma )}es una distribución de Rayleigh con parámetroσ{\displaystyle \sigma }, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Q1(0,incógnitaσ){\displaystyle F_{X}(x)=1-Q_{1}\left(0,{\frac {x}{\sigma }}\right)}
  • SiincógnitaMETROaincógnitawmill(σ){\displaystyle X\sim \mathrm {Maxwell} (\sigma )}es una distribución de Maxwell-Boltzmann con parámetroσ{\displaystyle \sigma }, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Q3/2(0,incógnitaσ){\displaystyle F_{X}(x)=1-Q_{3/2}\left(0,{\frac {x}{\sigma }}\right)}
  • Siincógnitaχk{\displaystyle X\sim \chi _{k}}es una distribución chi conk{\displaystyle k}grados de libertad, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qk/2(0,incógnita){\displaystyle F_{X}(x)=1-Q_{k/2}(0,x)}
  • Siincógnitanorteakagramoametroi(metro,Ω){\displaystyle X\sim \mathrm {Nakagami} (m,\Omega )}es una distribución de Nakagami conmetro{\displaystyle m}como parámetro de forma yΩ{\displaystyle \Omega }como parámetro de dispersión, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qmetro(0,2metroΩincógnita){\displaystyle F_{X}(x)=1-Q_{m}\left(0,{\sqrt {\frac {2m}{\Omega }}}x\right)}
  • SiincógnitaRidomi(ν,σ){\displaystyle X\sim \mathrm {Rice} (\nu ,\sigma )}es una distribución de Rice con parámetrosν{\displaystyle \nu }yσ{\displaystyle \sigma }, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Q1(νσ,incógnitaσ){\displaystyle F_{X}(x)=1-Q_{1}\left({\frac {\nu }{\sigma }},{\frac {x}{\sigma }}\right)}
  • Siincógnitaχk(λ){\displaystyle X\sim \chi _{k}(\lambda )}es una distribución chi no central con parámetro de no centralidadλ{\displaystyle \lambda }yk{\displaystyle k}grados de libertad, entonces su función de distribución acumulada viene dada porFincógnita(incógnita)=1Qk/2(λ,incógnita){\displaystyle F_{X}(x)=1-Q_{k/2}(\lambda ,x)}

Notas a pie de página

  1. JI Marcum (1960). Una teoría estadística de la detección de objetivos mediante radar pulsado: apéndice matemático, IRE Trans. Inform. Theory, vol. 6, 59-267.
  2. MK Simon y M.-S. Alouini (1998). Un enfoque unificado para el rendimiento de la comunicación digital sobre canales de desvanecimiento generalizados, Proceedings of the IEEE , 86(9), 1860-1877.
  3. 1 2 3 4 5 A. Annamalai y C. Tellambura (2001). Cota de Cauchy-Schwarz en la función generalizada de Marcum-Q con aplicaciones, Wireless Communications and Mobile Computing , 1(2), 243-253.
  4. 1 2 3 4 A. Annamalai y C. Tellambura (2008). Una representación integral exponencial simple de la función Q de Marcum generalizada Q M ( a , b ) para M de orden realcon aplicaciones. Conferencia de Comunicaciones Militares IEEE de 2008 , San Diego, CA, EE. UU.
  5. 1 2 3 4 5 6 7 Y. Sun, A. Baricz y S. Zhou (2010). Sobre la monotonicidad, la log-concavidad y las cotas ajustadas de las funciones Q generalizadas de Marcum y Nuttall. IEEE Transactions on Information Theory , 56(3), 1166–1186, ISSN 0018-9448 
  6. 1 2 Y. Sun y A. Baricz (2008). Desigualdades para la función Q de Marcum generalizada. Matemáticas Aplicadas y Computación 203(2008) 134-141.
  7. 1 2 3 4 5 6 N.M. Temme (1993). Aspectos asintóticos y numéricos de la distribución chi-cuadrado no central. Computers Math. Applic. , 25(5), 55-63.
  8. 1 2 3 4 5 6 A. Annamalai, C. Tellambura y John Matyjas (2009). "Una nueva variante de la función Q de Marcum generalizada Q M ( a , b ) con orden fraccional M y sus aplicaciones". 2009 6th IEEE Consumer Communications and Networking Conference , 1–5, ISBN  978-1-4244-2308-8
  9. 1 2 S. Andras, A. Baricz y Y. Sun (2011) La función Q de Marcum generalizada: un enfoque polinomial ortogonal. Acta Univ. Sapientiae Mathematica , 3(1), 60-76.
  10. 1 2 3 4 5 6 7 Y.A. Brychkov (2012). Sobre algunas propiedades de la función Q de Marcum. Transformadas Integrales y Funciones Especiales 23(3), 177-182.
  11. M. Abramowitz e IA Stegun (1972). Fórmula 10.2.12, Funciones de Bessel esféricas modificadas , Manual de funciones matemáticas , pág. 443
  12. WK Pratt (1968). Diferenciales parciales de la función Q de Marcum. Actas del IEEE , 56(7), 1220-1221.
  13. R. Esposito (1968). Comentario sobre diferenciales parciales de la función Q de Marcum. Actas del IEEE , 56(12), 2195-2195.
  14. VM Kapinas, SK Mihos, GK Karagiannidis (2009). Sobre la monotonicidad de las funciones Q generalizadas de Marcum y Nuttal. IEEE Transactions on Information Theory , 55(8), 3701-3710.
  15. 1 2 R. Li, PY Kam y H. Fu (2010). Nuevas representaciones y límites para la función Q de Marcum generalizada mediante un enfoque geométrico y una aplicación. IEEE Trans. Commun. , 58(1), 157-169.
  16. 1 2 M.K. Simon y M.-S. Alouini (2000). Límites de tipo exponencial en la función Q generalizada de Marcum con aplicación al análisis de probabilidad de error en canales con desvanecimiento. IEEE Trans. Commun. 48(3), 359-366.
  17. H. Guo, B. Makki, M. -S. Alouini y T. Svensson, "Una aproximación semilineal de la función Q de Marcum de primer orden con aplicación a sistemas de antenas predictoras", en IEEE Open Journal of the Communications Society , vol. 2, pp. 273-286, 2021, doi: 10.1109/OJCOMS.2021.3056393.
  18. 1 2 3 4 5 6 7 D.A. Shnidman (1989). El cálculo de la probabilidad de detección y la función Q de Marcum generalizada. IEEE Transactions on Information Theory, 35(2), 389-400.

Referencias

  • Marcum, JI (1950) "Tabla de funciones Q". Memorando de investigación RAND de la Fuerza Aérea de EE. UU. M-339 . Santa Mónica, CA: Rand Corporation, 1 de enero de 1950.
  • Nuttall, Albert H. (1975): Algunas integrales que involucran la función Q M , IEEE Transactions on Information Theory , 21(1), 95–96, ISSN 0018-9448 
  • Shnidman, David A. (1989): El cálculo de la probabilidad de detección y la función Q de Marcum generalizada, IEEE Transactions on Information Theory, 35(2), 389–400.
  • Weisstein, Eric W. Función Q de Marcum. De MathWorld—Un recurso web de Wolfram.