📚 الرئيسية

مسابقة دكتوراه 2012Université Abderrahmane Mira - Béjaïa — الموضوع 03

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

Concours d'entrée en Doctorat LMD, Option Analyse et Probabilité, Épreuve Probabilités et Statistique, Université de Béjaïa, Faculté des Sciences Exactes, 18/11/2012, durée 2h.

التمرين 1

Erreurs dans un livre et approximation de Poisson

#binomial-distribution#poisson-approximation#expectation#variance

Un livre de 1500 pages contient 1000 erreurs réparties au hasard. On ouvre le livre à une page quelconque et on désigne par XX le nombre d'erreurs rencontrées dans cette page. 1) Donner la distribution de XX ainsi que ses paramètres. 2) Calculer E(X)E(X) et var(X)\mathrm{var}(X). 3) On approxime la loi de XX par une loi de Poisson : est-ce légitime ? Déterminer son paramètre. 4) Calculer la probabilité de n'avoir aucune erreur. 5) Calculer la probabilité de rencontrer plus de deux erreurs.

الحل

1) Chaque erreur tombe sur la page choisie avec probabilité p=1/1500p=1/1500, indépendamment ; XB(1000,1/1500)X\sim\mathcal B(1000,1/1500). 2) E(X)=np=1000/1500=2/3E(X)=np=1000/1500=2/3 ; var(X)=np(1p)=23149915000.666\mathrm{var}(X)=np(1-p)=\tfrac23\cdot\tfrac{1499}{1500}\approx0.666. 3) nn grand, pp petit, np=2/3np=2/3 modéré : l'approximation de Poisson P(λ)\mathcal P(\lambda) avec λ=np=2/3\lambda=np=2/3 est légitime. 4) P(X=0)e2/30.513P(X=0)\approx e^{-2/3}\approx0.513. 5) P(X>2)=1P(X2)=1e2/3(1+23+(2/3)22)=1e2/3(1+0.667+0.222)10.5131.8890.030P(X>2)=1-P(X\le2)=1-e^{-2/3}(1+\tfrac23+\tfrac{(2/3)^2}2)=1-e^{-2/3}(1+0.667+0.222)\approx1-0.513\cdot1.889\approx0.030.

التمرين 1

Convergence de sin(nx)/(pi x) vers la masse de Dirac

#distributions#Riemann-Lebesgue#Dirac

On considère φD\varphi\in\mathcal{D} de support [a,b][a,b].

  1. Montrer que les intégrales
In=absin(nx)φ(x)dx,Jn=abcos(nx)φ(x)dxI_n=\int_a^b\sin(nx)\,\varphi(x)\,dx,\qquad J_n=\int_a^b\cos(nx)\,\varphi(x)\,dx

tendent vers 00 quand nn\to\infty. 2. En déduire que la distribution sin(nx)πx\dfrac{\sin(nx)}{\pi x} converge vers la distribution δ\delta de Dirac. (Indication : +sinttdt=π\int_{-\infty}^{+\infty}\frac{\sin t}{t}\,dt=\pi.)

الحل
  1. C'est le lemme de Riemann-Lebesgue : φ\varphi étant CC^\infty à support compact, une intégration par parties donne In,Jn=O(1/n)0I_n,J_n=O(1/n)\to0. 2. sin(nx)πx,φ=1πsin(nx)xφ(x)dx\langle\frac{\sin(nx)}{\pi x},\varphi\rangle=\frac1\pi\int\frac{\sin(nx)}{x}\varphi(x)\,dx ; par le changement t=nxt=nx et sinttdt=π\int\frac{\sin t}{t}dt=\pi, la partie φ(0)\varphi(0) donne φ(0)\varphi(0), le reste tend vers 00 (Riemann-Lebesgue sur φ(x)φ(0)x\frac{\varphi(x)-\varphi(0)}{x}). Donc sin(nx)πxδ\frac{\sin(nx)}{\pi x}\to\delta.

التمرين 2

Couple (U,V)=(X₁X₂, X₁/X₂) pour une densité en 1/x²

#joint-density#change-of-variables#conditional-density#independence

