📚 الرئيسية

مسابقة دكتوراه 2016Université des Sciences et de la Technologie Houari Boumediène (USTHB) — الموضوع 06

مسابقة تخصص · Probabilités & Statistiques · المدة: 2سا

Concours Doctorat LMD en Probabilités, Statistiques et Applications, épreuve de Probabilités (Sujet 3), Faculté des Mathématiques, Département de Probabilités et Statistique, USTHB, année universitaire 2016/2017, 3 novembre 2016, durée 2 heures.

التمرين 1

Exercice 1 — Vrai/Faux et espérance conditionnelle (loi exponentielle)

#markov-chain#stationary-distribution#continuous-time-markov#conditional-expectation#laplace-transform

A/ Les affirmations suivantes sont-elles vraies ou fausses ? Justifiez vos réponses et donnez l'affirmation correcte correspondante dans le cas où l'affirmation proposée est fausse.

  1. Une chaîne de Markov homogène irréductible possède une distribution stationnaire unique.
  2. Une chaîne de Markov homogène ayant une distribution stationnaire unique peut posséder des états de transition (transitoires).
  3. Le temps de séjour d'une chaîne de Markov à temps continu dans un état suit une distribution exponentielle.
  4. L'espérance d'un processus autorégressif AR(1) : Xt=φ0+φ1Xt1+εtX_t=\varphi_0+\varphi_1 X_{t-1}+\varepsilon_t, avec φ1<1|\varphi_1|<1 et εt\varepsilon_t un bruit blanc centré, est égale à φ0+φ1\varphi_0+\varphi_1.

B/ Soit (Ω,F,P)(\Omega,\mathcal F,P) un espace probabilisé, A\mathcal A une sous-σ\sigma-algèbre de F\mathcal F et XX une variable aléatoire intégrable sur (Ω,F,P)(\Omega,\mathcal F,P) de densité fXf_X.

  1. Définir la σ\sigma-algèbre engendrée par la v.a. XX.
  2. Soit YY une v.a. de loi exponentielle E(λ)\mathcal E(\lambda), indépendante de XX, et fX(s)f_X^{*}(s), s0s\ge0, la transformée de Laplace de XX. Calculer P(Y>X)P(Y>X).
الحل

A.

