Top Banner
1/5 2/5 3/5 4/5 5/5 7. Loi normale et th´ eor` eme central limite MTH2302D S. Le Digabel, ´ Ecole Polytechnique de Montr´ eal A2017 (v2) MTH2302D: loi normale 1/35
35

7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

Jun 26, 2020

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

7. Loi normale et theoreme central limite

MTH2302D

S. Le Digabel, Ecole Polytechnique de Montreal

A2017(v2)

MTH2302D: loi normale 1/35

Page 2: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Plan

1. Loi normale

2. Loi normale centree reduite

3. Approximation d’une binomiale

4. Loi lognormale

5. Theoremes limites

MTH2302D: loi normale 2/35

Page 3: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

1. Loi normale

2. Loi normale centree reduite

3. Approximation d’une binomiale

4. Loi lognormale

5. Theoremes limites

MTH2302D: loi normale 3/35

Page 4: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi normale

On dit qu’une variable aleatoire continue X suit une loi normale deparametres µ et σ2 si sa fonction de densite est

fX(x) =1

σ√

2πexp

(−(x− µ)2

2σ2

)pour tout x .

On denote ceci X ∼ N(µ, σ2).

MTH2302D: loi normale 4/35

Page 5: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi normale : proprietesProprietes de fX :

1. limx→±∞

fX(x) = 0.

2. fX(µ+ x) = fX(µ− x) (symetrie par rapport a l’axe x = µ).

3. fX atteint son maximum en x = µ (µ est le mode de X).

4. Les points d’inflexion du graphe de fX sont x = µ± σ.

MTH2302D: loi normale 5/35

Page 6: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi normale : proprietes (suite)

Si X ∼ N(µ, σ2) alors

1. P (X < µ− x) = P (X > µ+ x).

2. FX(µ− x) = 1− FX(µ+ x).

MTH2302D: loi normale 6/35

Page 7: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Moyenne et variance de la loi normale

Si X ∼ N(µ, σ2) alors

1. E(X) = µ.

2. V(X) = σ2.

MTH2302D: loi normale 7/35

Page 8: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

-4 -2 0 2 4

0.0

0.1

0.2

0.3

0.4

fonction de densité de X~N(0,1)

x

f(x)

MTH2302D: loi normale 8/35

Page 9: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

-4 -2 0 2 4

0.0

0.2

0.4

0.6

0.8

1.0

fonction de répartition de X~N(0,1)

x

F(x)

MTH2302D: loi normale 9/35

Page 10: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi normale : calcul avec des logiciels

I Excel :fX(x) = LOI.NORMALE(x, µ, σ, 0).FX(x) = LOI.NORMALE(x, µ, σ, 1).

I R :fX(x) = dnorm(x, mean=µ, sd=σ).FX(x) = pnorm(x, mean=µ, sd=σ).

MTH2302D: loi normale 10/35

Page 11: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

1. Loi normale

2. Loi normale centree reduite

3. Approximation d’une binomiale

4. Loi lognormale

5. Theoremes limites

MTH2302D: loi normale 11/35

Page 12: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi normale centree reduiteLorsque µ = 0 et σ2 = 1, la loi normale N(0, 1) est appelee centreereduite et on la denote par Z.

Sa fonction de densite est

φ(z) =1√2πe−

z2

2 .

Sa fonction de repartition est

Φ(z) =1√2π

∫ z

−∞e−

t2

2 dt .

Puisque cette integrale est difficile a evaluer, on a recours a unetable de loi normale pour calculer Φ(z). Voir livre page 476 (2emeedition) / page 512 (3eme edition) ou sur le site web du cours.

MTH2302D: loi normale 12/35

Page 13: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi normale centree reduite (suite)

Si X ∼ N(µ, σ2) alors

Z =X − µσ

∼ N(0, 1).

On peut donc ramener toute loi normale a une loi centree reduite.

MTH2302D: loi normale 13/35

Page 14: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Methodes de calculSi Z ∼ N(0, 1)

I Si b ≥ 0 alors P (Z ≤ b) = Φ(b).

I Si b < 0 alorsΦ(b) = P (Z ≤ b) = 1− P (Z ≤ −b) = 1− Φ(−b).

I P (Z ≥ b) = 1− P (Z ≤ b) = 1− Φ(b).

I P (a ≤ Z ≤ b) = P (Z ≤ b)− P (Z ≤ a) = Φ(b)− Φ(a).

Si X ∼ N(µ, σ2), alors P (X ≤ b) = Φ(b−µσ

)et

P (a ≤ X ≤ b) = P

(a− µσ≤ Z ≤ b− µ

σ

)= Φ

(b− µσ

)−Φ

(a− µσ

).

MTH2302D: loi normale 14/35

Page 15: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Exemples

I Exemple 1 : Verifier que µ = 0 et σ2 = 1 si X ∼ N(0, 1).

I Exemple 2 : Determiner Q1, Q2 et Q3 si X ∼ N(0, 1).

MTH2302D: loi normale 15/35

Page 16: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Exemple 3

Si Z ∼ N(0, 1), calculer

1. P (Z ≤ 1.25).

2. P (Z ≤ −0.52).

3. P (Z > −1).

4. Si X ∼ N(µ = 100, σ2 = 4), calculer P (98 < X ≤ 104).

5. Si P (Z ≤ b) = 0.6628, determiner b.

6. Si P (Z ≤ b) = 0.3446, determiner b.

MTH2302D: loi normale 16/35

Page 17: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Additivite

Soit X1, X2, . . . , Xn des variables aleatoires independantes avec

Xi ∼ N(µi, σ2i ) pour tout i.

Soit Y = a0 + a1X1 + a2X2 + · · ·+ anXn.

Alors Y ∼ N(µ, σ2), ou

µ = a0 + a1µ1 + a2µ2 + · · ·+ anµn ,

σ2 = a21σ

21 + a2

2σ22 + · · ·+ a2

nσ2n .

MTH2302D: loi normale 17/35

Page 18: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Exemple 4

Un assemblage consiste a inserer un arbre dans un palier selon leschema ci-dessous.

1 XX 2

Si X1 ∼ N(1.5, 0.0016) et X2 ∼ N(1.48, 0.0009) sont les deuxdiametres, le jeu entre les deux elements est Y = X1 −X2. Lesv.a. X1 et X2 sont independantes.

L’assemblage echoue si X1 < X2.

Dans quel pourcentage de cas l’assemblage echoue-t-il ?

MTH2302D: loi normale 18/35

Page 19: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

1. Loi normale

2. Loi normale centree reduite

3. Approximation d’une binomiale

4. Loi lognormale

5. Theoremes limites

MTH2302D: loi normale 19/35

Page 20: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Approx. d’une loi binomiale par une loi normale

Soit X ∼ B(n, p) une variable aleatoire suivant une loi binomiale.Alors X est la somme de variables de Bernoulli independantes deparametre p.

Si n est grand alors X suit approximativement une loi normaleN(µ = np, σ2 = np(1− p)).

Cette approximation est bonne si

I np > 5 lorsque p ≤ 12 .

I n(1− p) > 5 lorsque p > 12 .

MTH2302D: loi normale 20/35

Page 21: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Approx. d’une binomiale par une normale (suite)Puisque X ∼ B(n, p) est une variable discrete, on cherche acalculer des probabilites comme P (X = x).

Or ceci n’a pas de sens pour une v.a. continue et on doit corrigerla valeur cherchee pour pouvoir utiliser l’approximation de X parune loi normale.

Par exemple

Valeur cherchee Valeur corrigee

P (X = x) P (x− 12 ≤ X ≤ x+ 1

2)

P (a ≤ X ≤ b) P (a− 12 ≤ X ≤ b+ 1

2)

Cette correction est appelee correction pour la continuite.

MTH2302D: loi normale 21/35

Page 22: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Exemple 5

On lance une piece 200 fois. Quelle est la probabilite d’obtenir aumoins 110 piles ?

MTH2302D: loi normale 22/35

Page 23: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

1. Loi normale

2. Loi normale centree reduite

3. Approximation d’une binomiale

4. Loi lognormale

5. Theoremes limites

MTH2302D: loi normale 23/35

Page 24: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi lognormale

Une variable aleatoire X suit une loi lognormale de parametres µYet σ2

Y si

Y = ln(X) ∼ N(µY , σ2Y ) .

C’est equivalent a definir X = exp(Y ).

MTH2302D: loi normale 24/35

Page 25: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi lognormale : fonction de densite

La fonction de densite d’une variable aleatoire lognormale X deparametres µY , σ2

Y est

fX(x) =

1

xσY√

2πexp

(−(ln(x)− µY )2

2σ2Y

)si x > 0 ,

0 sinon.

MTH2302D: loi normale 25/35

Page 26: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi lognormale : fonction de repartition

La fonction de repartition d’une variable aleatoire lognormale X deparametres µY , σ2

Y est

FX(x) =

Φ(

ln(x)− µYσY

)si x > 0 ,

0 sinon.

MTH2302D: loi normale 26/35

Page 27: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi lognormale : moyenne et variance

Soit X une variable aleatoire lognormale de parametres µY , σ2Y .

Alors :

1. E(X) = exp(µY + 12σ

2Y ).

2. V(X) = exp(2µY + σ2Y )(exp(σ2

Y )− 1).

MTH2302D: loi normale 27/35

Page 28: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi lognormale : proprietes

Soit X1, X2, . . . , Xn des variables aleatoires lognormalesindependantes de parametres µYi , et σ2

Yipour i = 1, 2, . . . , n.

AlorsW = bXa1

1 Xa22 · · ·X

ann

suit une loi lognormale de parametres

µY = ln(b) +n∑i=1

aiµYi

et

σ2Y =

n∑i=1

a2iσ

2Yi.

MTH2302D: loi normale 28/35

Page 29: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Exemple 6

Soit

Y1 = ln(X1) ∼ N(4, 1)

Y2 = ln(X2) ∼ N(3, 0.5)

Y3 = ln(X3) ∼ N(2, 0.4)

Y4 = ln(X4) ∼ N(1, 0.01)

etW = e1.5

(X2.5

1 X0.22 X0.7

3 X3.14

).

Calculer P (2× 104 ≤W ≤ 6× 105).

MTH2302D: loi normale 29/35

Page 30: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

1. Loi normale

2. Loi normale centree reduite

3. Approximation d’une binomiale

4. Loi lognormale

5. Theoremes limites

MTH2302D: loi normale 30/35

Page 31: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Loi des grands nombres

Soient X1, X2, . . . , Xn des variables aleatoires independantes ayantla meme distribution, avec E(Xi) = µ pour i = 1, 2, . . . , n.

Alors pour tout ε > 0,

limn→∞

P

∣∣∣∣∣∣∣∣n∑i=1

Xi

n− µ

∣∣∣∣∣∣∣∣ > ε

= 0 .

MTH2302D: loi normale 31/35

Page 32: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Theoreme central limite (TCL)

Soit X1, X2, . . . , Xn une suite de variables aleatoiresindependantes, avec E(Xi) = µi et V(Xi) = σ2

i pouri = 1, 2, . . . , n.

Alors la variable aleatoire

Z =∑n

i=1 (Xi − µi)√∑ni=1 σ

2i

suit approximativement une loi normale N(0, 1) si n est grand.

MTH2302D: loi normale 32/35

Page 33: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

TCL : cas particulierSoit X1, X2, . . . , Xn une suite de variables aleatoiresindependantes et identiquement distribuees, avec E(Xi) = µ etV(Xi) = σ2 pour i = 1, 2, . . . , n.

Alors la variable aleatoire

Z =∑n

i=1Xi − nµσ√n

suit approximativement une loi normale N(0, 1) si n est grand.

Autres formulations (pour n →∞) :n∑i=1

Xi ∼ N(nµ, nσ2), ou

X =1n

n∑i=1

Xi ∼ N(µ, σ2/n) .

MTH2302D: loi normale 33/35

Page 34: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Exemple 7

On lance un de 100 fois. Quelle est la probabilite que la somme desresultats soit entre 340 et 360 ?

MTH2302D: loi normale 34/35

Page 35: 7. Loi normale et théorème central limite - GERAD · 1/52/53/54/55/5 7. Loi normale et th eor eme central limite MTH2302D S. Le Digabel, Ecole Polytechnique de Montr eal A2017 (v2)

1/5 2/5 3/5 4/5 5/5

Theoreme central limite : quelle valeur de n estassez grande ?

I Si les lois des Xi sont proches d’une loi normale alors pourn ≥ 4 l’approximation donnee par le theoreme central limiteest bonne.

I Si les lois des Xi sont moyennement proches d’une loi normale(p. ex. loi uniforme) alors pour n ≥ 12 l’approximation donneepar le theoreme central limite est bonne.

I Si les lois des Xi ne sont pas proches d’une loi normale alorspour n ≥ 100 l’approximation donnee par le theoreme centrallimite sera bonne (par exemple fonction de densite tresasymetrique).

MTH2302D: loi normale 35/35