3 Processus de branchement

3.1 Définitions

Soit (Xi,n)i (- N,n (- N une famille de v.a.i.i.d. à valeurs entières positives. On construit une suite (Zn )n (- N par :
   Z1  =  1
Zn+1   =  X1,n + ...+ XZn,n pour n > 1.
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.

  1. Notant pi = P(X1,1 = i) et m = E(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).
  2. Étudier la probabilité d’extinction de la population a = P(Z oo = 0). Pour cela réaliser une boucle Monte-Carlo de MC réalisations de ZN et estimer P(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?
  3. Étudier l’espérance E(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) = E(tX) pour t élément de [0, 1].
  1. Montrer que f prend ses valeurs dans [0, 1]. Que valent f(0), f(1), f'(0), f'(1) ?
  2. Soit fn(t) = E(tZn). Montrer que fn+1(t) = fn(f(t)) et que donc fn = fn (où le produit utilisé est la composition des applications).
  3. 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 ?
  4. Dans le cas étudié expérimentalement, calculer f et a.
  5. Calculer E(Zn) et E(Zn|Zn/=0).