En teoria de la probabilitat, la funció característica d'una variable aleatòria real és una eina matemàtica que proporciona informació completa sobre la distribució de probabilitat de la variable aleatòria i sovint en facilita l'estudi. A més, amb les funcions característiques es disposa, gràcies al teorema de continuïtat de Lévy, d'un mètode senzill i potent per estudiar la convergència en distribució d'una successió de variables aleatòries.
Donada una variable aleatòria real
definida sobre un espai de probabilitat
, la seva funció característica és la funció
(és a dir de valors complexos)
definida, per a tot real t, per la relació següent (on
,
i
denota l'operador esperança):
Expressions de la funció característica
[modifica]
Expressions integrals generals
[modifica]
Per definició de
:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\int _{\Omega }\mathrm {e} ^{i\,t\,X}\,d\mathbb {P} =\int _{\Omega }\cos(t\,X)\,d\mathbb {P} \,+\,i\,\int _{\Omega }\sin(t\,X)\,d\mathbb {P} \;\;\,\qquad \qquad \qquad \quad (1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f68f1c96c0b5e6796e51aa16406c333042a69f7)
Denotant per
la distribució de probabilitat de la variable aleatòria X:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\int _{\mathbb {R} }\mathrm {e} ^{i\,t\,x}\,d\mathbb {P} _{X}(x)=\int _{\mathbb {R} }\cos(t\,X)\,d\mathbb {P} _{X}(x)\,+\,i\,\int _{\mathbb {R} }\sin(t\,X)\,d\mathbb {P} _{X}(x)\qquad (2)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c80ba7d3276b9c42ba56eb0ed26845a7d4f72556)
- (segons el teorema de la mesura imatge)
Remarques:
- la definició (1) té sentit perquè per a tot real t, la variable aleatòria complexa
![{\displaystyle \mathrm {e} ^{i\,t\,X}=\cos \,(t\,X)+i\,\sin \,(t\,X)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc94d8db3a24def6b3acb741ff4f66fd8edc5b08)
- és fitada (té mòdul 1) i, per tant, és integrable respecte a la mesura de probabilitat
;
- l'equació (2) significa que la funció característica d'una variable aleatòria real X és la transformada de Fourier de la seva distribució de probabilitat
, mesura de probabilitat sobre l'espai mesurable (o probabilitzable)
, on
és la sigma-àlgebra de Borel de ![{\displaystyle \mathbb {R} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc9de9049e03e5e5a0cab57076dbe4a369c1e3a7)
- Quan X és discreta, amb valors
tals que per a tot k,
aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\sum _{k}\mathrm {e} ^{i\,t\,x_{k}}\,p_{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9e81f1a72a1270f194b2c3f7881f6547876c6e56)
- (suma finita o sèrie absolutament convergent)
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\int _{-\infty }^{+\infty }\mathrm {e} ^{i\,t\,x}\,f_{X}(x)\,dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/774b6fb6a00f407199ced5019525d1b71a2ed27a)
- (integral de Lebesgue; en els casos usuals coincideix amb la integral de Riemann)
La funció característica d'una variable aleatòria real X:
![{\displaystyle \varphi _{X}(0)=1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8498c848b1c0a1798a332adbfaa8b1dded342c14)
![{\displaystyle \forall \,t\in \mathbb {R} ,\,|\varphi _{X}(t)|\leq 1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/93207552e0e4f8f272d5eedcba5b234a0a643d91)
(on
és el conjugat del nombre complex z)
![{\displaystyle \forall \,a\in \mathbb {R} ,\,\forall \,b\in \mathbb {R} ,\,\forall \,t\in \mathbb {R} ,\,\varphi _{a\,X+b\,}(t)=\mathrm {e} ^{i\,b\,t}\,\varphi _{X}(a\,t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ac593cb7a5cbd6f3060bf6639433ae20af46678a)
- i en particular :
;
- per tant si
i
tenen la mateixa distribució (dita simètrica), la funció
és parella amb valors reals
(la tercera propietat es dedueix del teorema de convergència dominada; les altres són immediates)
Demostració de la continuïtat uniforme
Per
linealitat de l'esperança :
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\forall \,h\in \mathbb {R} ,\varphi _{X}(t+h)-\varphi _{X}(t)=\operatorname {E} \left(\mathrm {e} ^{i\,(t+h)\,X}-\mathrm {e} ^{i\,t\,X}\right)=\operatorname {E} \left(\mathrm {e} ^{i\,t\,X}\left(\mathrm {e} ^{i\,h\,X}-1\right)\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83462f72a923e3668bc65d4eefd07be377e3461e)
i per tant:
;
.
Ara bé:
![{\displaystyle \forall \,\omega \in \Omega ,\,\mathrm {e} ^{i\,h\,X(\omega )}-1\to 0{\text{ quan }}h\to 0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/641ba8710c449cfdd9d12ade251c2ba0ea9edb2c)
- i a més
![{\displaystyle \forall \,h\in \mathbb {R} ,\,\left|\mathrm {e} ^{i\,h\,X}-1\right|\leq 2}](https://wikimedia.org/api/rest_v1/media/math/render/svg/397ba8ff189d5f3b5e4f354088443bf2e98291f4)
Aleshores pel teorema de convergència dominada (atès que la variable aleatòria constant amb valor 2 és integrable):
.
Com que la funció majorant tendeix cap a 0 i no depén de
t, això acaba la demostració.
Si la variable aleatòria X segueix la distribució degenerada de valor
(és a dir:
; X és constant quasi segurament) aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\mathrm {e} ^{i\,t\,\mu }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94252a0954ea3f792395c200dc6b884f08b4a667)
Si la variable aleatòria X segueix la distribució binomial
(on
) aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\sum _{k=0}^{n}\mathrm {e} ^{i\,t\,k}\,{n \choose k}\,p^{k}\,(1-p)^{n-k}=\sum _{k=0}^{n}{n \choose k}\,\left(p\,\mathrm {e} ^{i\,t}\right)^{k}\,(1-p)^{n-k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2ba995b787357f8a42a13d3602b4648edc873e68)
d'on es dedueix (fórmula del binomi de Newton):
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\left(p\,\mathrm {e} ^{i\,t}+1-p\right)^{n}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/209eccdfe2e836ebe5d7e94cc29a8fd5ce71bd8f)
En particular, si la variable aleatòria X segueix la distribució de Bernoulli
(on
) aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=p\,\mathrm {e} ^{i\,t}+1-p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95a8ee72e986c2bbbe5e2817baa9bf854faba9f1)
Si la variable aleatòria X segueix la distribució geomètrica
(on
) aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\sum _{k=1}^{+\infty }\mathrm {e} ^{i\,t\,k}\,p\,(1-p)^{k-1}=p\,\mathrm {e} ^{i\,t}\,\sum _{k=1}^{+\infty }\left[(1-p)\,\mathrm {e} ^{i\,t}\right]^{k-1}={\frac {p\,\mathrm {e} ^{i\,t}}{1-(1-p)\,\mathrm {e} ^{i\,t}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/658f35e5bffe339fee58d1486ef633a58dddfa11)
Si la variable aleatòria X segueix la distribució de Poisson
(on
) aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\sum _{k=0}^{+\infty }\mathrm {e} ^{i\,t\,k}\,\mathrm {e} ^{-\lambda }\,{\frac {\lambda ^{k}}{k\,!}}=\mathrm {e} ^{-\lambda }\,\sum _{k=0}^{+\infty }{\frac {\left(\lambda \,\mathrm {e} ^{i\,t}\right)^{k}}{k\,!}}=\mathrm {e} ^{-\lambda }\,\mathrm {e} ^{\lambda \,\mathrm {e} ^{i\,t}}=\mathrm {e} ^{\lambda \,\left(\mathrm {e} ^{i\,t}-1\right)}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2166af5ffcd5aad92b6ff2e797d47cc89ceaa793)
Si la variable aleatòria X segueix la distribució uniforme contínua
(on
i a < b) aleshores:
si
, i
.
En particular, si X segueix la distribució
(on
) aleshores:
si
, i
.
Si la variable aleatòria X segueix la distribució exponencial
(on
) aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\int _{0}^{+\infty }\mathrm {e} ^{i\,t\,x}\,\lambda \,\mathrm {e} ^{-\lambda \,x}\,dx=\lambda \,\int _{0}^{+\infty }\mathrm {e} ^{-(\lambda -i\,t)\,x}\,dx={\frac {\lambda }{\lambda -i\,t}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5010785bdf13da0f016baabe027860dc180c720f)
Si la variable aleatòria X segueix la distribució normal estàndard
aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\int _{-\infty }^{+\infty }\mathrm {e} ^{i\,t\,x}\,{\frac {1}{\sqrt {2\,\pi }}}\,\mathrm {e} ^{-{\frac {x^{2}}{2}}}\,dx=\mathrm {e} ^{-{\frac {t^{2}}{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f7fc2b57024e72cee912d113c8223456f300779c)
Prova
L'integrand és derivable respecte a
t, i la funció
![{\displaystyle \mathbb {R} \to \mathbb {R} ,\,x\mapsto x\,\mathrm {e} ^{-{\frac {x^{2}}{2}}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3bbd186fb306886cbf038322b91c1efbb1650da3)
és integrable sobre
![{\displaystyle \mathbb {R} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/786849c765da7a84dbc3cce43e96aad58a5868dc)
. Així, doncs, es pot
derivar sota el signe integral; la funció
![{\displaystyle \ \varphi _{X}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2504d25fcfba857b3a5c609338fbc8e703fbddc0)
és derivable i compleix la relació següent:
;
en integrar per parts:
.
Per tant, en resoldre aquesta equació diferencial:
![{\displaystyle \exists \,C\in \mathbb {C} ,\,\forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=C\,\mathrm {e} ^{-{\frac {t^{2}}{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8ba5a582f973f877b0aea9123b4c70390b94d09e)
Per fi :
, Q.E.D.
(també es pot utilitzar el
teorema de Cauchy per a
funcions holomorfes)
Distribució de Laplace centrada en el 0
[modifica]
Si la variable X segueix la distribució de Laplace centrada en 0,
, amb
aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)={\frac {1}{2b}}\int _{-\infty }^{+\infty }\mathrm {e} ^{i\,t\,x}\mathrm {e} ^{-{\frac {\vert x\vert }{b}}}\,dx={\frac {1}{1+b^{2}t^{2}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/460e441dd76633f2350331dc314c87dbf1b959eb)
Prova
![{\displaystyle {\begin{aligned}{\frac {1}{2}}\int _{-\infty }^{\infty }e^{itx}e^{-{\frac {\vert x\vert }{b}}}\,dx&={\frac {1}{2}}\int _{-\infty }^{\infty }\cos(tx)\,e^{-{\frac {\vert x\vert }{b}}}\,dx+{\frac {i}{2}}\int _{-\infty }^{\infty }\sin(tx)\,e^{-{\frac {\vert x\vert }{b}}}\,dx=\int _{0}^{\infty }\cos(tx)\,e^{-{\frac {x}{b}}}\,dx\\&{\overset {(i)}{=}}\quad {\bigg [}{\frac {b^{2}t\sin(tx)-b\cos(tx)}{1+b^{2}t^{2}}}\,e^{-{\frac {x}{b}}}{\bigg ]}_{x=0}^{\infty }={\frac {1}{1+b^{2}t^{2}}},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/35ae2ba17010be37949cf79502fa8a3b74b6c6a7)
on a la igualtat (i) hem integrat per parts dues vegades.
Distribució de Cauchy simètrica
[modifica]
Si la variable aleatòria X segueix la distribució de Cauchy simètrica
(on
) aleshores:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X}(t)=\int _{-\infty }^{+\infty }\mathrm {e} ^{i\,t\,x}\,{\frac {1}{\pi }}\,{\frac {\gamma }{x^{2}+\gamma ^{2}}}\,dx=\mathrm {e} ^{-\gamma \,|\,t\,|}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92137959afaa8dfe0260a95738ce254803981e95)
Per demostrar-ho, es pot utilitzar el teorema dels residus (anàlisi complexa).
Prova
La funció
![{\displaystyle g:z\mapsto {\frac {\mathrm {e} ^{i\,t\,z}}{z^{2}+\gamma ^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/29920a9966e7bbd1de1efdf92b67742c9785c235)
(d'una variable complexa) és
holomorfa en l'obert
![{\displaystyle U=\mathbb {C} \setminus \{\gamma \,i,\,-\gamma \,i\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d0b572fdcd4611a3e7b51a7747865a15f3d40678)
(els dos punts singulars
![{\displaystyle \gamma \,i,\,-\gamma \,i}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f64ae2ee23a49d3a7c8bce1004c282a68a92ddb5)
són pols simples de
g).
Per a tot real R tal que R > 0, sigui el semidisc compacte
![{\displaystyle \Delta _{R}=\{z\in \mathbb {C} \,/\,|\,z\,|\leq R,\,\mathrm {Im} (z)\geq 0\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/06b1f4579af43d28269bed2284ab86cbc65d99e4)
la vora del qual és
, on
és el semicercle
.
Suposant
, el pol
és l'únic punt singular de la funció g en
; per tant, segons el teorema dels residus, si
:
.
Si
, aleshores
quan
.
En efecte : si
i
,
,
i per tant, si
:
.
Per a tot real positiu t, passant al límit quan
, s'obté:
.
Altrament dit:
.
Com que la distribució estudiada és simètrica, la funció
és parella, i se'n dedueix:
, Q.E.D.
Alternativament, Cramer [1] proposa la següent demostració utilitzant que la funció característica d'una distribució de Laplace centrada en el 0 és la densitat d'una distribució de Cauchy simètrica (a part d'una constant multiplicativa) i la fórmula d'inversió per a funcions característiques integrables.
Concretament, si
, hem calculat que la seva funció característica és
Aquesta funció característica és integrable:
Aleshores, per la fórmula d'inversió per a funcions característiques integrables (vegeu més avall les fórmules d'inversió), la funció de densitat de
,
, val:
D'aquí,
O, intercanviant els noms de les variables,
Però l'expressió de la dreta és exactament la funció característica de la variable aleatòria
amb distribució de Cauchy
avaluada en el punt
,
. Per tant,
![{\displaystyle \varphi _{X}(t)=\varphi _{X}{\big (}-(-t){\big )}=e^{-\gamma \,\vert t\vert }.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47cd1b80e24a28be889da5ee2f02d990da110d8a)
Cas de la distribució normal general
[modifica]
Sigui una variable aleatòria X amb distribució normal
(on
).
Aleshores:
.
Prova
En efecte, la variable aleatòria
![{\displaystyle X^{\ast }={\frac {X-\mu }{\sigma }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dea15ea03cc7357318c48c0552dbf09634939737)
segueix la distribució normal estàndard
. Per tant, (vegeu supra):
.
Com que
, se'n dedueix que:
.
Cas de la distribució de Cauchy general
[modifica]
Sigui una variable aleatòria X amb distribució de Cauchy
(on
).
Aleshores:
.
Prova
En efecte la variable aleatòria
![{\displaystyle \ Y=X-m}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bb0218ea12a033cb149925a2bf2b674694021e6a)
segueix la distribució de Cauchy simètrica
. Per tant, (vegeu supra):
.
Com que
, se'n dedueix que:
.
Perquè la funció característica és anomenada així
[modifica]
Com el seu nom ho indica, la funció característica d'una variable aleatòria (real) en caracteritza la distribució de probabilitat: dues variables aleatòries segueixen la mateixa distribució si i només si tenen la mateixa funció característica: és el teorema d'unicitat (vegeu infra).
Per aquesta raó, la funció característica d'una variable aleatòria X també és anomenada funció característica de la distribució d'X. Per exemple, es pot parlar de la funció característica de la distribució normal.
Fórmula d'inversió general. Donada una variable aleatòria real X, es denota per
la seva funció de distribució. Per a tot parell
de punts de continuïtat de
es compleix la relació següent:
![{\displaystyle F_{X}(b)-F_{X}(a)=\lim _{\tau \to +\infty }{\frac {1}{2\pi }}\int _{-\tau }^{+\tau }{\frac {\mathrm {e} ^{-i\,t\,a}-\mathrm {e} ^{-i\,t\,b}}{i\,t}}\,\varphi _{X}(t)\,dt}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69f405fe1ba92769952a5d407b89151579117603)
Això és una variant probabilista del teorema d'inversió de la transformació de Fourier.
Quan la funció característica és integrable no cal fer el pas límit en la fórmula d'inversió i, a més, la variable aleatòria té funció de densitat que s'obté com la transformada de Fourier de la funció característica. Concretament,
Fórmula de inversió per a funcions característiques integrables. Sigui
una variable aleatòria real amb funció característica
tal que
Aleshores
té funció de densitat
donada per
A més
és contínua i afitada.
Comentaris.
1. Sigui X una variable aleatòria amb distribució normal estàndard
. Segons hem vist, la seva funció característica és
, que és integrable ja que
Per tant, per la fórmula d'inversió,
Aleshores, si definim la transformada de Fourier amb la constant normalitzadora
,[2]
tindrem que la funció
és invariant per la transformada de Fourier:
2. Hi ha moltes funcions de densitat que no són contínues o no són afitades, i llavors la seva funció característica no és integrable. Per exemple, la distribució exponencial de paràmetre
,
, té funció de densitat
Aquesta funció té una discontinuïtat en el punt 0. La seva funció característica és
que té mòdul
Llavors,
El teorema d'inversió permet reconstruir (almenys en teoria) la funció de distribució d'una variable aleatòria a partir de la seva funció característica. Una conseqüència és l'important teorema d'unicitat:
Dues variables aleatòries reals són idènticament distribuïdes si i només si tenen la mateixa funció característica.
Prova del teorema d'unicitat
Es recorda que les funcions de distribució són creixents, contínues per la dreta en tot punt de
![{\displaystyle \mathbb {R} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/786849c765da7a84dbc3cce43e96aad58a5868dc)
, tendeixen a 0 en
![{\displaystyle -\infty }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ca2608c4b5fd3bffc73585f8c67e379b4e99b6f1)
i a 1 en
![{\displaystyle +\infty }](https://wikimedia.org/api/rest_v1/media/math/render/svg/bddbb0e4420a7e744cf71bd71216e11b0bf88831)
(però no són necessàriament contínues en
![{\displaystyle \mathbb {R} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/786849c765da7a84dbc3cce43e96aad58a5868dc)
: en particular, les funcions de distribució de les variables aleatòries discretes no ho són).
El teorema de la mesura imatge (vegeu aquí la relació (2)) té com a conseqüència immediata que si dues variables aleatòries X i Y són idènticament distribuïdes (és a dir
), aleshores
.
Recíprocament, siguin dues variables aleatòries X i Y tals que
. Aleshores, segons el teorema d'inversió, per a tot parell
de punts de continuïtat de
i
es compleix la relació següent:
![{\displaystyle \ F_{X}(b)-F_{X}(a)=F_{Y}(b)-F_{Y}(a)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/15563f9820daa1b5bb72a2ac638ad0ca6d457d52)
Ara bé, les funcions
,
són creixents i, per tant, el conjunt dels punts de discontinuïtat de cadascuna és finit o numerable; per consegüent:
- existeix una successió
de punts de continuïtat de
i
tal que
;
- per a tot real x existeix una successió
de punts de continuïtat de
i
que convergeix cap a x per la dreta.
Per a tot parell
d'enters naturals:
![{\displaystyle \ F_{X}(b_{m})-F_{X}(a_{n})=F_{Y}(b_{m})-F_{Y}(a_{n})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2be13b900b6d2276e999be7deae0ae272b3f8e3b)
Passant al límit quan
, com que
, se'n dedueix:
![{\displaystyle \forall \,m\in \mathbb {N} ,\,F_{X}(b_{m})=F_{Y}(b_{m})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7c84d42d561227003c05248252f693d16d48dd89)
Finalment, passant al límit quan
, com que
per la dreta i que
,
són contínues per la dreta en tot punt:
; altrament dit:
.
Se sap que la igualtat
![{\displaystyle \ F_{X}=F_{Y}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d95bcbcf9947e70fa61a428e38a7babf4fb7013e)
implica que
X i
Y tenen la mateixa distribució; això acaba la demostració.
El més sovint, el teorema d'unicitat s'utilitza de la manera següent per determinar la distribució de probabilitat d'una variable aleatòria real X: es calcula la funció característica
i es reconeix la funció característica d'una distribució clàssica que és, per tant, la distribució d'X (per exemple, vegeu infra la prova de l'estabilitat d'algunes distribucions de probabilitat).
Funció característica de la suma de variables aleatòries independents
[modifica]
Suma de dues variables aleatòries independents
[modifica]
Donades dues variables aleatòries reals independents X i Y (definides sobre el mateix espai de probabilitat), es compleix la relació següent:
.
En efecte,
.
Atès que X i Y són independents, també ho són, per a tot real t, les variables aleatòries
i
; per tant:
.
Remarca: el recíproc és fals. Existeixen variables aleatòries no independents les funcions característiques de les quals compleixen aquesta relació. Heus aquí un exemple ben conegut: donada una variable aleatòria X amb distribució de Cauchy simètrica
:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X+X}(t)=\varphi _{2X}(t)=\varphi _{X}(2\,t)=\mathrm {e} ^{-2\,\gamma \,|\,t\,|}=\mathrm {e} ^{-\gamma \,|\,t\,|}\,\mathrm {e} ^{-\gamma \,|\,t\,|}=\varphi _{X}(t)\,\varphi _{X}(t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e2c1853666cfaea59b8c58c55da8f95d7cb32fd)
Però és clar que X i X no són independents.
Donades n variables aleatòries reals independents
(definides sobre el mateix espai de probabilitat), es compleix la relació següent:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X_{1}+\,\cdots \,+X_{n}}(t)=\prod _{k=1}^{n}\varphi _{X_{k}}(t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5a8dee773d9616310f23992186e49398bdea90b)
- (per consegüent, el producte de funcions característiques també és una funció característica).
Se sap que la transformada de Fourier d'un producte de convolució és el producte ordinari de les transformades de Fourier.
Tenint en compte el teorema d'unicitat, la relació precedent s'interpreta així: si les variables aleatòries
són independents, aleshores:
: la distribució de probabilitat de la suma és el producte de convolució de les distribucions dels termes.
Per determinar la distribució de la suma, els dos punts de vista (producte de convolució de les distribucions de probabilitat, producte ordinari de les funcions característiques) són matemàticament equivalents. Tanmateix, el mètode de les funcions característiques és generalment més simple d'utilitzar.
Aplicació: estabilitat d'algunes distribucions de probabilitat
[modifica]
Siguin n variables aleatòries reals independents
.
- si per a tot k,
segueix la distribució binomial
, aleshores
segueix la distribució binomial ![{\displaystyle {\mathcal {B}}(m_{1}+\cdots +m_{n},\,p)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/54550d1c79a12f7b74ed3cb72a147e6e2308fbbe)
- si per a tot k,
segueix la distribució de Poisson
, aleshores
segueix la distribució de Poisson ![{\displaystyle {\mathcal {P}}(\lambda _{1}+\cdots +\lambda _{n})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/76c3127155cfc44f0c8f8861bcbace03f73f0513)
- si per a tot k,
segueix la distribució normal
, aleshores
segueix la distribució normal ![{\displaystyle {\mathcal {N}}(\mu _{1}+\cdots +\mu _{n},\,\sigma _{1}^{2}+\cdots +\sigma _{n}^{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/befb156db6a20db799a957de6bd57ca7f0e0d1ad)
- si per a tot k,
segueix la distribució de Cauchy
, aleshores
segueix la distribució de Cauchy ![{\displaystyle {\mathcal {C}}(m_{1}+\cdots +m_{n},\,\gamma _{1}+\cdots +\gamma _{n})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e9a028d7fad7b93e9e4c43a49a395eecbe8e3ab9)
Prova de l'estabilitat
Provem (per exemple) la segona afirmació (la prova de les altres és anàloga):
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\varphi _{X_{1}+\,\cdots \,+X_{n}}(t)=\prod _{k=1}^{n}\varphi _{X_{k}}(t)=\prod _{k=1}^{n}\mathrm {e} ^{\lambda _{k}\,\left(\mathrm {e} ^{i\,t}-1\right)}=\mathrm {e} ^{\lambda \,\left(\mathrm {e} ^{i\,t}-1\right)}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d6ba66d538ec2d2648cb3d0428526fcf7c4ea065)
- on
.
Es reconeix la funció característica de la distribució de Poisson
![{\displaystyle {\mathcal {P}}(\lambda )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/71176b16ebd66148ddd04db0540b51bc48b9a78c)
; segons el teorema d'unicitat, la variable aleatòria
![{\displaystyle X_{1}+\,\cdots \,+X_{n}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d60f2720a90f57e19ecefbe8a5b9a60b83fe21d3)
segueix aquesta distribució.
Funció característica i moments
[modifica]
Sigui una variable aleatòria real X.
Si el moment d'ordre m d'X existeix (finit), aleshores:
- la funció característica
és de classe
en ![{\displaystyle \mathbb {R} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/786849c765da7a84dbc3cce43e96aad58a5868dc)
, i per tant:
, on
.
Si
és m vegades derivable en el punt 0, aleshores:
- per a tot natural k tal que
el moment d'ordre k d'X existeix i:
![{\displaystyle \operatorname {E} \left(X^{k}\right)=(-i)^{k}\,\varphi _{X}^{(k)}(0)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f99fcdce7a29d8be9878a7364aa4014e67c96e98)
En particular, si
és infinitament derivable en el punt 0, aleshores tots els moments d'X existeixen.
Si la variable aleatòria X segueix la distribució de Poisson
, la seva funció característica és infinitament derivable en
: tots els moments d'X existeixen. Es comprova fàcilment que:
.
Per tant:
,
, i ![{\displaystyle \operatorname {Var} (X)=\operatorname {E} (X^{2})-\left[\operatorname {E} (X)\right]^{2}=\lambda }](https://wikimedia.org/api/rest_v1/media/math/render/svg/8035cfdf04fda95ffb1f7c42bd6dc5d2a7d79c0c)
(també es poden calcular directament com a sumes de sèries convergents).
Teorema de continuïtat de Lévy
[modifica]
Aquest teorema permet estudiar la convergència en distribució de les successions de variables aleatòries per mitjà de la convergència puntual de les seves funcions característiques.
Una successió
de variables aleatòries reals convergeix en distribució cap a una variable aleatòria real X si i només si:
quan
, on
és una funció contínua en el punt 0.
En aquest cas,
és la funció característica d'X.
Una successió
de variables aleatòries reals convergeix en distribució cap a una variable aleatòria real X si i només si:
quan
.
La segona versió exigeix que sigui coneguda per endavant la distribució límit.
Heus aquí unes quantes aplicacions clàssiques del teorema de continuïtat de Lévy.
Una aplicació clàssica del teorema de continuïtat de Lévy és la prova del teorema del límit central.
Teorema de convergència de Poisson
[modifica]
Una segona aplicació clàssica és la prova del teorema de convergència de Poisson:
- Sigui una successió real
tal que
(on
) i per a tot n,
.
- Si per a tot n, la variable aleatòria
segueix la distribució binomial
, aleshores la successió
convergeix en distribució cap a una variable aleatòria X amb distribució
.
Prova
En efecte:
, on
quan
; per tant:
quan
;
- tenint en compte el teorema de continuïtat de Lévy, això acaba la prova.
Una tercera aplicació clàssica és la prova de la llei feble dels grans nombres per a variables aleatòries integrables (és a dir amb esperança finita) i independents. S'enuncia així:
- Donada una successió
de variables aleatòries reals (definides sobre el mateix espai de probabilitat) independents i idènticament distribuïdes (abreujadament i.i.d), amb esperança finita, es posa:
.
- Si es defineix per a tot n:
, on
,
- aleshores la successió
convergeix en distribució cap a la constant
.
Prova
En efecte:
- les variables aleatòries
tenen la mateixa funció característica que denotem per
(sense índex). Com que per a tot n,
(moment d'ordre 1) existeix:
, on
quan
.
- Per independència:
.
- Aleshores:
,
- altrament dit:
![{\displaystyle \forall \,t\in \mathbb {R} ,\,\forall \,n\in \mathbb {N} ^{\ast },\,\varphi _{\;{\overline {X}}_{n}}(t)=\left[1+i\,{\frac {\mu \,t}{n}}+{\frac {t}{n}}\,\varepsilon _{1}\left({\frac {t}{n}}\right)\right]^{n}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f8d2a67a75b9539531f33df57a0435b65237a677)
- i se'n dedueix:
quan
.
- La funció
és la funció característica de la variable aleatòria constant amb valor
; això acaba la prova si es té en compte el teorema de continuïtat de Lévy.
Remarca: se sap que la convergència en distribució cap a una constant equival a la convergència en probabilitat cap a la mateixa constant.
Sigui
un vector aleatori de dimensió
, és a dir, una aplicació
tal que cada component
és una variable aleatòria. La seva funció característica és l'aplicació
definida per
Amb notació vectorial, si designem per
el producte escalar ordinari de dos vectors
,
Quan no hi hagi confusió, escriurem
en lloc de
.
Càlcul de la funció característica
[modifica]
Sigui
un vector aleatori discret amb funció de probabilitat
. Aleshores la seva funció característica és
Si
és un vector aleatori amb funció de densitat
. Aleshores la seva funció característica és
Les propietats de les funcions característiques unidimensionals es trasllades al cas vectorial. Les següents propietats es troben a Sato; per a les demostracions completes vegeu Cuppens.
, on
.
.
- la funció
és uniformement contínua.
- La funció
és hermítica: ![{\displaystyle \varphi (-{\boldsymbol {t}})={\overline {\varphi ({\boldsymbol {t}})}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3f5795fb7bb120d2bfb0181f820e7beb10ea1c0b)
- En aquesta propietat és convenient escriure tots els vectors en columna, tal com és habitual en Àlgebra lineal. Designarem per
la transposada d'una matriu (o vector)
. Sigui
un vector aleatori,
un vector d'escalars i
una matriu
. Definim
Aleshores, ![{\displaystyle \varphi _{\boldsymbol {Y}}({\boldsymbol {t}})=e^{i\,<{\boldsymbol {t}},{\boldsymbol {b}}>}\,\varphi _{\boldsymbol {X}}({\boldsymbol {A't}})=e^{i\,{\boldsymbol {t}}'{\boldsymbol {b}}}\,\varphi _{\boldsymbol {X}}({\boldsymbol {A't}}),\quad {\boldsymbol {t}}=(t_{1},\dots ,t_{k})'\in \mathbb {R} ^{k}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/426d0c47e9e602dc23eef385e1c16126c3681d8d)
- Teorema d'inversió. Necessitem algunes notacions: Recordem que un conjunt
, on
és la
-àlgebra de Borel sobre
, es diu que és un conjunt de continuïtat de (la distribució de)
si
, on
és la frontera de
. Donats dos vectors,
escriurem
(respectivament
) si
(respectivament
). Si
designarem per
el conjunt
; de manera anàloga es defineix
. Si
és un conjunt de continuïtat de
, aleshores
- Teorema d'unicitat. si
i
són dos vectors aleatoris, amb funcions característiques
i
respectivament, tals que
aleshores
i
tenen la mateixa distribució. Evidentment, el recíproc també és cert.
- Funció característica i independència. Les variables aleatòries
són independents si i només si ![{\displaystyle \varphi _{(X_{1},\dots ,X_{k})}(t_{1},\dots ,t_{k})=\varphi _{X_{1}}(t_{1})\cdots \varphi _{X_{k}}(t_{k}),\quad \forall t_{1},\dots ,t_{k}\in \mathbb {R} .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff2912049cd5fd4aff1a9fa9b502cd49167bd225)
Més generalment, els vectors aleatoris
-dimensionals
són independents si i només si
- Funció característica i suma de vectors aleatoris independents. Siguin
vectors aleatoris
-dimensionals independents i posem
Aleshores ![{\displaystyle \varphi _{\boldsymbol {Y}}({\boldsymbol {t}})=\varphi _{{\boldsymbol {X}}_{1}}({\boldsymbol {t}})\cdots \varphi _{{\boldsymbol {X}}_{k}}({\boldsymbol {t}}),\quad \forall {\boldsymbol {t}}\in \mathbb {R} ^{d}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/33619803a0b6d7ce521a6de8f303ded3c61432f1)
- Funció característica i moments. Recordem que es diu que un vector aleatori
té moment d'ordre
, on
, si
, i, en aquest cas, es defineix el moment d'ordre
per
Si el vector aleatori
compleix que
, on
és la norma d'un vector
, aleshores la funció característica
és de classe
i per a qualsevol
, amb
,
Recíprocament, si la funció característica
és de classe
per a
parell, aleshores el vector
té moments d'ordre
per qualsevol
, amb
.
- Funció característica i convergència en distribució. Sigui
una successió de vectors aleatoris
-dimensionals. Designem per
la funció característica del vector
. Aleshores la successió convergeix en distribució a un vector aleatori
si i només si
on
és una funció contínua en
. En aquest cas,
és la funció característica de ![{\displaystyle {\boldsymbol {X}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/899e933d518eefcbbd0c48512cc7887ee117d040)
Considerem un experiment que pot tenir
resultats diferents, que designarem per
, amb probabilitats
,
. Fem
repeticions independents i denotem per
el nombre de vegades que obtenim el resultat
, per
el nombre de vegades que obtenim el resultat
, i així successivament. Aleshores la probabilitat d'obtenir
vegades el resultat
,
vegades el resultat
, etc. amb
és
Es diu que el vector
segueix una distribució multinomial[5][6] de paràmetres
, i s'escriu
. Cal notar que cada component
té una distribució binomial de paràmetres
i
,
. De fet, una distribució multinomial és una extensió de la distribució binomial quan hi ha més de dos resultats possibles. La funció característica del vector
és
Càcul de la funció característica
Per a
![{\displaystyle t_{1},\dots ,t_{d}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/13ee448b81228fbea88bf97ec93b8b1039a9db89)
,
![{\displaystyle {\begin{aligned}\varphi (t_{1},\dots ,t_{d})&=E(e^{i\sum _{j=1}^{d}t_{j}X_{j}})=\sum _{x_{1},\dots ,x_{d}\in \{0,\dots ,n\}, \atop \sum _{j=1}^{d}x_{j}=n}{\frac {n!}{x_{1}!\cdots x_{d}!}}\,e^{i\sum _{j=1}^{d}t_{j}x_{j}}p_{1}^{x_{1}}\cdots p_{d}^{x_{d}}\\&=\sum _{x_{1},\dots ,x_{d}\in \{0,\dots ,n\}, \atop \sum _{j=1}^{d}x_{j}=n}{\frac {n!}{x_{1}!\cdots x_{d}!}}\,{\big (}p_{1}e^{it_{1}}{\big )}^{x_{1}}\cdots {\big (}p_{d}e^{it_{d}})^{x_{d}}={\big (}p_{1}e^{it_{1}}+\dots +p_{d}e^{it_{d}}{\big )}^{n},\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1b5687082db7937b964593281b36d311529b382a)
on a l'última igualtat hem aplicat la
fórmula
![{\displaystyle (a_{1}+a_{2}+\cdots +a_{d})^{n}=\sum {\binom {n}{x_{1},\dots ,x_{d}}}a_{1}^{x_{1}}\cdots a_{d}^{x_{d}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/20d4d8b80d999677a9f8204db39b8d374a0bafdb)
on la suma es fa sobre totes les
![{\displaystyle d}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e85ff03cbe0c7341af6b982e47e9f90d235c66ab)
-
ples ![{\displaystyle (x_{1},\dots ,x_{d})\in \mathbb {\{} 0,1,\dots ,n\}^{d}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f7adb4f73802b66ad40660bdbd437e1330cd699b)
tals que
![{\displaystyle x_{1}+\cdots +x_{d}=n}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80d9e0feb7cf7fe06713680648d6863d01e42e84)
.
A partir d'aquesta funció característica podem calcular de manera senzilla
:
d'on
Distribució normal multivariant
[modifica]
Vegeu Anderson.[7] En aquest exemple escriurem tots els vectors en columna. Un vector aleatori
es diu que segueix una distribució normal
-dimensional
on
és la matriu identitat, si té funció de densitat
Cal notar que les components del vector són independents, cadascuna amb una distribució normal estàndard
. La seva funció característica és
Càcul de la funció característica
Les variables
![{\displaystyle X_{1},\dots ,X_{d}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b0316987de06d0f2686274256e042705c267e180)
són independents i totes tenen distribució
![{\displaystyle {\mathcal {N}}(0,1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a3eeb356405c0b33b680b5caa425ada4e9f53e8b)
. En efecte, per exemple, la funció de densitat marginal de
![{\displaystyle X_{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f70b2694445a5901b24338a2e7a7e58f02a72a32)
és
![{\displaystyle {\begin{aligned}f_{X_{1}}(x_{1})&=\int _{-\infty }^{\infty }\cdots \int _{-\infty }^{\infty }f_{(X_{1},\dots ,X_{d})}(x_{1},\dots ,x_{d})\,dx_{2}\cdots dx_{d}\\&={\frac {1}{\sqrt {2\pi }}}\,e^{-x_{1}^{2}/2}\,{\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }e^{-x_{2}^{2}/2}\,dx_{2}\cdots {\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }e^{-x_{d}^{2}/2}\,dx_{d}={\frac {1}{\sqrt {2\pi }}}\,e^{-x_{1}^{2}/2}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a0140efc49b04b7c3f0901527fc1a69609234a84)
Per tant, d'una banda
![{\displaystyle X_{1}\sim {\mathcal {N}}(0,1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/150fa5b8e754a96d4807f705a5aef1ab9542153b)
. I de l'altra, tenim que
d'on
són independents. Aleshores, utilitzant la relació entre variables independents i funcions característiques i l'expressió de la funció característica de la distribució normal
que hem calculat abans, tenim que per qualsevol
,
![{\displaystyle \varphi _{(X_{1},\dots ,X_{d})}(t_{1},\dots ,t_{d})=\varphi _{X_{1}}(t_{1})\cdots \varphi _{X_{d}}(t_{d})=e^{-t_{1}^{2}/2}\cdots e^{-t_{d}^{2}/2}=e^{-(t_{1}^{2}+\cdots +t_{d}^{2})/2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3269482c66ea41eb4e3fe62e82ee2f1cca167d0f)
Sigui
una matriu
definida positiva[8] i
un vector d'escalars. La matriu
té una matriu arrel quadrada
definida positiva (i per tant simètrica) única,[9] que compleix
. Definim
Per la fórmula que hem vist abans, la funció característica de
serà, per
,
D'altra banda, atès que
d'on
I per les propietats de la matriu de variàncies-covariàncies, la matriu de variàncies-covariàncies del vector
serà:
S'escriu
. Utilitzant la fórmula del canvi de variables per a vectors aleatoris amb densitat, podem calcular la funció de densitat de
, que és:
on
és el determinant de la matriu
.
En el cas que hem vist fins ara, la matriu de variàncies-covariàncies del vector normal multidimensional
era no singular, és a dir,
. Utilitzant la funció característica es pot definir un vector normal multidimensional de manera que inclogui el cas que la matriu de variàncies covariàncies sigui singular i que s'anomena vector normal multidimensional singular o degenerat ;[10][11] aquest vector està concentrat en una varietat lineal (estricte) de
i no té funció de densitat. Específicament, sigui
una matriu
definida no negativa i
un vector d'escalars; un vector aleatori
, es diu que és normal multidimensional, i s'escriu
si té funció característica
Quan
es diu que és un vector normal multidimensional singular; en aquest cas, també el vector d'esperances és
i la matriu de variàncies és
, però si el rang de
és
, aleshores la distribució de
està concentrada en una varietat lineal de dimensió
i, per tant, no té funció de densitat.
- ↑ Cramer, Harald. Métodos matemáticos de Estadística. 4a.. Madrid: Aguilar, 1970, p. 114 i 283.
- ↑ {{{títol}}}. 2a edició. títol=Real analysis|editorial=Addison-Wesley, Advanced Book Program/World Science Division|data=1983|lloc=Reading, Mass|isbn=978-0-201-14179-5|nom=Serge|cognom=Lang|pàgines=363}}
- ↑ Johnson, N. L.; Kotz, S.; Balakrihsnan, N. Discrete Multivariate Distributions. Nova York: Wiley, 1997. ISBN 0-471-12844-1.
- ↑ Forbes, C.; Evans, M.; Hastings, N.; Peacock, B. Statistical distributions.. 4th ed.. Oxford: Wiley-Blackwell, 2010, pp.135-136. ISBN 978-0-470-62724-2.
- ↑ Anderson, T. W.. An introduction to multivariate statistical analysis. 3a edició. Hoboken, N.J.: Wiley-Interscience, 2003. ISBN 0-471-36091-0.
- ↑ Per definició, una matriu definida positiva és simètrica
- ↑ Seber, G. A. F.. A matrix handbook for statisticians. Hoboken, N.J.: Wiley-Interscience, 2008, p. 225, propietat 10.3.2. ISBN 978-0-470-22678-0.
- ↑ Bryc, Wlodzimierz. The normal distribution : characterizations with applications. Nova York: Springer-Verlag, 1995. ISBN 0-387-97990-5.
- ↑ Per altres definicions alternatives, vegeu Seber, G. A. F.. A matrix handbook for statisticians. Hoboken, N.J.: Wiley-Interscience, 2008, p. 436. ISBN 978-0-470-22678-0.
- Cuppens, Roger. Decomposition of multivariate probabilities. Nova York: Academic Press, 1975. ISBN 0-12-199450-3.
- Feller, William, Introducción a la teoría de las probabilidades y sus aplicaciones (vol. 2), Mèxico : Edit. Limusa, 1978.
- Lukacs, Eugen, Characteristic Functions. London: Griffin,, 1960 (primera edició); 1970 (segona edició revisada i ampliada).
- Rényi, Alfred, Wahrscheinlichkeitsrechnung, mit einem Anhang über Informations-theorie. V. E. B. Deutscher Verlag der Wissenschaften, Berlin, 1962. Traducció al francès: Calcul des probabilités avec un appendice sur la théorie de l'information. Paris: Dunod, 1966.
- Sato, Ken-iti. Lévy processes and infinitely divisible distributions. Cambridge, U.K.: Cambridge University Press, 1999, p. 9. ISBN 0-521-55302-4.