3 Processus de branchement
3.1 Définitions
Soit (Xi,n)i,n une famille de v.a.i.i.d. à valeurs entières positives. On construit une suite
(Zn )n par : Le processus est appelé processus de branchement. Par exemple, Zn peut représenter le nombre d’individus dans
une population à la n-ème génération et Xi,n le nombre d’enfants du ième individu.
3.2 Partie expérimentale
On se limitera ici à étudier des v.a. prenant les valeurs 0, 1 et 2.
- Notant pi = (X1,1 = i) et m = (X1,1), exprimez p1 et p2 en fonction de m et p0. Assurez-vous par
un test dans votre programme que le choix de m et p0 est permise. Écrire un programme donnant une
réalisation de la suite (Zn)1....,N avec m, p0 et N en entrées. Réfléchir aux débordements possibles
et imposer une limite supérieure MAXZ à Zn. (On rappelle que 231 - 1 = 2147483647).
- Étudier la probabilité d’extinction de la population a = (Z = 0). Pour cela réaliser une boucle
Monte-Carlo de MC réalisations de ZN et estimer (ZN = 0). Introduire un compteur d’échecs
NBECH pour compter les cas où il y a débordement. Considérer différentes valeurs de p0 et m. (
p0 = 0, p0 > 0, m = 1, m > 1). Prendre N = 100, 1000, 10000. Qu’observez-vous?
- Étudier l’espérance (ZN |ZN 0) et comparer la avec mN . Quelle conjecture formulez vous?
L’animation en ligne donne un exemple de simulation des premières générations du processus.
3.3 Partie théorique
Soit X = X1,1 et f la fonction génératrice de X : f(t) = (tX) pour t élément de
[0, 1].
- Montrer que f prend ses valeurs dans [0, 1]. Que valent f(0), f(1), f'(0), f'(1) ?
- Soit fn(t) = (tZn). Montrer que fn+1(t) = fn(f(t)) et que donc fn = fn (où le produit utilisé
est la composition des applications).
- Comparer a et fn(0). Montrer que a vérifie f(a) = a. Étudier les variations de f. Que se passe-t-il
si p0 = 0 , si m < 1, si m > 1 ?
- Dans le cas étudié expérimentalement, calculer f et a.
- Calculer (Zn) et (Zn|Zn0).