Introduction Syst Lineaires C7 PDF

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 22

INTRODUCTION AUX SYSTEMES LINEAIRES

Cours 7. Résolution des équations d’état


La solution des équations d’état permet de calculer la valeur de chaque état à tout moment dans
le temps, c´est donc la réponse du système.

Cette solution obéit à la solution classique des équations différentielles, dans un contexte
multivariable.

REPONSE DES SYSTEMES, pour EQUATIONS DIFFERENTIELLES MONOVARIABLES.

La réponse d’un système correspond à la solution de l’équation différentielle du modèle, qui


pour les systèmes linéaires invariants dans le temps peu se diviser en deux parties :
• La solution en régime transitoire, ou solution homogène (réponse libre) : qui représente la
transition d’un état initial a un état final du système, et ce calcule pour une entrée égale a
zéro.
• La solution en régime permanent, ou solution particulière (réponse forcée) : qui représente
la réponse du système pour un temps infini, et qui se calcule pour une entrée spécifique.

1
INTRODUCTION AUX SYSTEMES LINEAIRES 2

Réponses des systèmes de premier ordre.


Modèle mathématique de la forme:
τDy + y = u
Avec:
𝑦 : Sortie ; 𝑢 : entrée ; 𝜏 : constante de temps.

Solution en Régime Transitoire

On considère une entrée nulle:


τDy + y = 0
L’équation caractéristique est :
τD + 1 = 0
Ou sa seule racine est :
−1
D=
τ
La réponse libre est de la forme :
yT = Ce −t τ

Ou 𝐶 est un paramètre qui dépend des conditions initiales.

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 3

Solution en Régime Permanent


Celle si dépend du type d’entrée.

Pour une entrée en forme d’échelon, valeur constante à changement brusque appliqué a un
moment donné, pour tout 𝑡 > 0 l’entrée 𝑢 = 𝐻, ou 𝐻 est une constante, l’équation du système
prend la forme:
𝜏𝐷𝑦 + 𝑦 = 𝐻
Pour calculer la réponse forcée, on suppose une réponse de la forme
yEs = A

Ou 𝐴 est une constante.


On obtient la dérivé : DyEas = 0

On substitue la solution supposée et sa dérivé dans l’équation, et on obtient :


A=H
yE = H

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 4

Réponse du Système
La réponse complète du système, avec une entrée en forme d’échelon, est la somme des deux
réponses, vu que le système est linéaire et donc satisfait le principe de superposition :
y = yT + y E
y = Ce −t τ + H

Le coefficient 𝐶 se détermine avec les valeurs de conditions initiales, su on suppose que pour
𝑡 = 0 → 𝑦0 = 0 :
On a :
0 = Ce −0 τ + H
C = −H

La réponse du système est de la forme :


y = − He −t τ + H u
98.17%
Cette solution est représente graphiquement 63.21% y
H
sur la figure :

Ou pour 𝑡 = 𝜏 → 𝑢 = 63,21%𝐻 0 τ 4τ t

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 5

Pour une entrée en forme de rampe, valeur constante à changement brusque appliqué a un
moment donné, pour tout 𝑡 > 0 l’entrée 𝑢 = 𝐻𝑡, ou 𝐻 est une constante, l’équation du système
prend la forme:
𝜏𝐷𝑦 + 𝑦 = 𝐻𝑡
Pour calculer la réponse forcée, on suppose une réponse de la forme
𝑌𝐸𝑠 = 𝐴 + 𝐵𝑡

Ou 𝐴 et 𝐵 sont des constantes.


On obtient la dérivée : DyEs = B

On substitue la solution supposée et sa dérivée dans l’équation, et on obtient :


τB + A + Bt = Ht
On obtient deux équations, une pour les coefficients de t et une autre pour les indépendants :
B=H
τB + A = 0 ⇔ A = −τH
Donc la réponse forcée sera :
y E = −τH + Ht

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 6

Réponse du Système
La réponse complète du système, avec une entrée en forme de rampe, est la somme des deux
réponses, vu que le système est linéaire et donc satisfait le principe de superposition :
y = yT + y E
y = Ce −t τ − τH + Ht

