الرئيسية

مسابقة دكتوراه 2015Université Djilali Liabès - Sidi Bel Abbès — الموضوع 02

مسابقة تخصص · Probabilités & Statistiques

Concours d'accès à l'École doctorale « Calcul Stochastique, Statistique et Applications », 1ère Épreuve : Recherche opérationnelle et Équations différentielles stochastiques et Processus stochastique, Sujet 2, Faculté des Sciences Exactes, Université Djillali Liabès de Sidi Bel Abbès, année universitaire 2014/2015.

التمرين 1

Exercice 01 — Optimisation d'un profit de vente par la méthode du simplexe

#operations-research#linear-programming#simplex-method#duality#optimization

Avant l'arrivage massif de nouveaux modèles, un vendeur de téléphones portables veut écouler son stock composé de huit appareils, quatre kits « mains libres » et dix-neuf cartes prépayées. Il peut proposer un téléphone avec deux cartes (profit net de 7 dinars), ou préparer un coffret composé d'un téléphone, d'un kit « mains libres » et de trois cartes (profit net de 9 dinars). Il peut vendre n'importe quelle quantité de ces offres dans la limite du stock. Quelle quantité de chaque offre doit-il préparer pour maximiser son profit net ? (Méthode du simplexe.)

الحل

Exercice 01

Modélisation. Soit xx le nombre d'offres « téléphone ++ 2 cartes » et yy le nombre de coffrets « téléphone ++ 1 kit ++ 3 cartes ». Les stocks donnent

