📚 الرئيسية

مسابقة دكتوراه 2016Université de Ghardaïa — الموضوع 01

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

Concours d'accès à la formation de troisième cycle, année universitaire 2016-2017, Spécialité : Équations différentielles et aux dérivées partielles, Épreuve 2 : EDO-EDP, Université de Ghardaïa, durée 1h30.

التمرين 1

Exercice 1 — Orbites fermées dans le plan

#dynamical-systems#closed-orbits#poincare-bendixson#planar-systems

Les systèmes suivants admettent-ils des orbites fermées non triviales dans le plan ? Justifier.

(A){x˙=x2+yy˙=xy1(B){x˙=6x+πyy˙=πx+y(A) \begin{cases} \dot{x} = -x^2 + y \\ \dot{y} = x - y - 1 \end{cases} \qquad (B) \begin{cases} \dot{x} = 6x + \pi y \\ \dot{y} = -\pi x + y \end{cases}

الحل

Système (A).

Critère de Bendixson : div(F)=(x2+y)x+(xy1)y=2x1\mathrm{div}(F) = \frac{\partial(-x^2+y)}{\partial x} + \frac{\partial(x-y-1)}{\partial y} = -2x - 1. Cette expression n'est pas de signe constant sur R2\mathbb{R}^2, donc Bendixson ne conclut pas directement. On cherche les points fixes : x2+y=0-x^2+y=0 et xy1=0x-y-1=0 donnent xx21=0x-x^2-1=0, discriminant négatif — pas de points fixes réels. Par Poincaré-Bendixson, sans point fixe et avec une orbite bornée, il pourrait y avoir une orbite fermée. Cependant, l'absence de point fixe ne garantit pas l'existence. Analyse plus fine nécessaire.

Système (B).

La matrice est (6ππ1)\begin{pmatrix} 6 & \pi \\ -\pi & 1 \end{pmatrix}. tr(A)=7>0\mathrm{tr}(A) = 7 > 0, donc le point fixe (0,0)(0,0) est instable. Par Bendixson-Dulac, div(F)=6+1=7>0\mathrm{div}(F) = 6+1 = 7 > 0 : constant positif, donc il n'existe aucune orbite fermée dans R2\mathbb{R}^2.

(A):analyse neˊcessaire;(B):aucune orbite fermeˊe (div>0).\boxed{(A) : \text{analyse nécessaire}; \quad (B) : \text{aucune orbite fermée (div}>0).}

التمرين 2

Exercice 2 — Équation différentielle à retard

#delay-differential-equation#piecewise-solution#retarded-ode

En considérant les cas a=0a = 0 et a0a \neq 0, résoudre sur l'intervalle [0,2][0, 2] l'équation différentielle à retard suivante :

x˙(t)=ax(t)+bx(t1)\dot{x}(t) = ax(t) + bx(t-1)

x(t)=θ(t)=1+tx(t) = \theta(t) = 1 + t pour t[1,0]t \in [-1, 0], a=ctea = \mathrm{cte} et b=cte0b = \mathrm{cte} \neq 0.

Indication : considérer l'équation sur chacun des intervalles [0,1][0,1] et [1,2][1,2].

الحل

Cas a=0a = 0.

Sur [0,1][0,1] : x˙(t)=bx(t1)=b(1+(t1))=bt\dot{x}(t) = bx(t-1) = b(1+(t-1)) = bt. Donc x(t)=x(0)+bt22=1+b2t2x(t) = x(0) + b\frac{t^2}{2} = 1 + \frac{b}{2}t^2.

Sur [1,2][1,2] : x˙(t)=bx(t1)=b(1+b2(t1)2)\dot{x}(t) = bx(t-1) = b(1 + \frac{b}{2}(t-1)^2). On intègre : x(t)=x(1)+b1t(1+b2(s1)2)ds=(1+b2)+b(t1)+b26(t1)3x(t) = x(1) + b\int_1^t (1+\frac{b}{2}(s-1)^2)ds = (1+\frac{b}{2}) + b(t-1) + \frac{b^2}{6}(t-1)^3.

Cas a0a \neq 0.

Sur [0,1][0,1] : x˙(t)ax(t)=b(1+t1)=bt\dot{x}(t) - ax(t) = b(1+t-1) = bt. EDO linéaire du premier ordre : x(t)=eat(1+b0tseasds)x(t) = e^{at}\left(1 + b\int_0^t s e^{-as}ds\right). On intègre par parties : 0tseasds=taeat+1a2(1eat)\int_0^t se^{-as}ds = \frac{-t}{a}e^{-at} + \frac{1}{a^2}(1-e^{-at}).

Sur [1,2][1,2] : x˙(t)ax(t)=bx(t1)\dot{x}(t) - ax(t) = bx(t-1)x(t1)x(t-1) est la solution trouvée sur [0,1][0,1]. On résout comme EDO linéaire d'ordre 1 avec second membre connu.

x(t)=eat(1+b0tθ(s1)easds) sur [0,1].\boxed{x(t) = e^{at}\left(1 + b\int_0^t \theta(s-1)e^{-as}ds\right) \text{ sur } [0,1].}

التمرين 3

Exercice 3 — Fonction de Lyapunov et stabilité asymptotique

#lyapunov-function#asymptotic-stability#lotka-volterra#dynamical-systems

Pour le système :

{x˙=x(yb)y˙=y(xa),ab0\begin{cases} \dot{x} = x(y - b) \\ \dot{y} = y(x - a) \end{cases}, \quad a \cdot b \neq 0

trouver une fonction de Lyapunov définie au voisinage de (0,0)(0,0) qui montre que les solutions issues du domaine (x/a)2+(y/b)2<1(x/a)^2 + (y/b)^2 < 1 convergent vers l'origine quand t+t \to +\infty (stabilité asymptotique locale).

Rappel : une fonction V(x,y)V(x,y) définie au voisinage de (0,0)(0,0) est dite de Lyapunov si VV est définie positive et V˙\dot{V} est semi-définie négative.

الحل

Solution.

Considérons V(x,y)=x2a2+y2b2V(x,y) = \frac{x^2}{a^2} + \frac{y^2}{b^2} (définie positive au voisinage de (0,0)(0,0)).

V˙=2xa2x˙+2yb2y˙=2xa2x(yb)+2yb2y(xa)\dot{V} = \frac{2x}{a^2}\dot{x} + \frac{2y}{b^2}\dot{y} = \frac{2x}{a^2}x(y-b) + \frac{2y}{b^2}y(x-a)

=2x2(yb)a2+2y2(xa)b2= \frac{2x^2(y-b)}{a^2} + \frac{2y^2(x-a)}{b^2}

Au voisinage de (0,0)(0,0) avec (x/a)2+(y/b)2<1(x/a)^2+(y/b)^2 < 1, pour a,b>0a,b > 0 : x<ax < a et y<by < b, donc yb<0y-b < 0 et xa<0x-a < 0, ce qui donne V˙<0\dot{V} < 0 (sauf en (0,0)(0,0)). Par le théorème de Lyapunov, l'origine est asymptotiquement stable, et les solutions issues du domaine {(x/a)2+(y/b)2<1}\{(x/a)^2+(y/b)^2 < 1\} convergent vers (0,0)(0,0).

V(x,y)=x2a2+y2b2,V˙<0 sur le domaine donneˊ.\boxed{V(x,y) = \frac{x^2}{a^2} + \frac{y^2}{b^2}, \quad \dot{V} < 0 \text{ sur le domaine donné.}}