Le coefficient 𝐶 se détermine avec les valeurs de conditions initiales, su on suppose que pour
𝑡 = 0 → 𝑦0 = 0 :
On a :
0 = Ce −0 τ − τH + H 0
C = τH
La réponse du système est de la forme :
y = τHe −t τ − τH + Ht
u = Ht

a
Cette solution est représente y = -τH+Ht
graphiquement sur la figure : b

0 t

-τH

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 7

Pour une entrée en forme sinusoïdale, valeur constante à changement brusque appliqué a un
moment donné, pour tout 𝑡 > 0 l’entrée 𝑢 = 𝐻 sin 𝜔𝑡 ou 𝐻 est une constante, l’équation du
système prend la forme:
τDy + y = H sin ωt

Pour calculer la réponse forcée, on suppose une réponse de la forme


𝑌𝐸𝑠 = 𝐴 sin 𝜔𝑡 + 𝐵 cos 𝜔𝑡

Ou 𝐴 et 𝐵 sont des constantes.


On obtient la dérivée : 𝐷𝑦𝐸𝑠 = 𝐴𝜔 cos 𝜔𝑡 + B𝜔sin 𝜔𝑡

On substitue la solution supposée et sa dérivée dans l’équation, et on obtient :


τ ( Aω cos ωt − Bω sin ωt ) + A sin ωt + B cos ωt = H sin ωt

On obtient deux équations, une pour les coefficients de sin et une autre pour les coefficients de
cos, et on obtient:
H − τωH
A= ; B =
1 + τ 2ω 2 1 + τ 2ω 2
Donc la réponse forcée sera :
H − τωH
yE = sin ω t + cos ωt
1 + τ 2ω 2 1 + τ 2ω 2

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 8

Réponse du Système
La réponse complète du système, avec une entrée en forme sinusoïdale, est la somme des deux
réponses, vu que le système est linéaire et donc satisfait le principe de superposition :
y = yT + y E
H − τωH
y = Ce −t τ + sin ω t + cos ωt
1 + τ 2ω 2 1 + τ 2ω 2

Le coefficient 𝐶 se détermine avec les valeurs de conditions initiales, su on suppose que pour
𝑡 = 0 → 𝑦0 = 0 :
On a :
H − τωH τωH
0 = Ce −0 τ + sin ω 0 + cos ω 0 C =
1 + τ 2ω 2 1 + τ 2ω 2 ; 1 + τ 2ω 2

La réponse du système est de la forme : u


τωH − t τ
sin (ωt + φ )
H
y= e + H
1 + τ 2ω 2 1+τ ω2 2 H y
1 + τ 2ω 2

Avec : φ = tan (τω )


−1
φ
1/ω
0 t

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 9

Réponses des systèmes de deuxième ordre.


Modèle mathématique de la forme:
D 2 y + 2ξω n Dy + ω n y = ω n u
2 2

Avec :
𝑦 : Sortie ; 𝑢 : Entrée ; 𝜉 : Coefficient d’amortissement ; 𝜔𝑛 : Fréquence naturelle

Solution en Régime Transitoire

On considère une entrée nulle:


D 2 y + 2ξωn Dy + ω n y = 0
2

L’équation caractéristique est :


D 2 + 2ξω n D + ω n = 0
2

Les deux racines sont :


D1, 2 = −ξω n ± ω n ξ 2 − 1

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 10

Il existe donc 4 possibles formes de la réponse libre, suivant la valeur du coefficient 𝜉 :

Si 𝜉 > 1 : deux racines réelles positives, solution de la forme :


 −ξω +ω ξ 2 −1  t  −ξω −ω ξ 2 −1  t
yT = C1e  
+ C2 e  
n n n n

Si 𝜉 = 1 : deux racines réelles égales, solution de la forme :


yT = C1e −ξωnt + C2te −ξωnt

Si 0 < 𝜉 < 1 : deux racines complexes conjuguées, solution de la forme :