max P=7x+9ys.c.{x+y8(teˊleˊphones)y4(kits)2x+3y19(cartes)x,y0.\max\ P=7x+9y\quad\text{s.c.}\quad\begin{cases}x+y\leq 8 &(\text{téléphones})\\ y\leq 4 &(\text{kits})\\ 2x+3y\leq 19 &(\text{cartes})\\ x,y\geq 0.\end{cases}

Simplexe (variables d'écart s1,s2,s3s_1,s_2,s_3) :

  • yy entre (coût réduit 9-9), s2s_2 sort (min(8,4,19/3)=4\min(8,4,19/3)=4) : y=4y=4.
  • xx entre (coût réduit 7-7), s3s_3 sort : x=72x=\tfrac72.
  • s2s_2 ré-entre (coût réduit 32>0\tfrac32\gt 0), s1s_1 sort (min=1\min=1) : on obtient x=5, y=3x=5,\ y=3.

Tous les coûts réduits sont alors 0\geq 0 : optimum atteint.

x=5 offres,  y=3 coffrets,Pmax=75+93=62 dinars.\boxed{x=5\ \text{offres, }\ y=3\ \text{coffrets},\qquad P_{\max}=7\cdot 5+9\cdot 3=62\ \text{dinars}.}

Vérification des stocks : téléphones 5+3=85+3=8 (épuisés), cartes 25+33=192\cdot5+3\cdot3=19 (épuisées), kits 343\leq 4.

التمرين 2

Exercice 02 — EDS linéaire générale : existence, moments et variance

#stochastic-calculus#stochastic-differential-equation#ito-formula#existence-uniqueness#linear-sde

Soient a,α,b,βa,\alpha,b,\beta quatre constantes réelles. On considère l'EDS

dXt=(a+αXt)dt+(b+βXt)dBt,X0=x.dX_t=(a+\alpha X_t)\,dt+(b+\beta X_t)\,dB_t,\qquad X_0=x.

  1. Montrer que cette équation admet une solution unique.
  2. On note m(t)=E(Xt)m(t)=E(X_t) et M(t)=E(Xt2)M(t)=E(X_t^2). (i) Montrer que m(t)m(t) est l'unique solution de yαy=a, y(0)=xy'-\alpha y=a,\ y(0)=x. (ii) Écrire la formule d'Itô pour Xt2X_t^2. En déduire que M(t)M(t) est l'unique solution de y(2α+β2)y=2(a+bβ)m(t)+b2, y(0)=x2y'-(2\alpha+\beta^2)y=2(a+b\beta)m(t)+b^2,\ y(0)=x^2. (iii) En déduire E(Xt)E(X_t) et Var(Xt)\mathrm{Var}(X_t) dans le cas a=β=0a=\beta=0.
الحل

Exercice 02

1. Existence et unicité

Les coefficients μ(x)=a+αx\mu(x)=a+\alpha x et σ(x)=b+βx\sigma(x)=b+\beta x sont affines en xx. Ils sont donc globalement lipschitziens,

μ(x)μ(y)=αxy,σ(x)σ(y)=βxy,|\mu(x)-\mu(y)|=|\alpha|\,|x-y|,\qquad |\sigma(x)-\sigma(y)|=|\beta|\,|x-y|,

et vérifient la condition de croissance linéaire μ(x)+σ(x)C(1+x)|\mu(x)|+|\sigma(x)|\leq C(1+|x|). Le théorème d'existence et d'unicité des solutions fortes d'EDS s'applique :

l’EDS admet une unique solution forte.\boxed{\text{l'EDS admet une unique solution forte.}}

2. (i)

Sous forme intégrale, Xt=x+0t(a+αXs)ds+0t(b+βXs)dBsX_t=x+\int_0^t(a+\alpha X_s)ds+\int_0^t(b+\beta X_s)dB_s. En prenant l'espérance (l'intégrale d'Itô est centrée) :

m(t)=x+0t(a+αm(s))ds  m(t)=αm(t)+a,m(0)=x,m(t)=x+\int_0^t\big(a+\alpha m(s)\big)ds\ \Longrightarrow\ \boxed{m'(t)=\alpha m(t)+a,\quad m(0)=x,}

soit yαy=a, y(0)=xy'-\alpha y=a,\ y(0)=x (équation linéaire, solution unique).

2. (ii)

Formule d'Itô pour Xt2X_t^2 :

d(Xt2)=2XtdXt+dXt=[2Xt(a+αXt)+(b+βXt)2]dt+2Xt(b+βXt)dBt.d(X_t^2)=2X_t\,dX_t+d\langle X\rangle_t=\Big[2X_t(a+\alpha X_t)+(b+\beta X_t)^2\Big]dt+2X_t(b+\beta X_t)\,dB_t.

Le terme de dérive se développe en (2α+β2)Xt2+2(a+bβ)Xt+b2(2\alpha+\beta^2)X_t^2+2(a+b\beta)X_t+b^2. En prenant l'espérance :

M(t)=(2α+β2)M(t)+2(a+bβ)m(t)+b2,M(0)=x2,\boxed{M'(t)=(2\alpha+\beta^2)M(t)+2(a+b\beta)m(t)+b^2,\quad M(0)=x^2,}

soit y(2α+β2)y=2(a+bβ)m(t)+b2, y(0)=x2y'-(2\alpha+\beta^2)y=2(a+b\beta)m(t)+b^2,\ y(0)=x^2.

2. (iii) Cas a=β=0a=\beta=0

Alors m(t)=αm(t), m(0)=xm(t)=xeαtm'(t)=\alpha m(t),\ m(0)=x\Rightarrow m(t)=x\,e^{\alpha t}, donc E(Xt)=xeαtE(X_t)=x\,e^{\alpha t}.

Pour MM : M2αM=b2, M(0)=x2M'-2\alpha M=b^2,\ M(0)=x^2, d'où

M(t)=e2αt(x2+b20te2αsds)=x2e2αt+b22α(e2αt1).M(t)=e^{2\alpha t}\Big(x^2+b^2\int_0^t e^{-2\alpha s}ds\Big)=x^2 e^{2\alpha t}+\frac{b^2}{2\alpha}\big(e^{2\alpha t}-1\big).

Donc

E(Xt)=xeαt,Var(Xt)=M(t)m(t)2=b22α(e2αt1).\boxed{E(X_t)=x\,e^{\alpha t},\qquad \mathrm{Var}(X_t)=M(t)-m(t)^2=\frac{b^2}{2\alpha}\big(e^{2\alpha t}-1\big).}