📚 الرئيسية

مسابقة دكتوراه 2021Université Mohamed El Bachir El Ibrahimi - Bordj Bou Arréridj — الموضوع 01

مسابقة تخصص · Systèmes Dynamiques · المعامل: 1 · المدة: 1سا 30د

Concours d'accès au Doctorat LMD 2020-2021, filière Mathématiques Appliquées, épreuve 1 « Équations Différentielles Ordinaires » (coefficient 1), Université Mohamed El Bachir El Ibrahimi de Bordj Bou Arréridj, Faculté des Mathématiques et de l'Informatique, le 06 mars 2021, durée 1h30.

التمرين 1

Exercice 1 — Non-unicité pour y' = ∛(y²)

#cauchy-problem#non-uniqueness#separable-ode#peano-existence

(04 pts) On considère l'équation différentielle

y=y23y'=\sqrt[3]{y^{2}}

  1. Montrer que l'unicité des solutions n'est pas assurée sur l'axe OxOx.
  2. Déterminer la solution de cette équation pour y0y\neq 0.
الحل

1.

La fonction f(y)=y2/3f(y)=y^{2/3} est continue sur R\mathbb{R} (existence de solutions par le théorème de Peano), mais elle n'est pas lipschitzienne au voisinage de y=0y=0 : en effet

f(y)f(0)y0=y1/3y0+\frac{|f(y)-f(0)|}{|y-0|}=|y|^{-1/3}\xrightarrow[y\to 0]{}+\infty

donc le théorème de Cauchy–Lipschitz ne s'applique pas aux points (x0,0)(x_{0},0) de l'axe OxOx. Concrètement, par le point (x0,0)(x_{0},0) passent au moins deux solutions distinctes :

y10ety2(x)=(xx0)327y_{1}\equiv 0\qquad\text{et}\qquad y_{2}(x)=\frac{(x-x_{0})^{3}}{27}

(on vérifie que y2=(xx0)29=(y22)1/3y_{2}'=\frac{(x-x_{0})^{2}}{9}=\bigl(y_{2}^{2}\bigr)^{1/3}).

L’uniciteˊ est violeˊe en tout point de l’axe Ox\boxed{\text{L'unicité est violée en tout point de l'axe }Ox}

2.

Pour y0y\neq 0, l'équation est à variables séparables :

y2/3dy=dx3y1/3=x+Cy^{-2/3}\,dy=dx\quad\Longrightarrow\quad 3\,y^{1/3}=x+C

y(x)=(x+C3)3,CR\boxed{y(x)=\Bigl(\frac{x+C}{3}\Bigr)^{3},\qquad C\in\mathbb{R}}

التمرين 2

Exercice 2 — Problème de Cauchy à variables séparables

#separable-ode#initial-value-problem#integration

(04 pts) Résoudre le problème de Cauchy suivant :

yy+ex(yy1)=0,y(0)=1y\,y'+e^{x}\bigl(y\,y'-1\bigr)=0,\qquad y(0)=1

الحل

1.

L'équation s'écrit yy(1+ex)=exy\,y'\,(1+e^{x})=e^{x}, soit

yy=ex1+exddx(y22)=ex1+exy\,y'=\frac{e^{x}}{1+e^{x}}\quad\Longleftrightarrow\quad\frac{d}{dx}\Bigl(\frac{y^{2}}{2}\Bigr)=\frac{e^{x}}{1+e^{x}}

En intégrant :

y22=ln(1+ex)+C\frac{y^{2}}{2}=\ln\bigl(1+e^{x}\bigr)+C

La condition y(0)=1y(0)=1 donne 12=ln2+C\frac{1}{2}=\ln 2+C, donc C=12ln2C=\frac{1}{2}-\ln 2 et

y2=1+2ln1+ex2y^{2}=1+2\ln\frac{1+e^{x}}{2}

Comme y(0)=1>0y(0)=1\gt 0, on prend la racine positive :

y(x)=1+2ln1+ex2\boxed{y(x)=\sqrt{1+2\ln\frac{1+e^{x}}{2}}}

التمرين 3

Exercice 3 — Système linéaire X' = AX avec paramètre α

#linear-systems#eigenvalues#parameter-discussion#matrix-exponential