( ( )
yT = e −ξωnt C1 sin ω n 1 − ξ 2 t − C 2 cos ω n 1 − ξ 2 t ( ))
Si 𝜉 = 0 : deux racines imaginaires pures (sans partie réelle), solution de la forme :
yT = C1 sin ω n t − C 2 cos ω n t

Avec 𝐶1 et 𝐶2 des coefficients dépendants des conditions initiales.

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 11

Solution en Régime Permanent

Pour une entrée en forme d’échelon, valeur constante à changement brusque appliqué a un
moment donné, pour tout 𝑡 > 0 l’entrée 𝑢 = 𝐻, ou 𝐻 est une constante, l’équation du système
prend la forme:
D 2 y + 2ξωn Dy + ω n y = ω n H
2 2

Pour calculer la réponse forcée, on suppose une réponse de la forme


yEs = A

Ou 𝐴 est une constante.


On obtient la première et la deuxième dérivé : DyEs = 0 ; D yEs = 0
2

On substitue la solution supposée et es dérivées dans l’équation, et on obtient :


ωn 2 A = ωn 2 H
yE = H

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 12

Réponse du Système
La réponse complète du système, avec une entrée en forme d’échelon, est la somme des deux
réponses, vu que le système est linéaire et donc satisfait le principe de superposition :
y = yT + y E
Avec 𝑦𝑇 qui peut prendre les quatre formes possibles, et les coefficients calculées avec les
conditions initiales égales a zéro t = 0 ⇒ y0 = 0; Dy0 = 0 :

ξ >1 
yT =  − H −
H ξω n − ω(n ξ 2
−1 )e  −ξω + ω ξ 2 −1  t
 n n

+
(
H ξωn − ωn ξ 2 − 1 )e  −ξω −ω ξ 2 −1  t
 n n

+H
 2ωn ξ 2 − 1  2ωn ξ − 1 2
 

ξ =1 yT = − He −ξωnt − ξωn Hte −ξωnt + H

0 < ξ <1
yT = e −ξω n t

 n( ω − ξ
) (
 H sin ω 1 − ξ 2 t + ξωn H cos ω 1 − ξ 2
2 n
) ( ) 
t+ H

 n 1 

ξ =0 yT = − H cos ω n t + H

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 13

La représentation de cette réponse est présentée sur la figure:

y → 0 < ξ <1
y →ξ = 0

H y →ξ >1

y →ξ =1

0 t

ξ > 1 : Réponse sur-amortie.


ξ = 1 : Réponse avec amortissement critique.
0 < ξ < 1 : Réponse sub-amortie.
ξ = 0 : Réponse sans amortissement.

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 14

SOLUTION DES EQUATIONS EN REPRESENTATION D’ETAT.

Nous nous intéressons plus particulièrement à la solution de l’équation différentielle d’état


(l’équation dynamique) :
𝑥̇ (𝑡) = 𝐴(𝑡)𝑥 (𝑡) + 𝐵(𝑡)𝑢(𝑡)

Car une fois obtenu sa solution, le calcul de la sortie du système revient simplement a une
opération algébrique.

La réponse du système de se type sera composé de deux partie : la réponse libre et la réponse
forcée.

Réponse Libre : Dans le cas général où le modèle est temps-variant, il est nécessaire de définir
la matrice de transition d’état qui intervient dans la solution de l’équation d’état.

Matrice de transition : La solution de l’équation d’état homogène (non commandée):


𝑥̇ (𝑡) = 𝐴(𝑡)𝑥 (𝑡) 𝑥 (𝑡0 ) = 𝑥0

où 𝐴(𝑡) est continue par rapport à 𝑡, est donnée par la réponse libre d’un système, qui est la
réponse du système à ses seules conditions initiales :
𝑥(𝑡) = Φ(𝑡, 𝑡0 )𝑥0 ∀𝑡 ≥ 𝑡0

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 15

Où Φ(𝑡, 𝑡0) est la matrice de transition, solution de l’équation différentielle :