Soit XX de densité fX(x)=1x21[1,+[(x)f_X(x)=\frac1{x^2}\mathbf1_{[1,+\infty[}(x). Soient X1,X2X_1,X_2 indépendantes de même loi que XX, U=X1X2U=X_1X_2, V=X1/X2V=X_1/X_2. 1) Déterminer la densité du couple (U,V)(U,V) et tracer son domaine. 2) Déterminer les densités de UU et de VV. 3) Déterminer la densité conditionnelle de VV sachant U=uU=u. 4) UU et VV sont-elles indépendantes ? 5) On pose Z=UZ=\sqrt U ; calculer la densité de ZZ.

الحل

La densité jointe de (X1,X2)(X_1,X_2) est 1x12x22\frac1{x_1^2x_2^2} sur [1,)2[1,\infty)^2. Le changement u=x1x2u=x_1x_2, v=x1/x2v=x_1/x_2 donne x1=uvx_1=\sqrt{uv}, x2=u/vx_2=\sqrt{u/v}, jacobien (x1,x2)/(u,v)=12v|\partial(x_1,x_2)/\partial(u,v)|=\frac1{2v}. Le domaine est u1u\ge1, 1/uvu1/u\le v\le u (car x1,x21x_1,x_2\ge1). 1) fU,V(u,v)=1(uv)(u/v)12v=12u2vf_{U,V}(u,v)=\frac1{(uv)(u/v)}\cdot\frac1{2v}=\frac1{2u^2v} sur ce domaine. 2) fU(u)=1/uu12u2vdv=12u2[lnv]1/uu=lnuu2f_U(u)=\int_{1/u}^u\frac1{2u^2v}dv=\frac1{2u^2}[\ln v]_{1/u}^u=\frac{\ln u}{u^2}, u1u\ge1. La densité de VV s'obtient par symétrie : fV(v)=12v2f_V(v)=\frac1{2v^2} pour v1v\ge1 et 12\frac12 pour 0<v10<v\le1 (loi de X1/X2X_1/X_2). 3) fVU=u(v)=fU,V(u,v)fU(u)=12vlnuf_{V|U=u}(v)=\frac{f_{U,V}(u,v)}{f_U(u)}=\frac1{2v\ln u} sur [1/u,u][1/u,u]. 4) Elle dépend de uu via le support, donc UU et VV ne sont pas indépendantes. 5) Z=UZ=\sqrt U, u=z2u=z^2, fZ(z)=fU(z2)2z=ln(z2)z42z=4lnzz3f_Z(z)=f_U(z^2)\cdot2z=\frac{\ln(z^2)}{z^4}2z=\frac{4\ln z}{z^3} pour z1z\ge1.

التمرين 2

Équation intégrale de Volterra résolue par une EDO linéaire

#équation intégrale#EDO linéaire#Volterra

A) Soit fC1([a,b])f\in C^1([a,b]) avec 0[a,b]0\in[a,b] et λR\lambda\in\mathbb{R}^*. On considère

g(x)λ0xg(t)dt=f(x).(A)g(x)-\lambda\int_0^x g(t)\,dt=f(x).\tag{A}
  1. Montrer que gg vérifie une équation différentielle linéaire du premier ordre que l'on intégrera.
  2. En déduire que (A) admet une solution unique g(x)=f(x)+λ0xeλ(xt)f(t)dtg(x)=f(x)+\lambda\int_0^x e^{\lambda(x-t)}f(t)\,dt.

B) 1. Déduire que g(x)λ0xeμ(xt)g(t)dt=f(x)g(x)-\lambda\int_0^x e^{\mu(x-t)}g(t)\,dt=f(x) admet une solution unique. 2. Calculer la solution pour f(x)=eμxf(x)=e^{\mu x}.

الحل

A1. En dérivant (A) : g(x)λg(x)=f(x)g'(x)-\lambda g(x)=f'(x), avec g(0)=f(0)g(0)=f(0). A2. Facteur intégrant eλxe^{-\lambda x} : g(x)=f(x)+λ0xeλ(xt)f(t)dtg(x)=f(x)+\lambda\int_0^x e^{\lambda(x-t)}f(t)\,dt. B1. Même méthode : la solution unique est g(x)=f(x)+λ0xe(μ+λ)(xt)f(t)dtg(x)=f(x)+\lambda\int_0^x e^{(\mu+\lambda)(x-t)}f(t)\,dt. B2. Pour f=eμxf=e^{\mu x} : 0xe(μ+λ)(xt)eμtdt=e(μ+λ)xeμxλ\int_0^x e^{(\mu+\lambda)(x-t)}e^{\mu t}dt=\frac{e^{(\mu+\lambda)x}-e^{\mu x}}{\lambda}, d'où g(x)=eμx+(e(μ+λ)xeμx)=e(μ+λ)xg(x)=e^{\mu x}+\bigl(e^{(\mu+\lambda)x}-e^{\mu x}\bigr)=e^{(\mu+\lambda)x}.

