BCE Maths appliquées HEC ECE 2015Sujet, corrigé et rapport du jury

Mercredi 29 avril 2015, de 8 h. à 12 h.

La présentation, la lisibilité, l'orthographe, la qualité de la rédaction, la clarté et la précision des raisonnements entreront pour une part importante dans l'appréciation des copies.
Les candidats sont invités à encadrer dans la mesure du possible les résultats de leurs calculs.
Ils ne doivent faire usage d'aucun document : I'utilisation de toute calculatrice et de tout matériel électronique est interdite. Seule l'utilisation d'une règle graduée est autorisée.
Si au cours de l'épreuve, un candidat repère ce qui lui semble être une erreur d'énoncé, il la signalera sur sa copie et poursuivra sa composition en expliquant les raisons des initiatives qu'il sera amené à prendre

Soit n un entier supérieur ou égal à 2 et B = (e_1, e_2, …, e_n) la base canonique de ℝ^n.
Soit v un vecteur donné de ℝ^n de coordonnées v_1, v_2, …, v_n dans la base B et qui vérifie ∑_(i = 1)^n v_i = 1.
Soit f l'application définie sur ℝ^n qui à tout vecteur x = (x_1, x_2, …, x_n) ∈ ℝ^n, associe le vecteur f(x) défini par: f(x) = x − (∑_(i = 1)^n x_i)v.
1.a) Montrer que f est un endomorphisme de ℝ^n.
b) Montrer que f ∘ f = f.
2. Déterminer le spectre de f.
3.a) Montrer que le vecteur y appartient à l'image de f, notée Imf, si et seulement si f(y) = y.
b) Montrer que la dimension de Imf est inférieure ou égale à n − 1.
c) Montrer que pour tout i ∈ [ [1, n − 1] ], on a (e_i − e_(i + 1)) ∈ Imf.
d) En déduire une base et la dimension de Imf. Quel est le rang de f ?
4.a) Déterminer une base du noyau de f.
b) Quels sont les sous-espaces propres de f ?
c) L'endomorphisme f est-il diagonalisable ?
5. Écrire la matrice M de l'endomorphisme f dans la base canonique de ℝ^n et la matrice M^′ de f dans une base de vecteurs propres.

L'objectif de cette partie est de comparer deux estimateurs sans biais et convergents du paramètre inconnu 1/λ. Pour n entier de ℕ^∗, soit (X_1, X_2, …, X_n) un n-échantillon de variables aléatoires à valeurs strictement positives, indépendantes et de même loi que X.
On pose pour tout n ∈ ℕ^∗ : Y_n = ∑_(i = 1)^n X_i, X¯_n = 1/nY_n et M_n = max(X_1, X_2, …, X_n).

b) Établir pour tout n ∈ ℕ^∗, l'existence de l'espérance E(M_n) de la variable aléatoire M_n.
c) En posant z = 1 − e^(− λx), justifier pour tout a > 0, l'égalité :

d) En déduire que l'on a : E(M_n) = − n/λ∫_0^1 z^(n − 1)ln(1 − z)dz.
e) Montrer que la fonction z ⟼ (1 − z)(1 − ln(1 − z)) définie sur l'intervalle [ 0,1 [, est une primitive de la fonction z ⟼ ln(1 − z).
À l'aide d'une intégration par parties, en déduire pour tout n ∈ ℕ^∗ une relation entre E(M_(n + 1)) et E(M_n).
f) On pose pour tout n ∈ ℕ^∗ : u_n = ∑_(j = 1)^n 1/j. Déduire de la question précédente que E(M_n) = 1/λu_n.
4. On pose pour tout n ∈ ℕ^∗ : M_n^′ = (M_n)/(u_n) et v_n = ∑_(j = 1)^n 1/(j^2). On admet que V(M_n) = 1/(λ^2)v_n.
a) Calculer E(M_n^′) et V(M_n^′).
b) Justifier la convergence de la série de terme général v_n. Déterminer lim_(n → + ∞)u_n.
c) Déduire des questions 4.a) et 4.b) que M_n^′ est un estimateur sans biais et convergent du paramètre 1/λ.
d) Soit Q la fonction polynomiale définie sur ℝ par Q(x) = ∑_(j = 1)^n(x − 1/j)^2.

À l'aide de l'étude de Q, établir l'inégalité : u_n^2 ⩽ nv_n.
e) Comparer alors V(M_n^′) et V(X¯_n). Conclure.

Les notations et le contexte sont ceux de la partie I.
Dans cette partie, on suppose que la durée de vie (en heures) d'un composant électronique est une variable aléatoire X à valeurs strictement positives suivant la loi exponentielle de paramètre λ (avec λ > 0 ).
On suppose qu'en cas de panne, le composant électronique est immédiatement remplacé par un composant neuf dont la durée de vie est indépendante et de même loi que celle des composants précédents.
Pour j ∈ ℕ^∗, on note X_j la variable aléatoire égale à la durée de vie du j-ème composant.
Pour n entier de ℕ^∗, on constitue ainsi un n-échantillon ( X_1, X_2, …, X_n ) de variables aléatoires à valeurs strictement positives, indépendantes et de même loi que X.
On note (x_1, x_2, …, x_n) la réalisation du n-échantillon (X_1, X_2, …, X_n) et on pose : x¯ = 1/n∑_(i = 1)^n x_i.
5.a) Donner une interprétation de 1/λ. Dans quelle unité s'exprime 1/λ ?
b) Soit p ∈ ]0, 1[. Déterminer en fonction de p et λ, l'unique réel h_p pour lequel on a P([X > h_p]) = p.
c) Proposer un estimateur H_n sans biais et convergent pour le paramètre h_p.
d) On suppose que sur un échantillon de 100 composants, on a obtenu : ∑_(i = 1)^(100)x_i = 10^5 heures.