(12 pts) Soit le système différentiel X(t)=AX(t)X'(t)=A\,X(t), t0t\geq 0, où αR\alpha\in\mathbb{R} et AA est la matrice de lignes L1=(α, 1)L_{1}=(\alpha,\ 1) et L2=(α, α4)L_{2}=(\alpha,\ \alpha-4) :

x=αx+y,y=αx+(α4)yx'=\alpha x+y,\qquad y'=\alpha x+(\alpha-4)y

Étudier et résoudre le système suivant les valeurs du paramètre α\alpha (valeurs propres, diagonalisation éventuelle, solution générale). (La page des questions détaillées étant absente du document scanné, on résout complètement le système selon α\alpha.)

الحل

1.

Polynôme caractéristique :

χA(λ)=λ2(2α4)λ+(α(α4)α)=λ2(2α4)λ+α25α\chi_{A}(\lambda)=\lambda^{2}-(2\alpha-4)\lambda+\bigl(\alpha(\alpha-4)-\alpha\bigr)=\lambda^{2}-(2\alpha-4)\lambda+\alpha^{2}-5\alpha

Discriminant :

Δ=(2α4)24(α25α)=4α+16=4(α+4)\Delta=(2\alpha-4)^{2}-4(\alpha^{2}-5\alpha)=4\alpha+16=4(\alpha+4)

Valeurs propres :

λ±=(α2)±α+4(α4)\lambda_{\pm}=(\alpha-2)\pm\sqrt{\alpha+4}\quad(\alpha\geq -4)

2.

Cas α>4\alpha\gt -4 : deux valeurs propres réelles distinctes, AA est diagonalisable. Un vecteur propre associé à λ\lambda est vλ=(1, λα)v_{\lambda}=(1,\ \lambda-\alpha) (première ligne : (αλ)+(λα)=0(\alpha-\lambda)+(\lambda-\alpha)=0). Solution générale :

X(t)=C1eλ+t(1λ+α)+C2eλt(1λα),λ±=α2±α+4\boxed{X(t)=C_{1}e^{\lambda_{+}t}\begin{pmatrix}1\\ \lambda_{+}-\alpha\end{pmatrix}+C_{2}e^{\lambda_{-}t}\begin{pmatrix}1\\ \lambda_{-}-\alpha\end{pmatrix},\quad\lambda_{\pm}=\alpha-2\pm\sqrt{\alpha+4}}

3.

Cas α=4\alpha=-4 : valeur propre double λ=α2=6\lambda=\alpha-2=-6, et A+6IA+6I (lignes (2, 1)(2,\ 1) et (4, 2)(-4,\ -2)) est de rang 1 : AA n'est pas diagonalisable. Vecteur propre v=(1, 2)v=(1,\ -2) ; vecteur de Jordan w=(0, 1)w=(0,\ 1) vérifiant (A+6I)w=v(A+6I)w=v. Solution générale :

X(t)=e6t[C1(12)+C2(t(12)+(01))]\boxed{X(t)=e^{-6t}\Bigl[C_{1}\begin{pmatrix}1\\ -2\end{pmatrix}+C_{2}\Bigl(t\begin{pmatrix}1\\ -2\end{pmatrix}+\begin{pmatrix}0\\ 1\end{pmatrix}\Bigr)\Bigr]}

4.

Cas α<4\alpha\lt -4 : valeurs propres complexes conjuguées λ=(α2)±iω\lambda=(\alpha-2)\pm i\omega avec ω=(α+4)\omega=\sqrt{-(\alpha+4)}. En prenant les parties réelle et imaginaire de eλtvλe^{\lambda t}v_{\lambda} avec vλ=(1, λα)=(1, 2+iω)v_{\lambda}=(1,\ \lambda-\alpha)=(1,\ -2+i\omega) :

X(t)=e(α2)t[C1(cosωt2cosωtωsinωt)+C2(sinωt2sinωt+ωcosωt)]\boxed{X(t)=e^{(\alpha-2)t}\Bigl[C_{1}\begin{pmatrix}\cos\omega t\\ -2\cos\omega t-\omega\sin\omega t\end{pmatrix}+C_{2}\begin{pmatrix}\sin\omega t\\ -2\sin\omega t+\omega\cos\omega t\end{pmatrix}\Bigr]}

Dans tous les cas, C1,C2C_{1},C_{2} sont déterminées par la condition initiale X(0)X(0).