Φ̇ (𝑡, 𝑡0 ) = 𝐴(𝑡)Φ(𝑡, 𝑡0 ) ∀ 𝑡 ≥ 𝑡0

La matrice de transition possède les propriétés suivantes.


- Φ(𝑡2 , 𝑡2 )Φ(𝑡1 , 𝑡0 ) = Φ(𝑡2 , 𝑡0 )
- Φ(𝑡, 𝑡0 ) est inversible pour tout 𝑡 et tout 𝑡0
- Φ−1 (𝑡, 𝑡0 ) = Φ(𝑡0 , 𝑡) ∀ 𝑡, 𝑡0

Dans le cas général d’un système temps-variant, la matrice de transition est rarement accessible
littéralement, elle doit être souvent obtenue par intégration numérique. Pour les systèmes LTI, la
matrice de transition peut être calculée de manière explicite.

La réponse forcée d’un système est la réponse du système au seul signal d’entrée et pour des
conditions initiales nulles. Dans le cas d’un système représenté par son modèle d’état, elle est
donnée par :
𝑡
𝑥(𝑡) = � Φ(𝑡, 𝜏)𝐵 (𝜏)𝑢(𝜏)𝑑𝜏
𝑡0

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 16

A partir de la donnée de la matrice de transition, il est possible de déterminer la solution de


l’équation dynamique d’état qui correspond à la réponse complète du système par application du
principe de superposition. C'est-à-dire la somme de la réponse libre plus la réponse forcée.

La solution de l’équation dynamique d’état :


𝑥̇ (𝑡) = 𝐴(𝑡)𝑥(𝑡) + 𝐵(𝑡)𝑢(𝑡) 𝑥(𝑡0 ) = 𝑥0
est donnée par :
𝑡
𝑥(𝑡) = Φ(𝑡, 𝑡0 )𝑥0 + � Φ(𝑡, 𝜏)𝐵 (𝜏)𝑢(𝜏)𝑑𝜏
𝑡0
et la sortie s’écrit alors :
𝑡
𝑦(𝑡) = 𝐶 (𝑡)Φ(𝑡, 𝑡0 )𝑥0 + 𝐶(𝑡) � Φ(𝑡, 𝜏 )𝐵(𝜏 )𝑢(𝜏 )𝑑𝜏 + 𝐷(𝑡)𝑢(𝑡)
𝑡0

CAS PARTICULIER DES SYSTEMES LTI. Solution dans le domaine temporel.


Soit le modèle LTI :
𝑥̇ (𝑡) = 𝐴𝑥(𝑡) + 𝐵𝑢(𝑡)
𝑦(𝑡) = 𝐶𝑥(𝑡) + 𝐷𝑢(𝑡)

complètement caractérisé par le triplet (𝐴, 𝐵, 𝐶, 𝐷), où 𝑥 ∈ ℝ𝑛 est le vecteur d’état, 𝑢 ∈ ℝ𝑚 est
le vecteur de commande et 𝑦 ∈ ℝ𝑟 est le vecteur de sortie.

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 17

On définie l’ordre (𝑛) du modèle comme la dimension du vecteur d’état.

Si la matrice de transmission directe 𝐷 = 0 alors le modèle est dit strictement propre dans un
autre cas il est dit propre.

Le modèle LTI homogène (non commandé) de condition initiale 𝑥(𝑡0 ) = 𝑥0 ,


𝑥̇ (𝑡) = 𝐴𝑥(𝑡)
𝑦(𝑡) = 𝐶𝑥(𝑡)
constitué par une matrice d’équations de premier ordre, a pour matrice de transition :
Φ(𝑡, 𝑡0 ) = 𝑒 𝐴(𝑡−𝑡0 )
où il est possible d’obtenir la solution par le développement en série de Taylor :

2 2 𝑛 𝑛
𝐴𝑡
𝐴 𝑡 𝐴 𝑡 𝐴𝑖 𝑡 𝑖
𝑒 = 1 + 𝐴𝑡 + + ⋯+ +⋯=�
2! 𝑛! 𝑖!
𝑖=0
Sa réponse libre est donc donnée par :
𝑦(𝑡) = 𝐶𝑒 𝐴(𝑡−𝑡0 ) 𝑥0

