📚 الرئيسية

مسابقة دكتوراه 2018Université Badji Mokhtar - Annaba — الموضوع 02

مسابقة تخصص · Systèmes Dynamiques

Concours de la formation doctorale « Mathématiques appliquées », Épreuve d'EDO, Université Badji Mokhtar - Annaba, 20/10/2018. (Copie manuscrite ; variante distincte de l'épreuve d'EDO tapée.)

التمرين 1

Stabilité de la solution nulle d'un système linéaire à coefficients variables

#linear-system#stability#fundamental-matrix#bounded-solutions#non-autonomous

Soit le système différentiel dxdt=A(t)x(t),\frac{dx}{dt}=A(t)\,x(t),x(t)Rnx(t)\in\mathbb{R}^n et A(t)A(t) est une matrice n×nn\times n continue. Montrer que la solution nulle de ce système est stable si et seulement si toutes les solutions sont bornées sur [0,+[[0,+\infty[.

الحل

Matrice fondamentale

Soit R(t)R(t) la matrice fondamentale (résolvante) vérifiant R(t)=A(t)R(t)R'(t)=A(t)R(t), R(0)=IdR(0)=\mathrm{Id}. Toute solution s'écrit x(t)=R(t)x0x(t)=R(t)x_0.

(\Leftarrow) Solutions bornées \Rightarrow stabilité

Si toutes les solutions sont bornées, alors chaque colonne de R(t)R(t) est bornée, donc supt0R(t)=:M<\sup_{t\ge0}\|R(t)\|=:M<\infty (principe de la borne uniforme, ou simplement finitude en dimension finie). Pour ε>0\varepsilon>0, en prenant x0<δ=ε/M\|x_0\|<\delta=\varepsilon/M : x(t)=R(t)x0Mx0<ε,t0.\|x(t)\|=\|R(t)x_0\|\le M\|x_0\|<\varepsilon,\quad\forall t\ge0. Donc la solution nulle est stable.

(\Rightarrow) Stabilité \Rightarrow solutions bornées

Si la solution nulle est stable, il existe δ>0\delta>0 tel que x0δx(t)1\|x_0\|\le\delta\Rightarrow\|x(t)\|\le1 pour tout t0t\ge0. Par linéarité, pour un x0x_0 quelconque non nul, x~0=δx0/x0\tilde x_0=\delta x_0/\|x_0\| donne une solution bornée par 11, d'où R(t)x0=x0δR(t)x~0x0δ.\|R(t)x_0\|=\frac{\|x_0\|}{\delta}\,\|R(t)\tilde x_0\|\le\frac{\|x_0\|}{\delta}. Toute solution est donc bornée. Solution nulle stable    toutes les solutions borneˊes.\boxed{\text{Solution nulle stable}\iff\text{toutes les solutions bornées.}} (Spécificité du cas linéaire : stabilité = bornitude de la résolvante.)

التمرين 2

Équation de Clairaut y = x y' + √(b²+a² y'²)

#clairaut-equation#singular-solution#envelope#first-order-ode

Résoudre l'équation différentielle y=xdydx+b2+a2(dydx)2.y=x\,\frac{dy}{dx}+\sqrt{b^2+a^2\Big(\frac{dy}{dx}\Big)^2}.

الحل

Forme de Clairaut

Avec p=yp=y', l'équation est y=xp+b2+a2p2y=xp+\sqrt{b^2+a^2p^2}, de type Clairaut y=xp+f(p)y=xp+f(p) avec f(p)=b2+a2p2f(p)=\sqrt{b^2+a^2p^2}.

Dérivation

En dérivant par rapport à xx : p=p+xp+f(p)p  p(x+f(p))=0.p=p+xp'+f'(p)p'\ \Rightarrow\ p'\big(x+f'(p)\big)=0.

Solutions générales (droites)

p=0p=Cp'=0\Rightarrow p=C, d'où la famille de droites y=Cx+b2+a2C2,CR.\boxed{y=Cx+\sqrt{b^2+a^2C^2},\qquad C\in\mathbb{R}.}

Solution singulière (enveloppe)

x+f(p)=0x+f'(p)=0 avec f(p)=a2pb2+a2p2f'(p)=\dfrac{a^2p}{\sqrt{b^2+a^2p^2}}, soit x=a2pb2+a2p2x=-\dfrac{a^2p}{\sqrt{b^2+a^2p^2}}. En reportant y=xp+b2+a2p2=b2b2+a2p2y=xp+\sqrt{b^2+a^2p^2}=\dfrac{b^2}{\sqrt{b^2+a^2p^2}}. On élimine pp : de xx on tire x2a2=a2p2b2+a2p2\dfrac{x^2}{a^2}=\dfrac{a^2p^2}{b^2+a^2p^2} et y2b2=b2b2+a2p2\dfrac{y^2}{b^2}=\dfrac{b^2}{b^2+a^2p^2}, dont la somme vaut 11 : x2a2+y2b2=1.\boxed{\frac{x^2}{a^2}+\frac{y^2}{b^2}=1.} L'enveloppe est une ellipse (solution singulière).

التمرين 3

Solution polynomiale de (1−2x²)y''+2y'+4y=0

#linear-ode#polynomial-solution#power-series#reduction-of-order

Chercher une solution particulière sous forme d'un polynôme de l'équation (12x2)d2ydx2+2dydx+4y=0.(1-2x^2)\,\frac{d^2y}{dx^2}+2\,\frac{dy}{dx}+4y=0.

الحل

Recherche polynomiale

On cherche y=k=0Nakxky=\sum_{k=0}^{N}a_kx^k. Injectons y=a0+a1x+a2x2+y=a_0+a_1x+a_2x^2+\cdots : y=kakxk1,y=k(k1)akxk2.y'=\sum k a_kx^{k-1},\qquad y''=\sum k(k-1)a_kx^{k-2}. Le terme (12x2)y(1-2x^2)y'' produit k(k1)akxk22k(k1)akxk\sum k(k-1)a_kx^{k-2}-2\sum k(k-1)a_kx^{k}.

Relation de récurrence

En regroupant le coefficient de xkx^k : (k+2)(k+1)ak+2+2(k+1)ak+1+(42k(k1))ak=0.(k+2)(k+1)a_{k+2}+2(k+1)a_{k+1}+\big(4-2k(k-1)\big)a_k=0. Cherchons un polynôme de bas degré. Le facteur 42k(k1)4-2k(k-1) s'annule pour k(k1)=2k(k-1)=2, i.e. k=2k=2. Essayons y=a0+a1x+a2x2y=a_0+a_1x+a_2x^2 (degré 2) :

  • coeff de x2x^2 dans l'équation : 221a2+4a2=4a2+4a2=0 -2\cdot2\cdot1\,a_2+4a_2=-4a_2+4a_2=0\ \checkmark (automatique) ;
  • coeff de x0x^0 : 2a2+2a1+4a0=02a_2+2a_1+4a_0=0 ;
  • coeff de x1x^1 : 22a2 (de y)?2\cdot2\,a_2\ (\text{de }y'')? Reprenons proprement avec y=a0+a1x+a2x2y=a_0+a_1x+a_2x^2 : y=a1+2a2xy'=a_1+2a_2x, y=2a2y''=2a_2. (12x2)(2a2)+2(a1+2a2x)+4(a0+a1x+a2x2)=0.(1-2x^2)(2a_2)+2(a_1+2a_2x)+4(a_0+a_1x+a_2x^2)=0. Développons : 2a24a2x2+2a1+4a2x+4a0+4a1x+4a2x22a_2-4a_2x^2+2a_1+4a_2x+4a_0+4a_1x+4a_2x^2.
  • x2x^2 : 4a2+4a2=0 -4a_2+4a_2=0\ \checkmark
  • x1x^1 : 4a2+4a1=0a1=a24a_2+4a_1=0\Rightarrow a_1=-a_2
  • x0x^0 : 2a2+2a1+4a0=02a22a2+4a0=0a0=0.2a_2+2a_1+4a_0=0\Rightarrow 2a_2-2a_2+4a_0=0\Rightarrow a_0=0.

Solution

Avec a2=1a_2=1 : a1=1a_1=-1, a0=0a_0=0, donc y(x)=x2x\boxed{y(x)=x^2-x} est une solution polynomiale. (Vérification : y=2y''=2, (12x2)2+2(2x1)+4(x2x)=24x2+4x2+4x24x=0 (1-2x^2)\cdot2+2(2x-1)+4(x^2-x)=2-4x^2+4x-2+4x^2-4x=0\ \checkmark.)