BCE Maths appliquées ESSEC ECE 2018, épreuve 2Sujet, corrigé et rapport du jury

Lundi 7 mai 2018, de 14 h. à 18 h.

On s'intéresse à l'évolution d'une population de petits organismes (typiquement des insectes) pendant une "saison" reproductrice de durée maximale T où T ∈ 𝕀ℕ^∗. Les insectes sont supposés vivre une unité de temps, au bout de laquelle ils meurent en pondant un certain nombre d'œufs. Au moment du dépôt d'un œuf, un processus chimique, la diapause est susceptible de se mettre en marche qui entraîne l'arrêt de maturation de l'œuf jusqu'à la saison suivante. Ainsi, à chaque date t de la saison, une génération d'insectes s'éteint, en déposant des œufs. Immédiatement, une proportion p(t) de ces œufs se mettent en diapause. Les œufs qui ne sont pas entrés en diapause éclosent avant la date t + 1, donnant naissance à une nouvelle génération d'insectes, qui s'éteindra à la date t + 1 en déposant des œufs, etc... Comme, à la fin de la saison, tous les organismes vivants de la population meurent, hormis les œufs qui sont en diapause, ce sont ces derniers qui seront à l'origine d'une nouvelle population qui éclora à la saison suivante. Il est donc fondamental pour la survie de la lignée que les organismes adoptent une stratégie maximisant le nombre d'œufs en diapause accumulés jusqu'à la date où la saison s'achève.
Au cours du problème, on s'intéressera plus particulièrement au cas où la durée de la saison est une variable aléatoire τ pouvant prendre des valeurs entières entre 1 et T. Pour t ∈ {1, 2, …, T}, l'événement [τ = t] signifiera donc que la saison s' arrête à la date t.
Toutes les variables aléatoires intervenant dans le problème sont définies sur un espace probabilisé ( Ω, A, P ). Pour toute variable aléatoire Y, on notera E(Y) son espérance lorsqu'elle existe.

Dans cette question, on définit l'évolution formelle du nombre d'œufs en diapause entre les dates 0 et T.
On note D(t) = nombre d'œufs en diapause à la date t. Les œufs pondus à la date t qui entrent en diapause sont comptabilisés à la date t + 1.
N(t) = nombre moyen d'œufs produits à la date t
p(t) = proportion des œufs produits à la date t qui entrent en diapause
Par convention, la date 0 d'une saison est celle où les insectes nés des œufs en diapause de la saison précédente pondent N(0) œufs. On suppose pour simplifier :

Enfin, on suppose qu'à chaque date t de la saison, un individu produit en moyenne α œufs ( α étant un réel strictement positif). Par simplicité, on supposera que α reste constant pendant toute la saison.
1)
(a) Montrer que D(t + 1) = D(t) + p(t)N(t) pour tout entier t tel que 0 ≤ t ≤ T − 1.
(b) Montrer que N(t + 1) = α(1 − p(t))N(t) pour tout entier t tel que 0 ≤ t ≤ T − 1.
2) On suppose dans cette question que α ≤ 1.
(a) Montrer que pour tout entier t tel que 0 ≤ t ≤ T − 1, N(t + 1) ≤ N(t).
(b) Montrer que pour tout entier t tel que 0 ≤ t ≤ T − 1,

(c) Montrer que pour tout entier t tel que 0 ≤ t ≤ T,

(d) Montrer que pour tout entier t tel que 0 ≤ t ≤ T, D(t) ≤ N(0).
(e) On suppose que p(0) = 1.
i) Montrer que pour tout entier t tel que 1 ≤ t ≤ T, N(t) = 0.
ii) Montrer que pour tout entier t tel que 1 ≤ t ≤ T, D(t) = N(0).
iii) En déduire que si α ≤ 1, la meilleure stratégie adaptée à la saison est que les N(0) œufs produits à la date 0 entrent en diapause immédiatement.
3) On suppose désormais α > 1 jusqu'à la fin du problème.

On introduit maintenant τ une variable aléatoire à valeurs dans {1, 2, …, T} qui représente la date où s'achève la saison. On suppose que pour tout t ∈ {1, 2, …, T}, P(τ = t) > 0.
(a) Montrer que pour tout t ∈ {1, 2, …, T}, P(τ ≥ t) > 0. On définit alors H(t) = P_([τ ≥ t])(τ = t).
(b) Montrer que H(t) = (P(τ = t))/(P(τ ≥ t)).
(c) Montrer que H(T) = 1.
(d) Calculer H(t) pour t ∈ {1, 2, …, T} si τ suit une loi uniforme sur {1, 2, …, T}.
(e)
i) Soient T réels λ_1, λ_2, …, λ_T tels que 0 < λ_1 < λ_2 < ⋯ < λ_T = 1. Par convention, on pose λ_0 = 0. Soient q_1 = λ_1, q_2 = λ_2 − λ_1, …, q_T = λ_T − λ_(T − 1).
Montrer que (q_i)_(1 ≤ i ≤ T) définit une loi de probabilité sur {1, 2, …, T}.
ii) Calculer H(t) si τ suit la loi précédente.
iii) On suppose que T ≥ 2 et de plus que pour tout entier n tel que 1 ≤ n ≤ T − 1, on a λ_(n + 1) − λ_n ≥ λ_n − λ_(n − 1). Montrer que t ↦ H(t) est croissante sur {1, 2, …, T}.
On suppose désormais que t ↦ H(t) est croissante. Le but est maintenant de trouver une stratégie adéquate pour maximiser la quantité E(lnD(τ)). On va commencer par regarder un exemple simple.
4) On suppose ici que T = 2, que H est donnée par H(1) = 1/2 et H(2) = 1 et que α = 4.
(a)
i) Déterminer P(τ = 1).
ii) Quelle est la loi de τ ?
(b) Montrer que pour D(1) et N(1) donnés, D(2) est maximum pour p(1) = 1.
(c) On suppose p(1) = 1. Montrer que