1. Fausse. Une chaîne irréductible ne possède pas toujours de distribution stationnaire (cas transitoire ou récurrent nul sur un espace d'états infini). Affirmation correcte : une chaîne irréductible récurrente positive possède une unique distribution stationnaire (en particulier, toute chaîne finie irréductible).

2. Vraie. Si la chaîne ne possède qu'une seule classe récurrente (positive), la distribution stationnaire est unique et attribue une masse nulle aux états transitoires ; la chaîne peut donc parfaitement posséder des états transitoires. (Exemple : 121\leftrightarrow2 récurrents et 313\to1 transitoire.)

3. Vraie. Par la propriété d'absence de mémoire, le temps de séjour dans un état ii d'une chaîne de Markov à temps continu est exponentiel, de paramètre qii=jiqij-q_{ii}=\sum_{j\ne i}q_{ij}.

4. Fausse. En prenant l'espérance : E[Xt]=φ0+φ1E[Xt1]E[X_t]=\varphi_0+\varphi_1 E[X_{t-1}], d'où par stationnarité E[Xt]=φ01φ1\boxed{E[X_t]=\frac{\varphi_0}{1-\varphi_1}} (et non φ0+φ1\varphi_0+\varphi_1).

B.

1. La σ\sigma-algèbre engendrée par XX est σ(X)={X1(B):BB(R)}={{ω:X(ω)B}:BB(R)},\sigma(X)=\big\{X^{-1}(B):B\in\mathcal B(\mathbb R)\big\}=\big\{\{\omega:X(\omega)\in B\}:B\in\mathcal B(\mathbb R)\big\}, c'est-à-dire la plus petite sous-σ\sigma-algèbre de F\mathcal F rendant XX mesurable.

2. Comme YE(λ)Y\sim\mathcal E(\lambda) est indépendante de XX, on a P(Y>x)=eλxP(Y>x)=e^{-\lambda x} pour x0x\ge0. En conditionnant par XX : P(Y>X)=E[P(Y>XX)]=E[eλX]=0eλxfX(x)dx.P(Y>X)=E\big[P(Y>X\mid X)\big]=E\big[e^{-\lambda X}\big]=\int_0^{\infty}e^{-\lambda x}f_X(x)\,dx. P(Y>X)=fX(λ)\boxed{P(Y>X)=f_X^{*}(\lambda)} c'est-à-dire la transformée de Laplace de XX évaluée au point λ\lambda.

التمرين 2

Exercice 2 — File M/M/1 avec découragement (processus de naissance et de mort)

#birth-death-process#mm1-queue#balking#stationary-distribution#positive-recurrence

On considère un système d'attente du type M/M/1M/M/1. Les clients se présentent devant la file selon un processus de Poisson d'intensité λ>0\lambda>0. Les temps de service sont exponentiels de paramètre μ>0\mu>0. On fait l'hypothèse supplémentaire suivante : un client qui arrive devant la file choisit de rester ou non en fonction du nombre de clients déjà présents dans le système. Si ce nombre est nn, le client choisit de rester avec probabilité qnq_n ou de quitter instantanément le système avec la probabilité 1qn1-q_n. On suppose que (qn)nN(q_n)_{n\in\mathbb N} est une suite décroissante de nombres positifs, et on note qq sa limite. Les temps d'interarrivée, les temps de service et les choix des clients sont des variables aléatoires indépendantes dans leur ensemble. On note XtX_t le nombre de clients présents dans le système (attente et service) à l'instant tt.

  1. Montrer que (Xt)t0(X_t)_{t\ge0} est un processus de naissance et de mort et donner ses taux de transition.
  2. À quelle condition ce processus est-il récurrent positif ?
  3. On suppose que qn=1n+1q_n=\dfrac{1}{n+1}, pour tout nNn\in\mathbb N. a. Déterminer la distribution stationnaire. b. Quel est le nombre moyen de clients dans le système en régime stationnaire ?
الحل

1.

Le nombre de clients XtX_t ne varie que par sauts de ±1\pm1 : une arrivée (taux λ\lambda) qui décide de rester (probabilité qnq_n dans l'état nn) fait passer nn+1n\to n+1, et un service achevé (taux μ\mu) fait passer nn1n\to n-1. Grâce à l'absence de mémoire des lois exponentielles et à l'indépendance globale, (Xt)(X_t) est un processus de naissance et de mort de taux λn=λqn  (n0),μn=μ  (n1).\boxed{\lambda_n=\lambda\,q_n\ \ (n\ge0),\qquad \mu_n=\mu\ \ (n\ge1).}

2.

Pour un processus de naissance et de mort, on pose ρn=k=0n1λkμk+1=(λμ)nk=0n1qk.\rho_n=\prod_{k=0}^{n-1}\frac{\lambda_k}{\mu_{k+1}}=\Big(\frac{\lambda}{\mu}\Big)^{n}\prod_{k=0}^{n-1}q_k. La chaîne est récurrente positive si et seulement si n0ρn<\sum_{n\ge0}\rho_n<\infty. Comme qnqq_n\to q, le critère de d'Alembert donne ρn+1ρn=λμqnλqμ\dfrac{\rho_{n+1}}{\rho_n}=\dfrac{\lambda}{\mu}\,q_n\to\dfrac{\lambda q}{\mu}, d'où la condition λqμ<1(i.e. λq<μ).\boxed{\frac{\lambda q}{\mu}<1\quad\text{(i.e. } \lambda q<\mu\text{).}}

3.a.

Avec qn=1n+1q_n=\dfrac{1}{n+1} : k=0n1qk=k=0n11k+1=1n!\displaystyle\prod_{k=0}^{n-1}q_k=\prod_{k=0}^{n-1}\frac{1}{k+1}=\frac{1}{n!}. En posant r=λμr=\dfrac{\lambda}{\mu}, ρn=1n!(λμ)n=rnn!,n0ρn=er<.\rho_n=\frac{1}{n!}\Big(\frac{\lambda}{\mu}\Big)^{n}=\frac{r^{n}}{n!},\qquad \sum_{n\ge0}\rho_n=e^{r}<\infty. (Ici q=limqn=0q=\lim q_n=0, la chaîne est donc toujours récurrente positive.) La distribution stationnaire πn=ρn/kρk\pi_n=\rho_n/\sum_k\rho_k vaut πn=eλ/μ(λ/μ)nn!,n0,\boxed{\pi_n=e^{-\lambda/\mu}\,\frac{(\lambda/\mu)^{n}}{n!},\quad n\ge0,} c'est-à-dire une loi de Poisson de paramètre r=λ/μr=\lambda/\mu.

3.b.

Le nombre moyen de clients en régime stationnaire est l'espérance d'une loi de Poisson(λ/μ)(\lambda/\mu) : E[X]=λμ.\boxed{E[X]=\frac{\lambda}{\mu}.}

التمرين 3

Exercice 3 — Processus de Poisson et processus de Poisson composé

#poisson-process#compound-poisson#gamma-distribution#generating-function#superposition

A/

  1. Donner une définition formelle d'un processus de Poisson {N(t), t0}\{N(t),\ t\ge0\} de taux λ\lambda.
  2. Soient {N1(t), t0}\{N_1(t),\ t\ge0\} et {N2(t), t0}\{N_2(t),\ t\ge0\} deux processus de Poisson indépendants de taux respectifs λ1\lambda_1 et λ2\lambda_2, et N(t)=N1(t)+N2(t)N(t)=N_1(t)+N_2(t), t0t\ge0, le processus mélange. Montrer que {N(t), t0}\{N(t),\ t\ge0\} est un processus de Poisson de taux λ1+λ2\lambda_1+\lambda_2.
  3. Un magasin possède un stock de Q1Q_1 unités d'un produit 1 et Q2Q_2 unités d'un produit 2. Les deux produits sont retirés de la production et ne peuvent plus fournir d'autres unités au magasin. Les clients qui demandent le produit 1 arrivent selon un processus de Poisson de taux λ1\lambda_1. Indépendamment, les clients qui demandent le produit 2 arrivent selon un processus de Poisson de taux λ2\lambda_2. Chaque client demande une unité du produit concerné. Les deux produits servent de substitut l'un à l'autre : si un client demande un produit épuisé, il est satisfait avec l'autre produit lorsque celui-ci est disponible. a. Quelle est la distribution de probabilité du temps nécessaire à l'épuisement total du stock des deux produits ? b. Quelle est la probabilité de l'épuisement du stock du produit 1 avant le produit 2 ?

B/ Soit le processus stochastique {X(t), t0}\{X(t),\ t\ge0\} défini par

X(t)=\sum_{i=1}^{N(t)}D_i \tag{1}

{N(t), t0}\{N(t),\ t\ge0\} est un processus de Poisson de taux λ\lambda, et {Di}i1\{D_i\}_{i\ge1} une suite de v.a. i.i.d. positives, indépendante de {N(t)}\{N(t)\} et de moyenne finie μ\mu.

  1. Calculer E[X(t)]E[X(t)] en fonction de λ\lambda et μ\mu, en supposant les DiD_i discrètes.
  2. On suppose les DiD_i discrètes, avec aj=P(D1=j)a_j=P(D_1=j), j=0,1,j=0,1,\dots, et de fonction génératrice A(z)A(z). Pour tout t0t\ge0, soit rj(t)=P(X(t)=j)r_j(t)=P(X(t)=j) et R(z,t)=j=0rj(t)zjR(z,t)=\sum_{j=0}^{\infty}r_j(t)z^{j}, z1|z|\le1. Montrer que R(z,t)=eλt[1A(z)]R(z,t)=e^{-\lambda t[1-A(z)]}.
الحل

A.1.

Un processus de comptage {N(t), t0}\{N(t),\ t\ge0\} est un processus de Poisson de taux λ>0\lambda>0 si :

  • N(0)=0N(0)=0 ;
  • il est à accroissements indépendants ;
  • pour tous 0s<t0\le s<t, N(t)N(s)P(λ(ts))N(t)-N(s)\sim\mathcal P\big(\lambda(t-s)\big) (accroissements stationnaires suivant une loi de Poisson).

A.2.

Pour s<ts<t, N(t)N(s)=[N1(t)N1(s)]+[N2(t)N2(s)]N(t)-N(s)=[N_1(t)-N_1(s)]+[N_2(t)-N_2(s)] est la somme de deux v.a. indépendantes de lois P(λ1(ts))\mathcal P(\lambda_1(t-s)) et P(λ2(ts))\mathcal P(\lambda_2(t-s)), donc N(t)N(s)P((λ1+λ2)(ts))N(t)-N(s)\sim\mathcal P\big((\lambda_1+\lambda_2)(t-s)\big). Comme N(0)=0N(0)=0 et que l'indépendance et la stationnarité des accroissements se conservent par addition de processus indépendants, {N(t)}\{N(t)\} est un processus de Poisson de taux λ1+λ2\lambda_1+\lambda_2. (On peut aussi le voir par les fonctions génératrices : E[zN(t)]=e(λ1+λ2)t(1z)E[z^{N(t)}]=e^{-(\lambda_1+\lambda_2)t(1-z)}.)

A.3.a.

Les deux produits étant substituables, le stock total Q1+Q2Q_1+Q_2 diminue d'une unité à chaque arrivée, quel que soit le produit demandé. Les arrivées forment le processus mélange, de taux λ1+λ2\lambda_1+\lambda_2. Le temps d'épuisement total TT est donc la somme de Q1+Q2Q_1+Q_2 interarrivées exponentielles i.i.d. de paramètre λ1+λ2\lambda_1+\lambda_2 : il suit une loi Gamma (Erlang) fT(t)=(λ1+λ2)Q1+Q2(Q1+Q21)!tQ1+Q21e(λ1+λ2)t,t0.\boxed{f_T(t)=\frac{(\lambda_1+\lambda_2)^{Q_1+Q_2}}{(Q_1+Q_2-1)!}\,t^{\,Q_1+Q_2-1}\,e^{-(\lambda_1+\lambda_2)t},\quad t\ge0.}

A.3.b.

À chaque arrivée, la demande concerne le produit 1 avec probabilité p=λ1λ1+λ2p=\dfrac{\lambda_1}{\lambda_1+\lambda_2} et le produit 2 avec Q=λ2λ1+λ2Q=\dfrac{\lambda_2}{\lambda_1+\lambda_2}. Le produit 1 s'épuise avant le produit 2 si, parmi les Q1+Q21Q_1+Q_2-1 premières arrivées, on compte au moins Q1Q_1 demandes du produit 1 : P(produit 1 0˘0e9puis0˘0e9 avant 2)=k=Q1Q1+Q21(Q1+Q21k)pkQQ1+Q21k.\boxed{P(\text{produit 1 \u00e9puis\u00e9 avant 2})=\sum_{k=Q_1}^{Q_1+Q_2-1}\binom{Q_1+Q_2-1}{k}p^{k}\,Q^{\,Q_1+Q_2-1-k}.}

B.1.

Par la formule de Wald (conditionnement sur N(t)P(λt)N(t)\sim\mathcal P(\lambda t), les DiD_i i.i.d. de moyenne μ\mu et indépendantes de NN) : E[X(t)]=E[N(t)]E[D1]=λμt.E[X(t)]=E[N(t)]\,E[D_1]=\boxed{\lambda\mu t.}

B.2.

En conditionnant sur N(t)=nN(t)=n et par indépendance des DiD_i, E[zX(t)N(t)=n]=A(z)nE\big[z^{X(t)}\mid N(t)=n\big]=A(z)^{n}. D'où R(z,t)=n0A(z)neλt(λt)nn!=eλteλtA(z)=eλt[1A(z)].R(z,t)=\sum_{n\ge0}A(z)^{n}\,e^{-\lambda t}\frac{(\lambda t)^{n}}{n!}=e^{-\lambda t}\,e^{\lambda t A(z)}=\boxed{e^{-\lambda t[1-A(z)]}.}