Attention : A différence des scalaires (𝑎. 𝑏 = 𝑏. 𝑎) les opérations matricielles ne sont pas
commutatives (𝐴. 𝐵 ≠ 𝐵. 𝐴), donc :
𝑒 (𝑎+𝑏)𝑡 ≠ 𝑒 𝑎𝑡 +𝑒 𝑏𝑡

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 18

Theoème de Cayley-Hamilton.
Permet de calculer l’exponentiel matriciel à partir d’un nombre fini d’opérations. Ce théorème
exprime que toute matrice carrée 𝐴 est solution de son équation caractéristique. On note 𝑄𝐴(𝑠)
le polynôme caractéristique de 𝐴 ∶ 𝑄𝐴(𝑠) = 𝑑𝑒𝑡(𝑠𝐼 − 𝐴). Si 𝐴 est une matrice carrée de taille
𝑛, 𝑄𝐴(𝑠) est un polynôme de degré 𝑛.
𝑄𝐴(𝑠) = 𝑑𝑒𝑡(𝑠𝐼 − 𝐴) = 𝑠 𝑛 + 𝑎𝑛−1 𝑠 𝑛−1 + . . . + 𝑎1 𝑠 + 𝑎0

Le théorème de Cayley-Hamilton assure que 𝐴 vérifie :


𝐴𝑛 + 𝑎𝑛−1 𝐴𝑛−1 + . . . + 𝑎1 𝐴 + 𝑎0 𝐼 = 0

En d’autres termes, 𝐴𝑛 s’exprime comme une combinaison linéaire de puissances inférieures


de 𝐴 :
𝐴𝑛 = −𝑎𝑛−1 𝐴𝑛−1 − . . . − 𝑎1 𝐴 – 𝑎0 𝐼

Par conséquent, il est possible d’exprimer



𝐴𝑖 𝑡 𝑖
� = 𝑒 𝐴𝑡 = Φ(𝑡)
𝑖!
𝑖=0
en ne faisant intervenir que des puissances de A inférieures strictement à 𝑛, c’est-à-dire, qu’il
existe un jeu de coefficients (𝛼0, 𝛼1. . . 𝛼𝑛 − 1) tels que :
Φ(𝑡) = 𝑒 𝐴𝑡 = 𝛼𝑛−1 (𝑡)𝐴𝑛−1 + . . . + 𝛼1 (𝑡)𝐴 + 𝛼0 (𝑡)𝐼

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 19

La réponse forcée d’un système LTI est donnée par :


𝑡
𝑥(𝑡) = � 𝑒 𝐴(𝑡−𝜏) 𝐵𝑢(𝜏)𝑑𝜏
𝑡0
Le modèle LTI,
𝑥̇ (𝑡) = 𝐴𝑥(𝑡) + 𝐵𝑢(𝑡)
𝑦(𝑡) = 𝐶𝑥(𝑡) + 𝐷𝑢(𝑡)

A pour solution à l’Équation dynamique d’état :


𝑡
𝐴(𝑡−𝑡0 )
𝑥(𝑡) = 𝑒 𝑥0 + � 𝑒 𝐴(𝑡−𝜏) 𝐵𝑢 (𝜏)𝑑𝜏
𝑡0
Et la réponse du système est :
𝑡
𝐴(𝑡−𝑡0 )
𝑦(𝑡) = 𝐶𝑒 𝑥0 + 𝐶 � 𝑒 𝐴(𝑡−𝜏) 𝐵𝑢(𝜏 )𝑑𝜏 + 𝐷𝑢(𝑡)
𝑡0

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 20

CAS PARTICULIER DES SYSTEMES LTI. Solution dans le domaine de Laplace.

Il est aussi possible de déterminer la matrice de transition en passant par la transformée de


Laplace du système.

Soit le modèle LTI :


