Nous allons prouver que si X1 et X2 sont deux variables exponentielles indépendantes E(λ), alors les deux variables suivantes :
Z1=2λX1⋅cos(2πe−λX2)Z2=2λX1⋅sin(2πe−λX2)
...sont deux variables aléatoires indépendantes suivant la loi normale centrée réduite N(0,1).
Votre variable Z est simplement Z1.
La preuve est plus simple si on utilise les résultats des questions 1 et 2.
Le Rayon R :
Posons R=2λX1.
D’après votre question 2, nous savons que R suit une loi de Rayleigh(1).
Sa fonction de densité (PDF) est fR(r)=re−r2/2 pour r≥0.
L’Angle θ :
Posons U2=e−λX2. D’après votre question 1, U2 suit une loi uniforme U([0,1]).
Posons θ=2πU2=2πe−λX2.
Puisque U2 est uniforme sur [0,1], θ est une variable uniforme sur l’intervalle [0,2π].
Sa fonction de densité (PDF) est fθ(θ)=2π1 pour θ∈[0,2π].
Puisque X1 et X2 sont indépendantes, R et θ sont également indépendantes.
Étape 2 : Changement de variables (Polaire → Cartésien)¶
Nous avons maintenant deux variables indépendantes R et θ, et nous cherchons la loi de :
Z1=Rcos(θ)Z2=Rsin(θ)
C’est la transformation standard des coordonnées polaires (R,θ) vers les coordonnées cartésiennes (Z1,Z2).
Nous allons trouver la densité de probabilité jointe fZ1,Z2(z1,z2) en utilisant la formule de changement de variable.
Densité jointe de (R,θ) :
Puisqu’elles sont indépendantes, la densité jointe est le produit des densités marginales :
fR,θ(r,θ)=fR(r)⋅fθ(θ)=(re−r2/2)⋅(2π1)fR,θ(r,θ)=2πre−r2/2
(Valable pour r≥0 et θ∈[0,2π])
Formule de changement de variable :
La formule générale est fZ1,Z2(z1,z2)=fR,θ(r,θ)⋅∣J∣−1 où J est le Jacobien de la transformation (r,θ)→(z1,z2).
Une façon plus intuitive de l’écrire est d’égaliser les éléments de probabilité :
fZ1,Z2(z1,z2)dz1dz2=fR,θ(r,θ)drdθ
Nous savons que l’élément de surface cartésien dz1dz2 est lié à l’élément de surface polaire drdθ par dz1dz2=rdrdθ. (Le r est le Jacobien de la transformation).
Donc :
fZ1,Z2(z1,z2)⋅(rdrdθ)=fR,θ(r,θ)drdθfZ1,Z2(z1,z2)⋅r=fR,θ(r,θ)
En isolant fZ1,Z2(z1,z2), on obtient :
fZ1,Z2(z1,z2)=rfR,θ(r,θ)
Calcul de la densité de (Z1,Z2) :fZ1,Z2(z1,z2)=r2πre−r2/2fZ1,Z2(z1,z2)=2π1e−r2/2
Nous avons la densité jointe, mais elle est exprimée en fonction de r. Il suffit de remplacer r par son expression en fonction de z1 et z2.
De la transformation polaire, nous savons que r2=z12+z22.
En substituant :
fZ1,Z2(z1,z2)=2π1exp(−2z12+z22)
Maintenant, factorisons cette expression :
fZ1,Z2(z1,z2)=[2π1exp(−2z12)]⋅[2π1exp(−2z22)]
Nous reconnaissons que :
fZ1,Z2(z1,z2)=fZ1(z1)⋅fZ2(z2)
où fZ(z)=2π1e−z2/2 est la fonction de densité d’une loi normale centrée réduite N(0,1).
Puisque la densité jointe est le produit des densités marginales, cela prouve que :
Z1 suit la loi N(0,1).
Z2 suit la loi N(0,1).
Z1 et Z2 sont indépendantes.
La variable Z=Z1 de votre question suit donc bien la loi normale centrée réduite.