Centrale Mathématiques 1 TSI 2008

Un système différentiel autonome de dimension 2 est un système différentiel qui se présente sous la forme

où et sont des fonctions de classe sur un ouvert de . Les fonctions inconnues et sont des fonctions de la variable définies sur un intervalle de à valeurs dans . On utilise aussi la notation matricielle

où le vecteur représente le couple et où est de classe .
Un point critique de ce système est un point de tel que

Étant donné un intervalle borné ou non, le couple de fonctions est une solution de ( ) sur si les applications et sont de classe et vérifient pour tout . Une valeur initiale est un triplet ( ) de nombres réels. Une solution ( ) vérifie cette valeur initiale si , .
Si est une solution de sur , l'ensemble des points , est la trajectoire de cette solution. Cette trajectoire est maximale si pour toute trajectoire ayant un point commun avec , on a .
On s'intéresse au comportement des trajectoires de certains systèmes autonomes. On rappelle le résultat suivant, qui est une conséquence du théorème de CauchyLipschitz :
Tout point de appartient à une trajectoire maximale unique.
Remarque sur la rédaction Il sera demandé de tracer des courbes. Il est entendu qu'il faut le faire dans un repère cartésien, mais que la précision souhaitée n'est pas
celle d'un dessin industriel. Il suffira de tracer à main levée dans un carré de moins de 12 cm de côté, en s'attachant à bien mettre en valeur les éléments mathématiques de la figure, selon le cas : points particuliers, tangentes etc. Les courbes tracées point par point ne seront pas prises en compte.
Notation Dans le plan euclidien orienté, on note la matrice de la rotation de centre 0 , d'angle .

Les résultats obtenus dans cette partie seront utilisés par la suite.

I.A.1) Soit .
a) Donner la solution générale de l'équation différentielle sur .
b) Donner la solution particulière vérifiant la condition initiale .
I.A.2) Soit deux réels. Donner la solution générale réelle de l'équation différentielle

dans les deux cas : .
I.B - On considère l'équation différentielle autonome .
I.B.1) Déterminer les solutions constantes.
I.B.2) Trouver une primitive de .
I.B.3)
a) Montrer que les solutions qui ne sont pas constantes peuvent se mettre sous la forme

où est une constante non nulle.
b) Tracer, dans un repère orthonormal, la courbe représentative de la solution vérifiant la condition initiale .

On va étudier le système linéaire :

où est une matrice , ayant deux valeurs propres distinctes, réelles ou complexes, et non nulles. On admet que la solution générale du système ( ) est de la forme

où et sont deux solutions telles que le déterminant, dans une base orthormale directe, de la matrice ne s'annule jamais.
II.A - On considère le cas où est diagonale, donc de la forme . II.A.1) Utiliser I.A.1) pour obtenir la solution générale de , sous la forme où est une matrice et est un vecteur dont les composantes sont des constantes arbitraires.
II.A.2) Déterminer les trajectoires maximales passant par , .
II.A.3) Soit une valeur initiale avec . Démontrer que la solution a une trajectoire dont l'équation cartésienne est de la forme

On déterminera en fonction de , et en fonction de .
II.A.4) Décrire les trajectoires maximales passant par le point dans chacun des cas suivants :
(i) ;
(ii) ;
(iii) ,
et les tracer dans un même repère.
II.B - Maintenant n'est plus nécessairement diagonale, mais elle a toujours des valeurs propres réelles . Soit la matrice .
II.B.1) Montrer qu'il existe une matrice inversible qui s'écrit en colonne : telle que si , le nouveau vecteur vérifie .
II.B.2) En déduire que la solution générale de peut se mettre sous la forme où est un vecteur de constantes arbitraires et une matrice dont le déterminant ne s'annule pas.
II.B.3) Donner la condition sur les pour que toutes les solutions tendent vers lorsque .
II.B.4) On étudie le système différentiel : .
a) Calculer et (notations reprises du II.B.1)). Vérifier que pour on peut prendre une matrice de rotation où est un angle à déterminer. Que vaut ?
b) Trouver la solution générale.
c) On donne la valeur initiale . Calculer la solution du système ci-dessus, et la solution du système .
d) Tracer la trajectoire maximale de passant par le point (on pourra d'abord tracer celle correspondant à ).
II.B.5) Mêmes questions avec le système

II.C - On considère le cas où , où et sont deux réels, .
II.C.1) Calculer les valeurs propres de .
II.C.2) Montrer que si est solution, alors et sont toutes deux solutions de l'équation

En déduire que la solution générale de est de la forme , où est un vecteur de constantes arbitraires.
II.C.3) Donner la condition sur pour que toutes les solutions tendent vers 0 lorsque .
II.C.4) On étudie le système

a) Trouver la solution vérifiant la valeur initiale . Soit la trajectoire correspondante.
b) Déterminer l'équation polaire de , sous la forme .
c) Tracer .

On va étudier le cas non linéaire

où et est définie par .

III.A.1) Calculer la matrice jacobienne de et donner le développement de Taylor de au voisinage d'un point quelconque ( ).
III.A.2) Trouver les points critiques de ( ). Montrer que chaque point critique constitue une trajectoire maximale.
III.A.3) Décrire la trajectoire maximale passant par (0.5, 0.5) et celle passant par (on pourra utiliser la question I.B).
III.A.4) Trouver l'ensemble des points ( ) où la tangente à la trajectoire passant par ( ) est
a) horizontale ;
b) verticale ;
c) de coefficient directeur 1 ;
d) de coefficient directeur -1 .
III.B - On admet que pour un tel système différentiel, l'allure des trajectoires au voisinage de chaque point critique est la même que pour le système dit linéarisé au voisinage de 0 .
La figure ci-dessous représente quelques trajectoires de .

En utilisant la partie II, justifier l'allure de ces courbes au voisinage des points fixes.

III.C - La figure précédente a été tracée en appliquant la méthode d'Euler. C'est un algorithme qui permet l'approximation d'une solution au moyen d'une suite , construite par récurrence de la manière suivante :

Une interpolation des points permet de tracer une courbe qui constitue l'approximation cherchée de la trajectoire.
III.C.1) Soit la fonction définie par

Montrer que la méthode d'Euler consiste à construire la suite ( ) définie par

On appelle aussi ( ) une orbite de .
III.C.2) Déterminer les points fixes de , c'est-à-dire les solutions de l'équation

III.D - On veut montrer que si est pris dans un voisinage de , toutes les orbites de convergent vers . On prend dans et on pose .
III.D.1) Commençons par le vérifier pour toutes les orbites partant d'un point de la première bissectrice. Soit .
a) Montrer que si , alors .
b) Soit tel que . Soit la suite définie par récurrence en posant . Montrer que .
c) En déduire que si et , alors et .
III.D.2) Pour tout on appelle norme de le réel .
a) Vérifier qu'une suite de vecteurs ( ) converge vers une limite si, et seulement si, la norme tend vers 0 .
b) Pour chaque vecteur on définit

Montrer que si , alors .
c) Partant de , on définit la suite ( ) comme l'orbite de , c'est-à-dire . Montrer que .
d) En déduire que si , alors .