Variance et lois des variables aléatoires
Variance et lois des variables aléatoires
Exercice 1
ex
Soit F : R → R la fonction définie pour x ∈ R par F (x) = 2 1(−∞,0) (x) + 1[0,∞) (x).
1. Montrer que F est une fonction de répartition d’une mesure de probabilité.
2. En cas de réponse positive à la question précédante, on se donne X une v.a.r. de fonction de répartition F . La
variable X admet-elle une densité ?
Exercice 2
2
Soit X une v.a.r. de densité f (x) = xe−x /2
1(0,∞) (x).
1. Vérifier que f est une densité de probabilité.
2. La variable aléatoire Y = X 2 est-elle à densité ? Reconnaı̂tre la loi de Y .
3. Calculer l’espérance et la variance de Y .
Exercice 3
Soit X ∼ N (0, 1).
1. Déterminer les lois de eX (cette loi est appelée loi log-normale), |X| et X 2 (loi du χ2 (1)). Pour ce faire, on déterminera
les fonctions de répartition et les densités.
2. Calculer l’espérance et la variance de ces v.a.r..
3. Calculer le moment d’ordre k (k ∈ N∗ ) de X.
Exercice 4
Une v.a.r. X suit une loi de Weibull W(α, β) si sa densité est égale à :
β
fX (x) = αβxβ−1 e−αx 1R+ (x), α, β > 0.
fX (x)
hX (x) = , x ≥ 0.
1 − FX (x)
Expliciter hX .
3. Calculer le moment d’ordre k (k ∈ N∗ ) de X. En déduire l’espérance et la variance de X.
Exercice 5
Montrer qu’une v.a. X non identiquement nulle à valeurs positives suit une loi exponentielle si et seulement si elle est
diffuse et sans mémoire, i.e. pour tous s, t > 0 : P(X > s + t|X > s) = P(X > t).
Exercice 6
La loi bêta de 1ère espèce βI (a, b), a, b > 0, admet pour densité :
Z 1
1
f (x) = xa−1 (1 − x)b−1 1[0,1] (x), où B(a, b) = xa−1 (1 − x)b−1 dx.
B(a, b) 0
Γ(a)Γ(b)
1. Montrer que B(a, b) = Γ(a+b) pour tous a, b > 0.
2. Soit U ∼ U[0,1] . Montrer que U et U 2 suivent des lois bêta particulières dont on précisera les paramètres.
3. Calculer le moment d’ordre k (k ∈ N∗ ) de X ∼ βI (a, b).
4. En déduire l’espérance et la variance de X.
Exercice 7
Soit X ∼ N (0, 1). On s’intéresse à la v.a. X + = max(X, 0).
R
1. Que vaut R
h(x) dδ0 (x) ?
2. À l’aide de la caractérisation par des fonctions tests, montrer que la loi de X + est un mélange d’une masse de Dirac
en 0 et d’une loi à densité qu’on déterminera.
3. Déterminer l’espérance et la variance de X + .
Exercice 8
a
Soit f : R → R+ la fonction définie pour x ∈ R par f (x) = π(a2 +x2 ) où a > 0 est un paramètre.
1. Montrer que f est une densité de probabilité. La loi associée est appelée loi de Cauchy (centrée) de paramètre a > 0
et est notée C(a).
2. Soit X ∼ C(1). Montrer que X n’admet pas de premier moment.
3. Déterminer la loi de log |X|.
Exercice 9
Soit X la durée de vie, en années, d’un composant électronique. On suppose que X ∼ E(θ) (loi exponentielle), θ > 0.
En outre, l’exploitant a une politique le conduisant à changer systématiquement tout composant dont la durée de vie a
atteint 5 ans. Soit Y la durée de vie d’un composant. Déterminer la loi de Y appelée loi tronquée.
Exercice 10
Soit X ∼ E(θ). Déterminer la loi de Y = bXc où b·c est la partie entière.
Exercice 11
Soit (X, Y ) et (U, V ) deux vecteurs aléatoires de densités respectives
1 1
f (x, y) = (1 + xy)1[−1,1]2 (x, y) et g(u, v) = 1[−1,1]2 (u, v).
4 4
1. Vérifier que f et g sont bien des densités sur R2 .
2. Déterminer les densités marginales de (X, Y ) et (U, V ) notées respectivement fX , fY , gU et gV .
3. Justifier que X et U d’une part, et que Y et V d’autre part ont même loi mais que, cependant, (X, Y ) et (U, V ) ne
suivent pas la même loi.
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
Exercice 1
1. La fonction F est croissante de R dans [0, 1], continue à droite, limitée à gauche (càdlàg) — elle est en fait continue
partout sauf en 0, l’ouverture/fermeture des bornes sur les indicatrices est alors primordiale —, en limx→−∞ F (x) = 0
et limx→∞ F (x) = 1. Donc F est une fonction de répartition.
2. La fonction F est dérivable sauf en 0. Un candidat
R 0 naturelle pour la densité de X est sa dérivée : F 0 (x) = 12 ex 1(−∞,0)
1
presque partout. Par contre, on vérifie que F (x) dx = 2 6= 1, donc X n’admet pas de densité. On aurait pu aussi
remarquer que P(X = 0) = 21 contredisant ainsi le fait que X admette une densité.
Exercice 2
1. La fonction f est mesurable positive et
Z Z ∞
2 2
f (x) dx = xe−x /2
dx = [−e−x /2 ∞
]0 = 1.
R 0
2. Utilisons la méthode des fonctions tests. Soit g une fonction mesurable bornée et calculons
Z Z ∞
2
g(y)fY (y) dy = E(g(Y )) = E(g(x2 )) = g(x2 )xe−x /2 dx.
0
√ √
On fait le changement de variables y = x2 , d’où x = y et dx = dy/2 y :
Z ∞
1
g(y) e−y/2 dy.
0 2
De même, Z ∞ Z ∞
2 1 −x/2
2
E(Y ) = x e dx = [−x2 e−x/2 ]∞
0 + 2xe−x/2 dx = 4E(Y ) = 8.
0 2 0
D’où V(Y ) = 8 − 22 = 4.
Exercice 3
1. Soit g une fonction mesurable bornée.
— 2
e−x /2
Z
X x
E(g(e )) = g(e ) √ dx.
R 2π
On pose y = ex , x = ln y et dx = dy/y :
∞ 2
e−(ln y)
Z
E(g(eX )) = g(y) √ dy.
0 y 2π
— 2
∞
e−x /2
Z
E(g(|X|)) = g(x)2 √ dx.
0 2π
— 2
∞ Z ∞
e−x /2 e−y/2
Z
E(g(X 2 )) = 2 g(x2 ) √ dx = g(y) √ dy.
0 2π 0 2πy
Nous pouvons également tenter de caractériser ces lois à la de la fonction de répartition, même si ce n’est la méthode
la plus facile à mettre en oeuvre.
— Soit t ∈ R alors
F (t) = P(eX ≤ t) = P(X ≤ ln t) = Φ(ln t),
où Φ est la fonction de répartition d’une loi gaussienne centrée et réduite. Pour obtenir la densité, on calcule la
dérivée : 2
e−(ln t) /2
F 0 (t) = √ .
t 2π
R∞
On vérifie facilement que 0 F 0 (t) dt = F (∞) − F (−∞) = 1.
— Pour t ≤ 0, P(|X| ≤ t) = 0, sinon, pour t > 0, nous avons
F (t) = P(|X| ≤ t) = Φ(t) − Φ(−t).
Comme précédemment,
2 ∞
e−t /2
Z
F 0 (t) = 2 √ 1[0,∞) , et F 0 (t) dt = F (∞) − F (0) = 1.
2π 0
Exercice 4
1. Pour t ≤ 0, P(X ≤ 0) = 0. Si t > 0, alors
Z t
β β β
P(X ≤ t) = αβxβ−1 e−αx dx = [−e−αx ]t0 = 1 − e−αt .
0
Réciproquement, supposons que X soit sans mémoire. Alors, pour tout entier n ≥ 0 et t > 0 :
Ceci montre que P(X > nt) = (P(X > t))n . Maintenant, pour tout p ∈ Z et q ∈ N∗ nous obtenons
P(X > p/q) = (P(X > 1/q))p et P(X > 1) = P(X > q/q) = (P(X > 1/q))q .
p
Finalement, P(X > p/q) = (P(X > 1))p/q = e q ln P(X>1) .
Supposons un instant que P(X > 1) = 0 alors pour tout n ≥ 1, 0 = P(X > 1) = P(X > 1/n)n donc P(X > 1/n) = 0 ; il
vient alors par postivité de X que
[ X
P(X 6= 0) = P(X < 0) + P(X > 0) = P {X > 1/n} ≤ P(X > 1/n) = 0
n≥1 n≥1
impliquant P(X = 0) = 1. Par conséquent si X est supposée positive et non identiquement nulle alors P(X > 1) > 0.
Enfin, X est diffuse donc sa fonction de répartition et sa queue de distribution est continue. Par densité de Q dans
R, on obtient par unicité du prolongement par continuité, que P(X > t)et ln P(X>1) pour tout t ≥ 0 et X suit une loi
exponentielle.
En fait, si X n’est pas supposée diffuse, elle suit une loi géométrique ou exponentielle suivant que P(X = 0) > 0 ou non.
Exercice 6
R∞
1. Rappelons que Γ(x) = 0
tx−1 e−t dx. Alors
Z ∞ Z ∞
Γ(a)Γ(b) = ta−1 sb−1 e−(t+s) dsdt.
0 0
En faisant le changement de variable (u, v) = (t + s, t/(t + s)) ou de manière équivalente (t, s) = (uv, u(1 − v)). Il est
clair que (t, s) = φ(u, v) = (uv, u(1 − v)) est un C 1 -difféomorphisme de (0, ∞) × (0, 1) dans (0, ∞)2 . La jacobienne
de φ est donnée par
v u
Jacφ(u, v) = ,
1 − v −u
et son déterminant est −u. D’où, par la formule du changement de variables
Z ∞Z 1 Z ∞ Z 1
Γ(a)Γ(b) = (uv)a−1 ub−1 (1 − v)b−1 e−u u dvdu = ua+b−1 e−u du v a−1 (1 − v)b−1 dv = Γ(a + b)B(a, b).
0 0 0 0
2. Il est clair que U ∼ βI (1, 1). Pour U 2 , considérons g une fonction mesurable bornée
Z 1 Z 1
2 2 dy
E(g(U )) = g(x ) dx = g(y) √ .
0 0 2 y
a a(a+1) ab
4. Ainsi, E(X) = a+b et E(X 2 ) = (a+b)(a+b+1) si bien que V(X) = (a+b)2 (a+b+1) .
Exercice 7
R
1. Par définition si h = 1A alors h dδ0 = δ0 (A) = 1A (0) = h(0). Cette égalité est encore vraie si h est une fonction
étagée et par convergence monotone si h est une fonction mesurableRpositive. Si h est intégrable, en décomposant
h = h+ − h− en parties positives et négatives, nous obtenons encore h dδ0 = h(0).
2. Soit g une fonction mesurable bornée et calculons
2 r
e−x /2
Z Z
+ 1 1 2 −x2 /2
E(g(X )) = g(x ∨ 0) √ dx = g(0) + g(x) e 1R+ (x) dx.
R 2π 2 2 R π
q 2
Ceci montre que la loi de X + est une combinaison convexe de δ0 et de la loi de densité π2 e−x /2 1R+ (x).
3. Il suffit de remplacer g par la fonction identité :
1 ∞
r
1 h −x2 /2 i∞
Z
+ 2 −x2 /2 1
E(X ) = x e 1R+ (x) dx = √ −e =√ .
2 0 π 2π 0 2π
Exercice 8
1. La fonction f est clairement positive et
Z ∞
1
f (x) dx = arctan(x/a) = 1.
R π −∞
2. Par définition ∞
|x|
Z
E|X| = dx = ∞.
−∞ π(1 + x2 )
Pour la dernière égalité, il suffit de minorer à l’infini par une en 1/|x| qui n’est pas intégrable.
3. Soit g une fonction mesurable bornée. Alors
Z Z ∞
1 1
E(g(log |X|)) = g(log |x|) 2
dx = 2 g(log x) dx
R π(1 + x ) 0 π(1 + x2 )
Z ∞
ey
Z
1
=2 g(y) dy = g(y) 1R (y) dy.
−∞
2y
π(1 + e ) R π cosh(y) +
Exercice 9
Par définition, Y = X ∧ 5. Soit g une fonction mesurable bornée, alors
Z ∞ Z
E(g(Y )) = g(x ∧ 5)θe−θx dx = g(x)θe−θx 1[0,5] (x) dx + g(5)e−5θ .
0 R
θe−5θ
Autrement dit, Y ∼ (1 − e−5θ )h + e−5θ δ5 où h est la densité h(x) = 1
1−e−5θ [0,5]
(x).
Exercice 10
Soit g une fonction mesurable bornée. Calculons
Z ∞ ∞
X Z k+1 ∞
X
E(g(bxc)) = g(bxc)θe−θx dx = g(k) θe−θx dx = (1 − e−θ ) g(k)e−kθ .
0 k=0 k k=0
Exercice 11
1. Les fonctions f et g sont positives. De plus :
Z 1Z 1 1
1 1 x2 y 1 1
Z Z
1
(1 + xy) dxdy = x+ dy = 2 dy = 1,
−1 −1 4 4 −1 2 −1 4 −1
et Z
1 1
dudv = λ2 ([−1, 1]2 ) = 1.
−1,1 4 4
2. Pour calculer la loi de X il suffit d’intégrer la densité jointe par rapport à sa seconde variable :
Z 1
1 1
fX (x) = (1 + xy) dy1[−1,1] (x) = 1[−1,1] (x)
−1 4 2
Remarquons que la densité f est symétrique : f (x, y) = f (y, x). Ainsi, fY = fX .
1 1
D’un autre côté g(u, v) = 1[−1,1] (u) 1[−1,1] (v).
2
| {z 2
} | {z }
gU (u) gV (v)
3. On remarque que fU = fV = gU = gV si bien que U, V, X, Y sont identiquement distribuées. Cependant, f 6= g et
donc (X, Y ) et (U, V ) ne sont pas identiquement distribuées. Ceci montre un point important, de la loi du couple
nous pouvons toujours déduire les lois marginales ; par contre, les lois marginales ne suffisent pas en général à calculer
la loi du couple.
Si E(Y 2 ) = 0 alors Y = 0 presque sûrement et l’inégalité de Cauchy-Schwartz est trivialement vérifiée. Supposons
donc E(Y 2 ) > 0. Puisque la fonction g est de signe constant, nous obtenons que le discriminant est négatif ou nul :
4E(XY )2 ≤ 4E(X 2 )E(Y 2 ), d’où l’inégalité annoncée. Notons que le cas d’égalité signifie que g s’annule en un seul
point λ0 et dans ce cas X + λY = 0 presque sûrement, c’est à dire X et Y sont colinéaires dans L2 .
3. Soit α ∈ (0, 1) alors puisque X est positive
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
TD2 : Indépendance
Exercice 1
Soit X une v.a.r. de densité
fθ (x) = cθ x1[0,θ] (x), x ∈ R,
où θ > 0 est fixé. Soient X1 , . . . , Xn des variables aléatoires i.i.d. de même loi que X. On pose :
n
1X
Xn = Xi et Yn = max(X1 , . . . , Xn ).
n i=1
Exercice 2
Soient X1 , . . . , Xn des v.a.r. i.i.d.. On note µ = E(X1 ) et V(X1 ) = σ 2 . La variance empirique est donnée par l’expression
n
1X
Sn2 = (Xi − X n )2 .
n i=1
Montrer que
n
1X 2 2
Sn2 = X − Xn
n i=1 i
et exprimer l’espérance de Sn2 en fonction de σ 2 . (On pourra considérer dans un premier temps µ = 0).
Exercice 3
Soient X et Y deux variables aléatoires indépendantes de loi uniforme standard U[0,1] . Déterminer les densités des variables
aléatoires suivantes U = X + Y , V = X − Y , Z = XY , T = X/Y (suggestion : on pourra commencer par déterminer les
densités des couples (U,V) et (Z,T)).
Exercice 4
Soient X et Y deux v.a.r. indépendantes de densités respectives
1 2
f (x) = √ 1(−1,1) (x) et g(y) = ye−y /2
1R∗+ (y).
π 1 − x2
Montrer que U = XY suit une loi centrée réduite. (On pourra considérer le changement de variables (x, y) → (u, v) =
(xy, y)).
On cherche à montrer que U = XY suit une loi βI (a, b + c). On rappelle également que
Z 1
Γ(a)Γ(b)
B(a, b) = xa−1 (1 − x)b−1 dx et B(a, b) = .
0 Γ(a + b)
Exercice 7
Soit X = (U, V ) un vecteur aléatoire de R2 de densité fX . On considère Y = (R, Θ) ses coordonnées polaires.
1. Montrer que Y admet une densité fY que l’on explicitera en fonction fX .
2. Application : Les véhicules spatiaux désirant s’arrimer à la Station Spatiale Internationale (ISS) s’appuient sur le
système de guidage GPS pour la phase d’approche de la station. Cependant, à faible distance de l’ISS, les signaux
émits par la constellation des satellites qui consistuent le système GPS sont fortement perturbés par des phénomènes
de réflexions multiples sur la structure métallique de la station. L’onde électromagnétique reçue par le récepteur GPS
du véhicule spatial se présente donc comme la superposition de deux ondes en quadrature dont les amplitudes U et
V sont des variables√aléatoires gaussienne N (0, σ 2 ) supposées indépendantes (pour des raisons d’isotropie). L’étude
de l’amplitude R = U 2 + V 2 de l’onde reçue est de première importance pour assurer un guidage fiable lors de la
manoeuvre d’arrimage.
(a) À l’aide de la question 1, expliciter la densité du couple (R, Θ).
(b) Montrer que R et Θ sont indépendantes.
(c) Reconnaı̂tre la loi de Θ.
(d) Donner la densité de R. Cette loi s’appelle loi de Rayleigh.
(e) Calculer l’espérance et la variance de R.
(U, V ).
1. Soit g : R2 → R une fonction borélienne positive. Exprimer E(g ◦ φ−1 (X, Y )) en fonction des densités de X et Y
notées fX et fY respectivement.
2. À l’aide du changement de variable (x, y) = φ(u, v), montrer que
e−u/2
Z
1
E(g(U, V )) = g(u, v) 1R+ (u) 1(0,2π) (v) dudv.
R2 2 2π
1. Déterminer K.
2. Déterminer les lois marginales.
3. Les variables X et Y sont-elles indépendantes ?
Exercice 10
Soit Y ∼ E(θ) et ε une variable aléatoire indépendante de loi de Rademacher symétrique, i.e. P(ε = 1) = 1−P(ε = −1) = p
avec p = 1/2. Montrer que la variable aléatoire Z = εY est à densité et la déterminer. Cette loi est appelée loi exponentielle
symétrique.
Exercice 11
Soit X une variable aléatoire de loi de Cauchy C(1) dont la densité est donnée par
1
f (x) = , x ∈ R.
π(1 + x2 )
Exercice 12
Soient X1 et X2 deux v.a.r. indépendantes de lois exponentielles de paramètres λ1 , λ2 > 0.
1. Déterminer les lois de m := min(X1 , X2 ) et M := max(X1 , X2 ).
2. Supposant que λ1 = λ2 montrer que m et M − m sont indépendantes.
Exercice 13
Un singe tape au hasard sur un clavier une lettre ou un symbole et recommence une infinité de fois de façon indépendante.
Montrer que, presque sûrement, il aura tapé en entier Les Misérables de Victor Hugo sans aucune faute de frappe. N’a-t-on
pas un résultat plus fort encore ?
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
TD2 : Indépendance
Exercice 1
R
1. La fonction fθ est positive, il s’agit de trouver cθ telle que f dλ = 1. D’où
Z θ
1 = cθ x dx = cθ θ2 /2 =⇒ cθ = 2/θ2 .
0
Pour t < 0, Gn (t) = 0 et pour t ≥ θ, Gn (t) = 1. Un bon candidat pour la densité de Yn est G0n (t) = nfθ F (t)n−1 =
2n−1
2n xθ2n 1[0,θ] (x). Cette fonction est positive et
θ 2n θ
x2n−1
Z
x
2 = 2n = 1.
0 θ2n 2nθ2n 0
Exercice 2
Il suffit de développer le carré :
n n n
1X 2 2 X 2 1X 2 2
Sn2 = Xi − X n Xi + X n = Xi − X n .
n i=1 n i=1
n i=1
On peut montrer (et on montrera, c’est une application de la LGN) que Sn2 converge presque sûrement vers σ 2 , cependant,
en moyenne Sn2 sous-estime la variance théorique : on dit que Sn2 est un estimateur consistant biaisé. Les logiciels donnent
n
en général la variance non biaisée n−1 Sn2 .
Exercice 3
Soit g une fonction mesurable bornée sur R2 . Notons que par indépendance, f(X,Y ) (x, y) = 1[0,1]2 (x, y). Calculons
Z 1 Z 1
E(g(X + Y, X − Y )) = g(x + y, x − y) dxdy.
0 0
On pose (u, v) = φ(x, y) = (x + y, x − y). C’est un changement de variable linéaire, inversible. L’application φ est en
particulier un C 1 -difféomorphisme de R2 . L’image de [0, 1]2 par φ est le carré ABCD où A = (0, 0), B = (1, −1),
C = (2, 0) et D = (1, 1). On notera ∆ ce carré. La formule du changement de variable donne :
Z Z
g(x + y, x − y) dxdy = g(u, v)|det Dφ−1 (u, v)| dudv.
[0,1]2 ∆
Or 1 1
−1 1
Dφ (u, v) = 2
1
2 =⇒ |det Dφ−1 (u, v)| =
2 − 21 2
La densité jointe de (U, V ) est donc f(U,V ) (u, v) = 21 1∆ .
On calcule la densité de U en intégrant f(U,V ) par rapport à v :
Z Z u Z 2−u
1 1 1
fU (u) = 1∆ (u, v)dv = 1[0,1] (u) dv + 1[1,2] (u) dv = u1[0,1] (u) + (2 − u)1[1,2] (u).
R 2 2 −u 2 −2+u
De façon similaire,
Z Z 2+v Z 2−v
1 1 1
fV (v) = 1∆ (u, v)du = 1[−1,0] (v) du + 1[0,1] (v) du = (1 + v)1[−1,0] (v) + (1 − v)1[0,1] (v).
R 2 2 −v 2 v
Exercice 4
Soit h : R2 → R+ mesurable positive, par indépendance de X et Y
Z
1 2
E(h(XY, Y )) = h(xy, y)1(−1,1) (x) √ 1R+ (y)ye−y /2 dxdy
2 π 1 − x2
ZR Z v
v 2
= h(u, v) √ e−v /2 dudv
2
π v −u 2
R∗
+ −v
Z Z ∞ 2
ve−v /2
= h(u, v) √ dvdu
R |u| π v 2 − u2
2 ∞ 2
e−u /2 e−u /2
Z
2
L(u) = e−t /2
dt = √ .
π 0 2π
Ainsi, XY suit une loi normale N (0, 1).
θa+b
f(X,Y ) (x, y) = xa−1 y b−1 e−θ(x+y) 1R∗+ ×R∗+ (x, y).
Γ(a)Γ(b)
On fait le changement de variable (u, v) = φ(x, y) = (x + y, x/(x + y)). On calcule φ−1 (u, v) = (uv, u(1 − v)) = (x, y).
On note que φ est un C 1 -difféomorphisme de R∗+ × R∗+ dans R∗+ × (0, 1) puis on calcule
v u
Dφ−1 (u, v) = =⇒ |det Dφ−1 (u, v)| = | − uv − u + uv| = u.
(1 − v) −u
Finalement,
1 ∞
θa+b
Z Z
E(g(U, V )) = g(u, v) (uv)a−1 ub−1 (1 − v)b−1 e−θu u dudv =
0 0 Γ(a)Γ(b)
Z 1Z ∞
θa+b
ua+b−1 e−θu v a−1 (1 − v)b−1 dudv
0 0 Γ(a)Γ(b)
Z 1Z ∞
θa+b 1
= ua+b−1 e−θu v a−1 (1 − v)b−1 dudv
0 0 Γ(a + b) B(a, b)
2. De la question précédente, on observe que la densité de (U, V ) est à variable séparées, elle s’écrit comme le produit
de deux densités unidimensionelles.
3. On reconnaı̂t la loi bêta pour V et la loi Γ(a + b, θ) pour U .
X Z V
4. On remarque que V = X+Y = Z+1 où de manière équivalente Z = 1−V . Il est clair que puisque U et V sont
indépendentes, U est indépendante de toute transformation de V d’où Z et U sont indépendantes. Si g est mesurable
bornée, nous avons en posant z = v/(1 − v), dv = dz/(z + 1)2
1 ∞
v a−1 (1 − v)b−1 z a−1
Z Z
E(g(Z)) = g(v/(1 − v)) dv = g(z) dz.
0 B(a, b) 0 (z + 1)a+b
Exercice 7
1. On fait le changement de variable (u, v) = (r cos(θ), r sin(θ)).
2. (a) C’est l’application de la question 1 pour (U, V ) un couple de gaussiennes indépendantes.
(b) Remarquer que la densié de (R, θ) est à variables séparées (s’écrit comme un produit de densité).
(c) C’est une loi uniforme sur [0, 2π).
(d) Soit on intègre la marginale, soit on utilise l’indépendance pour identifier la densité marginale.
(e) C’est des IPP en faisant apparaı̂tre à certains endroit l’intégrale d’une densité gaussienne.
2. On calcule √ !
cos(v)
√
2 u
− u sin(v) 1
Dφ(u, v) = sin √ =⇒ |det Dφ(u, v)| = .
√v u cos(v) 2
2 u
2
On obtient 1 − e−U /2
∼ U([0, 1]).
6. On est capable de simuler W, Z indépendantes et uniformes sur [0, 1], alors φ(−2 ln(1 − W ), 2πZ) est une gaussienne
bivariée.
D’où K = 1/π.
2. Nous avons : 2 2
e−y /2 e−y /2
Z Z
2
fY (y) = f (x, y) dx = e−2(x−3y/2) dx = √ .
R π R 2π
Un calcul très simllaire en choisissant l’autre décomposition de Gauss pour la forme quadratique conduit à la densité
de X.
3. Les variables X et Y sont clairement dépendantes car la densités du couple ne s’écrit pas comme le produit des
densités.
Exercice 10
Soit g une fonction mesurable bornée, alors
Z ∞ Z 0 Z
1
E(g(εY )) = P(ε = 1) g(y)θe−θy dy + P(ε = −1) g(y)θe−θ|y| dy = g(y) θe−θ|y| dy.
0 −∞ R 2
Exercice 11
1. Il s’agit de faire le changement de variable y = 1/x qui est un C 1 -difféormisphe de R∗ dans lui-même :
Z Z
1 1 1
g(1/x) dx = g(y) 2 dy.
R
2
π(1 + x ) R y π(1 + y −2
Ainsi, 1/X suit encore une loi de Cauchy.
2. On fait le changement de variable (u, v) = (y, y/z) qui est un C 1 -difféormorphisme de R × R∗ dans lui-même. On
calcule
−1 1 0
Dφ (u, v) = =⇒ |det Dφ−1 (u, v)| = |u|/v 2 .
1/v −u/v 2
La formule du changement de variable donne
2
/2−u2 /2v 2
e−u |u|
Z
E(g(Y, Y /Z)) = g(u, v) dudv.
R2 2π v2
Pour calculer la loi de Y /Z il ne nous reste plus qu’à intégrer la densité de la loi jointe par rapport à u :
Z ∞
1 u2 2 1
fY /Z = 2 ue− 2 (1−1/v ) du = .
v π 0 π(1 + v 2 )
Exercice 12
1. Calculons la loi du couple (X ∧ Y, X ∨ Y ). Soit g une fonction mesurable bornée :
Z ∞Z ∞
E(g(X ∧ Y, X ∨ Y )) = g(x ∧ y, x ∨ y)λ1 λ2 e−λ1 x−λ2 y dxdy
0 0
Z Z
−λ1 x−λ2 y
= g(x, y)λ1 λ2 e 10≤x≤y dxdy + g(y, x)λ1 λ2 e−λ1 x−λ2 y 10≤y≤x dxdy.
R2+ R2+
2. Calculons la loi de (m, M − m) en repartant de la densité jointe de (X, Y ) car nous ne connaissons pas la densité de
(m, M ) mais seulement celles de ses marginales :
Z Z
2 −λ(s+t)
E(g(m, M − m)) = g(s, t − s)210≤s≤t λ e dsdt = g(u, v)2λ2 e−2λu−λv dudv
R2 R2+
En posant
B = {∀k ∈ N∗ , (Xkn0 +1 , Xkn0 +2 , . . . , X(k+1)n0 ) 6= m}
et
B1 = {∀k ∈ {1, . . . , n}, (Xkn0 +1 , Xkn0 +2 , . . . , X(k+1)n0 ) 6= m},
nous obtenons, A{ ⊂ B ⊂ Bn . Nous avons évidemment P(A{ ) ≤ P(B) ≤ P(Bn ). Par indépendance,
n n0 n
Y
n 1
P(Bn ) = P((Xkn0 +1 , Xkn0 +2 , . . . , X(k+1)n0 ) 6= m) = P((Xkn0 +1 , Xkn0 +2 , . . . , X(k+1)n0 ) 6= m) = 1 − .
N
k=1
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
Exercice 1
Déterminer la fonction caractéristique de X dans les cas suivants :
1. PX = δa , a ∈ R ;
2. PX = 21 (δ−1 + δ1 ) ;
3. X ∼ B(n, p) ;
4. PX = n≥1 2−n δn ;
P
5. X ∼ P(Θ) ;
6. X ∼ U[−1,1] ;
7. X admet pour densité fX (x) = (1 − |x|)1[−1,1] , x ∈ R ;
8. pour tout B ∈ B(R), Z
1 1
PX (B) = 1B (0) + 1B (x) dx;
2 4 [−1,1]
9. X ∼ E(Θ).
Exercice 2
Soit X ∼ N (0, 1). Déterminer la fonction caractéristique de Y = X + = max(X, 0).
Exercice 3
Soient X1 et X2 deux v.a. indépendantes de loi B(n, p) et B(m, p) où p ∈ (0, 1). Montrer à l’aide des fonctions ca-
ractéristiques que X1 + X2 suit une binomiale B(n + m, p).
Exercice 7
Soient (Xk )k≥1 une suite de v.a. i.i.d. et N une v.a. à valeurs
Pdans N indépendante de la suite (Xk )k≥1 . On suppose que
n
E(X12 ) < ∞ et E(N 2 ) < ∞. Pour tout n ≥ 1, on note Sn = k=1 Xk .
1. Calculer E(SN ) et V(SN ).
2. Calculer la foncton caractéristique de SN et retrouver les résultats de la questions précédantes.
Exercice 8 (Extrait de l’examen de janvier 2017)
On considère X1 , X2 , X3 et X4 quatre variables aléatoires réelles indépendantes de même loi N (0, 1). On pose Y =
X1 X2 + X3 X4 .
1. Montrer que la fonction caractéristique de X1 X2 vaut
1
t → φX1 X2 (t) = E(eitX1 X2 ) = √ .
1 + t2
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
Exercice 1
On calcule facilement :
1. φX (t) = eta ;
2. φX (t) = cos(t) ;
3. φX (t) = (peit + (1 − p))n ;
eit
4. φX (t) = 2−eit ;
−λ(1−eit )
5. φX (t) = e ;
sin(t)
6. φX (t) = t ;
2(1−cos(t))
7. φX (t) = t2 ;
1 1 sin(t)
8. φX (t) = 2 + 2 t ;
θ
9. φX (t) = it−θ .
Exercice 2
On calcule
∞
Z r
1 1 2 itx−x2 /2
φX (t) = + e dx.
02 2 π
R∞q2 itx −x2 /2
R∞q2 itx −x2 /2
Notons φ(t) = 0 πe e dx. Alors φ est C 1 et sa dérivée est φ0 (t) = 0 π ixe e dx. Une intégration par
partie donne " #∞
∞
r Z r r
0 2 itx −x2 /2 2 itx −x2 /2 2
φ (t) = −i e e −t e e dx = i − tφ(t).
π 0 π π
0
Il s’agit donc de résoudre
Rt une EDO linéaire du premier ordre avec second membre constant. L’équation homogène a pour
2
solution φ(t) = ke− 0 s ds = ke−t /2 .
2
Pour l’équation avec second membre, on applique la méthode de la variation de la constante. On pose y(t) = k(t)e−t /2
puis on calcule
r r
2 0 0 −t2 /2 −t2 /2 −t2 /2 0 2 t2 /2
i = y (t) + ty(t) = k (t)e − tk(t)e + tk(t)e =⇒ k (t) = i e .
π π
Ainsi, la solution de l’EDO avec secod membre est de la forme
" r Z t #
2 2 2
φ(t) = e−t k + i ex /2 dx .
π 0
q R
2 2
2 t 2
Comme φ(0) = 1, on trouve k = 1 et donc φ(t) = e−t /2
+ ie−t π 0
ex /2
dx.
Exercice 3
On rappelle que φX1 (t) = (1 − p + peit )n et que les fonctions caractéristiques caractérisent la loi. Par indépendance,
où l’on reconnait la fonction caractéristique d’une B(n + m, p). Remarquons qu’il est primordiale que le paramètre p soit
identique pour les deux variables aléatoires.
Exercice 7
1. On calcule par Fubini (on utilise le fait que X1 , N ∈ L1 ) et en utilisant l’indépendance de N et (Xk )k≥1 .
∞ ∞ ∞
!
X X X
E(SN ) = E 1N =n Sn = E(1N =n Sn ) = P(N = n)nE(X1 ) = E(N )E(X1 ).
n=0 n=0 n=0
V(Sn ) = E(Sn2 ) − E(Sn )2 = E(N )E(X12 ) + E(X1 )2 E[N (N − 1)] − E(N )2 E(X1 )2
= E(N )V(X1 ) + E(X1 )2 V(N ).
2. On note φ la fonction caracteristique de X alors par Fubini et indépendance
∞
X
ψ(t) = E[eitSN ] = P(N = n)φ(t)n = GN ◦ φ(t),
n=0
où G(s) = E(s ) est la fonction génératrice de N . Remarquons que ψ est C 2 comme la composée de deux fonctions
N
3. On calcule
0 ∞
1 eitx+x 1 eitx−x
Z
1 1 1 1 1 1
φZ (t) = eitx e−|x| dx = + = − =
2 R 2 it + 1 −∞ 2 it − 1 0 2 it + 1 2 it − 1 1 + t2
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
Exercice 2
Soit (Un )n≥0 une suite de v.a.r. i.i.d. de loi normale centrée et de variance σ 2 > 0. Pour tout θ ∈ R, on définit la suite
(Xn )n≥1 par Xn = θUn−1 + Un pour tout n ≥ 2 et X1 = U1 .
Montrer que pour tout n ≥ 1, X = (X1 , . . . , Xn )∗ est une vecteur gaussien non dégénéré dont on précisera la densité,
l’espérance et la matrice de covariance.
Exercice 3
On considère X ∼ N (0, 1) et une v.a. ε indépendante de X et de loi 12 (δ−1 + δ1 ).
1. Montrer que la v.a. Y = εX est gaussienne et que les v.a. X et Y sont non corrélées.
2. Les variables X et Y sont-elles indépendantes ?
3. Le couple (X, Y ) est-il gaussien ?
Exercice 4
(Extrait du partiel de Décembre 2013)]
La fonction caractéristique d’une loi du Chi-deux à 1 degré de liberté est donnée pour t ∈ R par ψ(t) = 1/(1 − 2it)1/2 .
On considère trois v.a. indépendantes X, Y et Z de loi N (0, 1).
1. Donner la loi de la variable aléatoire U = X + Y + Z.
2. Montrer que (U, X − Y ) est un vecteur gaussien dont on précisera l’espérance et la matrice de covariance. En déduire
que U est indépendante des variables X − Y , Y − Z et Z − X.
3. Montrer que les variables aléatoires X + Y + Z, 2X − Y − Z et Y − Z sont indépendantes. Identifier les lois de
(2X − Y − Z)2 /6 et de (Y − Z)2 /2.
4. Vérifier l’égalité
1 3
(x − y)2 + (y − z)2 + (z − x)2 = (2x − y − z)2 + (y − z)2
2 2
pour tout x, y, z ∈ R et en déduire l’expression de la fonction caractéristique de V = (X − Y )2 + (Y − Z)2 + (Z − X)2 .
5. Les variables (X − Y )2 , (Y − Z)2 et (Z − X)2 sont-elles indépendantes ?
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
d’où Cov(X, Y ) = 1.
3. La matrice de covariance s’écrit donc
4 1
Σ= .
1 1
4. Le vecteur (X + Y, 2X − Y ) est gaussien par transformation linéaire du vecteur gaussien (X, Y ). Plus précisément :
X +Y 1 1 X
= .
2X − Y 2 −1 Y
| {z }
A
Exercice 2
En écrivant
··· ···
1 0 0
.. ..
X1 θ
1 . . U1
.. .. .. .. ..
=
. . . 1 . 0 . .
Xn . .. .. Un
.. . . 0
0 ··· 0 θ 1
Ceci montre que (X1 , . . . , Xn ) est un vecteur gaussien d’espérance nulle et de matrice de covariance :
On vérifie facilement que det Σ = det AA∗ σ 2 > 0 et donc que la loi gaussienne est non P dégénérée. De même, Σ−1 =
−2 −1 ∗ −1 − −1 −1 n−1
σ (A ) A . Pour une formule explicite de A 1 on peut écrire A = (I + θE) = k=0 (−1)k θk E k où E est la
matrice avec des 1 sous la diagonale et des zéros partout ailleurs. Nous vérifions que cette matrice rest nilpotente d’ordre
(k) (k)
n et plus exactement E k = (ei,j ) avec ei,j = 0 sauf pour j = k + i pour lequel cela vaut 1.
Exercice 3
1. Soit g une fonction continue bornée alors par indépendance
1 1
E(g(Y )) = E(1ε=−1 g(−X) + 1ε=1 g(X)) = E(g(−X)) + E(g(X)) = E(g(X)),
2 2
car X et −X ont même loi. Ainsi, Y ∼ N (0, 1). On calcule
1 1
Cov (X, Y ) = E(XεX) = E(−X 2 ) + E(X 2 ) = 0.
2 2
2. Si X et Y étaient indépendantes, alors le couple (X, Y ) serait gaussien de moyenne nulle et de matrice de covariance
l’identité. Par conséquent,
2 2
φ(X,Y ) (s, t) = e−t /2 e−s /2 .
Or, un calcul très similaire à celui de la question précédante implique
2 2
itX+isεY e−(t−s) /2
+ e−(t+s) /2
φ(X,Y ) (s, t) = E(e )= .
2
En spécifiant s = t = 1, on obtient e−1 d’une part et (1 + e−2 )/2 d’autre part. C’est une contradiction.
3. Si (X, Y ) est gaussien alors X + Y est une variable aléatoire gaussienne. Or X + Y = (1 + ε)X et on constate que
P(X(1 + ε) = 0) = 21 . Contradiction. Remarque importante : les notions de décorrélation et d’indépendance ne
sont équivalentes que dans le cas où (X, Y ) est un vecteur gaussien ; la normalité des marginales est insuffisante.
Calculons
1 1 1 1 2 0 3 0 0
BB ∗ = 2 −1 −1 1 −1 1 = 0 6 0 .
0 1 −1 1 −1 −1 0 0 2
Là encore, la matrice de covariance est diagonale et (X, Y, Z) est gaussien donc les trois variables considérées sont
indépendantes.
On remarque que en lisant la matrice de covariance que 2X−Y √ −Z ∼ N (0, 1) demême que Y√ −Z
∼ N (0, 1). Leurs
6 2
2
carrés suivent une loi du χ à un degré de liberté.
4. En utilisant l’égalité et l’indépendance de 2X − Y − Z avec Y − Z :
1
φV (t) = ψ(3t)ψ(3t) =
1 − 6it
2 2 2
5. Etant donné que X − Y suit une N (0, 2), on déduit que (X−Y 2
)
∼ χ2 (1). Il en va de même pour (Y −Z)
2 et (Z−X)
2 .
1 2
Supposons l’indépendance des trois variables aléatoires alors 2 V est la somme de trois χ (1) indépendantes si bien
que φ 12 V (t) = (1−2it)−3/2 or φ 12 V (t) = φ(t/2) = (1−3it)−1 par la question précédante. On obtient une contradiction.
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
TD5 : Convergences
Exercice 1
Soit (Xn )n≥1 une suite de v.a. de loi
1 1
PXn = 1− δ0 + (δ−n + δn ).
n2 2n2
1. Montrer que (Xn )n≥1 converge presque sûrement vers 0.
2. Que dire de la convergence en moyenne quadratique ?
Exercice 2
1. Soit (Xn )n≥1 une suite de v.a. i.i.d. positives admettant un premier moment. On note E(X1 ) = m < ∞ et on pose
Un = Xn /n2 , n ≥ 1. La suite (Un )n≥1 converge-t-elle presque sûrement ? dans L1 ? En probabilité ? Si oui, vers
quelle limite ?
2. On suppose désormais que les v.a. Xn admettent une densité
c
f (x) = 1x>0 , x ∈ R.
π(x2 + 1)
(a) Montrer que Un et Vn suivent des lois binomiales dont on précisera la valeur des paramètres. Les v.a. Un et Vn
sont-elles indépendantes.
(b) Montrer que (Un − Vn )/n tends presque sûrement vers une valeur déterministe que l’on précisera.
(c) Calculer E(Un ), E(Vn ) et en déduire la limite de E((Vn − Un )/n) quand n → ∞. Ce résultat pouvait-il se
déduire de la question précédente ?
(d) Calculer E(Un Vn ) et Cov(Un , Vn ).
(e) Montrer que V((Un − Vn )/n) tend vers 0 lorsque n → ∞. Ce résultat implique-t-il celui de la question 3b ? En
implique-t-il une forme plus faible (que l’on précisera alors) ?
(f) Pour n grand, construire un intervalle de confiance pour θ de niveau de risque 5% à partir des Xi , 1 ≤ i ≤ n,
et un autre à partir des Ui , 1 ≤ i ≤ n. Commenter. Quelle est l’utilité de supposer que n est grand ?
Exercice 5
Sur l’espace de Borel standard ([0, 1], B([0, 1]), λ), on considère la suite (Xn )n≥1 de v.a. définies pour ω ∈ [0, 1] par
1
Xn (ω) = √ 1(0,1/n) (ω).
ω
Montrer que (Xn )n≥1 converge en probabilité vers 0 mais ne tend pas vers 0 en moyenne quadratique.
3. (a) Montrer que, pour n donné, la famille d’événements At = {N (t) < n} est décroissante en t ≥ 0 pour l’inclusion,
i.e. s ≤ t implique At ⊂ As .
(b) En déduire que si P(At ) → 0 lorsque t → ∞ alors P(∩t≥0 At ) = 0.
4. On suppose désormais que E(X1 ) = m < ∞.
(a) Montrer que pour tout n entier, Tn admet un premier moment et en déduire que P(Tn > t) → 0 lorsque t → ∞.
Que peut-on en déduire sur le nombre de pannes quand t tend vers l’infini.
(b) Montrer que Tn /n converge presque-sûrement vers une valeur à déterminer.
(c) En déduire que la suite TN (t) /N (t) converge presque-sûrement vers m lorsque t tend vers l’infini (on pourra se
contenter d’une heuristique).
(d) Montrer que pour t ∈ R+
TN (t) t TN (t)+1 1
≤ < 1+ .
N (t) N (t) N (t) + 1 N (t)
(e) En déduire que t/N (t) converge presque sûrement vers m et interpéter ce résultat dans le contexte de l’énoncé.
Exercice 7
Soit (Xn )n≥1 une suite de v.a. i.i.d. de carré intégrable. On pose m = E(X1 ) et σ 2 = V(X1 ). Pour tout entier n ≥ 1, on
pose
n n
1X 1X
Xn = Xk et Zn = (Xk − X n )2 .
n n
k=1 k=1
1. Calculer E(Zn ) ((Ce calcul a déjà été effectué dans la feuille de TD 2.).
2. Montrer que Zn converge presque sûrement vers σ 2 .
Exercice 8
On considère une suite de v.a. indépendantes (Xn )n≥1 telles que, pour n ≥ 1, P(Xn = n) = 1/n = 1 − P(Xn = 0).
1. Montrer que (Xn )n≥1 converge vers 0 en probabilité.
Pn
2. On pose Yn = n−1 k=1 Xk . La suite (Yn )n≥1 converge-t-elle vers 0 en probabilité ?
3. En utilisant le lemme de Cesàro, montrer que (Xn )n≥1 ne peut converger presque sûrement.
4. Que dire de l’application du lemme de Cesàro à la convergence en probabilité ?
Exercice 9 (Extrait du rattrapage de février 2017)
Soit (Un )n≥1 (respectivement (Yn )n≥1 ) une suite de variables aléatoires i.i.d. de loi commune la loi uniforme sur [0, 1]
(respectivement de loi de Poisson de paramètre α > 0).
1. La loi faible des grands nombres L2 s’applique-t-elle à la suite (Un )n≥1 ? Justifier.
2. Soit f : R → R une fonction continue bornée. Calculer la limite
x1 + · · · + xn
Z
lim f dx1 . . . dxn .
n→∞ [0,1]n n
On pourra montrer que cette intégrale est l’espérance d’une fonctionnelle de la suite (Un )n≥1 .
Pn
3. Montrer que Zn = k=1 Yk suit une loi de Poisson de paramètre nα.
4. Montrer que (Zn /n)n≥1 converge en probabilité vers une limite que l’on précisera.
5. En déduire, pour toute fonction f : R → R continue et bornée, la limite
X (αn)k
lim e−αn f (k/n).
n→∞ k!
k≥0
Exercice 10
Soient X1 , . . . , Xn des variables aléatoires i.i.d. de loi commune µ. On note Fn la fonction de répartition empirique associée
et définie pour n ≥ 1 et t ∈ R par
n
1X
Fn (t) = 1[Xi ,∞) (t).
n i=1
Montrer que si F est la fonction de répartition de la loi µ, alors limn→∞ Fn (t) = F (t) presque sûrement pour tout t ∈ R.
Exercice 12
1. On considère deux v.a.r. i.i.d. X et Y définies sur (Ω, F, P). On pose Z = X − Y et on note φX (resp. φZ ) la
fonction caractéristique de X (resp. Z).
(a) Montrer que PZ = P−Z et que φZ (t) = |φX (t)|2 pour tout t ∈ R.
(b) Montrer que Z ∈ L2 si et seulement si X ∈ L2 .
2. Pour n ∈ N∗ , on définit la fonction fn par
(a) Montrer que fn est la densité d’une loi de probabilité Pn ; calculer l’espérance et la variance de Pn et calculer
φn sa fonction caractéristique.
(b) Montrer que si (Xn )n≥0 est une suite de v.a. de loi Pn , alors la suite (Xn )n≥0 convege en loi vers la loi uniforme
sur [−1, 1].
(c) A-t-on convergence de la suite (fn )n≥0 ? de la suite des variances des Xn ?
(d) Montrer que si on suppose en outre l’indépendance des Xn alors la moyenne empirique converge en probabilité
vers 0.
(e) Trouver un exemple de loi symétrique µ pour laquelle il n’existe pas de couple (X, Y ) de v.a. i.i.d. telles que
X − Y soit de loi µ.
est de loi µ.
(c) Montrer que la suite Yn converge en loi vers la loi normale centrée réduite.
(d) Identifier µ.
Partie I
Rappel :
Une fonction f : Rd → Rd est lipschitzienne si il existe une constante k ≥ 0 telle que kf (x) − f (y)k ≤ kkx − yk.
On munit Rd de la norme euclidienne notée k · k. Soient X un ensemble et (fε )ε∈X une famille de fonction lipschitzienne
de Rd dans Rd . Pour chaque ε ∈ X, on note cε la constante de Lipschitz, i.e. :
kfε (x) − fε (y)k
cε := sup = inf k ≥ 0 : ∀(x, y) ∈ Rd × Rd , kf (x) − f (y)k ≤ kkx − yk .
x,y∈Rd :x6=y kx − yk
Soit X0 ∈ Rd de loi µ et (εn )n≥0 une suite de variables aléatoires à valeurs dans X indépendantes et identiquement
distribuées. On suppose en outre que (εn )n≥0 est indépendante de X0 . On définit la suite de vecteurs aléatoires (Xn )n≥0
par récurrence :
X0 ∼ µ et pour tout n ≥ 0 : Xn+1 = fεn (Xn ).
1. Pourquoi Xn+1 = fε ◦ fε
n
◦ · · · ◦ fε (X0 ) est-elle de même loi que X
n−1 0
en+1 = fε ◦ fε ◦ · · · ◦ fε (X0 ) ?
0 1 n
(b) En utilisant le critère de Cauchy pour les séries numériques, montrer, sous la condition E(log(cε0 )) < 0, que
P∞ Qk
la série k=0 `=0 cε` est convergente presque-sûrement et en probabilité (on pourra utiliser le résultat de la
question 1. de la partie I).
(c) En déduire que (Xen )n≥0 converge presque-sûrement. On notera X e∞ la limite.
4. En supposant les conditions de la question 3 satisfaites, montrer que (Xn )n≥0 converge en loi vers la loi de Xe∞ .
5. Parmi les conditions E(cε0 ) < 1 et E(log(cε0 )) < 0, laquelle est la plus faible ?
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
TD5 : Convergences
Exercice 1
1
P
1. Soit ε > 0 alors P(|Xn | ≥ ε) = n2 ainsi n≥1 P(|Xn | > ε) < ∞ et Xn converge presque sûrement vers 0 par
Borel-Cantelli.
2. On aE(Xn2 ) = 1 donc Xn ne converge pas vers 0 dans L2 .
Exercice 2
1. Soit ε > 0, l’inégalité de Markov implique
X X m
P(|Un | > ε) ≤ <∞
n2 ε
n≥1 n≥1
donc Un converge presque sûrement vers 0, a fortiori elle converge en probabilité. Enfin, E|Un | = EX1 /n2 qui tends
vers 0 quand n tends vers l’infini. Ainsi Un converge dans L1 vers 0.
R ∞ dx 1
2. (a) On calcule facilement 0 π(1+x 2 ) = 2 d’où c = 2.
(b) Pas tout à fait car si X1 est effectivement positive nous avons par contre EX1 = ∞.
(c) Clairement E|Un | = ∞ pour tout n ≥ 1 donc Un ne converge pas vers 0 dans L1 . D’autre part,
Z ∞
2dx 2 2 1
P(|Un | ≥ ε) = P(X1 ≥ n2 ε) = 2
= 1 − arctan n2 ε = arctan 2 .
n2 ε π(1 + x ) π π n ε
3. (a) Les variables aléatoires Un et Vn sont des sommes de Bernoulli indépendantes aussi Un ∼ B(n, P(X ≤ 0)) et
Vn ∼ B(n, P(X > 0)). On calcule
1−θ 1+θ
p = P(X ≤ 0) = , et 1 − p = P(X > 0) = .
2 2
(b) Par la loi des grands nombres Un /n et Vn /n converge respectivement vers p et 1 − p si bien que (Un − Vn )/n
converge presque sûrement vers 2p − 1 = −θ.
(c) Nous avons E(Un ) = np et E(Vn ) = n(1 − p) si bien que E((Vn − Un )/n) tends vers 1 − 2p = θ lorsque n → ∞.
En appliquant le théorème de convergence dominée, la convergence considérée pouvait se déduire de la question
précédante.
(d) Comme Un + Vn = n, E(Un Vn ) = E(nUn ) − E(Un2 ) = n2 p − np(1 − p) − n2 p2 = n(n − 1)p(1 − p). On déduit
(e) On calcule
1
V((Un − Vn )/n) = [Cov(Un − Vn , Un − Vn ) = V(Un ) + V(Vn ) − 2Cov(Un , Vn )]
n2
2np(1 − p) + 2np(1 − p) 4p(1 − p)
= = −→n→∞ 0.
n2 n
Nous avons lim E(Vn − Un )/n = θ et lim V((Vn − Un )/n) = 0, on en déduit que (Vn − Un )/n converge en
probabilité vers θ.
(f) Nous avons calculé précédemment que E(X) = θ/4 or un estimateur consistant de la moyenne est la moyenne
empirique. On peut donc considérer
n
4X
θ̂ = Xi .
n
k=1
Alors θ̂ converge presque sûrement vers θ par la loi des grands nombres et le TCL implique :
√ θ̂ − θ
n =⇒ N (0, 1).
4σ
Le problème est qu’on ne connaı̂t pas σ, l’écart-type de X, c’est en réalité une fonction de θ le paramètre
inconnu. On remplace donc σ par son estimation σ̂ :
v
u n
uX
σ̂ = t (Xi − X n )2 .
k=1
Il est facile de voir en développant le carré que σ̂ converge presque sûrement vers σ. Puis le lemme de Slutsky
implique que
√ θ̂ − θ √ θ̂ − θ σ
n = n =⇒ N (0, 1).
4σ̂ 4σ σ̂
Ceci va permettre de construire un intervalle de confiance pour θ. On note t0 = 1, 96 alors
!
√ θ̂ − θ
4t0 σ̂ 4t0 σ̂
0.95 = P −t0 ≤ n ≤ t0 = P θ̂ − √ ≤ θ ≤ θ̂ + √ .
4σ̂ n n
1−θ
Dans les questions précédentes, nous avons montré que Un /n convergeait vers p = 2 . En ce sens,
2Un V n − Un
θ̃ = 1 − = .
n n
est un estimateur consistant de θ. Encore une fois,
Un
√ n −p
n =⇒ N (0, 1),
p(1 − p)
si bien que par la delta méthode
√ f (Un /n) − f (p) √ θ̃ − θ
n 0 = − n (1−θ)2 (1+θ) =⇒ N (0, 1),
f (p)p(1 − p)
4
où f (x) = 1 − 2x. Là encore, le lemme de Slutsky permet de remplacer θ par θ̃ au numérateur. D’où
!
t0 (1 − θ̃)2 (1 + θ̃) t0 (1 − θ̃)2 (1 + θ̃)
0.95 = 0.95 = P θ̃ − √ ≤ θ ≤ θ̃ + .
4 n 4
4. (a) Puisque Tn est positive, on peut calculer E(Tn ) = nm ce qui montre que Tn admet un premier moment. Par
l’inégalité de Markov,
nm
P(Tn > t) ≤ −→t→∞ .
t
De là, avec la question 1., on obtient P(N (t) < n) = P(Tn > t) d’où N (t) → ∞ lorsque t → ∞.
(b) La loi des grands nombres implique que Tn /m converge vers m presque sûrement.
TN (t)
(c) Comme N (t) tends vers l’infini, N (t) converge presque sûrement vers m par composition de limite.
(d) Ce n’est ni plus ni moins que l’inéquation de la question 2. renormalisée par N (t).
(e) Par encadrement, on obtient facilement que t/N (t) converge presque sûrement vers m lorsque t tends vers
l’infini. On a montré que le nombre moyen de pannes sur un long intervalle [0, t] est de l’ordre 1/m, c’est
inversement proportionnel au temps entre deux pannes.
Exercice 7
1. Voir le calcul effectué dans la feuille de TD 2.
2. De même, nous avons déjà montré que
n
1X 2
Zn = Xk − (X n )2 .
n
k=1
En appliquant la loi des grands nombres aux deux termes, on trouve le résultat souhaité.
Exercice 8
1. Soit ε > 0, il s’agit de calculer P(|Xn | > ε) = 1/n qui tends vers 0. D’où la convergence en probabilité.
2. Si Yn tend vers 0 en probabilité alors Yn converge vers δ0 en loi dont la fonction caractéristique est ϕδ0 (t) = 1
pour tout t ∈ R. Nous allons montrer qu’il existe t0 ∈ R tel que lim supn→∞ |ϕYn (t0 )| < 1 ce qui nous donnera la
contradiction recherchée.
À l’aide de l’indépendance des Xn , on calcule
n n n n
1 eitk/n 1 − eitk/n
h i Y Y Y
itYn Y itXk /n
ϕYn (t) = E e = E e = 1− + = 1− = zk,n (t).
k k k
k=1 k=1 k=1 k=1
Par ailleurs,
1 − eitk/n 1 − e−itk/n
2
|zk,n (t)| = 1 − 1−
k k
2(1 − cos(tk/n)) 2(1 − cos(tk/n)) 2(k − 1)
=1− + =1− (1 − cos(tk/n).
k k2 k2
De là, on considère t ∈ [−π/2; π/2] et on utilise les inégalités de convexités ln(1 + x) ≤ x pour x > −1 et cos(x) − 1 ≤
− π2 x pour x ∈ [0, π/2]. On obtient
( n 1/2 )
X 2(k − 1)
|ϕYn (t)| = exp log 1 − (1 − cos(tk/n)
k2
k=1
( n ) ( n
)
X k−1 2 1 X k−1
≤ exp (cos(tk/n) − 1) ≤ exp − |t| .
k2 π n k
k=1 k=1
Ceci montre que (Yn ) converge presque sûrement vers 0 donc en probabilité vers 0. Contradiction avec la question
précédente.
4. Il n’est plus valide.
Exercice 10
Il s’agit de la version simple de Glivenko-Cantelli, le vrai théorème de Glivenko-Cantelli est montrée dans le cours : pour
tout t ∈ R on veut montrer que Fn (t) converge presque sûrement vers F (t). Dans ce cas, Fn (t) converge presque-sûrement
par la loi des grands nombres vers E(1[X1 ,∞) (t)) = P(X1 ≤ t). Ceci signifie que pour tout t ∈ R, il existe un ensemble
négilgeable Nt tel que ω ∈/ Nt , limn→∞ Fn (t, ω) → F (t).
Nous avons σn Y + mn converge presque sûrement vers σY + m. Comme g est continue, g(σn Y + mn ) converge
presque sûrement vers g(σY + m) et supn≥0 |g(σn Y + mn )| ≤ kgk∞ . Le théorème de convergence dominée implique
Exercice 12
1. (a) Par hypothèse PX ⊗ PY = PY ⊗ PX si bien que L(X, Y ) = L(Y, X). En notant p : R2 → R définie par
p(x, y) = x − y, il vient alors L(p(X, Y )) = L(p(Y, X)). D’où PZ = P−Z .
Par indépendance de X et Y , nous avons φZ = φX φ−X = φX φX = |φX |2 .
(b) On suppose que X ∈ L2 alors par Cauchy-Schwartz X ∈ L1 et
Ainsi, Z ∈ L2 .
Réciproquement, supposons Z ∈ L2 alors par le théorème de transfert et le théorème de Fubini pour les fonctions
mesurables positives Z Z Z
Z 2 dP = (x − y)2 PX (dx)PY (dy) < ∞.
Par conséquent, l’inégalité de Markov implique que (X − y) ∈ L2 pour PY -presque tout y ∈ R. Il est alors clair
que X ∈ L2 — il suffit pour cela que X − y ∈ L2 pour un seul y ∈ R.
2. (a) Il est clair que fn est positive pour tout n ≥ 1. Puis, on calcule
Z 1 Z 1
1
sin2 (2πnx) dx = 1 − cos(4πnx) dx = 1.
−1 2 −1
Soit Xn une variable aléatoire de densité fn . Notons que |Xn | ≤ 1 presque sûrement, elle admet donc des
moments de tout ordre. De plus, par parité
Z 1
E(Xn ) = x sin2 (2πnx) dx = 0.
1
Finalement,
1 1
V(Xn ) = − .
3 8π 2 n2
Enfin, calculons la fonction caractéristique :
1 1 itx
Z
sin(t)
φ(t) = e dx = .
2 −1 t
On constate donc que limn→∞ φn (t) = φ(t) pour tout t ∈ R. Ceci montre que Xn converge en loi vers la loi
uniforme sur [−1, 1].
(c) La fonction y → sin2 (2πy) est 12 -périodique. Par ailleurs, la suite nx mod 12 est dense dans [0, 12 ) dès que x ∈
/ Q 1.
Ainsi, dès que x ∈/ Q, l’ensemble des valeurs d’adhérences de fn (x) est [0, 1]. Si x ∈ Q, la suite nx mod 12 est
périodique et l’ensemble des valeurs d’adhérence de fn (x) est fini. En fait on voit facilement que fn (x) converge
uniquement pour x = 0 mod 12 , seul cas où fn (x) ne possède qu’une valeur d’adhérence.
De plus limn→∞ V(Xn ) = 31 . Notons que cette variance limite n’est rien d’autre que la variance de la loi
uniforme sur [−1, 1].
(d) L’indépendance des Xn et le fait que Xn est de carré intégrable pour tout n implique par Bienaymé-Tchebichev
que, pour tout ε > 0, n ! P
n
1 X V(Xk ) 1
P Xk > ε ≤ k=12 2 ≤ → 0.
n n ε 3nε2
k=1
(e) Une telle loi µ symétrique admet une fonction caractéristique positive. Or, si µ est la loi uniforme sur [−1, 1],
alors µ̂(t) = sin(t)
t . En particulier, µ̂ est négative sur (−π, 0), c’est une contradiction avec la question 1.
(c) Dire le maximum dépasse un certain seuil signifie qu’il existe un Xi qui le dépasse d’où le résultat.
1 1 1
P max Xk > + δ ≤ nP X1 > + δ log n .
log n 1≤k≤n λ λ
2. (a) La fonction G est continue et croissante comme la composée de deux fonctions continues et décroissantes. De
plus, limt→∞ G(t) = 1 et limx→−∞ G(t) = 0.
(b) Soit t ∈ R alors
n
log n
= 1 − e−λ(t+log(n)/λ) ∼ exp −e−λt ne− log n → G(t).
lim P max Xk − ≤t
n→∞ 1≤k≤n λ
(c) Cela montre la convergence vers la la loi de Gambel puisque la suite de fonction de répartition convergence
ponctuellement vers G pour tout t ∈ R.
Des résultats de cet exercice, on obtient un nouvel estimateur du paramètre 1/λ de la loi exponentielle. Remarquons que
la largeur de l’IC associé est de l’ordre log n.
Partie I
1. Puisque les εn sont i.i.d. et indépendante de X0 , nous avons que L(ε0 , . . . , εn , X0 ) = L(εn , . . . , ε0 , X0 ). D’où L(Xn ) =
L(X̃n ).
2. (a) Montrer que pour tout k ≥ 0 :
n−1
Y
kX
ek+1 − X
ek k = kfε ◦ · · · ◦ fε (X0 ) − fε ◦ · · · ◦ fε (X0 )k ≤
0 k 0 k−1
εk kfεk (X0 ) − X0 k.
k=0
(b) Soient n, p ≥ 0
n+p−1
X n+p−1
X k−1Y
kX
en+p − X
en k ≤ kX̃k+1 − X̃k k ≤ cε` kfεk (X0 ) − X0 k.
k=n k=n `=0
(c) On calcule
n+p−1
X k−1Y n+p−1
X
EkX̃n+p − X̃n k ≤ cε` kfεk (X0 ) − X0 k = E(cε0 )k E[kfε0 (X0 ) − X0 k],
k=n `=0 k=n
en )) −→n→∞ E(h(X
E(h(Xn )) = E(h(X e∞ )),
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
Exercice 5
Soient X1 , . . . , Xn des variables aléatoires indépendantes de loi de Poisson. Déterminer la loi de Xj sachant X1 + · · · + Xn .
Exercice 6
Soit (X, Y ) une variable aléatoire à valeurs dans R2 de carré intégrable et satisfaisant
Exercice 7
Soit X une variable aléatoire exponentielle de paramètre θ > 0. Déterminer E(X|X > t) pour tout t ∈ R.
Exercice 8
Soient X, Y, Z trois v.a.r. telles que
— X est une loi uniforme sur l’intervalle [0, 1] ;
— la densité de Y conditionnellement à X est donnée par
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.
1ère année Probabilités générales - Cursus Mathématiques ENSAI
d f (y − 1) + f (y + 1)
g(y) = P(Y ≤ y) = .
dy 2
R
Il est immédiat que g est positive et que R
g(y) dy = 1 par invariance par translation de la mesure de Lebesgue.
D’où l’on déduite y = −z et x = 2z, ainsi Ker Σ = Vect (2, −1, 1).
(b) Comme Σ n’est pas inversible, X n’est pas à densité par rapport à la mesure de Lebesgue sur R3 .
(c) La variable X est supportée par m + (Ker Σ)⊥ .
2. On remarque que V(X1 ) = 2 6= 0 puis on calcule
Cov(X2 , X1 )
a= =1
V(X1 )
Ainsi,
1 1
f(X1 ,X2 ) (x, y) = √ exp{− [5(x − 1)2 + 2y 2 − 4xy + 4y]}.
2π 6 12
Finallement,
√
f(X,Y ) (x, y) 4π 1 1 1 [y−(x−1)]2
fX2 |X1 =x (y) = = √ exp{− [5(x − 1)2 + 2y 2 − 4xy + 4y] + (x − 1)2 ]} = √ e− 6 .
fX1 (x) 2π 6 12 4 6π
On reconnaı̂t la densité d’une N (x − 1, 3).
Une autre méthode, plus éléguante, qui s’inspire de la preuve du théorème sur le conditionnement des vecteurs
gaussiens : soit a ∈ R et posons Z = X2 − E(X2 ) − a(X1 − E(X1 )) = X2 − aX1 + a. La variable Z est gaussienne et
centrée. Il se trouve que l’on peut choisir a tel Z et X1 soit indépendante. En effet, 0 = Cov (Z, X1 ) = Cov (X2 , X1 )−
aV(X1 ). Ainsi la condition d’indépendance est valide pour a = CovV(X (X1 ,X2 )
1)
= 1. Nous obtenons donc X2 = X1 −1+Z.
Remarquons V(Z) = V(X1 ) + V(X2 ) − 2Cov (X1 , X2 ) = 3. Considérons g une fonction continue bornée, alors par
indépendance de Z et X1
avec σ 2 = 4/5. Ainsi, la loi conditionelle de X|X + 2Y est une gaussienne de moyenne X + 2Y et de variance
4/5.
Exercice 5
On pose
n
X X
Sn = Xj et Sn(j) = Xk .
j=1 k6=j
Pn (j)
λj , ainsi Sn ∼ P(λ) et λ(j) = λk ainsi Sn ∼ P(λ(j) ). On cherche à calculer
P
De même on pose λ = j=1 k6=j
(j)
P(Xj = `, Sn = k) P(Sn = k − `)P(Xj = `)
P(Xj = `|Sn = k) = = .
P(Sn = k) P(Sn = k)
A l’aide de la question 1., on remplace
(j)
e−λ e−λj (λ(j) )k−` λ`j k
P(Xj = `|Sn = k) = .
e−λ λk `
Il suffit de remarquer que λk = λk−` λ` , de poser p = λj /λ et de constater que 1 − p = λ(j) /λ. Ainsi :
k `
P(Xj = `|Sn = k) = p (1 − p)k−` .
`
On conclut que
λj
L(Xj |X1 + · · · + Xn ) = B X1 + · · · + Xn , .
λ1 + · · · + λn
Exercice 6
1. La formule des probabilités totales donnent
Exercice 7
Comme P(X > t) = 1t<0 + 1t≥0 e−θt pour tout t > 0. Il vient pour t < 0, E(X|X > t) = E(X) = 1/θ et pour t ≥ 0
Z ∞ Z ∞
E(X1(t,∞) (X)) −θx −θx ∞
E(X|X > t) = =e θt
xθe θt
dx = e [−xe ]t + e θt
e−θx dx = t + 1/θ.
e−θt t t
Exercice 8
1. La densité de la loi jointe de (X, Y, Z) est donnée par
f(X,Y,Z) (x, y, z) = fZ|X=x,y=y (z)fY |X=x (y)fX (x) = 1D (x, y)1R∗+ (z)(y − x)2 e−(z+1)(y−x)
2. À l’aide de la densité conditionelle, on remarque que L(Z|X, Y ) = E(Y − X). On calcule la loi de Z en intégrant la
loi jointe par rapport aux marginales X et Y :
Z Z 1 Z ∞
fZ (z) = f(X,Y,Z) (x, y, z) dxdy = 1R∗+ (z) (y − x)2 e−(z+1)(y−x) dydx
R2 0 x
Z ∞
1 2
= 1R+ (z)
∗ u2 (z + 1)e−(z+1)u du = 1R∗+ (z)
z+1 0 (z + 1)3
On fait le changement de variable (t, u, v) = φ(x, y, z) = (x, y −x, z(y −x)). Nous avons φ−1 (t, u, v) = (t, u+t, v/u) =
(x, y, z). On calcule
1 0 0
det φ−1 (t, u, v) = 1
1 0 = 1/u.
0 −v/u2 1/u
d’où
f(X,U,V ) (x, u, v) = 1[0,1] (x) 1R+ (u)ue−u 1[0,1] (v)e−v .
| {z } | {z } | {z }
fX (x) fU (u) fV (v)
Ceci montre que (U, V, W ) est un vecteur gaussien de moyenne nulle et de matrice de covariance Σ = AA∗ :
2 −1 −1 2 1 3 6 0 9
Σ = AA∗ = 1 1 1 −1 1 1 = 0 3 0
3 1 −4 −1 1 −4 9 0 26
Donc U ∼ N (0, 6), V ∼ N (0, 3) et W ∼ N (0, 26). Parmi ces trois couples, sauf (U, W ) n’est pas indépendant.
Cov (W,U ) 9
2. On veut trouver a ∈ R tel que Cov (W − aU, U ) = 0, c’est à dire a = V(U ) = 6 = 23 . On déduit donc
3 5 5
Z = W − U = X2 − X3 .
2 2 2
Ainsi, (U, Z) est gaussien et Cov (U, Z) = 0, ils sont indépendants. De plus, E(W |U ) = 32 U .
Ce document est mis à disposition selon les termes de la licence Creative Commons “Attribution
- Partage dans les mêmes conditions 4.0 International”.