Résolution d'un système X' = AX : Cas où la matrice A est diagonalisable, à valeurs propres toutes réelles
Définition :
Dire que \(A\) est diagonalisable signifie qu'il existe une matrice inversible \(P\) et une matrice diagonale \(D\) telle que l'on ait \(A = P D P^{-1}\)
Rappelons que les colonnes de \(P\) sont formées des composantes d'une base de vecteurs propres, les éléments diagonaux correspondants de \(D\) étant les valeurs propres associées.
L'idée est de ramener la résolution du système \(X' = AX\) à celle d'un système \(U' = D\)U, où \(D\) est diagonale.
Posons \(U(t) = P^{-1} X(t)\).
Puisque \(P\) est une matrice constante, on a \(U'(t) = P^{-1} X'(t)\).
On trouve alors \(U'(t) = P^{-1} X'(t) = P^{-1} A X(t) = P^{-1} PD P^{-1} X(t) = D P^{-1} X(t) = D U(t)\).
Or, puisque la matrice \(D\) est diagonale, on sait résoudre le système \(U'(t) = D U(t)\) (voir page précédente).
On revient aux solutions du système initial en utilisant la relation \(X(t) = P U(t)\).
Règle : Pour résumer
Si \(U(t)\) est la solution générale de \(U'(t) = D U(t)\), la solution générale de \(X' = AX\) est \(X(t) = P U(t)\).
Remarque :
Pour trouver la solution générale de \(X' = AX\), il faut trouver les matrices \(D\) et \(P\), mais il n'est pas nécessaire de calculer la matrice \(P^{-1}\).
Autre présentation de la solution générale
Il résulte facilement de ce qui précède que, si \((V_1, ... , V_n)\) est une base de vecteurs propres de \(A\), le vecteur \(V_i\) étant associé à la valeur propre \(\lambda_i\), les fonctions \(\textrm{exp}(\lambda_it)V_i(i=1,...,n)\) forment une base de l'espace des solutions.
La solution générale de \(X' = AX\) s'écrit donc
\(X(t)=\alpha_1\textrm{exp}(\lambda_1t)V_1+...+\alpha_n\textrm{exp}(\lambda_nt)V_n\)
où les \(a_i (i = 1 ... n)\) sont des constantes réelles.
Conditions initiales
Pour \(t = 0\), la solution prend la valeur initiale \(X_0 = a_1V_1+...+a_nV_n\). Autrement dit les constantes \((a_1, ... , a_n)\) sont les composantes de la valeur initiale \(X_0\) dans la base de vecteurs propres \(( V_1, ... , V_n )\).
En fait, la condition initiale \(X_0\) est en général donnée dans la base canonique, sous la forme
\((2)\displaystyle{\left\{\begin{array}{ccc}x_1(0)=c_1 \\ \vdots \\ x_n(0)=c_n\end{array}\right.}\)
Les constantes \(a_1,..., a_n\) peuvent se calculer par la formule de changement de base
\(\displaystyle{\left[\begin{array}{ccc}\alpha_1 \\ \vdots \\ \alpha_n \end{array}\right]} = P^{-1}\displaystyle{\left[\begin{array}{ccc} c_1 \\ \vdots \\ c_n\end{array}\right]}\)
Le calcul de \(P^{-1}\) peut donc être utile quand on veut imposer à la solution des conditions initiales données.
Dans la pratique, si on doit trouver une seule solution avec les conditions initiales (2) imposées, chercher les constantes \(a_1,...,a_n\) telles que
\(\alpha_1V_1 + ... + \alpha_nV_n = \displaystyle{\left[\begin{array}{ccc}c_1 \\ \vdots \\ c_n\end{array}\right]}\)
revient à résoudre un système d'équations linéaires d'ordre \(n\) (voir exemple).
Exemple : Exemple de résolution d'un système X' = AX, où la matrice A est diagonalisable
Soit à résoudre le système
\(\displaystyle{\left\{\begin{array}{lll}x'=5x-y+9z \\ y'=3x+4y \\ z'=x+y+z\end{array}\right.}\)
La matrice \(A = \displaystyle{\left(\begin{array}{ccc}5 & -1 & 9 \\ 3 & 4 & 0 \\ 1 & 1 & 1\end{array}\right)}\) a pour valeurs propres \(\lambda_1 = 1\), \(\lambda_2 = 2\), \(\lambda_3 = 7\).
Des vecteurs propres correspondants sont \(V_1 = \displaystyle{\left[\begin{array}{ccc} 9 \\ -9 \\5\end{array}\right]},V_2 = \displaystyle{\left[\begin{array}{ccc} -2 \\ 3 \\ 1\end{array}\right]},V_3 = \displaystyle{\left[\begin{array}{ccc} 3 \\ 3 \\ 1\end{array}\right]}\)
La solution générale du système s'écrit donc
\(X(t) = \displaystyle{\left[\begin{array}{ccc}x(t) \\ y(t) \\ z(t)\end{array}\right]}=\alpha\textrm{e}^tV_1+b\textrm{e}^{2t}V_2+c\textrm{e}^{7t}V_3=\alpha\textrm{e}^t\displaystyle{\left[\begin{array}{ccc} 9 \\ -9 \\5\end{array}\right]}+b\textrm{e}^{2t}\displaystyle{\left[\begin{array}{ccc} -2 \\ 3 \\ 1\end{array}\right]}+c\textrm{e}^{7t}\displaystyle{\left[\begin{array}{ccc} 3 \\ 3 \\ 1\end{array}\right]}\)
On l'écrit plutot sous la forme :
\(x(t) = 9a e^t - 2 b e^{2t} + 3 c e^{7t}\)
\(y(t) = - 9a e^t + 3 b e^{2t} + 3 c e^{7t}\)
\(z(t) = 5a e^t + b e^{2t} + c e^{7t}\)
Cherchons la solution vérifiant \(x(0) = 1\), \(y(0) = 2\), \(z(0) = 0\).
En faisant \(t = 0\) dans le système ci-dessus, on obtient le système (non différentiel)
\(\displaystyle{\left\{\begin{array}{lll}9a-2b+3c=1 \\ -9a+3b+3c=2 \\ 5a+b+c=0\end{array}\right.}\)
En résolvant ce système linéaire, on trouve \(a = -1/12\), \(b = -1/10\), \(c = 31/60\).
La solution cherchée est donc
\(x(t) = - 3/4 e^t + 1/5 e^{2t} + 31/20 e^{7t}\)
\(y(t) = 3/4 e^t - 3/10 B e^{2t} + 31/20 e^{7t}\)
\(z(t) = - 5/12 e^t - 1/10 e^{2t} + 31/60 e^{7t}\)