𝑥̇ (𝑡) = 𝐴𝑥(𝑡) + 𝐵𝑢(𝑡)
𝑦(𝑡) = 𝐶𝑥(𝑡) + 𝐷𝑢(𝑡)

Si on considère les conditions initiales égales a zéro, ce qui est le cas pour les transformées
de Laplace, la solution de l’équation d’état est de la forme :
𝑡
𝑥(𝑡) = Φ(𝑡)𝑥(0) + � Φ(𝑡 − 𝜏)𝐵𝑢 (𝜏)𝑑𝜏
0
Il est possible de démontrer que :
Φ(𝑠) = ℒ{Φ(𝑡)} = [𝑠𝐼 − 𝐴]−1

Il est donc possible d’obtenir la solution de l’équation d’état avec la transformée inverse de
Laplace de Φ(𝑠) :
Φ(𝑡) = ℒ −1 {Φ(𝑠)}

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 21

Exemple : On veut obtenir la matrice de transition du système linéaire représenté par la


matrice :
0 −2
𝐴=� �
1 −3
Il est donc possible d’obtenir :
−1
−1 1 0 0 −2 𝑠 2 −1
Φ(𝑠) = [𝑠𝐼 − 𝐴] = �𝑠 � �−� �� = � �
0 1 1 −3 −1 𝑠 + 3
1 𝑠 + 3 −2 1 𝑠 + 3 −2
Φ(𝑠) = � �= 2 � �
Δ𝑠 1 𝑠 (s + 3s + 2) 1 𝑠
(𝑠 + 3) 2
⎡ − ⎤
( )(
𝑠+1 𝑠+2 ) ( )(
𝑠+1 𝑠+2 ⎥ )
Φ(𝑠) = ⎢
⎢ 1 𝑠 ⎥
⎣(𝑠 + 1)(𝑠 + 2) (𝑠 + 1)(𝑠 + 2) ⎦

Pour obtenir la transformée inverse il faut décomposer les éléments en fractions, par
exemple pour le premier élément de la matrice (𝐴11 ) :
(𝑠 + 3) 𝐾1 𝐾2
= +
(𝑠 + 1)(𝑠 + 2) (𝑠 + 1) (𝑠 + 2)
(𝑠−𝑠𝑖 )𝑁(𝑠)
Avec : 𝐾𝑖 = ( )
, donc :
𝐷 𝑠

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.
INTRODUCTION AUX SYSTEMES LINEAIRES 22

(𝑠 + 1)(𝑠 + 3)
𝐾1 = � =2
(𝑠 + 1)(𝑠 + 2) 𝑠=−1
(𝑠 + 2)(𝑠 + 3)
𝐾21 = � = −1
(𝑠 + 1)(𝑠 + 2) 𝑠=−2

La matrice s’écrit :
2 −1 −2 2
⎡ + + ⎤
(𝑠 + 1 ) (𝑠 + 2) (𝑠 + 1) (𝑠 + 2)
Φ(𝑠) = ⎢ ⎥
⎢ 1 + −1 −1
+
2 ⎥
⎣(𝑠 + 1) (𝑠 + 2) (𝑠 + 1) (𝑠 + 2)⎦
La matrice de transition du système est :
−𝑡 −2𝑡 −𝑡 −2𝑡
( ) −1 ( )
Φ t = ℒ {Φ 𝑠 } = � −𝑡 2𝑒 − 𝑒 −2𝑒 + 2𝑒 �
𝑒 − 𝑒 −2𝑡 −𝑒 −𝑡 + 2𝑒 −2𝑡

Pour des conditions initiales connues, par exemple 𝑥1 (0) = 𝑥2 (0) = 1, il est possible de
déterminer la valeur des états en tout temps :
𝑥1 (𝑡) 1 𝑒 −2𝑡
� � = Φ(t) � � = � −2𝑡 �
𝑥2 (𝑡) 1 𝑒

Dulhoste Jean-François.- Universidad de Los Andes. Esc. Ing. Mecánica. Dpto. Ciencias Térmicas - Grenoble INP. ENSE3.Gipsa-lab.

Vous aimerez peut-être aussi