Automatique Continue: Bernard BAYLE

Automatique Continue
Bernard BAYLE
Télécom Physique Strasbourg

Cours Automatique Continue
Enseignant responsable
Bernard Bayle : Pr. Télécom Physique Strasbourg.
Chercheur dans l’équipe AVR d’ICube.
http://eavr.u-strasbg.fr/~bernard
Organisation
CM : 12 × 1h45 (' 50 min.+10 min. de pause+50 min.)
TD : 10 × 1h45
TP : 1 × 4h Matlab + 3 × 4h Automatique
Enseignant responsable
Bernard Bayle : Pr. Télécom Physique Strasbourg.
Chercheur dans l’équipe AVR d’ICube.
http://eavr.u-strasbg.fr/~bernard
Evaluation (documents autorisés)

CCF : 1h45
Quelques règles :
Portable éteint.
Mains sur la table.
Je détruis :
Les 20 minutes (et je sors les étudiants les lisant).
Les téléphones qui sonnent en cours.
Notion de système
Système
Etymologiquement : ensemble organisé.
Notion de système
Système
Système et Automatique
Procédé de nature quelconque : électrique, mécanique,
chimique, économique, . . . d’entrée u et de sortie y .
Notion de système
Système
Système et Automatique
Procédé de nature quelconque : électrique, mécanique,
chimique, économique, . . . d’entrée u et de sortie y .
Système dynamique
Procédé évoluant au cours du temps sous l’action de son
entrée u.
Notion de système asservi
Système asservi
Contre-réaction : entrée du système ajustée en réaction aux
informations de sortie.
Système asservi
A la douche . . .
Douche à un bouton. . .
Système asservi
A la douche . . .
. . . après la douche brûlante, la douche à deux robinets. . .
Système asservi
A la douche . . .
. . . un thermostat et c’est réglé !
Système asservi
A la douche . . .
A quoi ça sert ?
Avoir la température désirée. . .
Système asservi
A la douche . . .
A quoi ça sert ?
Avoir la température désirée. . .
. . . le débit d’eau chaude baisse ?
Systèmes asservis au cours du temps (1)
Clepsydre - Ktesibios d’Alexandrie (270 avant JC) - . . .

Antique horloge à eau.
Machine de Watt 1769

Régulation de vitesse d’une machine à vapeur.
Théories classiques
Stabilité - Maxwell, Routh, Lyapunov, 1868-1893.
Servomécanismes - Nyquist, Bode, Evans, 1932-1947.
Théories modernes
Représentation d’état, Bellman, Kalman, Pontryagin, 1950.
Ce que vous allez apprendre. . .
Modélisation
Ecriture des lois de la physique, lorsque les paramètres du
système sont relativement bien connus.
Modélisation
Identification d’un modèle dans le cas contraire.
Modélisation
Analyse
Caractérisation des différentes propriétés du système.
Modélisation
Analyse
Caractérisation des différentes propriétés du système.
Commande
Modification de l’entrée du système par un correcteur. Mise en
œuvre de ce correcteur.
Objectifs et cadre de l’étude
Systèmes étudiés
Cas des systèmes monovariables : c’est-à-dire possédant une
seule entrée u et une seule sortie y .
Systèmes étudiés
Analogique ou numérique ?
Cas des systèmes à temps continu : u(t) et y (t) fonctions d’une
variable continue.
Systèmes étudiés
Analogique ou numérique ?
Cas des systèmes à temps continu : u(t) et y (t) fonctions d’une
variable continue.
Linéaire ou non ?
Etude des systèmes au comportement linéaire : soit linéaires,
soit linéarisés autour d’un point de fonctionnement.
Première partie I
Systèmes à temps continu

Plan
1 Modélisation
Mise en équation d’un système physique
Propriétés des systèmes à temps continu
Représentations temporelles des systèmes
Représentation par fonction de transfert
Correspondances entre représentations
2 Réponses des systèmes à temps continu

Réponse temporelle
Réponse temporelle des systèmes élémentaires
Réponse harmonique
Réponse harmonique des systèmes élémentaires
Simplification de la représentation d’un système
3 Analyse des systèmes à temps continu

Commandabilité et observabilité
Stabilité
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Premier exemple : moteur à courant continu (1)
Description
Energie électrique : alimentation de l’induit mobile (rotor)
Champ magnétique de l’inducteur (séparé et constant)
Forces électriques (addition=balais+collecteur) : rotation
Système : entrée=alimentation induit, sortie=rotation rotor
Modélisation
Mécanique : principe fondamental de la dynamique
Electrique : loi des mailles
Entrée=tension d’induit, sortie=vitesse de rotation du rotor
i(t)
R L
i(t)
γ(t) e(t)
u(t) f u(t)

ω(t)

Modélisation
Modélisation mécanique : γ − f ω = J dω
dt
di
Modélisation électrique : Ri + L dt +e =u
Loi couple (par construction) : γ = Km i
Loi vitesse (par construction) : e = Ke ω
i(t)
R L
i(t)
γ(t) e(t)
u(t) f u(t)

ω(t)

Equations mécaniques : γ − f ω = J dωdt et γ = Km i
di
Equations électriques : Ri + L dt + e = u et e = Ke ω
Relation entrée u - sortie ω avec (Ke = Km = K ) :

di

Ki = f ω + J dω
dt
di

Ki = f ω + J dω
dt
di 2
K dt = f dω d ω
dt + J dt 2
di

Ki = f ω + J dω
dt
di 2
K dt = f dω d ω
dt + J dt 2
R dω
L dω d 2ω

K f ω + J dt + K f dt + J dt 2 +K ω =u
di

Ki = f ω + J dω
dt
di 2
K dt = f dω d ω
dt + J dt 2
R dω
L dω d 2ω

d 2ω RJ+Lf dω Rf +K 2 K
dt 2
+ LJ dt + LJ ω= LJ u
di

Ki = f ω + J dω
dt
di 2
K dt = f dω d ω
dt + J dt 2
R dω
L dω d 2ω

d 2ω RJ+Lf dω Rf +K 2 K
dt 2
+ LJ dt + LJ ω= LJ u
Relation entrée-sortie
Equation différentielle d’ordre 2, linéaire à coef. constants.
Deuxième exemple : suspension (1)
Description
Véhicule : masse mv
Modèle suspension : ressort amorti (ks , fs )
Roue : masse mr , modèle roue + pneu : ressort idéal (kr )
Système : entrée=profil route, sortie=altitude véhicule
système à l’équilibre
mv zv
4
ks fs
zr
mr
caillou !
kr
~z zc
~x
Modélisation
Deux solides à l’équilibre
Mécanique : principe fondamental de la statique
Equilibre
−kr ∆lr
−ks ∆ls
ks fs
mv
véhicule 4
mr roue
kr
ks fs −mr g
ks ∆ls
− m4v g
Modélisation
Deux solides en mouvement
Mécanique : principe fondamental de la dynamique
Variations −kr (zr − zc )
−ks (zv − zr )
−fs dtd (zv − zr )
ks fs
véhicule mv
zv 4
mr roue
zr
kr
ks fs fs dtd (zv − zr )
ks (zv − zr )
zc
Modélisation

mv d 2 zv dzv dzr
Véhicule : 4 dt 2 = −ks (zv − zr ) − fs dt − dt
2

Roue : mr ddtz2r = −kr (zr − zc ) + ks (zv − zr ) + fs dzv
dt − dzr
dt
−kr (zr − zc )
−ks (zv − zr )
−fs dtd (zv − zr )
véhicule mv
zv 4
mr roue
zr
fs dtd (zv − zr )
ks (zv − zr )
Relation entrée zc - sortie zv ?
Troisième exemple : régulateur de niveau (1)
Description
Réservoir alimenté en liquide
Niveau réglé par une vanne d’entrée
Système : entrée=débit vanne, sortie=niveau liquide
vanne entrée
q0 + qe
h0 + h
vanne sortie
q0 + qs
R
Modélisation
Point de fonctionnement nominal (q0 , h0 )
Résistance du tuyau de sortie : variations hauteur de
liquide/débit de sortie (variations notées h et qs )
Hauteur de liquide : fonction de la section C et de la
différence de débit entrée - sortie
vanne entrée
q0 + qe
h0 + h
vanne sortie
q0 + qs
R
Modélisation
Résistance (selon nature de l’écoulement)
hauteur de liquide
h0 + h
α
h0
q0 q0 + qs débit de sortie
h
R = tan α(q0 ) = qs
Section : Cdh = (qe − qs )dt

Equation résistance : R = qhs
Equation écoulement : Cdh = (qe − qs )dt
Relation entrée qe - sortie h :


C dh
dt = qe −
h
R

C dh
dt = qe −
h
R
RC dh
dt + h = Rqe

C dh
dt = qe −
h
R
RC dh
dt + h = Rqe
Equation différentielle d’ordre 1, linéaire à coef. constants.
Quatrième exemple : échangeur thermique (1)
Description
Echange thermique entre circuit de vapeur et circuit d’eau
Débit de vapeur réglé par une vanne d’entrée
Système : entrée=débit vapeur, sortie=température eau
wv
θve
θee
θe θm
θv
Modélisation
Chaleur amenée par la vapeur : qv = wv cv (θve − θv )
(wv : débit massique et cv : chaleur massique de la vapeur)
Bilan de chaleur : qe = θv −θ
R
e
(R : résistance thermique moyenne de l’échangeur)

Flux de chaleur net : Cv dθ
dt = qv − qe
v
(Cv : capacité thermique de la vapeur)
wv
θve
θee
θe θm
θv
Chaleur amenée : qv = wv cv (θve − θv )

θv −θe
Bilan de chaleur : qe = R
Flux de chaleur net : Cv dθ v
dt = qv − qe
Equations du système :

θv −θe
dt = qv − qe
Flux de chaleur net, vapeur :
θv −θe
Cv dθ
dt = wv cv (θve − θv ) − R
v

θv −θe
dt = qv − qe
θv −θe
Cv dθ
v
θv −θe
Flux de chaleur net, eau : Ce dθ
dt = we ce (θee − θe ) +
e
R

θv −θe
dt = qv − qe
θv −θe
Cv dθ
v
θv −θe
e
R
Mesure : θm (t) = θe (t − τ )

θv −θe
dt = qv − qe
θv −θe
Cv dθ
v
θv −θe
e
R
Mesure : θm (t) = θe (t − τ )
Système non linéaire avec retard.
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Linéarité
Définition
Soit y1 (t) et y2 (t) les réponses d’un système Σ excité
séparément par les entrées u1 (t) et u2 (t) et α ∈ R.
Le système est linéaire si sa sortie vaut αy1 (t) + y2 (t) en

réponse à l’entrée αu1 (t) + u2 (t).
u1 (t) y1 (t)
Σ
linéarité αu1 (t) + u2 (t) αy1 (t) + y2 (t)

Σ
u2 (t) y2 (t)
Σ
Principe de superposition.
Invariance
Définition
Soit y (t) la réponse d’un système Σ d’entrée u(t).
Le système est invariant si une même commande, appliquée à

deux instants différents produit la même sortie aux instants
considérés.
u(t) y (t) invariance u(t + τ ) y (t + τ )

Σ Σ
Principe de causalité
Définition
Un signal f (t) à temps continu est causal si f (t) = 0, ∀t < 0.
Définition
Soit y (t) la réponse d’un système d’entrée u(t).
Le système est causal si, ∀t < 0, u(t) = 0 ⇒ y (t) = 0.

La réponse du système ne précède pas son excitation.
Tout système physiquement réalisable est causal.
Hypothèse
Tous les signaux et systèmes étudiés sont causaux.
Linéarité et invariance des systèmes étudiés
Invariance
Commande appliquée à l’instant t + τ : même sortie qu’en
appliquant cette commande à l’instant t, décalée de τ .
MCC : échauffement du circuit = variation de la résistance
Linéarité
Linéarité : hypothèses de modélisation.
Suspension : pas de frottement sec
Niveau de liquide : étude en un point de fonctionnement
Echangeur thermique : non linéaire ?
Retard
Sans influence sur la linéarité et l’invariance.
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Représentation externe par une équation différentielle
Propriété
Un système à temps continu linéaire et invariant sans retard est
décrit par une équation différentielle linéaire à coefficients
constants.
Notations
n m
X d i y (t) X d i u(t)
ai = bi
dt i dt i
i=c i=0
Propriété
constants.
Notations
n m
ai = bi
dt i dt i
i=c i=0
ai et bi ∈ R tel que ac , an , b0 et bm 6= 0
Propriété
constants.
Notations
n m
ai = bi
dt i dt i
i=c i=0
n, m ∈ N tel que m 6 n pour un système causal
Propriété
constants.
Notations
n m
ai = bi
dt i dt i
i=c i=0
n : ordre du système, c 6 n : classe du système
Propriété
constants.
Notations
n m
ai = bi
dt i dt i
i=c i=0
n : ordre du système, c 6 n : classe du système
y (t) : n CI pour y et m CI pour u
Relation entrée-sortie du MCC
Entrée : tension d’induit u
Sortie : vitesse de rotation du rotor ω
d 2 ω RJ + Lf dω Rf + K 2 K
2
+ + ω= u
dt LJ dt LJ LJ
2
+ + ω= u
dt LJ dt LJ LJ
n = 2 : système d’ordre 2
2
+ + ω= u
dt LJ dt LJ LJ
c = 0 6 n : système de classe 0
2
+ + ω= u
dt LJ dt LJ LJ
c = 0 6 n : système de classe 0
dω
ω(t) : 2 CI pour ω, dt (0) et ω(0)
Représentation d’état du système (1)
Propriété
décrit par une infinité de représentations sous forme de
systèmes différentiels d’ordre un, à coefficients constants.
Représentation d’état
dx
dt = Ax + Bu
y = Cx + Du
Notations
dx
dt = Ax + Bu, équation d’évolution, d’état
y = Cx + Du, équation de sortie, de mesure
Notations
dx
x de dimension n × 1 : vecteur d’état, représentation

interne du système
Notations
dx

interne du système
A, B, C et D : matrices constantes
Notations
dx

interne du système
A de dimension n × n : matrice d’évolution (ou d’état)
Notations
dx

interne du système
B de dimension n × 1 : matrice de commande (ou d’entrée)
Notations
dx

interne du système
C de dimension 1 × n : matrice d’observation
Notations
dx

interne du système
D scalaire : coefficient de transmission directe
Notations
dx

interne du système
D scalaire : coefficient de transmission directe
x(t) et y (t) : CI x(t0 )
Modèle d’état du MCC
Entrée u, sortie ω
Variables d’état indépendantes x1 = ω et x2 = i
Réécriture des équations :
di
Ri + L + Kω = u
dt
dω
Ki − f ω = J
dt
dx2
Rx2 + L + Kx1 = u
dt
dx1
Kx2 − fx1 = J
dt
dx2
Rx2 + L + Kx1 = u
dt
dx1
Kx2 − fx1 = J
dt
! ! !
− Jf K

d x1 J x1 0
= + u,
dt x2 − KL − RL x2 1
L

x1
y = 1 0
x2
dω
Variables d’état indépendantes x1 = ω et x2 = dt
dω
+ + ω = u
dt 2 LJ dt LJ LJ
dω
dx2 RJ + Lf Rf + K 2 K
+ x2 + x1 = u
dt LJ LJ LJ
dx1
x2 =
dt
dω
dx2 RJ + Lf Rf + K 2 K
+ x2 + x1 = u
dt LJ LJ LJ
dx1
x2 =
dt
! ! !
d x1 0 1 x1 0
= 2 + u,
dt x2 − Rf LJ
+K
− RJ+Lf x2 K
LJ LJ
x1
y = 1 0
x2
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Transformée de Laplace : définition
Définition
Soit f (t) un signal à temps continu, prenant la valeur f (t) à
l’instant t. La transformée de Laplace de f (t) est définie par :
Z +∞
F (s) = L{f (t)} = f (t)e−st dt.
0
Propriété
Soit s = σ + jω. La transformée de Laplace est généralement
définie sur un demi-plan complexe pour lequel σ ∈ ]σ0 , +∞[.
La valeur σ0 définissant la limite de convergence est appelée
abscisse de convergence de la transformée.
Transformée de Laplace : calcul
Tables de transformées
Autant que possible, on utilise des tables de transformées
pré-calculées :
δ(t) → 1
U(t) → 1s
...
Calcul
Calcul direct de l’intégrale
Exemple de calcul
Calcul de la transformée de Laplace de f (t) = e−at U(t).
Exemple de calcul
Définition (s = σ + jω) :
Z +∞
F (s) = e−at e−st dt,
0
Z +∞
= e−at e−(σ+jω)t dt,
0
Z +∞
= e−(σ+a)t e−jωt dt.
0
Exemple de calcul
Convergence en valeur absolue de l’intégrale :

Z +∞ Z +∞
−(σ+a)t −jωt
|e e |dt = |e−(σ+a)t ||e−jωt |dt,
0 0
Z +∞
= |e−(σ+a)t |dt,
0
Z +∞
= e−(σ+a)t dt,
0
" #t→+∞
e−(σ+a)t
= − .
σ+a
t=0
Converge en valeur absolue ssi σ + a > 0, donc σ > σ0 = −a.

Exemple de calcul
Calcul :
Z +∞
F (s) = e−(s+a)t dt,
0
" #t→+∞
e−(s+a)t
= −
s+a
t=0
1
=
s+a
Exemple de calcul
Transformée de e−at U(t)

Finalement :
1
F (s) = L{e−at U(t)} =
s+a
pour tout s = σ + jω si et seulement si σ > −a.
Transformée de Laplace : propriétés principales
linéarité L{f (t) + αg(t)} = F (s) + αG(s), ∀α ∈ R
retard L{f (t − τ )} = e−τ s F (s), ∀τ ∈ R

nR o
t
intégration L 0 f (τ )dτ = F (s) s
n o
dérivation en t L dfdt(t) = sF (s) − f (0)
valeur finale limt→∞ f (t) = lims→0 sF (s)

nR o
+∞
convolution L{f (t) ∗ g(t)} = L −∞ f (τ )g(t − τ )dτ = F (s)G(s)
Transformée de Laplace : inversion
Définition
Soit F (s) la transformée de Laplace de f (t). La transformée de
Laplace inverse de F (s) s’écrit :
Z
−1 1
f (t) = L {F (s)} = F (s)est ds.
2πj Γ
Tables de transformées
Décomposition en éléments simples pour utiliser les tables de
transformées.
. . .sinon formule d’inversion et calcul des résidus.

Transformée de Laplace : intérêt
Réécriture du modèle du système

Pour les systèmes linéaires à temps continu : possibilité de
transformer les équations différentielles décrivant l’évolution
dynamique du système en équations algébriques en s.
Conditions initiales
n o
Dérivation en t : L dfdt(t) = sF (s) − f (0)
Termes liés aux CI6= 0

Relation entrée-sortie du MCC en t
Variable t :
RJ + Lf dω(t) LJ d 2 ω(t) K
ω(t) + 2
+ 2 2
= u(t),
Rf + K dt Rf + K dt Rf + K 2
Variable t :
ω(t) + 2
+ 2 2
= u(t),
Transformée de Laplace à CI nulles :
RJ + Lf LJ K
Ω(s) + 2
sΩ(s) + 2
s2 Ω(s) = U(s),
Rf + K Rf + K Rf + K 2
Variable t :
ω(t) + 2
+ 2 2
= u(t),
Transformée de Laplace à CI nulles :
RJ + Lf LJ K
Ω(s) + 2
sΩ(s) + 2
s2 Ω(s) = U(s),
Relation entrée-sortie du MCC en s

RJ + Lf LJ 2 K
1+ 2
s+ 2
s Ω(s) = U(s),
Fonction de transfert : définition
Soit un système linéaire invariant sans retard d’entrée u(t) et
de sortie y (t).
Relation entrée-sortie :
n m
ai = bi .
dt i dt i
i=c i=0
de sortie y (t).
En appliquant la transformée de Laplace à CI nulles :

n
X m
X
ai si Y (s) = bi si U(s),
i=c i=0
soit :
Y (s) = G(s)U(s).
de sortie y (t).
Définition
On appelle fonction de transfert (FT) du système la fraction
rationnelle :
Y (s)
G(s) = .
U(s)
Le terme synonyme transmittance est souvent utilisé.
Fonction de transfert : propriétés
Forme de la fonction de transfert

Dans le cas des systèmes linéaires invariants sans retard la FT
prend la forme d’une fraction rationnelle :
Pm
N(s) bi si
G(s) = = Pi=0
n i
D(s) i=c ai s
Caractéristiques :
racines de N(s) : m zéros
racines de D(s) : n pôles
zéros et les pôles ∈ C

b0 + b1 s + . . . bm sm
G(s) =
ac sc + ac+1 sc+1 + . . . + an sn
Caractéristiques :
b0
K = ac : gain statique

Qm
bm (s − zi )
G(s) = Qni=1
an i=1 (s − pi )
Caractéristiques :
bm
an : coefficient de gain
Relation entrée-sortie du MCC en s :

RJ + Lf LJ 2 K
1+ 2
s+ 2
s Ω(s) = U(s),
Fonction de transfert du MCC

K
Ω(s) LJ
G(s) = = RJ+Lf Rf +K 2
.
U(s) s2 + s +
LJ LJ
donc :
K RJ + Lf Rf + K 2
N(s) = et D(s) = s2 + s+ .
LJ LJ LJ
Caractéristiques :
pas de zéro
pôles (tels que D(s) = 0) ?
On montre (. . .) :
KG
Ω(s) τel τem
G(s) = = ,
U(s) s+ 1
s+ 1
τel τem
L
τel = ,
R
RJ
τem = ,
Rf + K 2
K
et KG = .
Rf + K 2
Caractéristiques :
pas de zéro
pôles (tels que D(s) = 0) ?
On montre (. . .) :
KG
Ω(s) τel τem
G(s) = = ,
U(s) s+ 1
s+ 1
τel τem
donc deux pôles :
1
p1 = −
τel
1
et p2 = − .
τem
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Passage représentation d’état/fonction de transfert
Transformée de Laplace de la représentation d’état :
dx
dt = Ax + Bu sX (s) = AX (s) + BU(s)
y = Cx + Du Y (s) = CX (s) + DU(s)
dx
y = Cx + Du Y (s) = CX (s) + DU(s)
soit :
X (s) = (sI − A)−1 BU(s)

Y (s) = CX (s) + DU(s)
dx
y = Cx + Du Y (s) = CX (s) + DU(s)
soit :
X (s) = (sI − A)−1 BU(s)

Y (s) = CX (s) + DU(s)

Y (s) = C(sI − A)−1 B + D U(s)
Calcul des pôles
Fonction de transfert :
Y (s)
= C(sI − A)−1 B + D,
U(s)
Calcul des pôles
Y (s)
= C(sI − A)−1 B + D,
U(s)
coT (sI − A)
= C B + D,
det(sI − A)
Calcul des pôles
Y (s)
= C(sI − A)−1 B + D,
U(s)
coT (sI − A)
= C B + D,
det(sI − A)
C coT (sI − A) B

= + D.
det(sI − A)
Calcul des pôles
Y (s)
= C(sI − A)−1 B + D,
U(s)
coT (sI − A)
= C B + D,
det(sI − A)
C coT (sI − A) B

= + D.
det(sI − A)
Equation caractéristique du système = équation des pôles :
det(sI − A) = 0.
Calcul des pôles
Y (s)
= C(sI − A)−1 B + D,
U(s)
coT (sI − A)
= C B + D,
det(sI − A)
C coT (sI − A) B

= + D.
det(sI − A)
Equation caractéristique du système = équation des pôles :
det(sI − A) = 0.
Pôles et valeurs propres

Pôles du système = valeurs propres de A
Passage fonction de transfert/représentation d’état (1)
Non-unicité de la représentation d’état

∃ une infinité de représentations équivalentes.

Changement de variable régulier z = P −1 x :
d dz
(Pz) = P = APz + Bu,
dt dt
y = CPz + Du,

d dz
(Pz) = P = APz + Bu,
dt dt
y = CPz + Du,
soit :
dz
= P −1 APz + P −1 Bu,
dt
y = CPz + Du.
Type de représentation d’état

Obtention d’un modèle d’état avec A diagonale, pôles distincts

Décomposition en éléments simples :

n
Y (s) X αi
= + D.
U(s) s − λi
i=1


n
Y (s) X αi
= + D.
U(s) s − λi
i=1
Choix des variables d’état telles que :
1
Xi (s) = U(s),
s − λi


n
Y (s) X αi
= + D.
U(s) s − λi
i=1
sXi (s) = λi Xi (s) + U(s),



n
Y (s) X αi
= + D.
U(s) s − λi
i=1
dxi
= λi xi + u.
dt
Obtention d’un modèle d’état avec A diagonale

dxi 1 Pn
Pôles distincts : dt
= λi xi + u, Xi (s) = s−λi
U(s) et Y (s) = i=1 αi Xi (s) + DU(s)

dxi 1 Pn
Pôles distincts : dt
= λi xi + u, Xi (s) = s−λi
U(s) et Y (s) = i=1 αi Xi (s) + DU(s)
   
λ1 0 ... ... ... 0 1
0 λ2 0 ... ... 0 1
   
dx . . . ... ... . . . . . . . . .
 x + . . . u,
 
= 
. . .
dt  ... ... . . . . . . . . .
 . . . 
 
0 ... ... 0 λn−1 0  1
0 ... ... ... 0 λn 1

y = α1 α2 . . . . . . . . . αn x + Du.
Représentation d’état sous forme modale


Modèle d’état avec A diagonale, pôles pas tous distincts


n
Y (s) α1 α2 αp X αi
= + 2
+ ... + p
+ +D
U(s) s − λ1 (s − λ1 ) (s − λ1 ) s − λi
i=p+1


n
= + 2
+ ... + p
+ +D
U(s) s − λ1 (s − λ1 ) (s − λ1 ) s − λi
i=p+1
Variables d’état telles que :

1
X1 (s) = U(s),
s − λ1
...
1
Xp (s) = U(s),
(s − λ1 )p
1
Xi (s) = U(s), pour i = p + 1, . . . , n.
s − λi


n
= + 2
+ ... + p
+ +D
U(s) s − λ1 (s − λ1 ) (s − λ1 ) s − λi
i=p+1

1
X1 (s) = U(s),
s − λ1
...
1
Xp (s) = Xp−1 (s),
(s − λ1 )
1
Xi (s) = U(s), pour i = p + 1, . . . , n.
s − λi


n
= + 2
+ ... + p
+ +D
U(s) s − λ1 (s − λ1 ) (s − λ1 ) s − λi
i=p+1

sX1 (s) = λ1 X1 (s) + U(s),
sX2 (s) = λ1 X2 (s) + X1 (s),
...
sXp (s) = λ1 Xp (s) + Xp−1 (s),
sXi (s) = λi Xi (s) + U(s), pour i = p + 1, . . . , n.

Pôles pas tous distincts : sX1 (s) = λ1 X1 (s) + U(s), sX2 (s) = λ1 X2 (s) + X1 (s), . . .
sXp (s) = λ1 Xp (s) + Xp−1 (s) et sXi (s) = λi Xi (s) + U(s), i = p + 1, . . . , n

Pôles pas tous distincts : sX1 (s) = λ1 X1 (s) + U(s), sX2 (s) = λ1 X2 (s) + X1 (s), . . .
sXp (s) = λ1 Xp (s) + Xp−1 (s) et sXi (s) = λi Xi (s) + U(s), i = p + 1, . . . , n
λ1 0 ... ... ... ... ... ... 0 x1 1

    
1 λ1 0 ... ... ... ... ... 0   x2   0 
0 1 λ1 0 ... ... ... ... 0   x3   0 
    
. . . ... ... ... ... ... ... ... . . .   . . .  . . . 
    
dx
= . . . ... 0 1 λ1 0 ... ... 0   xp  +  0  u,
    
dt
0 ... ... ... 0 λp+1 0 ... 0   xp+1   1 
    
. . . ... ... ... ... ... ... ... . . .   . . .  . . . 
    
0 ... ... ... ... ... 0 λn−1 0  x  1
n−1
0 ... ... ... ... ... ... 0 λn xn 1

y = α1 α2 ... ... αn x + Du.
Matrice d’évolution sous forme réduite de Jordan


Obtention d’un modèle d’état avec A compagne

Fonction de transfert (an = 1) :
Y (s) N(s) b0 + b1 s + . . . bm sm
G(s) = = = ,
U(s) D(s) a0 + a1 s + . . . + sn

Y (s) N(s) b0 + b1 s + . . . bm sm
G(s) = = = ,
U(s) D(s) a0 + a1 s + . . . + sn
Variable d’état x1 telle que :
U(s) = D(s)X1 (s),

Y (s) = N(s)X1 (s),

Y (s) N(s) b0 + b1 s + . . . bm sm
G(s) = = = ,
U(s) D(s) a0 + a1 s + . . . + sn
Variable d’état x1 telle que :
U(s) = (sn + an−1 sn−1 + . . . + a0 )X1 (s),

Y (s) = (bm sm + bm−1 sm−1 + . . . + b0 )X1 (s).

Y (s) N(s) b0 + b1 s + . . . bm sm
G(s) = = = ,
U(s) D(s) a0 + a1 s + . . . + sn
Variables d’état x2 , . . . , xn telles que :
dx1
= x2 ,
dt
... ...
dxn−1
= xn .
dt

Y (s) N(s) b0 + b1 s + . . . bm sm
G(s) = = = ,
U(s) D(s) a0 + a1 s + . . . + sn
Finalement, variables d’état x1 , x2 , . . . , xn telles que :
dxi d i x1
= .
dt dt i

dxi d i x1 dxi−1
U(s) = (sn + an−1 sn−1 + . . . + a0 )X1 (s), dt
= dt i
et dt
= xi .

dxi d i x1 dxi−1
U(s) = (sn + an−1 sn−1 + . . . + a0 )X1 (s), dt
= dt i
et dt
= xi .
Il vient :
d n x1 d n−1 x1
= −a n−1 − . . . − a0 x1 + u,
dt n dt n−1

dxi d i x1 dxi−1
U(s) = (sn + an−1 sn−1 + . . . + a0 )X1 (s), dt
= dt i
et dt
= xi .
Il vient :
dxn dxn−1
= −an−1 − . . . − a0 x1 + u,
dt dt

dxi d i x1 dxi−1
U(s) = (sn + an−1 sn−1 + . . . + a0 )X1 (s), dt
= dt i
et dt
= xi .
Il vient :
dxn
= −an−1 xn − . . . − a0 x1 + u.
dt

dxi d i x1 dxi−1
U(s) = (sn + an−1 sn−1 + . . . + a0 )X1 (s), dt
= dt i
et dt
= xi .
soit :
   
0 1 0 ... ... ... 0 0
 0 0 1 0 0 
... ... 0
  
dx 
=  . . . . . . . .. ... ... ...  x + . . . u.
  
dt 
 0 ... ... ... ... 0 1  . . . 
−a0 −a1 . . . ... ... . . . −an−1 1
Matrice d’évolution sous forme compagne horizontale.


Il reste alors à déterminer la sortie.

Sortie :
Y (s) = bm sm X1 (s) + (bm−1 sm−1 + . . . + b0 )X1 (s).


Sortie :
Fonction de transfert strictement propre (m < n) :

y = b0 b1 . . . bm 0 . . . 0 x.

Sortie :

y = b0 b1 . . . bm 0 . . . 0 x.
Fonction de transfert non strictement propre (n = m) :
Y (s) = bn sn X1 (s) + (bn−1 sn−1 + . . . + b0 )X1 (s),


Sortie :

y = b0 b1 . . . bm 0 . . . 0 x.

Y (s) = bn (−an−1 sn−1 . . . − a0 )X1 (s) + U(s) + (bn−1 sn−1 . . . + b0 )X1 (s).

Sortie :

y = b0 b1 . . . bm 0 . . . 0 x.

y = b0 − a0 bn b1 − a1 bn . . . bn−1 − an−1 bn x + bn u.

Récapitulatif
   
0 1 0 ... ... ... 0 0
 0 0 1 0 0 
... ... 0
  
dx 
 x + . . .  u
=  . . . . . . . .. ... ... ...
dt 
 0
  
... ... ... ... 0 1  . . . 
−a0 −a1 . . . ... ... . . . −an−1 1

y = b0 − a0 bn b1 − a1 bn . . . bn−1 − an−1 bn x + bn u.
Forme canonique commandable.

Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Représentation d’état : calcul de la réponse (1)
Problème
1 Résolution de l’équation d’état :
dx
= Ax + Bu.
dt
2 Obtention de la réponse du système :
y = Cx + Du.
Réponse libre de l’état

Réponse libre : solution de l’équation d’état homogène :
dx
= Ax.
dt
Matrice de transition
On montre :
x(t) = Φ(t, t0 )x(t0 ).
avec Φ(t, t0 ) matrice de transition , d’ordre n tq :
dΦ(t, t0 )
= AΦ(t, t0 ), ∀t > t0 .
dt
Réponse forcée de l’état

Réponse forcée : réponse au seul signal d’entrée, CI nulles :
dx
= Ax + Bu.
dt
On montre : Z t
x(t) = Φ(t, τ )Bu(τ )dτ.
t0
Réponse complète
Réponse complète du système, par superposition :
Z t
x(t) = Φ(t, t0 )x(t0 ) + Φ(t, τ )Bu(τ )dτ,
t0
Réponse complète
Z t
t0
y (t) = Cx(t) + Du(t),

Réponse complète
Z t
t0
y (t) = Cx(t) + Du(t),

Z t
y (t) = C Φ(t, t0 )x(t0 ) + Φ(t, τ )Bu(τ )dτ + Du(t).
t0
Matrice de transition
Pour un système linéaire invariant :
Φ(t, t0 ) = eA(t−t0 )
où :
A2 t 2 An t n
eAt = 1 + At + + ... + + ...
2! n!
Réponse complète
Z t
y (t) = C eA(t−t0 ) x(t0 ) + eA(t−τ ) Bu(τ )dτ + Du(t).
t0
Calcul de la matrice de transition

Dans le cas où A est diagonale, chacun des termes de la
diagonale de eAt est l’exponentielle du terme correspondant de
A. En effet, ∀ai , i = 1, 2, . . . , n :
ai2 t 2 an t n
eai t = 1 + ai t + + ... + i + ...
2! n!

A. En effet, ∀ai , i = 1, 2, . . . , n :
ai2 t 2 an t n
eai t = 1 + ai t + + ... + i + ...
2! n!
Méthode générale
Dans le cas où A est non diagonale : diagonalisation préalable
Ad = P −1 AP, puis changement de variable d’état ou :
eAt = PeAd t P −1 .

A. En effet, ∀ai , i = 1, 2, . . . , n :
ai2 t 2 an t n
eai t = 1 + ai t + + ... + i + ...
2! n!
Autres méthodes
Transformée de Laplace, théorème de Sylvester.
Moteur à courant continu Maxon F2260, bobinage 885
Problème : calculer la sortie y (t) = ω(t) en réponse à une
tension constante de 12 V pour les CI en t = 0 :
ω = 2230 tours/min = 233, 5 rad/s,

i = 0, 168 A,
Représentation d’état :
! ! !
− Jf K

d x1 J x1 0
= + u,
L

x1
y = 1 0
x2
avec x1 = ω et x2 = i.
Valeur numérique du modèle :
R = 1, 44 Ω
L = 5, 6 10−4 H
J = 1, 29 10−4 kg.m2
f = 7, 19 10−5 N.s.m−1
K = 0, 10 N.m.A−1
Expression numérique de la représentation d’état :
! ! !
d x1 −0, 55768 775, 19 x1 0
= + u,
dt x2 −178, 57 −2571, 4 x2 1785, 7

x1
y = 1 0
x2
Diagonalisation de la représentation d’état :

−55, 58 0
Ad =
0 −2516, 4

0, 9975 −0, 2945
P =
−0, 07080 0, 9557
dz
= P −1 APz + P −1 Bu,
dt
y = CPz + Du.
Diagonalisation de la représentation d’état :

239, 3745
z(0) =
17, 9099

dz −55, 58 0 563, 9
= z+ u,
dt 0 −2516, 4 1910, 3

y = 0, 9975 −0, 2945 z
Réponse :
Z t
−1
y (t) = CP Φ(t, 0)z(0) + Φ(t, τ )P Bu(τ )dτ ,
0
avec : −55,58(t−τ )
e 0
Φ(t, τ ) =
0 e−2516,4(t−τ )
soit :
y (t) = 238, 77 e−55,58t −5, 2738 e−2516,4t +121, 45(1−e−55,58t )−2, 6825(1−e−2516,4t )
ou :

y (t) = 118, 77 + 117, 32 e−55,58t − 2, 5913 e−2516,4t U(t).
Représentation externe : calcul de la réponse (1)
Méthode directe
Résolution de l’équation différentielle entrée-sortie :
n m
ai = bi
dt i dt i
i=c i=0
Problématique au delà de n = 2.
Transformée de Laplace
On utilise la transformée de Laplace et donc la FT du système,
si les CI sont nulles.
Représentation externe : calcul de la réponse (2)
D’après la définition de la FT, avec des CI sont nulles :
Y (s) = G(s)U(s).
Propriété de la transformée de Laplace vis-à-vis de la

convolution :
Z +∞
L{f (t) ∗ g(t)} = L f (τ )g(t − τ )dτ = F (s)G(s)
−∞
Théorème
La réponse d’un système linéaire invariant d’entrée u(t) et de
sortie y (t) peut s’écrire sous la forme :
y (t) = g(t) ∗ u(t).

Valeur numérique du modèle :
R = 1, 44 Ω
L = 5, 6 10−4 H
J = 1, 29 10−4 kg.m2
f = 7, 19 10−5 N.s.m−1
K = 0, 10 N.m.A−1
Problème : calculer la sortie y (t) = ω(t) en réponse à une

tension constante de 12 V pour des CI nulles.
9, 8975
G(s) = ,
(1 + 0, 0184s)(1 + 0, 0004s)
d’où :
9, 8975
Y (s) = U(s),
(1 + 0, 0184s)(1 + 0, 0004s)
9, 8975
G(s) = ,
(1 + 0, 0184s)(1 + 0, 0004s)
d’où :
9, 8975 12
Y (s) = ,
(1 + 0, 0184s)(1 + 0, 0004s) s
9, 8975
G(s) = ,
(1 + 0, 0184s)(1 + 0, 0004s)
d’où :
118, 77
Y (s) = .
s(1 + 0, 0184s)(1 + 0, 0004s)
118, 77
Y (s) = ,
s(1 + 0, 0184s)(1 + 0, 0004s)
devient :
118, 77 2, 2308 9, 9808 10−4
Y (s) = − + ,
s 1 + 0, 0184s 1 + 0, 0004s
On connait :
1
L{e−at U(t)} =
s+a
donc :
t τ
L{e− τ U(t)} = .
1 + τs
La transformée inverse de Laplace de :
118, 77 2, 2308 9, 9808 10−4

Y (s) = − + ,
s 1 + 0, 0184s 1 + 0, 0004s
donne :

y (t) = 118, 77 − 121, 45 e−55,58t + 2, 6825 e−2516,4t U(t).
Réponse impulsionnelle (1)
Definition (Impulsion de Dirac)

Soit f (t) une fonction continue en 0. Alors l’impulsion de Dirac
est la distribution δ(t) telle que :
Z +∞
f (τ )δ(τ )dτ = f (0).
−∞
Réponse impulsionnelle
Réponse impulsionnelle : réponse à une impulsion de Dirac.
Réponse impulsionnelle du système :
y (t) = g(t) ∗ δ(t)
soit, par définition de l’impulsion de Dirac :

Z t
y (t) = g(τ )δ(t − τ )dτ = g(t).
0
y (t) = g(t) ∗ δ(t)

Z t
y (t) = g(τ )δ(t − τ )dτ = g(t).
0
g(t) = L−1 {G(s)} : réponse impulsionnelle du système.
y (t) = g(t) ∗ δ(t)

Z t
y (t) = g(τ )δ(t − τ )dτ = g(t).
0
g(t) = L−1 {G(s)} : réponse impulsionnelle du système.
Remarque
Faible intérêt pratique.
Approximation de l’impulsion de Dirac
Z +∞
f (τ )δ(τ )dτ = f (0).
−∞
ϕ(t)
1
a
δ(t) = lima→0 ϕ(t)
0 a t
Réponse indicielle
Définition
On appelle réponse indicielle d’un système sa réponse à un
échelon unité : (
0, si t < 0,
U(t) =
1, si t > 0.
Cette réponse vaut :

Z t Z t
y (t) = g(t) ∗ U(t) = g(τ ) U(t − τ )dτ = g(τ )dτ.
0 0
Remarque
Intérêt pratique : caractérisation (identification) du système.
Réponse indicielle
D1
105 %
100 %
95 %
90 %
Amplitude
10 %
tm t1 Temps
t5%
Systèmes du 1er ordre : réponse temporelle (1)
Définition
Un système linéaire invariant à temps continu d’ordre un est
décrit par une équation différentielle d’ordre un à coefficients
constants reliant son entrée u(t) et sa sortie y (t) :
dy (t)
y (t) + τ = Ku(t)
dt
où τ et K sont des constantes réelles non nulles ; τ est la
constante de temps du système et K son gain statique.
t
La réponse indicielle est y (t) = α + βe− τ avec α et β deux
constantes réelles dépendant des CI.
Détermination des paramètres de :
t
y (t) = α + βe− τ .
terme constant : régime permanent de la sortie

une partie variable : régime transitoire
cas τ > 0 : stable
t
y (t) = α + βe− τ .

cas τ > 0 : stable
A l’instant t = 0 :
y (0) = α + β.
t
y (t) = α + βe− τ .

cas τ > 0 : stable
dy (t)
Quand t → ∞, dt =0:
lim y (t) = K ,
t→∞
et lim y (t) = α,
t→∞
t
y (t) = α + βe− τ .

cas τ > 0 : stable
soit les paramètres :
α = K,
β = y (0) − K .
t
y (t) = α + βe− τ .

cas τ > 0 : stable
Finalement :
t t
y (t) = K (1 − e− τ ) + y (0)e− τ .
Réponse indicielle
16
14
12
10
Amplitude
0
0 0.01 0.02 0.03 0.04 0.05 0.06
Temps (s)
Réponse indicielle d’un système du premier ordre de constante de

temps τ = 0, 01 s et de gain statique K = 10 pour différentes CI
Réponse indicielle
100 %
95 %
63 %
Amplitude
τ 3τ
Temps
Caractéristiques de la réponse indicielle d’un système d’ordre 1

Systèmes du 2nd ordre : réponse temporelle (1)
Définition
Un système linéaire invariant à temps continu d’ordre deux est
décrit par une équation différentielle d’ordre deux à coefficients
constants reliant son entrée u(t) et sa sortie y (t).
Définition
On considère des systèmes dont l’équation différentielle se met
sous la forme canonique :
dy (t) d 2 y (t)
ωn2 y (t) + 2ξωn + = K ωn2 u(t),
dt dt 2
où ξ, ωn et K sont des constantes réelles strictement positives
Définition
On considère des systèmes dont l’équation différentielle se met
sous la forme canonique :
dy (t) d 2 y (t)
ωn2 y (t) + 2ξωn + = K ωn2 u(t),
dt dt 2
où ξ, ωn et K sont des constantes réelles strictement positives
ξ : coefficient d’amortissement
ωn : pulsation propre non amortie ou pulsation naturelle
K : gain statique
Réponse indicielle fonction de la valeur de ξ :
p
−ξωn t sin( 1 − ξ 2 ω t + ϕ),


 α + βe n si 0 < ξ < 1,
α + (β + γt)e −ξω n t , si ξ = 1,
 √ √
−(ξ+ ξ 2 −1)ω t −(ξ− ξ 2 −1)ω t
α + βe + γe , si ξ > 1,
 n n
avec α, β, γ ∈ R dépendant des CI.

p
−ξωn t sin( 1 − ξ 2 ω t + ϕ),


 α + βe n si 0 < ξ < 1,
α + (β + γt)e −ξω n t , si ξ = 1,
 √ √
−(ξ+ ξ 2 −1)ω t −(ξ− ξ 2 −1)ω t
α + βe + γe , si ξ > 1,
 n n
ξ > 1 : aucune oscillation

p
−ξωn t sin( 1 − ξ 2 ω t + ϕ),


 α + βe n si 0 < ξ < 1,
α + (β + γt)e −ξω n t , si ξ = 1,
 √ √
−(ξ+ ξ 2 −1)ω t −(ξ− ξ 2 −1)ω t
α + βe + γe , si ξ > 1,
 n n
ξ > 1 : aucune oscillation

ξ < 1 : pseudo-oscillations
p i. e. oscillations de pulsation
fixe ωp = 1 − ξ 2 ωn , dont l’amplitude décroît
exponentiellement vers zéro. On appelle ωp
pseudo-pulsation ou pulsation amortie.
Réponse indicielle
16
14
ξ=0,2
12
ξ=0,42
10
ξ=0,71
Amplitude
ξ=1
8
ξ=1,5
ξ=3,48
0
0 0.05 0.1 0.15 0.2 0.25
Temps (s)
Réponses indicielles d’un système du second ordre pour différentes

valeurs du coefficient d’amortissement
Temps de réponse
Pas de loi simple : abaques ou simulation.
Courbe ωn t5% fonction de ξ

Temps de réponse
Pas de loi simple : approximation.
Approximation :
deux pôles réels, associés à deux constantes de temps
(τ1 >> τ2 )
− τt − τt
α + βe 1 + γe 2 ⇒ t5% ' 3τ1
deux pôles complexes : la réponse indicielle est comprise
à l’intérieur d’une enveloppe exponentielle connue :
p 3
α + βe−ξωn t sin( 1 − ξ 2 ωn t + ϕ) ⇒ t5 % '
ξωn
Réponse indicielle
20
18
16
14
12
Amplitude
10
0
0 0.05 0.1 0.15 0.2 0.25
Temps (s)
Réponse indicielle d’un système du second ordre pseudo-oscillant et

enveloppe des pseudo-oscillations
Réponse indicielle
12
10
8
Amplitude
0
0 1 2 3 4 5 6 7 8 9 10
Temps (sec)
Temps de réponse à 5 % d’un système du second ordre de coefficient

d’amortissement 0, 6 et d’un premier ordre de constante de temps ξω1 n
Premier dépassement (cas pseudo-oscillant)

Formes analytiques :
π − √ ξπ
t1 = p , D1 = e 1−ξ2 .
1 − ξ 2 ωn
Compromis optimal amortissement-rapidité ?

100
90
80
Amplitude du premier dépassement (%)
70
60
50
40
30
20
10
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Amortissement
Correspondance entre premier dépassement (D1% ) et coefficient

d’amortissement, pour un système du second ordre
Premier dépassement (cas pseudo-oscillant)

Compromis optimal amortissement-rapidité est obtenu pour :
√
2
ξ= ' 0, 7
2
Temps de réponse : D1 % = 5 % de la valeur finale : t1 = t5 % .

Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Analyse harmonique
Définition
On considère le cas d’un système linéaire invariant de FT G(s),
en régime permanent sinusoïdal de pulsation ω. On appelle
G(s = jω) réponse harmonique.
Propriété
La réponse du système à une entrée sinusoïdale A sin ωt est :
y (t) = A |G(jω)| sin (ωt + Arg{G(jω)}) .

Analyse harmonique
Définition
On considère le cas d’un système linéaire invariant de FT G(s),
en régime permanent sinusoïdal de pulsation ω. On appelle
G(s = jω) réponse harmonique.
Analyse harmonique
Analyse harmonique : étude de la fonction G(jω) :
comportement fréquentiel du système (signal périodique)
diagrammes mettant en correspondance module et
argument
Diagrammes harmoniques (1)
Diagramme de Bode
Le diagramme de Bode est constitué de deux courbes :
module en décibels (dB) : GdB (ω) = 20 log10 |G(jω)|
argument en degrés (deg) : ϕ(ω) = Arg{G(jω)}
Vocabulaire
On utilise traditionnellement les termes de gain et de phase,
plutôt que les termes modules et argument.
Intérêt du diagramme de Bode

Module d’un produit de nombres complexes = produit de
leurs modules. Module en dB d’un produit de nombres
complexes = somme de leurs modules en dB :
20 log10 |G1 (jω)G2 (jω)| = 20 log10 |G1 (jω)| + 20 log10 |G2 (jω)|,
= G1dB (jω) + G2dB (jω).
Argument d’un produit de nombres complexes = somme

des arguments :
Arg{G1 (jω)G2 (jω)} = ArgG1 (jω) + ArgG2 (jω).

Diagramme de Nyquist
Le diagramme de Nyquist est le lieu de G(jω) dans le plan
complexe, lorsque ω varie de −∞ à +∞.
Remarque
Ce diagramme est orienté selon les ω croissants. En général
on choisi l’échelle du diagramme de Nyquist de sorte que le
point complexe d’abscisse −1, dit point critique apparaisse et
puisse être situé par rapport au lieu de G(jω).
Lieu de Black-Nichols
Le lieu de Black-Nichols est le lieu orienté des points de
coordonnées (ϕ(ω), GdB (ω)) lorsque ω varie de −∞ à +∞. On
tache aussi de faire apparaître le point critique de coordonnées
(−180, 0) sur ce lieu.
Remarque
Ce diagramme est orienté selon les ω croissants. En général
on choisi l’échelle du diagramme de Nyquist de sorte que le
point complexe d’abscisse −1, dit point critique apparaisse et
puisse être situé par rapport au lieu de G(jω).
Valeur numérique des paramètres du modèle :
KG
Ω(s) τel τem
G(s) = = ,
U(s) s+ 1
s+ 1
τel τem
L
τel = ,
R
RJ
τem = ,
Rf + K 2
K
et KG = .
Rf + K 2
9, 8975
G(s) = .
(1 + 0, 0184s)(1 + 0, 0004s)
Diagramme de Bode
20
0
Amplitude(dB)
−20
−40
−60
−80
0 1 2 3 4 5
10 10 10 10 10 10
−20
−40
−60
Phase (deg)
−80
−100
−120
−140
−160
0 1 2 3 4 5
10 10 10 10 10 10
Pulsation (rad/s)
1
−1
Axe imaginaire
−2
−3
−4
−5
−6
−2 0 2 4 6 8 10
Axe réel
Lieu de Black−Nichols
20
10
−10
−20
Gain (dB)
−30
−40
−50
−60
−70
−180 −160 −140 −120 −100 −80 −60 −40 −20

Phase (deg)
Systèmes du 1er ordre : réponse harmonique (1)
Fonction de transfert
La FT d’un système du premier ordre est donc :
Y (s) K
G(s) = = .
U(s) 1 + τs
Réponse harmonique
La réponse harmonique associée est :
K
G(jω) = .
1 + jτ ω
Description de la réponse harmonique :

étude du comportement asymptotique du régime
permanent sinusoïdal
extrapolation par des valeurs choisies
K
Description de la réponse harmonique : G(jω) = 1+jτ ω .
ω G(jω) gain gain phase
équivalent (dB) (deg)

1
ω K K KdB = 20 log10 K 0
τ
1 K K
√ KdB − 3 −45
τ 1+j 2
1 −jK K
ω KdB − 20 log10 τ − 20 log10 ω −90
τ τω τω
Diagramme de Bode
K
dB
20
3 dB
Amplitude(dB)
10
−2
0 0d
B/d
éca
de
−10
−20
0 1 2 3 4
10 10 10 10 10
−10
Phase (deg)
−20
−30
−40 −45
−50
−60
−70
−80 ωc
−90
0 1 2 3 4
10 10 10 10 10
Pulsation (rad/s)
Diagramme de Bode d’un système du premier ordre de constante de

temps τ = 0, 01 s et de gain statique K = 10
20 K
dB
15
10
5
Gain (dB)
−5
−10
−15
−90
−20
−180 −160 −140 −120 −100 −80 −60 −40 −20
Phase (deg)
Lieu de Black-Nichols d’un système du premier ordre de constante de

temps τ = 0, 01 s et de gain statique K = 10
1
−1
Axe imaginaire
−2
−3
−4
−5
K
−2 0 2 4 6 8 10
−1 Axe réel K/2 = 5
Diagramme de Nyquist d’un système du premier ordre de constante

de temps τ = 0, 01 s et de gain statique K = 10
Systèmes du 2nd ordre : réponse harmonique (1)
La FT du système du second ordre est :
Y (s) K ωn2
G(s) = = 2 .
U(s) ωn + 2ξωn s + s2
Y (s) K ωn2
G(s) = = 2 .
Pôles = solutions de ωn2 + 2ξωn s + s2 = 0 :

p
p1,2 = −(ξ ± j 1 − ξ 2 )ωn si 0 < ξ 6 1
p
et p1,2 = −(ξ ± ξ 2 − 1)ωn si ξ > 1.
Y (s) K ωn2
G(s) = = 2 .
Axe imaginaire
ωn
p1 p
ωn 1 − ξ2
cos(ψ) = ξ ψ 0 Axe réel

−ξωn
p
-ωn 1 − ξ2
p2
K ωn2
Description de la réponse harmonique : G(jω) = 2
ωn −ω 2 +2jξωn ω
.
ω G(jω) gain gain phase
équivalent (dB) (deg)
ω ωn K K KdB = 20 log10 K 0
K K
ωn 2jξ 2ξ KdB − 6 − 20 log10 ξ −90
−K ωn2 K ωn2
ω ωn ω2 ω2
KdB + 40 log10 ωn − 40 log10 ω −180
Résonance
√
Si ξ 6 22 , il peut y avoir un phénomène de résonance,
c’est-à-dire que le gain passe par un maximum en :
p
ωr = 1 − 2ξ 2 ωn .
Le dépassement en gain correspondant vaut :
1
Gr = p
2ξ 1 − ξ 2
Diagramme de Bode
K
20 dB
Amplitude(dB)
0
−4
ξ= 0,2 0d
−20 B/d
0,42 éca
0,71 de
−40 1
1,5
3,48
−60
0 1 2 3 4
10 10 10 10 10
−20
−40
−60
Phase (deg)
−80 −90
−100
−120
−140
−160
−180 ωc
0 1 2 3 4
10 10 10 10 10
Pulsation (rad/s)
Diagramme de Bode d’un système du second ordre K = 10,

ωn = 100, ξ variable
30
K ξ= 0,2
20 dB
0,42
10 0,71
1
0
1,5
Gain (dB)
−10
3,48
−20
−30
−40
−50
−180 −160 −140 −120 −100 −80 −60 −40 −20

Phase (deg)
Lieu de Black-Nichols d’un système du second ordre K = 10,

0
−5 1,5 3,48
0,71 1
0,42
ξ= 0,2
−10
Axe imaginaire
−15
−20
−25
−10 −5 0 5 K= 10 15
Axe réel
Diagramme de Nyquist d’un système du second ordre K = 10,

Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Simplification de modèles (1)
Forme générique d’une FT :
Qp Qq ξi 1 2
K i=1 (1 + τi s)
i=p+1 (1 + 2 ωn i s + ωn 2 s )
i
G(s) = c Qp0 Qq 0 ξj
.
s 1 2
j=1 (1 + τj s) j=p0 +1 (1 + 2 ωn j s + ω 2 s )
nj
Qp Qq ξi 1 2
K i=1 (1 + τi s)
i=p+1 (1 + 2 ωn i s + ωn 2 s )
i
.
s 1 2
j=1 (1 + τj s) j=p0 +1 (1 + 2 ωn j s + ω 2 s )
nj
Termes du premier ordre : pôles et zéros réels

Qp Qq ξi 1 2
K i=1 (1 + τi s)
i=p+1 (1 + 2 ωn i s + ωn 2 s )
i
.
s 1 2
j=1 (1 + τj s) j=p0 +1 (1 + 2 ωn j s + ω 2 s )
nj
Termes du premier ordre : pôles et zéros réels

Termes du second ordre : pôles et zéros complexes
conjugués
Simplification du modèle du MCC
Modèle du second ordre :
9, 8975
G(s) = .
(1 + 0, 0184s)(1 + 0, 0004s)
Modèle du premier ordre :
9, 8975
G(s) = .
1 + 0, 0184s
Réponse harmonique relative au pôle dominant ?
1
G1 (jω) = ,
1 + 0, 0184jω
Réponse harmonique relative au pôle secondaire ?
1
G2 (jω) =
1 + 0, 0004jω
KdB = 20
(−1)
ω
0
(−1)
(KG1 )dB
1
τem
ω
0
−90
Arg{KG1 }
−180
KdB = 20
(−1)
ω
0
(−1)
(−1) (G2 )dB
(KG1 )dB
1 1
τem τel
ω
0
Arg{G2 }
−90
Arg{KG1 }
−180
KdB = 20
(−1)
ω
0
(−1)
(−1) (G2 )dB
(KG1 )dB
(−
2)
GdB
1 1
τem τel
ω
0
Arg{G2 }
−90
Arg{KG1 }
−180 Arg{G}
KdB = 20
(−1)
ω
0
(−1)
(KG1 )dB
(−
2)
GdB
1
τem
ω
0
−90
Arg{KG1 }
−180 Arg{G}
0
−0.02
−0.04
−0.06
−0.08
Amplitude
−0.1
−0.12
−0.14
−0.16
−0.18
−0.2
0 0.02 0.04 0.06 0.08 0.1 0.12
Temps (s)
Erreur de modélisation dans l’hypothèse d’un modèle d’ordre un du

MCC : erreur sur la réponse indicielle pour un moteur Maxon F2260
Règles générales : pôles réels

Lorsque deux pôles sont suffisamment distincts le pôle le plus
près de l’axe des imaginaires, c’est-à-dire le plus petit en valeur
absolue, associé à la constante de temps la plus lente, est
prépondérant.
Si l’on doit faire une approximation pour simplifier l’étude d’un

système, dont le modèle est d’ordre élevé, on négligera donc
les pôles les plus rapides.
Si les pôles sont proches, il peut devenir plus hasardeux

d’effectuer une telle simplification.
Règles générales : pôles complexes conjugués

On pourra, de même considérer que la dynamique liée à une
paire de pôles complexes conjugués est négligeable devant
celle liée à un pôle simple ou à une autre paire de pôles
complexes conjugués si la pulsation naturelle associée à cette
paire est grande devant la pulsation naturelle de l’autre paire,
ou devant la pulsation associée au pôle simple.
Règles générales : zéros

Cas similaire : on simplifiera les zéros entre eux de la même
manière.
En revanche, on procédera avec prudence pour ce qui est de

négliger un zéro prépondérant au vu de la valeur des pôles.
Exemple 1
K (1 + τ2 s)
G(s) =
(1 + τ1 s)(1 + τ3 s)
avec τ1 > τ2 > τ3 .
A. N. :
τ1 = 1 s
τ2 = 0, 1 s
τ3 = 0, 01 s
Caractéristiques
Tracé du diagramme de Bode
Tracé du lieu des pôles et zéros
Exemple 2
K (1 + τ1 s)
G(s) = 2

1 + ωs2 + ωs 2 (1 + τ3 s)
2
1
avec τ1 > ω2 > τ3 .
A. N. :
τ1 = 1 s
ω2 = 10 rad/s
τ3 = 0, 01 s
Caractéristiques
Tracé du diagramme de Bode
Tracé du lieu des pôles et zéros
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Commandabilité et observabilité : généralités
Contexte
Notions fondamentales pour l’étude des systèmes en vue de
leur commande, définies à partir du modèle d’état du système
considéré
Objectif
Caractérisation d’après les critères de Kalman.
Commandabilité (1)
Définition
Commandabilité d’un système = capacité à voir son
comportement évoluer sous l’action de sa commande (son
entrée).
Commandabilité d’un système

Une variable d’état xi est commandable s’il est possible de
déterminer une commande u(t) sur un intervalle [t0 , tf ]
conduisant tout état initial xi (t0 ) en 0 en un temps t1 avec
t0 6 t1 6 tf .
Si cette propriété est vraie pour tout t0 et pour toute variable du
vecteur d’état, alors le système est dit complètement
commandable.
Commandabilité (2)
Critère de commandabilité (de Kalman)

Un système linéaire invariant d’équation dynamique d’état :
dx
= Ax + Bu
dt
est commandable si et seulement si la matrice de
commandabilité :
C = B AB . . . An−1 B

est de rang n.
On parle couramment de commandabilité de la paire (A, B).

Système de niveau de liquide
Système de niveau de liquide régi par la loi d’évolution :
dh
RC + h = Rqe .
dt
vanne entrée
q0 + qe
h0 + h
vanne sortie
q0 + qs
R
Modèle d’état
Modèle physique :
dh
RC + h = Rqe .
dt
Système d’ordre 1 : variable d’état = variation h de la hauteur

de liquide, entrée (commande) = variations u = qe du débit
d’entrée autour du débit nominal q0 .
Equation d’état :
dx 1 1
=− x+ u
dt RC C
Commandabilité
1
Système commandable car B = C 6= 0
hr
vanne entrée
q0 + qe
qr
h0 + h
vanne sortie
q0 + qs
R
Système de niveau de liquide modifié

Système de niveau de liquide régi par les lois d’évolution :
Cdh = (qe + qr − qs )dt Cr dhr = −qr dt,

qs = Rh qr = Rhrr .
Modèle d’état
Système d’ordre 2 : variables d’état = h et hr , commande = u.
Equation d’état :
! ! !
1 1 1
d x1 − RC Rr C x1 C
= + u.
dt x2 0 − Rr1Cr x2 0
Commandabilité
1 1

− RC
C = (B AB) = C 2
.
0 0
Système non complètement commandable, car C de rang

n − 1 = 1.
x2 = hr non commandable par la vanne d’entrée. . .
Observabilité (1)
Définition
Observabilité d’un système = possibilité de déterminer son état
à partir des mesures de sa sortie.
Observabilité d’un système

Une variable d’état xi est observable s’il est possible de
déterminer xi (t0 ) à partir de la connaissance de la sortie y(t)
sur un intervalle [t0 , tf ].
Si cette propriété est vraie pour tout t0 et pour toute variable du
vecteur d’état, alors le système est dit complètement
observable.
Observabilité (2)
Critère d’observabilité (de Kalman)

Un système linéaire invariant de représentation d’état :
dx
= Ax + Bu,
dt
y = Cx + Du
est observable si et seulement si la matrice d’observabilité :

 
C
 CA 
O= ... 

CAn−1
est de rang n.
On parle couramment d’observabilité de la paire (A, C).

vanne entrée
q0 + qe
h0 + h
vanne sortie
q0 + qs
R
Système de niveau de liquide

Représentation d’état du système :
dx 1 1
= − x+ u,
dt RC C
y = x,
Observabilité
Système observable car C = 1.
hr
vanne entrée
q0 + qe
qr
h0 + h
vanne sortie
q0 + qs
R
Système de niveau de liquide avec réservoir auxiliaire

!
1 1 1
dx − RC Rr C x1 C u,
= +
dt 0 − Rr1Cr x2 0

x1
y = 1 0
x2
Système de niveau de liquide avec réservoir auxiliaire
!
1 1 1
dx − RC Rr C x1 C u,
= +
dt 0 − Rr1Cr x2 0

x1
y = 1 0
x2
Observabilité

C 1 0
O= = 1 1 .
CA − RC Rr C
Système complètement observable car la matrice

d’observabilité est de rang n = 2.
Plan
1 Modélisation

Réponse temporelle
Réponse harmonique

Stabilité
Stabilité et état d’équilibre (1)
Etat d’équilibre
Un état d’équilibre est un état du système qui reste non modifié
lorsque le système est abandonné à lui-même.
Première définition de la stabilité : stabilité interne

De manière naturelle, on dira qu’un système est stable si,
écarté de sa position d’équilibre, il y revient.
Caractérisation d’un état d’équilibre

dx
Dans un état d’équilibre dt = Ax + Bu s’écrit :
Ax = 0.
A = application linéaire sur l’espace d’état :
dim(Im A) + dim(Ker A) = rangA + dim(Ker A) = n.
Conclusions :
si KerA = 0 : x = 0 est unique solution ;
si KerA 6= 0, alors rang A 6= n et donc det A = 0 : ∃x 6= 0
Cas où det A 6= 0
Cas du MCC, sortie=vitesse de rotation, entrée=tension
d’induit : !
− Jf K
J
A= .
− KL − RL
fR K2
Alors det A = JL + JL > 0.
Le seul état d’équilibre possible du système est obtenu quand
x1 = ω = 0 et x2 = i = 0. . .
Remarque
Il n’existe qu’une seule solution au système en x = 0. Le seul
état d’équilibre est donc l’origine de l’espace d’état.
Cas où det A = 0
Cas du MCC, sortie=position angulaire θ :
 
0 1 0
 
A=
0 − Jf K 
J 
0 − KL − RL
Alors det A = 0. Résolution de Ax = 0 : tous les triplets t.q.

x2 = x3 = 0 sont des états d’équilibre . . .
Remarque
Dans ce cas KerA 6= 0 et ∃ une infinité d’états d’équilibre tous
compris dans un espace de dimension n−rang A, que l’on
obtient en résolvant Ax = 0.
Vocabulaire
Différents types d’équilibre :
équilibre instable
équilibre asymptotiquement stable
équilibre simplement stable
équilibre instable équilibre asymptotiquement stable équilibre simplement stable

le
tab
ins
si
n
m able
io
pl
at
em
rb
st
r tu
en
pe
t
t
en
m
le e
ab qu
st toti
p
ym
as
nouvel équilibre
équilibre initial
Stabilité BIBO
Deuxième définition de la stabilité : stabilité BIBO

Un système est stable si toute entrée bornée produit une sortie
bornée.
Cette définition caractérise la stabilité entrée bornée-sortie
bornée
ESBS ou plus couramment BIBO en anglais.
Critères algébriques de stabilité (1)
Théorème
Un système linéaire invariant à temps continu est stable si tous
ses pôles sont à partie réelle strictement négative.
Ce théorème est valable pour tout système, qu’il soit en boucle

ouverte ou fermée. Pour un système d’ordre élevé, il faut
généralement recourir à une résolution numérique pour
déterminer les pôles du système.
Théorème
Un système linéaire invariant à temps continu est stable si tous
ses pôles sont à partie réelle strictement négative.
Ce théorème est valable pour tout système, qu’il soit en boucle

ouverte ou fermée. Pour un système d’ordre élevé, il faut
généralement recourir à une résolution numérique pour
déterminer les pôles du système.
Critère de Routh-Hurwitz
Soit D(s) = an sn + an−1 sn−1 + · · · + ac sc le dénominateur de la FT considérée. Le
système est stable si les ai , ∀i = c, c + 1, . . . , n sont de
même signe et du même signe que les éléments de la première
colonne du tableau suivant (dit tableau de Routh) :
an an−2 an−4 ...
an−1 an−3 an−5 ...
an−1 an−2 −an an−3 an−1 an−4 −an an−5

bn−1 = an−1
bn−3 = an−1
bn−5 ...
bn−1 an−3 −an−1 bn−3 bn−1 an−5 −an−1 bn−5

cn−1 = bn−1
cn−3 = bn−1
... ...
c b −bn−1 cn−3
dn−1 = n−1 n−3
c
... ... ...
n−1
... ... ... ...

Critère de Routh-Hurwitz
Le nombre de changements de signe dans la première colonne
indique le nombre de pôles instables du système.
Cas où apparaît un zéro dans la première colonne du tableau

On le remplace par un scalaire ε > 0, et on poursuit le calcul.
Si le scalaire directement sous ε est positif, il existe un pôle à
partie réelle nulle. Sinon, il existe un pôle à partie réelle
positive.
Stabilité des systèmes d’ordre un, deux et trois (1)
Système d’ordre un
Y (s) K
G(s) = =
U(s) 1 + τs
Unique pôle p = − τ1 .
Condition de stabilité
Un tel système est stable si τ > 0.
Système d’ordre deux

K ωn2
G(s) =
ωn2 + 2ξωn s + s2
Une paire de pôles :

p
p1,2 = −(ξ ± j 1 − ξ 2 )ωn , si 0 < ξ 6 1,
p
p1,2 = −(ξ ± ξ 2 − 1)ωn , si ξ > 1.
Condition de stabilité
Pôles complexes : système stable ssi ξωn > 0.
Pôlespréels : le système est stable système stable ssi
(ξ − ξ 2 − 1)ωn > 0.
Système toujours stable si ξ > 0 et ωn > 0.

Système d’ordre trois

K
G(s) = ,
1 + a1 s + a2 s2 + a3 s3
avec a1 , a2 et a3 tous positifs.

Pôles ? −→ Routh.
Tableau de Routh
a3 a1 0
a2 1 0
a1 a2 −a3
a2 0 0
1 0 0
0 0 0
Conclusions :
système stable ssi a1 a2 > a3
sinon, il existe deux pôles à partie réelle négative (deux
changements de signe)
Cas où a1 a2 = a3 : nouveau tableau de Routh

a3 a1 0
a2 1 0
ε 0 0
1 0 0
0 0 0
Conclusion : le système possède alors un pôle à partie réelle

nulle.
Système d’ordre trois : condition de stabilité

système stable si a1 a2 > a3
système instable si a1 a2 < a3
système stable au sens large
Exemples numériques :
éq. car. : s3 + 2s2 + s + 1 = 0, racines : −1, 755 et −0, 123 ± 0, 745j : stable ;
éq. car. : 2s3 + s2 + s + 1 = 0, racines : −0, 739 et 0, 119 ± 0, 814j : instable ;
éq. car. : s3 + s2 + s + 1 = 0, racines : −1 et ±j : limite de stabilité.
Etude d’un cas
On s’intéresse à un procédé de fonction de transfert :
5000
F1 (s) = .
s3 + 61s2 + 560s + 500
1 Calculer les pôles du système. Le système est-il stable ?

2 Appliquer le critère de Routh-Hurwitz pour résoudre le
même problème.
3 Quelle est l’influence d’un zéro ? Pour cela ajouter un zéro,
de sorte que le système ait pour fonction de transfert :
0, 2s + 1
F2 (s) = 5000 .
s3 + 61s2 + 560s + 500
Pôles :
système du troisième ordre sans zéro.
on remarque que p1 = −1 est racine du polynôme
dénominateur, ce qui nous permet de factoriser :
s3 + 61s2 + 560s + 500 = (s + 1)(s2 + 60s + 500).
racines de s2 + 60s + 500 :

√
p2/3 = −30 ± ∆0 ,
avec ∆0 = 302 − 500 = 202 . On obtient :
p2 = −10
p3 = −50
Conclusion
Pôles du système tous réels et négatifs : système stable.
Critère de Routh-Hurwitz :
1 560 0
61 500 0
61×560−1×500
551, 8 = 61 0 0
500 0 0
0 0 0
Conclusion
Tous les coefficients du polynôme dénominateur et tous les
éléments de la première colonne du tableau de Routh sont
positifs : système stable.
Etude d’un nouveau cas
On s’intéresse maintenant à un procédé de fonction de
transfert :
0, 2s + 1
F2 (s) = 5000 .
s3 + 61s2 + 560s + 500
Critère de Routh-Hurwitz ?
Conclusion
Un zéro n’introduit pas d’exception : le système reste stable.
Deuxième partie II
Asservissements à temps continu

Plan
4 Systèmes asservis à temps continu

Les différents types d’asservissement
Asservissement classique
Asservissement par retour d’état
5 Analyse des systèmes asservis à temps continu

Stabilité
Précision des systèmes à temps continu
6 Commande des systèmes à temps continu

Cahier des charges
Correcteurs PID
Méthodes harmoniques de synthèse de correcteur
Méthode du lieu des racines
Synthèse de correcteur par retour d’état
Plan


Stabilité

Cahier des charges
Correcteurs PID
Plan


Stabilité

Cahier des charges
Correcteurs PID
Notion de système asservi (1)
Schéma général
Schéma bouclé : entrée du procédé corrigée selon la différence
entre grandeur de référence et grandeur de mesure
Correcteur : produit la commande du procédé.
perturbation
référence + erreur commande

correcteur procédé
−
mesure
capteur
bruit
Notion de système asservi (2)
Techniques d’asservissement envisagées. . .

. . . pour conférer au système le comportement dynamique et le
régime permanent souhaités
Si l’état n’est pas complètement mesurable

Représentation externe du système par sa fonction de transfert
pour effectuer la correction : techniques fréquentielles,
placement de pôle, ou méthodes empiriques.
Si l’état est complètement mesurable

Utilisation prioritaire de la représentation d’état du système et
correction de type retour d’état.
Plan


Stabilité

Cahier des charges
Correcteurs PID
Schéma d’un système asservi (1)
Schéma général d’un asservissement à temps continu

perturbation
référence + erreur commande sortie

−
mesure
capteur
bruit
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
Vocabulaire et notations
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
Yr (s) : référence (ou grandeur de consigne)
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
Y (s) : sortie (ou grandeur réglée)
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
Ym (s) : mesure
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
E(s) : erreur (ou écart) de l’asservissement
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
C(s) : correcteur
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
U(s) : commande du procédé
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
P(s) : perturbation
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
B(s) : bruit de mesure
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
CG(s) = C(s)G1 (s)G2 (s) : fonction de transfert de la chaîne
directe (ou chaîne d’action)
Schéma bloc
P(s)
+
Yr (s) + E(s) U(s) + Y (s)
C(s) G1 (s) G2 (s)
−
Ym (s) +
H2 (s) H1 (s)
B(s)
H(s) = H1 (s) H2 (s) : fonction de transfert de la chaîne de
retour (ou chaîne de contre-réaction)
Fonction de transfert d’un système asservi (1)
Yr (s) + E(s) U(s) Y (s)

C(s) G(s)
−
Ym (s)
H(s)
Définitions
Yr (s) + E(s) U(s) Y (s)

C(s) G(s)
−
Ym (s)
H(s)
Définitions
Fonction de transfert en boucle ouverte du système :
FTBO : C(s)G(s)H(s),
notée CGH(s)
Yr (s) + E(s) U(s) Y (s)

C(s) G(s)
−
Ym (s)
H(s)
Définitions
Fonction de transfert en boucle fermée du système :
Y (s)
FTBF :
Yr (s)
Yr (s) + E(s) U(s) Y (s)
C(s) G(s)
−
Ym (s) U(s) = C(s)E(s),

H(s) E(s) = Yr (s) − Ym (s),
Ym (s) = H(s)Y (s),
Y (s) = G(s)U(s).
d’où :
U(s) = C(s) (Yr (s) − H(s)Y (s)) ,

Y (s) = CG(s) (Yr (s) − H(s)Y (s)) ,
et finalement :
Y (s) CG(s)
= .
Yr (s) 1 + CGH(s)
Equation caractéristique
Equation des pôles : 1 + CGH(s) = 0
Yr (s) + E(s) U(s) Y (s)

C(s) G(s)
−
En notant :
N(s)
CG(s) = ,
D(s)
on a :
Y (s) CG(s) N(s)
= = .
Yr (s) 1 + CG(s) N(s) + D(s)
Yr (s) + E(s) U(s) Y (s)

C(s) G(s)
−
Equation des pôles : N(s) + D(s) = 0.
Influence des pôles et zéros de la BO sur les pôles de la BF.
Régulation de la vitesse du rotor d’un MCC
Commande = tension d’entrée yr pour une vitesse de rotation
désirée ωr . Mesure = vitesse du rotor est mesurée par une
génératrice tachymétrique placée sur l’axe et qui délivre une

tension ym proportionnelle à la vitesse de rotation de l’axe.

Yr (s) + U(s) K Ω(s)

C(s) LJ
+K 2
s2 + RJ+Lf s+ Rf LJ
− LJ
Ym (s)
Kω

Yr (s) + U(s) K Ω(s) Ym (s)

C(s) LJ
Kω
+K 2
s2 + RJ+Lf s+ Rf LJ
− LJ

Yr (s) Ωr (s) + U(s) K Ω(s)

1
Kω
Kω C(s) LJ
+K 2
s + LJ s+ Rf LJ
2 RJ+Lf
−
Comportement fréquentiel d’un système asservi
(
Y (jω) CG(jω) ' CG(jω), si CGH(jω) 1,
= → 1
Yr (jω) 1 + CGH(jω) ' H(jω) , si CGH(jω) 1.
BO : CG(s)
(−1)
ω
0
(−
2)
ϕ(ω)
1 1
τ1 τ2 ω
0
−90
−180
(
= → 1
H(s) = 1 ω
0
ϕ(ω)
ω
0
−90
−180
(
= → 1
H(s) = 1 (−1) ω
0
(−
2)
BO : CG(s)
ϕ(ω)
1 1
τ1 ωc τ2 ω
0
−90
−180
(
= → 1
H(s) = 1 (−1) ω
0
BF
(−
2)
BO : CG(s)
ϕ(ω)
1 1
τ1 ωc τ2 ω
0
−90
−180
Plan


Stabilité

Cahier des charges
Correcteurs PID
Schéma d’un système asservi par retour d’état
Schéma général d’un asservissement par retour d’état

perturbation
référence commande sortie

état
mesure
capteur
bruit
Représentation d’état du système corrigé (1)
Retour d’état
Loi de commande u proportionnelle à l’état :
u = −Kx,

où : K = K1 K2 . . . Kn en monovariable.
Retour d’état
u = −Kx,

Retour d’état et réglage dynamique

Nouveau modèle d’évolution :
dx
= Ax − BKx
dt
Retour d’état
u = −Kx,

Retour d’état et réglage dynamique

Nouveau modèle d’évolution :
dx
= (A − BK )x.
dt
Choix de K : réglage du comportement dynamique du système.
Retour d’état et référence
Problème : introduire la consigne dans la loi de commande.
Objectif : y∞ = yr (sortie en régime permanent).
Première modification de la loi initiale

Soient x∞ et u∞ les régimes permanents de x et u. La
commande par retour d’état :
u − u∞ = −K (x − x∞ )
est telle que u = u∞ dès lors que x = x∞ .

Première modification de la loi initiale

Soient x∞ et u∞ les régimes permanents de x et u. La
commande par retour d’état :
u − u∞ = −K (x − x∞ )
est telle que u = u∞ dès lors que x = x∞ .
Reformulation du problème
Calculer x∞ et u∞ tels que y∞ = yr .
Equations d’état du système en régime permanent :
0 = Ax∞ + Bu∞ ,
y∞ = Cx∞ + Du∞ .
0 = Ax∞ + Bu∞ ,
y∞ = Cx∞ + Du∞ .
En posant :
x∞ = N x y∞ ,
u∞ = Nu y∞ ,
0 = Ax∞ + Bu∞ ,
y∞ = Cx∞ + Du∞ .
En posant :
x∞ = N x y∞ ,
u∞ = Nu y∞ ,
il vient :
0 A B Nx
y∞ = y
1 C D Nu ∞
0 = Ax∞ + Bu∞ ,
y∞ = Cx∞ + Du∞ .
En posant :
x∞ = N x y∞ ,
u∞ = Nu y∞ ,
il vient :
0 A B Nx
y∞ = y
1 C D Nu ∞
soit :
0 A B Nx
= .
1 C D Nu
Après inversion :
−1
Nx A B 0
=
Nu C D 1
Après inversion :
−1
Nx A B 0
=
Nu C D 1
l’équation initiale :
u = u∞ − K (x − x∞ )
Après inversion :
−1
Nx A B 0
=
Nu C D 1
u = u∞ − K (x − x∞ )
conduit donc à :
u = Nu y∞ − K (x − Nx y∞ ).
Après inversion :
−1
Nx A B 0
=
Nu C D 1
u = u∞ − K (x − x∞ )
conduit donc à :
u = Nu y∞ − K (x − Nx y∞ ).
On choisit donc la loi de commande :
u = Nu yr − K (x − Nx yr ).
Expression finale du retour d’état
On trouve :
u = −Kx + N̄yr ,
avec N̄ = Nu + KNx .
y
yr + u dx
dt
= Ax + Bu
N̄
y = Cx + Du
−
x
K
Expression finale du retour d’état
Modèle d’état du système corrigé en boucle fermée :
dx
= (A − BK )x + B N̄yr ,
dt
y = (C − DK )x + D N̄yr .
yr dx
= (A − BK )x + B N̄yr y
dt
y = (C − DK )x + D N̄yr
Exemple de commande par retour d’état (1)
Asservissement de la vitesse d’un MCC par retour d’état

Hypothèse : il faut disposer de mesure de son état.
Rappel sur les modèles du MCC pour y = ω :

x1 = ω et x2 = i
! ! !
− Jf K

d x1 J x1 0
= + u,
L

x1
y = 1 0
x2
dω
x1 = ω et x2 = dt
! ! !
d x1 0 1 x1 0
= 2 + u,
dt x2 − Rf LJ
+K
− RJ+Lf x2 K
LJ LJ

x1
y = 1 0
x2

Hypothèse : il faut disposer de mesure de son état.
Rappel sur les modèles du MCC pour y = ω :

x1 = ω et x2 = i
dω
x1 = ω et x2 = dt
Technologie :
mesure de la vitesse angulaire du rotor
mesure du courant d’induit = réalisation du variateur
mesure de l’accélération angulaire du rotor

y =ω
yr = ωr + u dx
dt
= Ax + Bu
N̄
y = Cx + Du
− −
x2 = i
K2
x1 = ω
K1
Retour d’état : u = N̄ωr − K1 ω − K2 i.

Asservissement de la vitesse d’un MCC par boucles

imbriquées
Y (s) = Ω(s)
Yr (s) = Ωr (s) + + U(s)
Cω (s) Ci (s) G(s)
− −
I(s)
Ω(s)
Boucles imbriquées : U(s) = Cω (s)Ci (s)Ωr (s) − Cω (s)Ci (s)Ω(s) − Ci (s)I(s).

Similitudes
Similitudes
Avec :
K1
Cω (s) = et Ci (s) = K2
K2
on a le retour d’état U(s) = K1 Ωr (s) − K1 Ω(s) − K2 I(s).
Similitudes
Avec :
K1
Cω (s) = et Ci (s) = K2
K2
Similitudes
Avec :
K1
Cω (s) = et Ci (s) = K2
K2
Asservissement de position du MCC

Même analogie entre asservissement de position par retour
d’état et avec boucles imbriquées de courant, de position et
tachymétrique (voir TP sur le MCC).
Plan


Stabilité

Cahier des charges
Correcteurs PID
Plan


Stabilité

Cahier des charges
Correcteurs PID
Stabilité
Généralités
Notions et critères de stabilité définis dans le cas généralrstent
valables dans le cas d’un système asservi :
calcul des pôles du système
critère de Routh appliqué à la FTBF
Particularités
On dispose de plusieurs autres outils et résultats pour étudier
la stabilité d’un système en boucle fermée à partir du modèle
de la boucle ouverte.
Critère de Nyquist (1)
Principe
Critère de Nyquist : basé sur le tracé dans le diagramme de
Nyquist du lieu de la réponse harmonique en boucle ouverte ou
lieu de Nyquist.
Controur de Nyquist
Préalable au tracé du lieu de Nyquist : définition du contour de
Nyquist CN .
Axe imaginaire Axe imaginaire

+j∞ +j∞
jω1
+∞ +∞
0 Axe réel 0 Axe réel
−jω1
CN CN
Principe du critère
Stabilité en BF : dépend de la manière dont le lieu de Nyquist
entoure le point critique.
Critère de Nyquist
Soit un système possèdant P pôles à partie réelle > 0. Soit N
le nombre de tours fait par le lieu de Nyquist autour du point
critique, comptés positivement dans le sens trigonométrique.
Alors, le nombre de pôles instables du système en BF est :
Z = P − N.
Précaution
Attention : les tours sont comptés dans le sens trigo alors que
le CN est orienté dans le sens anti-trigo.
Corollaire
Système est stable en BF si le lieu fait P tours autour du point
critique.
Cas du MCC
On s’intéresse à la position θ du rotor. Par intégration de la
fonction de transfert en vitesse :
Θ(s) KG
= .
U(s) s (1 + τel s) (1 + τem s)
Cas du MCC : hypothèses
Modèle d’ordre deux :
Θ(s) KG
G(s) = =
U(s) s (1 + τem s)
et correction proportionnelle C(s) = Kp avec un capteur de

position assimilé à un gain pur Kθ .
Yr (s) Θ(s)
Kp G(s)
+
−
Y (s)
Kθ
Axe imaginaire
Paramétrage du contour de Nyquist :
+∞
partie du contour à partie imaginaire
positive
partie entourant l’origine : s = ρejα ,
+∞ avec α croissant de 0 à π2 et ρ tendant
0 Axe réel vers 0.
points de l’axe imaginaire : s = jω avec
ω ∈ [ρ, +∞[
CN fermeture du contour : s = ρejα avec α
décroissant de π2 à 0 et ρ tendant vers
l’infini.
Axe imaginaire
+∞
positive
0 Axe réel vers 0.
ω ∈ [ρ, +∞[
l’infini.
Axe imaginaire
+∞
positive
0 Axe réel vers 0.
ω ∈ [ρ, +∞[
l’infini.
Axe imaginaire
+∞
positive
0 Axe réel vers 0.
ω ∈ [ρ, +∞[
l’infini.
KCGH
CGH(s) = s(1+τem s) avec KCGH = Kp KG Kθ
Lieu de Nyquist
KCGH
Lieu de Nyquist
Contour : partie du contour de Nyquist entourant l’origine : s = ρejα , avec α croissant
de 0 à π2 et ρ tendant vers 0
Lieu, image du contour :
KCGH
CGH(s) = = r ejβ ,
ρejα
KCGH
avec r = ρ
tendant vers +∞ et β = −α décroissant de 0 à − π2 .
KCGH
Lieu de Nyquist
Contour : points de l’axe imaginaire s = jω avec ω ∈ [ρ, +∞[
Lieu image du contour :
KCGH
CGH(jω) = .
jω (1 + jτem ω)
Pas de particularité remarquable : tracé à partir du calcul de la partie réelle et de la
partie imaginaire de CGH(jω) lorsque ω varie de 0+ à +∞ :
KCGH τem KCGH

Re (CGH(jω)) = − 2 ω2
, et Im (CGH(jω)) = − 2 ω2
.
1 + τem ω 1 + τem
KCGH
Lieu de Nyquist
π
Contour : fermeture du contour s = ρejα avec α décroissant de 2
à 0 et ρ tendant vers
l’infini.
Lieu, image du contour :
KCGH
CGH(s) = ,
τem ρ2 e2jα
réduit à l’origine car dans ce cas ρ tend vers +∞.
Lieu de Nyquist :
Première partie :
CGH(s) = r ejβ ,
Axe imaginaire K
avec r = CGH tendant vers +∞ et β = −α
ω → 0− ρ
décroissant de 0 à − π2 .
Deuxième partie : lorsque ω varie de 0+ à +∞ :
ω → −∞ Axe réel KCGH τem

ω → +∞ Re (CGH(jω)) = − ,
−1 +∞ 2 ω2
1 + τem
KCGH
Im (CGH(jω)) = − 2 ω2
.
ω 1 + τem
N
ω → 0+
Troisième partie : réduite à l’origine.
Axe imaginaire
ω → 0−
ω → −∞ Axe réel
−1 ω → +∞ +∞
N
ω → 0+
Application du critère de Nyquist

Pas de pôle instable : P = 0. Par ailleurs N = 0.
L’asservissement en position d’un MCC à l’aide d’un correcteur
purement proportionnel est toujours stable, ce qui est normal
puisque il s’agit d’un système du second ordre.
Attention
Si on ne néglige plus le pôle électrique du système, pour Kp
assez grand. . .
Axe imaginaire
ω → 0−
−1 ω → +∞ Axe réel
ω → −∞ +∞
N
ω → 0+
Attention
Si on ne néglige plus le pôle électrique du système, pour Kp
assez grand : N = −2 alors que l’on a toujours P = 0.
Alors Z = 2 : le système en BF possède une paire de pôles
instables (pôles complexes conjugués).
Axe imaginaire
ω → 0−
−1 ω → +∞ Axe réel
ω → −∞ +∞
N
ω → 0+
Critère du revers (1)
Hypothèses
Version simplifiée du critère de Nyquist, quand le système est à
minimum de phase (pôles et zéros à partie réelle négative
uniquement et pas de terme de retard).
Critère du revers
Soit un système à minimum de phase.
Le système est stable en BF si en parcourant le lieu de
transfert de la boucle ouverte selon les ω croissants on laisse le
point critique à gauche dans le plan de Nyquist (ou à droite
dans le plan de Black).
Critère du revers (2)
Cas du MCC
Système à phase minimale : on applique le critère du revers.
Axe imaginaire Axe imaginaire
Axe réel Axe réel

−1 −1
CGH(jω) CGH(jω)
Cas toujours stable Cas instable

(modèle simplifié 2nd ordre) (grand gain proportionnel
et pôle électrique pris en compte)
Marges de stabilité (1)
Marge de phase
Mϕ = 180 + Arg{CGH(jωc )},

où ωc est telle que |CGH(jωc )|dB = 0 dB.
Marge de gain
1
MG = ,
|CGH(jω−180 )|
ou :
MGdB = −20 log10 |CGH(jω−180 )|
avec ω−180 tel que Arg{CGH(jω−180 )} = −180 deg.
Théorème
Le système est stable en boucle fermée si la marge de phase
(la marge de gain) du système en boucle ouverte est positive.
Stabilité relative
La valeur de la marge de phase permet donc de donner une
image de la stabilité du système, comme l’amortissement.
1
0.9
0.8
0.7
Coefficient d’amortissement
0.6
0.5
0.4
0.3
0.2
0.1
10 20 30 40 50 60 70 80
Marge de phase (deg)
Coefficient d’amortissement ξ vs marge de phase Mϕ

Marges de stabilité et diagrammes harmoniques
Les marges de stabilité peuvent être observées indifféremment
dans l’un ou l’autre des diagrammes harmoniques.
Axe imaginaire
−1 ω−180 1/MG 0 Axe réel

G(ω−180 )
Mϕ ϕ(ωc )
ωc
CGH(jω)
GdB
ωc ω−180 ω
MGdB
GdB (ω−180 )
ϕ(deg )
ωc ω−180 ω
−90
ϕ(ωc )
Mϕ
−180
Lieu des racines (1)
Définition
Lieu des racines (lieu d’Evans) : lieu des pôles de la fonction de
transfert en BF lorsque le gain Kp de la chaîne directe varie de
0 à +∞.
Yr (s) + E(s) U(s) Y (s)

Kp G(s)
−
H(s)
Lieu des racines et stabilité

Système stable en BF : pôles dans le demi-plan complexe
gauche, pour Kp donné.
Lieu des racines : construction

Soit m le degré du numérateur de la FTBO du système et n son
ordre. Soit N(s) le numérateur de la FTBO du système et D(s)
son dénominateur. Les propriétés suivante du lieu des racines,
permettent sa construction :

Règle n. 1 : Symétrie par rapport à l’axe réel ;


Règle n. 2 : n branches ;
n points de départ, pour Kp = 0, confondus avec les pôles
de la FTBO ;
m points d’arrivée, pour Kp = +∞, confondus avec les
zéros de la FTBO.

Règle n. 3 : n − m branches infinies :

ces branches donnent n − m asymptotes faisant des
angles βλ = 2λ+1n−m π, avec λ = 0, 1 . . . , n − m − 1 avec
l’horizontale ;
intersection P
des asymptotes avec l’axe réel au point
n Pm
i=1 pi − i=1 zi
d’affixe δ = n−m .

Règle n. 4 : Branches du lieu appartenant à l’axe réel : un point

d’affixe réelle appartient au lieu si le nombre de pôles et zéros
réels de la FTBO situés à sa droite est impair ;

Règle n. 5 : Intersections du lieu avec l’axe réel : un point

d’affixe réelle x est un point potentiel de séparation si :

1 dN(s) 1 dD(s)
= .
N(s) ds s=x D(s) ds s=x
Lieu des racines d’un système du 2nd ordre, en boucle fermé
On s’intéresse au lieu des racines du procédé de fonction de
transfert :
1
G2 (s) =
(s − p1 )(s − p2 )
en boucle fermée, le retour étant unitaire.
Lieu des racines
1
0.8
0.6
0.4
0.2
Im
−0.2
−0.4
−0.6
−0.8
−1
−5 −4.5 −4 −3.5 −3 −2.5 −2 −1.5 −1 −0.5 0
Re
Lieu des racines d’un système du 3ème ordre, en boucle fermé
On s’intéresse au lieu des racines du procédé de fonction de
transfert :
1
G3 (s) =
s(s − p1 )(s − p2 )
en boucle fermée, le retour étant unitaire.
Lieu des racines
10
2
Im
−2
−4
−6
−8
−10
−14 −12 −10 −8 −6 −4 −2 0 2 4
Re
Plan


Stabilité

Cahier des charges
Correcteurs PID
Précision (1)
Précision statique
La précision statique du système est caractérisée par l’erreur
en régime permanent en réponse à un échelon.
Cette erreur est appelée erreur statique (ou erreur de position).
D’après le théorème de la valeur finale :
ε∞ = lim ε(t) = lim sE(s).

t→+∞ s→0
Précision dynamique
On parlera de précision dynamique dès que l’entrée du
système évolue de manière continue dans le temps : par
exemple on désigne par erreur de vitesse la valeur de l’erreur
quand l’entrée du système est une rampe.
Précision (2)
Sans perturbation ni bruit de mesure, l’expression de l’erreur
est :
E(s) = Yr (s) − H(s)Y (s),
et donc :
CGH(s)Yr (s)
E(s) = Yr (s) − .
1 + CGH(s)
Finalement :
Yr (s)
E(s) = .
1 + CGH(s)
Précision (3)
L’expression générale de l’erreur d’un système asservi est :
sYr (s)
ε∞ = lim .
s→0 1 + CGH(s)
Avec :
K 1 + β1 s + . . .
CGH(s) =
sc 1 + α1 s + . . .
On établit le tableau suivant :
E0 V0
entrée échelon Yr = rampe Yr =
s s2
E0
classe 0 ∞
1+K
V0
classe 1 0
K
classe 2 0 0
Précision (4)
Dualité stabilité-précision
20 log10 Kp
(Kp > 1)
ωc ωc0 > ωc ω
Kp GdB (ω)
GdB (ω)
ϕ(ω)
1 1
τ1 τ2 ω
−90
Mϕ
Mϕ0 < Mϕ
−180
Plan


Stabilité

Cahier des charges
Correcteurs PID
Plan


Stabilité

Cahier des charges
Correcteurs PID
Commande : généralités (1)
Problème posé ?
La commande u(t) d’un système à temps continu peut être
modifiée en asservissant ce système à l’aide d’un correcteur.
Les paramètres de réglages sont :
la forme du correcteur :
série
boucle interne
retour d’état
etc.
les paramètres du correcteur :
réglage des actions proportionnelle, intégrale et dérivée
réglage d’un correcteur par retour d’état
Cahier des charges

Problème posé sous forme d’un cahier des charges (faisable) =
ensemble de contraintes à satisfaire
Exemples :
dépassement maximal
précision statique ou dynamique
temps de réponse
bande passante
etc.
Caractéristiques classiques d’un bon asservissement

Stabilité :
système stable
amortissement maîtrisé
Mϕ ' 50 à 70 deg
ξ ' 0, 5 à 0, 8
Précision : selon la classe ou le gain du système (le cas

échéant ajout d’une ou plusieurs intégrations)
Rapidité : rapidité et bande passante sont liées (pulsation de
coupure élevée, système rapide)
Plan


Stabilité

Cahier des charges
Correcteurs PID
Correcteurs PID : introduction
Forme idéale du PID

Yr (s) + E(s) U(s) Y (s)
1
C(s) = Kp 1+ + τd s G(s)
τi s
−
H(s)

Yr (s) + E(s) U(s) Y (s)
1
C(s) = Kp 1+ + τd s G(s)
τi s
−
H(s)
Action proportionnelle :
précision améliorée / stabilité diminuée
temps de montée réduit et plus de dépassement
temps de réponse pas forcément diminué

Yr (s) + E(s) U(s) Y (s)
1
C(s) = Kp 1+ + τd s G(s)
τi s
−
H(s)
Action intégrale :
classe du système augmentée : précision améliorée
marge de phase diminuée de 90 deg par l’ajout d’une
intégration pure
saturation : dispositif d’anti-saturation

Yr (s) + E(s) U(s) Y (s)
1
C(s) = Kp 1+ + τd s G(s)
τi s
−
H(s)
Action dérivée :
augmentation de la bande passante du système
augmentation de la stabilité, à bande passante égale
correcteur dérivé non causal : correction approchée, par
avance de phase

Yr (s) + E(s) U(s) Y (s)
1
C(s) = Kp 1+ + τd s G(s)
τi s
−
H(s)
Actions combinées :
P = proportionnel
PI = proportionnel et intégral (ou retard de phase)
PD = proportionnel et dérivé (ou avance de phase)
PID = proportionnel, intégral et dérivé
Adéquation correcteurs/systèmes à asservir (1)
Correcteur PI

1 Kp (1 + τi s)
C(s) = Kp 1 + = .
τi s τi s
Avantages :
intégration : convient donc bien lorsque l’on souhaite
annuler l’erreur statique d’un système de classe 0
correcteur le plus utilisé
Inconvénient :
action intégrale : saturation éventuelle de la commande
Correcteur PI
1+τi s
C(s) = Kp τi s
CdB (ω)
PI
20 log10 Kp
ϕ(ω)
1
τi ω
PI
−90
Correcteur à retard de phase
1+τi s
C(s) = aKp 1+aτi s , avec a > 1
CdB (ω)
20 log10 aKp
retard de phase ω
20 log10 Kp
ϕ(ω)
1 1
aτi τi ω
retard de phase
ϕm
−90
Comparaison correcteur PI et correcteur à retard de phase
1+τi s
C(s) = aKp 1+aτi s , avec a > 1 → PI approché, si a 1
CdB (ω)
20 log10 aKp
retard de phase ω
PI
20 log10 Kp
ϕ(ω)
1 1
aτi τi ω
retard de phase
PI
ϕm
−90
Correcteur à avance de phase et correcteur PD

1+τd s
C(s) = Kp 1+aτd s (a < 1) et C(s) = Kp (1 + τd s)
Avantages :
augmentation de la phase dans une certaine bande de
fréquence
système corrigé plus stable : convient donc bien pour la
correction des systèmes peu stables (systèmes de classe
supérieure ou égale à un)
Inconvénient :
PD idéal non causal : non réalisable physiquement
Correcteur PD
C(s) = Kp (1 + τd s)
CdB (ω) PD idéal
20 log10 Kp ω
ϕ(ω)
90 PD idéal
ω
1
τd
Correcteur à avance de phase
1+τd s
C(s) = Kp 1+aτd s , avec a < 1
CdB (ω)

Kp
20 log10 a
avance de phase
20 log10 Kp ω
ϕ(ω)
90
ϕM
avance de phase
ω
1 1
τd aτd
Comparaison correcteur à avance de phase et PD
1+τd s
C(s) = Kp 1+aτd s , avec a < 1 → PD approché, si a 1
CdB (ω) PD idéal

Kp
20 log10 a
avance de phase
20 log10 Kp
ω
ϕ(ω)
90 PD idéal
ϕM
avance de phase
ω
1 1
τd aτd
Correcteur à avance de phase : caractéristiques

Avance de phase maximale telle que :
1−a 1
sin ϕM = à ωM = √ .
1+a aτd
Exemple : a = 0, 1 → ϕM = 54, 9 deg
Correcteur à avance de phase : réglage

Calcul de a :
1 − sin ϕM
a=
1 + sin ϕM
Calcul de τd :
1
τd = √ .
aωM
Correcteur à avance et retard de phase

1 + τd s 1 + τi s
C(s) = bKp ,
1 + aτd s 1 + bτi s
avec a < 1 et b > 1.
Intérêt
Ce correcteur est bien évidemment plus général que les
correcteurs précédents. Il a vocation à corriger des systèmes
plus délicats à régler. Il n’est cependant pas nécessaire
d’utiliser ce type de correcteur si le cahier des charges peut
être rempli par un des correcteurs précédemment évoqués.
Approche le correcteur PID idéal.
Plan


Stabilité

Cahier des charges
Correcteurs PID
Méthodes harmoniques : correction PI
Correction PI
Réglage d’un correcteur PI :

1 Kp (1 + τi s)
C(s) = Kp 1 + = .
τi s τi s
en utilisant le diagramme de Bode.
Hypothèse
Système de classe zéro avec deux pôles réels, associés aux
constantes de temps τ1 et τ2 , τ1 τ2 .
GHdB (ω)
1
τ2 ω
1
τ1
ϕ(ω)
1 1
τ1 τ2 ω
−90
−180
CGHdB (ω) pour Kp = 1
(−
1)
GHdB (ω)
1
τ2 ω
1
τ1
(−
2)
ϕ(ω)
1 1
τ1 τ2 ω
Arg (CGH(ω))
−90
−180
CGHdB (ω)
(−
CGHdB (ω) pour Kp = 1 1)
GHdB (ω)
1
ωc τ2 ω
1
20 log Kp τ1
(−
2)
ϕ(ω)
1 1
τ1 ωc τ2 ω
Arg (CGH(ω))
−90
Mϕ désirée
Mϕ
−180
Synthèse de correcteur PI
1
G(s) =
(s + 1)(s + 5)
Cahier des charges :

erreur statique nulle
coefficient d’amortissement de 0, 7
Méthodes harmoniques : correction à avance de
phase
Correction à avance de phase

Réglage d’un correcteur PD approché :
1 + τd s
C(s) = Kp
1 + aτd s
en utilisant le diagramme de Bode.
Hypothèse
Système est de classe 1 avec un pôle réel, associé à la
constante de temps τ1 .
(−
1)
(−
2)
.
ϕ(ω)
GHdB (ω)
1
τ1 ω
−90
−180
(−
1)
(−
1)
(−
2)
ϕ(ω)
GHdB (ω)
1 1 1
τ1 τd
ωc aτd ω
−90
Mϕ désirée Arg (CGH(ω))
Mϕ
−180
(−
1)
(−
2)
(−
1) ω
20 log Kp
(−
2)
CGHdB (ω)
ϕ(ω)
GHdB (ω)
1 1 1
τ1 τd
ωc aτd ω
−90
Mϕ désirée Arg (CGH(ω))
−180
Synthèse de correcteur à avance de phase
150
G(s) =
s(s + 5)

marge de phase de 60 deg
bande passante 10 rad/s
Plan


Stabilité

Cahier des charges
Correcteurs PID
Méthode du lieu des racines (1)
Interprétation pour la correction des systèmes

Lieu des racines = lieu des pôles de la FTBO, lorsque le gain
de la chaîne directe varie.
Yr (s) + E(s) U(s) Y (s)

Kp CG(s)
−
H(s)
Synthèse par placement de pôles et zéros

Problème de synthèse = obtention des pôles et zéros conférant
au système le comportement dynamique souhaité,
généralement de type second ordre :
choix de pôles et zéros compensant les pôles présents
dans la FTBO
choix de pôles et zéros imposant la nouvelle dynamique du
système
réglage du gain
Axe imaginaire
ωn
p1 p
ωn 1 − ξ2
cos(ψ) = ξ ψ 0 Axe réel

−ξωn
p
-ωn 1 − ξ2
p2
Utilisation des courbes caractéristiques pour un système

équivalent ordre 2 pseudo-oscillant :
courbes de même amortissement
courbes de même temps de réponse
courbes de même bande passante
Lieu des racines 25
25
0.46 0.34 0.24 0.17 0.11 0.05
20
20
0.62
15
15
10
10 0.84
5
5
Axe imaginaire
−5
5
−10 0.84
10
−15
15
0.62
−20
20
0.46 0.34 0.24 0.17 0.11 0.05
−25
−15 −10 −5 25 0
Axe réel
Principe de compensation : précautions !

Tous les systèmes ne peuvent pas nécessairement être
transformés en système du second ordre.

Certaines compensations doivent être évitées :
s − (pi + ∆pi ) ∆pi

=1− .
s − pi s − pi

s − (pi + ∆pi ) ∆pi

=1− .
s − pi s − pi
on ne compense jamais un pôle à partie réelle positive


s − (pi + ∆pi ) ∆pi

=1− .
s − pi s − pi
on ne compense jamais un pôle à partie réelle positive

on ne compense pas autant que possible une paire de
pôles très oscillants
Synthèse de correcteur PI
1
G(s) =
(s + 1)(s + 5)

coefficient d’amortissement de 0, 7
Synthèse de correcteur à avance de phase
150
G(s) =
s(s + 5)

marge de phase de 60 deg
bande passante 10 rad/s
Méthodes empiriques de synthèse
Méthodes de Ziegler-Nichols
Essais :
identification de la réponse indicielle
identification du régime critique
Synthèse sous forme de tableaux.

K
pente α = τ
y (t)
K
t
tr tr + τ
coefficients du correcteur
type de correcteur Kp τi τd
1
P
αtr
0, 9 tr
PI
αtr 0, 3
1, 2
PID 2tr 0, 5tr
αtr
Plan


Stabilité

Cahier des charges
Correcteurs PID
Commande par retour d’état (1)
y
yr + u dx
dt
= Ax + Bu
N̄
y = Cx + Du
−
x
Commande par retour d’état

Représentation d’état en boucle fermée :
dx
= (A − BK )x + B N̄yr ,
dt
y = (C − DK )x + D N̄yr .
Commande par retour d’état

Représentation d’état en boucle fermée :
dx
= (A − BK )x + B N̄yr ,
dt
y = (C − DK )x + D N̄yr .
det (sI − (A − BK )) = 0 : degré n paramétrée par les n gains
K1 , K2 , . . . , Kn .
Problème : placement de n pôles p1 , p2 , . . . , pn donnant
l’équation caractéristique : (s − p1 )(s − p2 ) . . . (s − pn ) = 0.
Commandabilité et retour d’état

Si un système linéaire invariant est commandable, alors il
existe toujours un choix de gains permettant de placer
l’ensemble des pôles du système par un retour d’état.
Par ailleurs, dans le cas d’un système mono-entrée
mono-sortie, ce choix est unique.
Commande du MCC par retour d’état

Avec x1 = ω et x2 = i, représentation d’état du MCC :
! ! !
− Jf K

d x1 J x 1 0
= + 1 u,
dt x2 − KL − RL x2 L

x1
y = 1 0 .
x2
Soit :
f K

s 0 −J J
0
sI − (A − BK ) = − − 1 K1 K2 ,
0 s − KL − RL L
ou :
f
− KJ

s+ J
sI − (A − BK ) = K +K1
L s + R+K L
2
Equation caractéristique sI − (A − BK ) :

2 f R + K2 f (R + K2 ) + K (K + K1 )
s + + s+ = 0.
J L LJ
Identification des pôles à une paire de pôles complexes

conjugués :
(s − p1 )(s − p2 ) = s2 + 2ξωn s + ωn2 = 0.


2 f R + K2 f (R + K2 ) + K (K + K1 )
s + + s+ = 0.
J L LJ

conjugués :
(s − p1 )(s − p2 ) = s2 + 2ξωn s + ωn2 = 0.
On déduit :
f R + K2
+ = 2ξωn ,
J L
f (R + K2 ) + K (K + K1 )
= ωn2 ,
LJ

2 f R + K2 f (R + K2 ) + K (K + K1 )
s + + s+ = 0.
J L LJ

conjugués :
(s − p1 )(s − p2 ) = s2 + 2ξωn s + ωn2 = 0.
soit les gains du retour d’état :
f
K2 = L(2ξωn − ) − R,
J
LJωn2 − f (R + K2 )
K1 = − K.
K
Remarque
Dans le cas précédent l’identification reste faisable
analytiquement car le système n’est pas d’ordre trop élevé.
Toutefois, au delà du second ordre, la résolution devient vite
fastidieuse.
Problème
Trouver une méthode systématique pour un système d’ordre
quelconque.
Modèle dans sa forme canonique commandable :
0 1 0 ... ... ... 0
 
 0 0 1 0 ... ... 0  0
 
 ... ... ... ... ... ... ...  0
 
dx
=  ... ... ... ... ... ... . . .  x + . . . u.
   
dt 0 
 ... ... ... ... ... ...
  . . . 
 0 0 ... ... ... 0 1  1
−a0 −a1 ... ... ... ... −an−1
Retour d’état :
0 1 0 ... ... ... 0
 
 0 0 1 0 ... ... 0 
... ... ... ... ... ... ...
 
 
A − BK =  ... ... ... ... ... ... ... .
 
... ... ... ... ... ... 0
 
 
 0 0 ... ... ... 0 1 
−a0 − K1 −a1 − K2 ... ... ... ... −an−1 − Kn
Equation caractéristique det(sI − (A − BK )) :
sn + (an−1 + Kn )sn−1 + (an−2 + Kn−1 )sn−2 + . . . + a0 + K1 = 0.
Identification avec l’équation caractéristique sous forme

développée : sn + αn−1 sn−1 + αn−2 sn−2 + . . . + α0 = 0 :
Kn = αn−1 − an−1 ,
Kn−1 = αn−2 − an−2 ,
...
K1 = α0 − a0 .
Synthèse de correcteur par retour d’état, exemple (2)
! ! !
d x1 0 1 x1 0
= + u,
dt x2 −ω02 0 x2 1

x1
y = 1 0 .
x2

coefficient d’amortissement de 1
pulsation naturelle deux fois supérieure à la boucle ouverte
Remarque
Résolution possible car la représentation d’état du système est
sous forme canonique commandable.
Généralisation
Si ce n’est pas le cas, on recherche le changement de variable
x = Pz transformant le couple (A, B) en (Ac , Bc ), sous forme
canonique commandable.
Une fois obtenue la matrice de gain Kc du retour d’état
u = −Kc z, si l’on souhaite conserver la représentation d’état
initiale, on transformera la matrice de gain pour déduire les
gains nécessaires à la commande du système dans son état
initial : K = Kc P −1 .
Formule d’Ackerman
Permet de déterminer la matrice de gain du retour d’état :
calcul du changement de variable pour mettre le système
dans sa forme canonique commandable
calcul de la matrice de gain
transformation de cette matrice pour appliquer le retour
d’état dans la représentation initiale
Formule d’Ackerman
Soit C la matrice de commandabilité associée à la paire (A, B).
Alors, la matrice de gain :
K = 0 0 . . . 0 1 C−1 Pα (A),

où Pα (A) = An + αn−1 An−1 + . . . + α0 I est consitué à partir de

l’équation caractéristique issue du placement de pôles.
Limitations
Intérêt pratique/faible robustesse : non validité en cas de
mauvais conditionnement de la matrice de commandabilité.
Système d’ordre élevé : formule d’Ackerman inapplicable.
Synthèse de correcteur par retour d’état, exemple (2)
! ! !
d x1 0 1 x1 0
= + u,
dt x2 −ω02 0 x2 1

x1
y = 1 0 .
x2

coefficient d’amortissement de 1
pulsation naturelle deux fois supérieure à la boucle ouverte

Automatique Continue: Bernard BAYLE

Transféré par

Droits d'auteur :

Formats disponibles

Automatique Continue: Bernard BAYLE

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Automatique Continue: Bernard BAYLE

Transféré par

Droits d'auteur :

Formats disponibles

Automatique Continue

Télécom Physique Strasbourg

Evaluation (documents autorisés)

Clepsydre - Ktesibios d’Alexandrie (270 avant JC) - . . .

Machine de Watt 1769

Systèmes à temps continu

2 Réponses des systèmes à temps continu

3 Analyse des systèmes à temps continu

2 Réponses des systèmes à temps continu

3 Analyse des systèmes à temps continu

2 Réponses des systèmes à temps continu

3 Analyse des systèmes à temps continu

Relation entrée u - sortie ω avec (Ke = Km = K ) :

Relation entrée u - sortie ω avec (Ke = Km = K ) :

Relation entrée u - sortie ω avec (Ke = Km = K ) :

Relation entrée u - sortie ω avec (Ke = Km = K ) :

Relation entrée u - sortie ω avec (Ke = Km = K ) :

Relation entrée u - sortie ω avec (Ke = Km = K ) :

Variations −kr (zr − zc )

Section : Cdh = (qe − qs )dt

Relation entrée qe - sortie h :

Relation entrée qe - sortie h :

Relation entrée qe - sortie h :

Relation entrée qe - sortie h :

(R : résistance thermique moyenne de l’échangeur)

(Cv : capacité thermique de la vapeur)

Chaleur amenée : qv = wv cv (θve − θv )

Chaleur amenée : qv = wv cv (θve − θv )

Chaleur amenée : qv = wv cv (θve − θv )

Chaleur amenée : qv = wv cv (θve − θv )

Chaleur amenée : qv = wv cv (θve − θv )

2 Réponses des systèmes à temps continu

3 Analyse des systèmes à temps continu

Le système est linéaire si sa sortie vaut αy1 (t) + y2 (t) en

linéarité αu1 (t) + u2 (t) αy1 (t) + y2 (t)

Le système est invariant si une même commande, appliquée à

u(t) y (t) invariance u(t + τ ) y (t + τ )

Le système est causal si, ∀t < 0, u(t) = 0 ⇒ y (t) = 0.

2 Réponses des systèmes à temps continu

3 Analyse des systèmes à temps continu

x de dimension n × 1 : vecteur d’état, représentation

x de dimension n × 1 : vecteur d’état, représentation

x de dimension n × 1 : vecteur d’état, représentation

x de dimension n × 1 : vecteur d’état, représentation

x de dimension n × 1 : vecteur d’état, représentation

x de dimension n × 1 : vecteur d’état, représentation

x de dimension n × 1 : vecteur d’état, représentation

Réécriture des équations :

Réécriture des équations :

Réécriture des équations :

Réécriture des équations :

Réécriture des équations :

Réécriture des équations :

2 Réponses des systèmes à temps continu

3 Analyse des systèmes à temps continu

Convergence en valeur absolue de l’intégrale :

Converge en valeur absolue ssi σ + a > 0, donc σ > σ0 = −a.

Transformée de e−at U(t)

linéarité L{f (t) + αg(t)} = F (s) + αG(s), ∀α ∈ R