(d) Construire le tableau de variations sur ]0,1] de la fonction φ définie par :

(e) Déterminer p^∗(0) qui maximise E(lnD(τ)).

Par convention, on conviendra que si h est une fonction numérique définie sur {0, 1, 2, …, T}, on a

On pose

(a) Montrer que ψ_ξ est croissante sur [0, 1].
(b) Calculer ψ_ξ(1).
(c)
i) Calculer ψ_ξ(0). On pose désormais

ii) Montrer que que pour tout t ∈ {0, 1, 2, …, T − 1}, A(X(t)) ≤ X(t + 1) ≤ 1.
iii) Montrer que ξ ↦ A(ξ) est croissante sur [0, 1].
8) Justifier l'égalité de variables aléatoires :

On pose R^(0) = ln(D(1))/(N(0)) et R^(t) = ln(D(t + 1))/(D(t)) pour 1 ≤ t ≤ T − 1.
9)
(a) Montrer que

(b) Montrer que (D(1))/(N(0)) = (αX(1))/(1 + (α − 1)X(1)).
(c) Montrer que (D(t + 1))/(D(t)) = (αX(t + 1))/(X(t)[1 + (α − 1)X(t + 1)]) pour 1 ≤ t ≤ T − 1.

Pour x et y deux réels strictement positifs, on pose u(x, y) = lnα − lnx + lny − ln(1 + (α − 1)y).
(d) Montrer que R^(0) = u(1, X(1)).
(e) Montrer que R^(t) = u(X(t), X(t + 1)) pour 1 ≤ t ≤ T − 1.
(f) Conclure que

On voit donc que maximiser E(lnD(τ)) revient à choisir, à chaque date t telle que 1 ≤ t ≤ τ − 1, la valeur X(t + 1) vérifiant la contrainte A(X(t)) ≤ X(t + 1) ≤ 1 de façon à rendre maximale l'expression

On expose dans cette partie les deux premières étapes de la méthode de la programmation dynamique pour résoudre le problème.
10) Soit B un événement. On note 𝕀_B la variable aléatoire telle que

(a) Déterminer la loi de 𝕀_B.
(b) Soient B et C deux événements. Montrer l'égalité de variables aléatoires 𝕀_(B ∩ C) = 𝕀_B ⋅ 𝕀_C.
(c) On suppose que 0 < P(B) < 1. Si Y est une variable aléatoire prenant un nombre fini de valeurs, on définit la variable aléatoire notée E_B(Y) par

où B¯ désigne l'événement contraire de B.
i) Soient Y et Z deux variables aléatoires prenant un nombre fini de valeurs. Montrer que

ii) Montrer que E(E_B(Y)) = E(Y).
iii) Montrer que E_B(Y𝕀_B) = E_B(Y)𝕀_B.
11) On suppose dans cette question que, quand l'événement [τ = T] est réalisé, X(1), …, X(T − 1) sont connus. Comme on l'a vu précédemment, si on pose x = X(T − 1), le meilleur choix à faire est alors de prendre pour X(T) la valeur y^∗(x, T − 1) ∈ [A(x), 1] qui maximise u(x, y).
(a) Montrer que y^∗(x, T − 1) = 1.
(b) Montrer que u(x, 1) = − lnx.
12) On suppose maintenant que, quand l'événement (τ ≥ T − 1) est réalisé, X(1), …, X(T − 2) sont connus. La stratégie reste donc de choisir X(T − 1) et X(T) de façon à maximiser E(∑_(t = T − 2)^(τ − 1)u(X(t), X(t + 1))).
La variable aléatoire τ prend les deux valeurs T − 1 et T avec les probabilités respectives P_([τ ≥ T − 1])(τ = T − 1) et P_([τ ≥ T − 1])(τ = T).
(a) Montrer que

(b) Montrer que

(c) Montrer que

(d) On suppose que X(T − 1) est donné.
i) Montrer que le meilleur choix pour X(T) est 1 .
ii) Montrer que pour un tel choix u(X(T − 1), X(T)) = − lnX(T − 1).
(e) Montrer que P_([τ ≥ T − 1])(τ = T) = 1 − H(T − 1).

On veut maintenant choisir la stratégie optimale à la date T − 2.
(f) Montrer qu'on doit choisir pour X(T − 1) la valeur y^∗(X(T − 2), T − 2) ∈ [A(X(T − 2)), 1] de telle sorte que

soit maximal.
(g) Calculer φ^′(y),
(h) Construire le tableau de variation de φ dans le cas (H(T − 1))/((α − 1)(1 − H(T − 1))) ≤ 1.
(i) Construire le tableau de variation de φ dans le cas (H(T − 1))/((α − 1)(1 − H(T − 1))) ≥ 1.
(j) Donner la valeur de y^∗(X(T − 2), T − 2).