Introduction Syst Lineaires C7 PDF
Introduction Syst Lineaires C7 PDF
Introduction Syst Lineaires C7 PDF
Cette solution obéit à la solution classique des équations différentielles, dans un contexte
multivariable.
1
INTRODUCTION AUX SYSTEMES LINEAIRES 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 3
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
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
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
𝑌𝐸𝑠 = 𝐴 + 𝐵𝑡
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
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
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
Avec :
𝑦 : Sortie ; 𝑢 : Entrée ; 𝜉 : Coefficient d’amortissement ; 𝜔𝑛 : Fréquence naturelle
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
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
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
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
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
y → 0 < ξ <1
y →ξ = 0
H y →ξ >1
y →ξ =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 14
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.
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
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
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
Si la matrice de transmission directe 𝐷 = 0 alors le modèle est dit strictement propre dans un
autre cas il est dit propre.
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
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
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
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
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.