Système différentiel

Un système différentiel est un ensemble d'équations différentielles couplées, c'est-à-dire d'équations différentielles qui ne peuvent pas être résolues séparément.

Il s'agit en général d'équations différentielles ordinaires, mais un ensemble d'équations aux dérivées partielles peut aussi être qualifié de système différentiel.

Exemples

Équations différentielles couplées

Système de Lorenz :

\left\{{\begin{aligned}{\frac {\mathrm {d} x}{\mathrm {d} t}}&=\sigma \,[(y(t)-x(t)]\\{\frac {\mathrm {d} y}{\mathrm {d} t}}&=\rho \,x(t)-y(t)-x(t)\,z(t)\\{\frac {\mathrm {d} z}{\mathrm {d} t}}&=x(t)\,y(t)-\beta \,z(t)\end{aligned}}\right.

Ce système à seulement trois degrés de liberté est une simplification des équations de Navier-Stokes (voir ci-dessous), applicable à la convection de Rayleigh-Bénard pour des nombres de Rayleigh supérieurs à la valeur critique ( $\rho =\mathrm {Ra/Ra_{c}}$ ). C'est un des systèmes différentiels les plus simples conduisant à un comportement chaotique (ainsi qu'à des trajectoires périodiques).

Équations aux dérivées partielles couplées

Équations de Navier-Stokes :

\left\{{\begin{aligned}{\frac {\partial }{\partial t}}(\rho \,\mathbf {u} )+\nabla \cdot (\rho \,\mathbf {u} \otimes \mathbf {u} )&=-\nabla \cdot p\,\mathbf {I} +\nabla \cdot {\boldsymbol {\tau }}+\rho \,\mathbf {g} \\\rho \left({\frac {\partial \mathbf {u} }{\partial t}}+\mathbf {u} \cdot \nabla \mathbf {u} \right)&=-\nabla {\bar {p}}+\mu \,\nabla ^{2}\mathbf {u} +{\tfrac {1}{3}}mu\,\nabla (\nabla \cdot \mathbf {u} )+\rho \mathbf {g} \end{aligned}}\right.

où les grandeurs $p$ et ${\boldsymbol {\tau }}$ sont reliées à $\mathbf {u}$ par des relations non différentielles.

Les équations de Navier-Stokes décrivent le mouvement des fluides newtoniens, et constituent le cœur de la dynamique des fluides.

Système différentiel vs équation différentielle unique

La résolution d'un système d'équations différentielles peut être ramenée à celle d'une unique équation différentielle, d'ordre supérieur. Dans le système de Lorenz par exemple (cf. ci-dessus), on peut utiliser la 1^re équation pour exprimer $y(t)$ en fonction de $x(t)$ et ${\frac {\mathrm {d} x}{\mathrm {d} t}}$ , et reporter le résultat dans les deux autres équations. On peut ensuite extraire $z(t)$ de la 2^e équation pour l'exprimer en fonction de $x(t)$ , ${\frac {\mathrm {d} x}{\mathrm {d} t}}$ et ${\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}$ , et reporter le résultat dans la 3^e et dernière équation. On obtient ainsi une unique relation entre $x(t)$ , ${\frac {\mathrm {d} x}{\mathrm {d} t}}$ , ${\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}$ et ${\frac {\mathrm {d} ^{3}x}{\mathrm {d} t^{3}}}$ , c'est-à-dire une équation différentielle d'ordre 3.

Réciproquement, on peut transformer une équation différentielle d'ordre n en un système différentiel d'ordre 1 et de dimension n (c'est-à-dire un ensemble de n équations différentielles couplées, chacune d'ordre 1).

Cas particuliers

Systèmes différentiels linéaires

Un système différentiel linéaire est constitué d'équations différentielles linéaires (la linéarité porte sur les fonctions inconnues et sur leurs dérivées).

Un système différentiel d'ordre 1 et de dimension $n$ est équivalent à une unique équation différentielle d'ordre $n$ , et réciproquement. Plus généralement, tout système (de dimension quelconque) d'équations différentielles linéaires (d'ordre quelconque) peut s'écrire comme un système différentiel linéaire d'ordre 1, qu'on peut mettre sous la forme canonique :

{\frac {\mathrm {d} \mathbf {x} }{\mathrm {d} t}}=\mathbf {L} (t)\cdot \mathbf {x} (t)

où $\mathbf {x} (t)$ est un vecteur colonne rassemblant les $n$ fonctions inconnues $x_{1}(t),x_{2}(t),\ldots ,x_{n}(t)$ et $\mathbf {L} (t)$ une matrice carrée dont les éléments sont des fonctions connues de la variable $t$ . Le système a une solution unique si on lui adjoint $n$ conditions supplémentaires, en général sous la forme de conditions initiales :

\mathbf {x} (t_{0})=\mathbf {x} _{0}

où $t_{0}$ est un certain instant (« initial ») et $\mathbf {x} _{0}$ une colonne de $n$ constantes.

Comme tout système d'ordre 1 muni de conditions initiales, le système ci-dessus a une solution unique. On peut l'expliciter quand $\mathbf {L} (t)$ , la matrice des coefficients, commute avec sa dérivée $\mathrm {d} \mathbf {L} /\mathrm {d} t$ [alpha 1] :

\mathbf {x} (t)=\exp \left(\int _{t_{0}}^{t}\mathbf {L} (\tau )\,\mathrm {d} \tau \right)\cdot \mathbf {x} _{0}

où $\exp$ désigne l'opérateur d'exponentiation matricielle[1]^,[alpha 2]. Dans le cas général on ne sait exprimer la solution que sous la forme d'un développement en série[2]^,[alpha 3].

Systèmes différentiels autonomes

Quand on parle de systèmes autonomes la variable est en général le temps t. Un système différentiel est dit autonome si ses équations ne comportent aucune fonction de t autre que les fonctions inconnues et leurs dérivées (équations différentielles autonomes).

C'est notamment le cas du système de Lorenz ci-dessus, ainsi que des équations de Navier-Stokes si les paramètres ( $\rho$ , ${\vec {g}}$ , etc.) et les conditions aux limites ne dépendent pas explicitement du temps.

La particularité d'un système autonome est que par tout point de l'espace des solutions il passe une trajectoire et une seule. Dans le cas du système de Lorenz par exemple, par tout point A (de coordonnées $x_{\mathrm {A} },y_{\mathrm {A} },z_{\mathrm {A} }$ ) il passe une unique trajectoire $\{x(t),y(t),z(t)\}$ (au choix près de l'origine des temps).

Réciproque

Il est possible, dans une certaine mesure, de remonter des séries temporelles observées au système autonome qui les a engendrées, s'il est polynomial et suffisamment concis (jusqu'à 9 termes). La procédure, testée sur 28 cas théoriques impliquant jusqu'à 5 variables, est relativement robuste au bruit[3]^,[4].

Notes et références

Notes

C'est notamment le cas quand $\mathbf {L} (t)$ est une matrice constante ( $\mathbf {L} (t)=\mathbf {L} _{0}$ ), ou proportionnelle à une matrice constante ( $\mathbf {L} (t)=\mathbf {L} _{0}f(t)$ ), ou bien encore si elle est diagonale.
Pour vérifier que cette expression est bien une (la) solution du système différentiel et des conditions initiales ci-dessus, il suffit de calculer ${\frac {\mathrm {d} \mathbf {x} }{\mathrm {d} t}}$ en appliquant la définition de l'exponentielle d'une matrice : $\mathrm {e} ^{A}=I+A+{\frac {A^{2}}{2}}+\cdots +{\frac {A^{n}}{n!}}+\cdots$ .
On connaît une solution analytique fermée dans quelques rares cas où $\mathbf {L} (t)$ ne commute pas avec sa dérivée, notamment celui d'une matrice triangulaire[2].

Références

(en) Ariel Provost, Cécile Buisson et Olivier Merle, « From progressive to finite deformation and back », Journal of Geophysical Research: Solid Earth, vol. 109, n^o B2,‎ février 2004, p. 1-11, article n^o B02405 (DOI 10.1029/2001JB001734, lire en ligne, consulté le 10 juin 2018).
Daniel Pham, Techniques du Calcul Matriciel, Paris, Dunod, 1962, 387 p., p. 232-235.
« Est-il possible de retrouver les équations qui gouvernent la dynamique d’un système environnemental exclusivement à partir de séries de mesures ? », sur INSU, 1^er mars 2019 (consulté le 7 mars 2019).
(en) Sylvain Mangiarotti et Mireille Huc, « Can the original equations of a dynamical system be retrieved from observational time series? », Chaos (en), vol. 29,‎ 25 février 2019, p. 1-13, article n^o 023133 (DOI 10.1063/1.5081448).

Portail de l'analyse

Cet article est issu de Wikipedia. Le texte est sous licence Creative Commons - Attribution - Partage dans les Mêmes. Des conditions supplémentaires peuvent s'appliquer aux fichiers multimédias.

[1] C'est notamment le cas quand $\mathbf {L} (t)$ est une matrice constante ( $\mathbf {L} (t)=\mathbf {L} _{0}$ ), ou proportionnelle à une matrice constante ( $\mathbf {L} (t)=\mathbf {L} _{0}f(t)$ ), ou bien encore si elle est diagonale.

[3] Pour vérifier que cette expression est bien une (la) solution du système différentiel et des conditions initiales ci-dessus, il suffit de calculer ${\frac {\mathrm {d} \mathbf {x} }{\mathrm {d} t}}$ en appliquant la définition de l'exponentielle d'une matrice : $\mathrm {e} ^{A}=I+A+{\frac {A^{2}}{2}}+\cdots +{\frac {A^{n}}{n!}}+\cdots$ .

[5] On connaît une solution analytique fermée dans quelques rares cas où $\mathbf {L} (t)$ ne commute pas avec sa dérivée, notamment celui d'une matrice triangulaire[2].

[2] (en) Ariel Provost, Cécile Buisson et Olivier Merle, « From progressive to finite deformation and back », Journal of Geophysical Research: Solid Earth, vol. 109, n^o B2,‎ février 2004, p. 1-11, article n^o B02405 (DOI 10.1029/2001JB001734, lire en ligne, consulté le 10 juin 2018).

[Pham1962-4] Daniel Pham, Techniques du Calcul Matriciel, Paris, Dunod, 1962, 387 p., p. 232-235.

[6] « Est-il possible de retrouver les équations qui gouvernent la dynamique d’un système environnemental exclusivement à partir de séries de mesures ? », sur INSU, 1^er mars 2019 (consulté le 7 mars 2019).

[7] (en) Sylvain Mangiarotti et Mireille Huc, « Can the original equations of a dynamical system be retrieved from observational time series? », Chaos (en), vol. 29,‎ 25 février 2019, p. 1-13, article n^o 023133 (DOI 10.1063/1.5081448).