التمرين 3

Estimation du paramètre d'une loi de densité θx^(θ-1)

#parametric-estimation#method-of-moments#maximum-likelihood#efficiency

Soit XX de densité f(x,θ)=θxθ11[0,1](x)f(x,\theta)=\theta x^{\theta-1}\mathbf1_{[0,1]}(x), θ>0\theta>0. 1) Déterminer par la méthode des moments l'estimateur de θ\theta au vu de l'échantillon (X1,,Xn)(X_1,\dots,X_n). 2) Déterminer par le maximum de vraisemblance l'estimateur T1/θT_{1/\theta} de 1/θ1/\theta. 3) T1/θT_{1/\theta} est-il sans biais ? 4) Si oui, est-il efficace ?

الحل

1) E(X)=01xθxθ1dx=θθ+1E(X)=\int_0^1x\theta x^{\theta-1}dx=\frac\theta{\theta+1}. En égalant à Xˉ\bar X : θ^θ^+1=Xˉ\frac{\hat\theta}{\hat\theta+1}=\bar X, donc θ^MM=Xˉ1Xˉ\hat\theta_{MM}=\frac{\bar X}{1-\bar X}. 2) La log-vraisemblance (θ)=nlnθ+(θ1)lnXi\ell(\theta)=n\ln\theta+(\theta-1)\sum\ln X_i ; (θ)=nθ+lnXi=0\ell'(\theta)=\frac n\theta+\sum\ln X_i=0 donne θ^MV=nlnXi\hat\theta_{MV}=-\frac n{\sum\ln X_i}, donc l'estimateur de 1/θ1/\theta est T1/θ=1nlnXiT_{1/\theta}=-\frac1n\sum\ln X_i. 3) lnXiE(θ)-\ln X_i\sim\mathcal E(\theta) (loi exponentielle), d'espérance 1/θ1/\theta, donc E(T1/θ)=1/θE(T_{1/\theta})=1/\theta : sans biais. 4) var(T1/θ)=1n1θ2=1nθ2\mathrm{var}(T_{1/\theta})=\frac1{n}\cdot\frac1{\theta^2}=\frac1{n\theta^2}. L'information de Fisher pour 1/θ1/\theta atteint cette borne (modèle exponentiel), donc T1/θT_{1/\theta} est efficace.

التمرين 4

Intervalle de confiance pour la superficie moyenne des états

#confidence-interval#student-distribution#sampling#estimation

La population des USA compte 50 états. Un échantillon aléatoire de 5 états a fourni les superficies (en milliers de km²) : 147, 84, 24, 85, 159. 1) Trouver un intervalle de confiance à 95% pour la superficie moyenne. 2) Déterminer l'intervalle de confiance à 95% pour la superficie totale des USA. 3) La superficie totale réelle est 3620000 km² ; est-elle dans l'intervalle ?

الحل

1) Moyenne xˉ=147+84+24+85+1595=4995=99.8\bar x=\frac{147+84+24+85+159}5=\frac{499}5=99.8. Variance corrigée s2=14(xixˉ)2s^2=\frac1{4}\sum(x_i-\bar x)^2 ; les écarts donnent s52.5s\approx52.5. Avec t0.975,4=2.776t_{0.975,4}=2.776, l'IC est 99.8±2.77652.55=99.8±65.299.8\pm2.776\cdot\frac{52.5}{\sqrt5}=99.8\pm65.2, soit environ [34.6,165.0][34.6,165.0] milliers de km². 2) Superficie totale =50×=50\times moyenne, IC : [5034.6,50165.0]=[1730,8250][50\cdot34.6,50\cdot165.0]=[1730,8250] milliers de km², soit [1730000,8250000][1\,730\,000,8\,250\,000] km². 3) 36200003\,620\,000 km² appartient à cet intervalle : oui, l'aire réelle est comprise dans l'intervalle de confiance.