Donner une estimation de h_(1/2) (on donne ln2 ≃ 0, 7 ).
6. On admet sans démonstration que pour tout n ∈ ℕ^∗, la fonction de répartition F_(Y_n) de la variable aléatoire Y_n est donnée par :

Soit t > 0 un réel fixé et N(t) la variable aléatoire égale au nombre de pannes dans l'intervalle [ 0, t.
a) Déterminer la loi de la variable aléatoire N(t).
b) Donner une estimation " naturelle" du nombre moyen de pannes dans l'intervalle [ 0, t[.
7. L'objectif de cette question est de déterminer un intervalle de confiance asymptotique du paramètre inconnu 1/λ au niveau de confiance 1 − α(0 < α < 1).
a) Soit Φ la fonction de répartition de la loi normale centrée réduite et soit t_1 et t_2 deux réels vérifiant Φ(t_1) = 1 − α/2 et Φ(t_2) = α/2. Montrer que t_2 = − t_1.
b) On pose : R_n = √n(λX¯_n − 1). Justifier que la suite de variables aléatoires (R_n)_(n ∈ ℕ^∗) converge en loi vers une variable aléatoire suivant la loi normale centrée réduite.
c) En déduire que lim_(n → + ∞)P([ − t_1 ⩽ R_n ⩽ t_1]) = 1 − α.
d) Déterminer un intervalle de confiance asymptotique de 1/λ au niveau de confiance 1 − α.
e) Que se passe-t-il lorsque α est proche de 0 ou lorsque α est proche de 1 ?
f) Pour α = 0, 05, on donne t_1 ≃ 2. Montrer qu'avec l'échantillon de la question 5.d), la réalisation de l'intervalle de confiance asymptotique de 1/λ au niveau de confiance 0,95 est : [833,1250].

Les notations et le contexte sont ceux des parties précédentes.
Pour tout n ∈ ℕ^∗, on pose T_n = M_n − 1/λlnn et on note F_(T_n) la fonction de répartition de T_n.
8.a) Montrer que pour tout x réel, on a : F_(T_n)(x) = (F_X(x + 1/λlnn))^n.
b) Pour chaque réel x fixé, déterminer un entier naturel N_x (qui dépend de x ) tel que pour tout entier n ⩾ N_x, on a : x + 1/λlnn > 0.
En déduire que pour tout entier n ⩾ N_x, on a : F_X(x + 1/λlnn) = 1 − (e^(− λx))/n.
c) Montrer que pour tout x réel, on a : lim_(n → + ∞)F_(T_n)(x) = exp(− e^(− λx)) = e^(− e^(− λx)).
9. Soit F la fonction définie sur ℝ à valeurs réelles telle que : F(x) = exp(− e^(− λx)).
a) Justifier que F est de classe C^∞ sur ℝ et montrer que F réalise une bijection de ℝ sur l'intervalle ]0, 1[.
b) En déduire que F est la fonction de répartition d'une variable aléatoire T admettant une densité f_T continue sur ℝ que l'on déterminera; on dit que T suit la loi de Gumbel de paramètre λ.
c) Justifier la convergence en loi de la suite de variables aléatoires (T_n)_(n ∈ ℕ^∗) vers la variable aléatoire T.
10.a) Établir l'existence de l'espérance E(T) de la variable aléatoire T.
b) On pose : Z = − 1/λln(λX). Montrer que les variables aléatoires Z et T sont de même loi.
c) Justifier l'égalité : E(T) = − 1/λ∫_0^(+ ∞)e^(− t)lnt dt.
d) À l'aide de la concavité de la fonction ln sur ℝ_+^∗, établir l'inégalité : E(T) ⩾ 0.
11. On suppose dans cette question que λ = 1.
a) Expliciter la bijection réciproque G de la fonction F.
b) On considère le programme Scilab suivant :
x = linspace(− 2, 2, 400); y = (exp(− exp(− x))); plot(x, y), plot(y, x)
(i) Le réel 0 fait-il partie des nombres renvoyés par la commande x = linspace (− 2, 2, 400) ?
(ii) Quel sera le résultat de l'exécution de ce programme?
c) Soit U une variable aléatoire suivant la loi uniforme sur l'intervalle ]0, 1[.

Quelle est la loi de la variable aléatoire G(U) ?
d) Établir l'inégalité : E(T) ⩽ 1.
e) Par une méthode de votre choix, écrire en Scilab les commandes qui permettent de simuler la loi de T.
f) Écrire en Scilab les commandes qui permettent de renvoyer une valeur numérique approchée de E(T) en utilisant la méthode de Monte-Carlo.