Fonction de Transfert et Systèmes Linéaires
Fonction de Transfert et Systèmes Linéaires
22/02/2017 21:06 1
Soit :
d‘où :
22/02/2017 21:06 2
1
22/02/2017
22/02/2017 21:06 3
22/02/2017 21:06 4
2
22/02/2017
22/02/2017 21:06 6
3
22/02/2017
Exemples
a) Système du second ordre excité par un échelon
unitaire
Considérons un système régi par l’équation différentielle suivante :
On en déduit donc :
22/02/2017 21:06 7
22/02/2017 21:06 8
4
22/02/2017
Identifions :
d‘où :
22/02/2017 21:06 9
Au final :
22/02/2017 21:06 10
5
22/02/2017
Remarque :
Les expressions temporelles fournies par les tables ne sont valables
que pour t > 0 ; ces fonctions sont nulles pour t < 0. La présence de
u(t) dans ces expressions suffit à nous le rappeler.
22/02/2017 21:06 11
22/02/2017 21:06 12
6
22/02/2017
22/02/2017 21:06 13
Par ailleurs,
d’où :
22/02/2017 21:06 14
7
22/02/2017
15
dy t 1 de t
y t Conditions initiales : e 0 e0 ; y 0 y0
dt RC dt
L’analyse de la sortie y(t) passe par la résolution de cette équation.
16
8
22/02/2017
dy t 1 de t
y t
e(t) R y(t) dt RC dt
e 0 e0 ; y t y0
1
sY s y0 Y s s E s e0
soit : RC
1
s Y s s E s y0 e0
RC
17
s 1
Y s E s y0 e 0
s 1 RC s 1 RC
G s G0 s
soit :
Y s G s E s G0 s y0 e0
18
9
22/02/2017
Variables d’état :
x1(t)
dx t dx1 t 1
e t L 2 Rx2 t x1 t x2 t
dt dt C
dx1 t 1 dx2 t
x2 t 1 R 1
x1 t x2 t e t
dt C dt L L L
y t R x2 t y t R x2 t
19
x t
En utilisant la matrice x des variables d’état : x t 1
x2 t
dx1 t 1 1
x2 t dx 0 0
1 e t
dt C t C
x t
dx t dt 1 R
2 1 R 1
x1 t x2 t e t L
dt L L L L L
y t R x2 t y t 0 R x t
On obtient l’écriture matricielle : CI : x1 0 ; x2 0
10
22/02/2017
dx t
t A x t B e t
e A t B e d
At
x t e x 0
dt
y t C x t Résolution de Sans Second Membre 0 Avec Second Membre
l’équation d’état
Ce qui conduit au diagramme suivant :
dx
+ dt x
e B C y
+
21
dx
t A x t B e t s X s A X s B E s
dt
y t C x t Transformée Y s C X s
de Laplace
s I A X s B E s X s s I A1 B E s
Y s C s I A B E s
1
Y s C X s
La fonction de transfert s’écrit :
Y s
C s I A B
1
G s
E s
22
11
22/02/2017
Exemple
dx1 t
4 x1 t 0.5 x2 t 0.5 u t
dt
dx2 t
Équations d’état 8 x1 t Équations de sortie
dt
dx3 t
x2 t
dt y t 0.25 x2 t 0.05 x3 t
Déterminer :
• Le modèle d’état du système.
• La fonction de transfert :
Y s
G s
U s
23
Relations de passage
2
4
12
22/02/2017
Structure
Un système de commande comprend les éléments suivants :
Système :
Organe principal à commander, il dispose d’un ensemble entrées-sorties
permettant de le connecter aux éléments externes.
Actionneur :
Organe de puissance permettant de générer les actions motrices en
fonction des signaux de contrôle.
Correcteur :
Organe d’intelligence de la structure de commande. Sa fonction consiste
à déterminer l’erreur entre le signal à commander et la référence, puis à
agir en conséquence pour minimiser cette erreur.
25
Perturbation : w(t)
26
13
22/02/2017
Perturbation : w(t)
Capteur
27
xret(t)
Chaîne
Paramètres : de retour
yréf : grandeur de référence (consigne). (Capteur)
14
22/02/2017
Diagramme fonctionnel
e Chaîne y
directe (G)
Y s
G s Fonction de transfert en
E s boucle ouverte.
29
Diagramme fonctionnel
+ Chaîne
e y
directe (G)
- Y s
G Fonction de transfert Directe.
x s
Chaîne de X s
H Fonction de transfert de
retour (H) Y s retour.
30
15
22/02/2017
Diagramme fonctionnel
Fonctions de transfert
on distingue les fonctions de transfert suivantes :
+ Chaîne
e y
X s Fonction de transfert en
-
directe (G) GH
s boucle ouverte.
x
Y s G Fonction de transfert en
Chaîne de E s 1 GH boucle fermée.
retour (H)
s 1 Fonction de transfert
E s 1 GH de l’erreur.
31
Exemple
G(s)
k + + 1 1
U(s) Y(s)
T² - - s s
2m
T
1
T²
32
16
22/02/2017
Plan
Réponse temporelle.
Réponse harmonique.
34
17
22/02/2017
Comportement stable
Soit un système à entrée e0 et à sortie y0.
À l’instant t = t0, on fait varier l’entrée de e.
La sortie y passe de y0 à y1 et se stabilise à l’instant t1.
y1
t
t0 t1
35
Comportement instable
Soit un système à entrée e0 et à sortie y0.
À l’instant t = t0, on fait varier l’entrée de e.
La sortie y augmente indéfiniment et ne se stabilise pas.
y0
e(t) y(t)
Système
e1
e
e0
t
t0
36
18
22/02/2017
y0
e(t) y(t)
Système
e1
e
e0
t
t0
37
38
19
22/02/2017
Exemple
2s 2 5 s 2
E(s) Y(s)
G s
G(s) 3s 2 15 s 12
Numérateur :
2s2 5 s 2 0 2 s 0.5 s 2 0
2 zéros : z1 = -0.5 ; z2 = -2
Dénominateur :
3s2 15 s 12 0 3 s 4 s 1 0
2 (s + 0.5)(s + 2) 1 (1 2s)(1 0.5s) 2 pôles : p1 = -4 ; p2 = -1
G(s)
3 (s + 4)(s + 1) 6 (1 0.25s)(1 s)
39
• Un système est stable si la partie réelle de tous ses pôles est négative.
• Les pôles dominants sont ceux dont la partie réelle est la plus faible en valeur
absolue. Ces termes sont influents sur la réponse temporelle.
Im(s)
Pôles complexes
conjuguées à partie
réelle négative
Re(s)
40
20
22/02/2017
Définition
La caractéristique statique est définie par une relation établie (en régime
permanent) entre l’entrée e et la sortie y.
y (mA)
T(°C) i(mA)
Transmetteur 20 Gain statique : k
Y 16
k 0.32mA / C
Température : Signal électrique :
4 E 50
Emin = 50 °C Ymin = 4 mA
Emax = 100 °C Ymax = 20 mA e (°C)
50 100
41
Définition
Ce mode caractérise la phase transitoire d’un état d’équilibre à un autre
suite à la variation de la grandeur d’entrée.
Caractéristiques
• Le comportement dynamique dépend de la nature du système ainsi que
du type de l’entrée appliquée.
• Les excitations usuelles utilisées pour la caractérisation des
comportements dynamiques sont : l’impulsion de Dirac, l’échelon, la
rampe et l’excitation sinusoïdale.
42
21
22/02/2017
Objectif
Afin de connaître le comportement d’un système, on lui applique certaines
entrées particulières (appelées excitations) et on analyse les sorties
correspondantes (appelées réponses).
43
Types de réponses
Transformée de
Excitation Représentation Réponse
Laplace
(t)
Impulsion de
Dirac : 1 L t 1 Réponse impulsionnelle
e(t) = (t)
t
u(t)
Échelon
1 1
unitaire : L u t Réponse indicielle
e(t) = u(t) s
t
44
22
22/02/2017
Types de réponses
Transformée de
Excitation Représentation Réponse
Laplace
r(t)
Signal en rampe : m
m L r t Réponse en vitesse
e(t) = r(t) s2
t
e(t)
Excitation 0 Em
L e t
Em
sinusoïdale : Réponse harmonique
t
s2 02
e(t) = Em sin(0t)
45
Réponse temporelle
Régime transitoire
La réponse transitoire doit être
suffisamment rapide et amortie.
yind
Caractéristiques
Domaine (yf k% yf)
• Temps de réponse : tr
a
Durée au bout de laquelle la réponse yf
entre définitivement dans une bande
définie par k% de la valeur finale.
b
• Dépassement : D
Il est défini par :
a
D% 100 tr
t
b
46
23
22/02/2017
Réponse temporelle
Régime permanent
La réponse doit suivre l’entrée sans une
erreur considérable.
yind
Caractéristiques
• On applique au système une excitation
ef et on prélève la sortie yf en régime ef
permanent. f = yf – ef
• On détermine l’erreur statique par : yf
f yf ef
t
47
Réponse temporelle
Erreur de position
L’excitation appliquée à l’entrée est un
échelon de position : e t e u t
f
yind yind
yf ef ef
f = yf – ef
yf
t t
= 0 : Erreur de position nulle. 0 : Erreur de position non nulle.
48
24
22/02/2017
Réponse temporelle
Erreur de traînage
L’excitation appliquée à l’entrée est une
rampe : e t r t mt u t
e;y e;y
e e
y y
t t
= 0 : Erreur de traînage nulle. 0 : Erreur de traînage non nulle.
49
Débit d’entrée : Qe
Débit de sortie : Qs
Section transversale : S = 0.25 m2
Qe
Équation du réservoir :
dv dH
Qe Qs S
dt dt
Conditions initiales (t = 0) : H
H0 1 m
Qs0 0.02 m 3 .s 1
Pompe
Qe0 de
Échelon commande : volumétrique Qs
25
22/02/2017
Variation de niveau :
dH Qe Qs Qe0 Qe Qs Qe
dt S S S
dH 0.001 Qe
0.004 m.s1 H t H0 t 1 0.004 t
dt 0.25 S
Qe (m3.s-1) H (m)
1.04
0.021
1.00
0.02
t (s) t (s)
0 10 0 10
Réponse harmonique
e(t) y(t)
G(s)
E0 sin(t) Y0 sin(t+)
Transmittance : Y j
G j
E j
Module :
G dB 20 og G j
Phase :
arg G j argY j arg E j
52
26
22/02/2017
Réponse harmonique
53
Réponse harmonique
Plan de Bode
On représente, sur une échelle semi-logarithmique, les variations du gain
et de la phase en fonction de la pulsation :
abscisse : log 10
Magnitude (dB)
log 0
dB
ordonnée : G 20 log G s -10
-20
abscisse :
Phase (deg)
ordonnée :
arg G j -90
10
-1 0
10 10
1 2
10
3
10
Frequency (rad/sec)
54
27
22/02/2017
Réponse harmonique
Plan de Nyquist
C’est un lieu géométrique des points indicatifs de G(j) dans le plan
complexe, pour toutes les pulsations positives sous forme polaire :
Lieu de Nyquist Im
rayon vecteur :
G dB 20 log G s
angle polaire :
arg G j Re
55
Réponse harmonique
Plan de Black
C’est un lieu géométrique équivalent au plan de Bode mais dans un seul
plan :
GdB
56
28
22/02/2017
Systèmes primordiaux
Tout système linéaire peut être régi par des équations différentielles linéaires
(à coefficients constants).
L’ordre de ces équations définit le type du système :
Y s bm sm b1 s b0
G s
e(t)
Système y(t) E s an sn a1 s a0
d’ordre n
bm m
n m G j 0
an n
On s’intéresse particulièrement aux systèmes d’ordres 1 et 2 puisque :
• La plupart des systèmes physiques appartiennent à cette catégorie.
• Leur étude permet d’approcher pour les systèmes d’ordre supérieur.
57
Définition
Le système est représenté par une équation différentielle d’ordre 1 :
e(t) y(t) dy t
y t k e t
Système
d’ordre 1 T
dt
58
29
22/02/2017
Fonction de transfert
On exprime la transformée de Laplace de l’équation caractéristique :
k T
Y s E s y 0
E(s) Y(s)
Système
d’ordre 1 1 Ts 1 Ts
Régime forcé Régime libre
(permanent) (transitoire)
Y s k
Fonction de transfert : G s
E s 1 Ts
59
1. Réponse impulsionnelle
On applique à l’entrée une impulsion de Dirac : e(t) = (t)
k k k tT
Y s E s 1 y t e ut
1 Ts 1 Ts T
• Essai significatif pour la stabilité.
e,y
1/
• Système du 1er ordre :
e(t)
Naturellement stable k/T
y(t)
t
T
60
30
22/02/2017
2. Réponse indicielle
On applique à l’entrée un échelon : e(t) = e0 u(t) avec y(0) = 0
Y s
k
1 Ts
E s
k e
0
1 Ts s
y t k e0 1 e
tT
ut
• Temps de réponse à 5% :
e,y
y tr 0.95 k e0 tr 3T e0
• Gain statique : 1.05 ke0
y finale ke0 yf
k
e position :
• Erreur de
0.95 ke0
t
y fin e fin e0 k 1 dy t k e0
T tr
a
dt t 0 T
61
m (t-T)
tT
y t k m t T T e u t e,y
e(t) = mt k=1
km (t-T)
• Régime permanent (t) : e(t-T)
k1
y t k m t T
yf
k 1 y t m t T e t T
k 1 y t : sortie instable t
T
62
31
22/02/2017
Mise en équation
La réponse harmonique (ou fréquentielle) d’un système correspond à sa
sortie par rapport à une excitation sinusoïdale.
dy t
e(t)
Système y(t) T y t k e t
E0 sin(t) d’ordre 1 Y0 sin(t+) dt
k k E 2
Y s E s 20 2
1 Ts 1 Ts s
k E0
Par identification Y0
1 T
2
des termes
k k E0 2 arctg T
Y s E s
1 Ts 1 Ts s2 2
63
Plan de Bode
On représente, sur une échelle semi-logarithmique, les variations du gain
et de la phase en fonction de la pulsation .
Fonction de transfert : Y s k
G s
E s 1 Ts
Expression du gain : k
G dB 20 log G s 20 log
1 T
2
Expression de la phase :
arg G s arctg T
On utilise un axe logarithmique où l’on porte en abscisse log10 ().
64
32
22/02/2017
• En Basse Fréquence : 1
G dB 20 log k
T
Asymptote horizontale à 20 log(k)
• En Haute Fréquence :
1
G dB 20 log k 20 log T
T
Droite asymptotique de pente (-20)
Octave : 2 = 2 1 Décade : d = 10 1
G dB 2 G dB 1 6 dB G dB d G dB 1 20 dB
GdB
c 2 G dB c G0 1 dB
-5 Diagramme
- 10 asymptotique
c G dB c G0 3 dB -20 dB/décade
- 15 Courbe
2c G dB c G0 7 dB réelle
G0 – 20
c/100 c/10 c 10 c 100 c
66
33
22/02/2017
1
Point
d’inflexion
T 2 – /2
Asymptote horizontale : = -/2
67
Plan de Black
C’est une représentation de synthèse par rapport à la représentation de
Bode :
68
34
22/02/2017
k
G dB 20 log
1 T
2
GdB
arctg T – /2 – /4
(=0)
(k = 1)
0
• B.F. : G 20 log k
0 dB – 3 dB
0
(=c)
• H.F. :
G
dB
2
()
69
Plan de Nyquist
C’est une représentation dans le plan complexe où l’on associe à chaque
valeur de une forme polaire :
70
35
22/02/2017
k 1 jT k k T
G j j Re
1 jT 1 jT 1 T
2
1 T
2 0
– /4
() (=0)
• B.F. :
0 G j k 0 k j0 (=c)
• H.F. :
G j 0 90 0 j0
71
36