📚 الرئيسية

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

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

Concours d'accès au doctorat 2022/2023 — Épreuve 2 : Systèmes dynamiques / EDO — Sujet N°2 — 28 janvier 2023

التمرين 1

Champ de vecteurs à partir d'un flot ; réduction par changement de variable et point selle

#flot#champ de vecteurs#variétés invariantes#point selle

Partie A. On considère le flot défini sur R2\mathbb{R}^2 par φt(x0,y0)=(13y02e2t+(x013y021)et+1,  y0et).\varphi_t(x_0,y_0)=\Big(\tfrac13 y_0^2 e^{2t}+(x_0-\tfrac13 y_0^2-1)e^{-t}+1,\; y_0 e^{t}\Big). Déterminer le champ de vecteurs (système différentiel autonome) dont φt\varphi_t est le flot.

Partie B. On considère le système {x˙=4x2+4xyx+y2+1,y˙=8x28xy+4x2y2+y2.\begin{cases}\dot x=4x^2+4xy-x+y^2+1,\\ \dot y=-8x^2-8xy+4x-2y^2+y-2.\end{cases}

1. À l'aide du changement de variable u=2x+yu=2x+y, simplifier le système.

2. En déduire le flot du système.

3. Déterminer les points d'équilibre et les variétés invariantes.

4. Esquisser le portrait de phase.

الحل

Partie A — Champ de vecteurs

Posons X(t)=13y02e2t+(x013y021)et+1X(t)=\tfrac13 y_0^2 e^{2t}+(x_0-\tfrac13y_0^2-1)e^{-t}+1 et Y(t)=y0etY(t)=y_0e^t.

Composante yy : Y˙=y0et=Y\dot Y=y_0e^t=Y, donc y˙=y\dot y=y.

Composante xx : X˙=23y02e2t(x013y021)et\dot X=\tfrac23 y_0^2e^{2t}-(x_0-\tfrac13y_0^2-1)e^{-t}. Or y02e2t=Y2=y2y_0^2e^{2t}=Y^2=y^2 et (x013y021)et=X13y02e2t1=x13y21(x_0-\tfrac13y_0^2-1)e^{-t}=X-\tfrac13y_0^2e^{2t}-1=x-\tfrac13y^2-1. D'où X˙=23y2(x13y21)=y2x+1.\dot X=\tfrac23y^2-\big(x-\tfrac13y^2-1\big)=y^2-x+1. Le champ de vecteurs est donc   x˙=1x+y2,y˙=y.  \boxed{\;\dot x=1-x+y^2,\qquad \dot y=y.\;}

Partie B

1. Réduction. Comme 4x2+4xy+y2=(2x+y)2=u24x^2+4xy+y^2=(2x+y)^2=u^2, on a x˙=u2x+1\dot x=u^2-x+1 et y˙=2u2+4x+y2\dot y=-2u^2+4x+y-2. Alors u˙=2x˙+y˙=2(u2x+1)+(2u2+4x+y2)=2x+y=u.\dot u=2\dot x+\dot y=2(u^2-x+1)+(-2u^2+4x+y-2)=2x+y=u. Le système devient x˙=1x+u2, u˙=u\dot x=1-x+u^2,\ \dot u=u (identique à celui de la partie A, avec uu au lieu de yy).

2. Flot. u=u0etu=u_0e^{t} avec u0=2x0+y0u_0=2x_0+y_0, et comme en partie A x=13u02e2t+(x013u021)et+1,y=u2x.x=\tfrac13u_0^2e^{2t}+(x_0-\tfrac13u_0^2-1)e^{-t}+1,\qquad y=u-2x.

3. Équilibres et variétés invariantes. u˙=0u=0\dot u=0\Rightarrow u=0, puis x˙=1x=0x=1\dot x=1-x=0\Rightarrow x=1 : unique équilibre (x,u)=(1,0)(x,u)=(1,0), soit (x,y)=(1,2)(x,y)=(1,-2). La jacobienne en (1,0)(1,0) est (1001)\begin{pmatrix}-1&0\\0&1\end{pmatrix}, de valeurs propres 1-1 et 11 : c'est un point selle.

  • Variété stable : la droite u=0u=0, c'est-à-dire 2x+y=02x+y=0 (on y a x˙=1x\dot x=1-x, donc convergence vers x=1x=1).
  • Variété instable : la parabole x=1+13u2x=1+\tfrac13u^2, c'est-à-dire x=1+13(2x+y)2x=1+\tfrac13(2x+y)^2 (invariance : si x=1+u2/3x=1+u^2/3, alors x˙=23uu˙=23u2\dot x=\tfrac23u\dot u=\tfrac23u^2 et 1x+u2=23u21-x+u^2=\tfrac23u^2).

4. Portrait de phase. Point selle en (1,2)(1,-2) ; séparatrice stable 2x+y=02x+y=0, séparatrice instable la parabole x=1+13(2x+y)2x=1+\tfrac13(2x+y)^2. Les autres trajectoires suivent des branches d'allure hyperbolique entre ces deux séparatrices.

التمرين 2

Cycles limites en coordonnées polaires et ensembles limites α/ω

#coordonnées polaires#cycles limites#Poincaré-Bendixson#ensembles limites

Soit a1a\ge1. On considère le système {x˙=y+x(2a(x6+y6+3x4y2+3x2y4))(x2+y26),y˙=x+y(2a(x6+y6+3x4y2+3x2y4))(x2+y26).\begin{cases}\dot x=y+x\big(2-a(x^6+y^6+3x^4y^2+3x^2y^4)\big)(x^2+y^2-6),\\[3pt]\dot y=-x+y\big(2-a(x^6+y^6+3x^4y^2+3x^2y^4)\big)(x^2+y^2-6).\end{cases}

1. Écrire le système en coordonnées polaires.

2. Montrer que l'origine est le seul point d'équilibre.

3. Déterminer les cycles limites, préciser leur stabilité, et esquisser le portrait de phase.

4. Montrer, à l'aide du théorème de Poincaré–Bendixson, l'existence d'un cycle limite dans une couronne.

5. Déterminer les ensembles limites α\alpha et ω\omega de chaque trajectoire.

الحل

On remarque d'abord que x6+3x4y2+3x2y4+y6=(x2+y2)3x^6+3x^4y^2+3x^2y^4+y^6=(x^2+y^2)^3. En posant r2=x2+y2r^2=x^2+y^2, le facteur commun s'écrit F(r)=(2ar6)(r26)F(r)=(2-ar^6)(r^2-6).

1. Coordonnées polaires

rr˙=xx˙+yy˙=(x2+y2)F(r)=r2F(r),r2θ˙=xy˙yx˙=(x2+y2)=r2.r\dot r=x\dot x+y\dot y=(x^2+y^2)F(r)=r^2F(r),\qquad r^2\dot\theta=x\dot y-y\dot x=-(x^2+y^2)=-r^2. D'où r˙=r(2ar6)(r26),θ˙=1.\dot r=r\,(2-ar^6)(r^2-6),\qquad \dot\theta=-1.

2. Unique équilibre

Comme θ˙=10\dot\theta=-1\neq0, aucun point de rayon r>0r>0 ne peut être fixe. Le seul équilibre est l'origine r=0r=0.

3. Cycles limites et stabilité

Pour r>0r>0, r˙=0    2ar6=0\dot r=0\iff 2-ar^6=0 ou r26=0r^2-6=0, soit r1=(2a)1/6etr2=6.r_1=\Big(\tfrac2a\Big)^{1/6}\qquad\text{et}\qquad r_2=\sqrt6. Comme a1a\ge1, r121/6<6=r2r_1\le 2^{1/6}<\sqrt6=r_2. Le signe de r˙\dot r est :

  • 0<r<r10<r<r_1 : r˙<0\dot r<0 ;
  • r1<r<r2r_1<r<r_2 : r˙>0\dot r>0 ;
  • r>r2r>r_2 : r˙<0\dot r<0.

Donc le cercle r=r1r=r_1 est un cycle limite instable, le cercle r=r2=6r=r_2=\sqrt6 est un cycle limite stable, et l'origine est un foyer stable. Le portrait de phase : rotation (sens θ\theta décroissant), spirale vers l'origine pour r<r1r<r_1, spirale de r1r_1 vers r2r_2 dans la couronne intermédiaire, et spirale de l'extérieur vers r2r_2.

4. Théorème de Poincaré–Bendixson

Choisissons α,β\alpha,\beta avec r1<α<6<βr_1<\alpha<\sqrt6<\beta. Sur le cercle r=αr=\alpha on a r˙>0\dot r>0 (le flot entre dans la couronne) et sur r=βr=\beta on a r˙<0\dot r<0 (le flot entre aussi). La couronne K={αrβ}K=\{\alpha\le r\le\beta\} est donc positivement invariante et ne contient aucun équilibre. Par le théorème de Poincaré–Bendixson, KK contient une orbite périodique : c'est le cycle r=6r=\sqrt6.

5. Ensembles limites α\alpha et ω\omega

Selon le rayon initial r0r_0 :

  • r0=0r_0=0 : α=ω={0}\alpha=\omega=\{0\}.
  • 0<r0<r10<r_0<r_1 : ω={0}\omega=\{0\} et α={r=r1}\alpha=\{r=r_1\}.
  • r0=r1r_0=r_1 : α=ω={r=r1}\alpha=\omega=\{r=r_1\}.
  • r1<r0<r2r_1<r_0<r_2 : ω={r=6}\omega=\{r=\sqrt6\} et α={r=r1}\alpha=\{r=r_1\}.
  • r0=r2r_0=r_2 : α=ω={r=6}\alpha=\omega=\{r=\sqrt6\}.
  • r0>r2r_0>r_2 : ω={r=6}\omega=\{r=\sqrt6\} et α=\alpha=\varnothing (la trajectoire n'est pas bornée lorsque tt\to-\infty).