Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf ·...
Transcript of Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf ·...
![Page 1: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/1.jpg)
Regulation et Asservissement
Yassine Ariba
Dpt GEI - Icam, Toulouse.
version 1.0
Y. Ariba - Icam, Toulouse. Automatique Lineaire 1 / 131
![Page 2: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/2.jpg)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 2 / 131
![Page 3: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/3.jpg)
Sommaire
1 Modelisation
Exemples introductifs
Notion de systeme
Fonction de transfert
Schema fonctionnel2 Reponse d’un systeme
Reponse temporelle
Reponse frequentielle
Resume3 L’Automatique
Definition
Applications
Architecture de commande
Notion de stabilite
Objectif du cours
4 Stabilite et performances d’un
asservissement
Stabilite
Precision
Rapidite
Marges de stabilite5 Synthese de correcteurs
Introduction
Synthese directe
Action proportionnelle
Action integrale
Action derivee
Combinaisons d’actions
Proportionnel-Integral-Derive
(PID)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 3 / 131
![Page 4: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/4.jpg)
Modelisation Exemples introductifs
Qu’est ce que la modelisation ?
Dans toutes les disciplines, l’etude d’un processus physique/dispositif/mecanismenecessite d’etablir un modele :
pour comprendre et analyser le dispositif,
pour pouvoir predire son comportement,
pour utiliser des outils de simulation.
B un modele est toujours une representation approchee de la realite.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 4 / 131
![Page 5: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/5.jpg)
Modelisation Exemples introductifs
Qu’est ce que la modelisation ?
Dans toutes les disciplines, l’etude d’un processus physique/dispositif/mecanismenecessite d’etablir un modele :
pour comprendre et analyser le dispositif,
pour pouvoir predire son comportement,
pour utiliser des outils de simulation.
B un modele est toujours une representation approchee de la realite.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 4 / 131
![Page 6: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/6.jpg)
Modelisation Exemples introductifs
Qu’est ce que la modelisation ?
Dans toutes les disciplines, l’etude d’un processus physique/dispositif/mecanismenecessite d’etablir un modele :
pour comprendre et analyser le dispositif,
pour pouvoir predire son comportement,
pour utiliser des outils de simulation.
B un modele est toujours une representation approchee de la realite.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 4 / 131
![Page 7: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/7.jpg)
Modelisation Exemples introductifs
Exemple : le circuit RC
Considerons un circuit electronique
Appliquons la loi des mailles
e(t) = Ri(t) + v(t) sachant que i(t) = Cdv(t)
dt.
Nous avons le modele :RCv(t) + v(t) = e(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 5 / 131
![Page 8: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/8.jpg)
Modelisation Exemples introductifs
Exemple : le circuit RC
Considerons un circuit electronique
Appliquons la loi des mailles
e(t) = Ri(t) + v(t) sachant que i(t) = Cdv(t)
dt.
Nous avons le modele :RCv(t) + v(t) = e(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 5 / 131
![Page 9: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/9.jpg)
Modelisation Exemples introductifs
Exemple : le circuit RC
Considerons un circuit electronique
Appliquons la loi des mailles
e(t) = Ri(t) + v(t) sachant que i(t) = Cdv(t)
dt.
Nous avons le modele :RCv(t) + v(t) = e(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 5 / 131
![Page 10: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/10.jpg)
Modelisation Exemples introductifs
Exemple : structure masse-ressort
Considerons une structure mecanique
Appliquons le principe fondamental de la dynamique
mx(t) = f(t) + fressort(t) + famort(t).
Nous avons le modele :
mx(t) + cx(t) + kx(t) = f(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 6 / 131
![Page 11: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/11.jpg)
Modelisation Exemples introductifs
Exemple : structure masse-ressort
Considerons une structure mecanique
Appliquons le principe fondamental de la dynamique
mx(t) = f(t) + fressort(t) + famort(t).
Nous avons le modele :
mx(t) + cx(t) + kx(t) = f(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 6 / 131
![Page 12: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/12.jpg)
Modelisation Exemples introductifs
Exemple : structure masse-ressort
Considerons une structure mecanique
Appliquons le principe fondamental de la dynamique
mx(t) = f(t) + fressort(t) + famort(t).
Nous avons le modele :
mx(t) + cx(t) + kx(t) = f(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 6 / 131
![Page 13: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/13.jpg)
Modelisation Exemples introductifs
Exemple : base mobile
Considerons un robot a deux roues motrices independantes
On montre que les vitesses de translation et rotation du mobile sont
v = rωr + ωl
2et ω = r
ωr − ωl2l
Nous avons le modele : x = v cos θy = v sin θ
θ = ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 7 / 131
![Page 14: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/14.jpg)
Modelisation Exemples introductifs
Exemple : base mobile
Considerons un robot a deux roues motrices independantes
On montre que les vitesses de translation et rotation du mobile sont
v = rωr + ωl
2et ω = r
ωr − ωl2l
Nous avons le modele : x = v cos θy = v sin θ
θ = ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 7 / 131
![Page 15: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/15.jpg)
Modelisation Exemples introductifs
Exemple : base mobile
Considerons un robot a deux roues motrices independantes
On montre que les vitesses de translation et rotation du mobile sont
v = rωr + ωl
2et ω = r
ωr − ωl2l
Nous avons le modele : x = v cos θy = v sin θ
θ = ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 7 / 131
![Page 16: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/16.jpg)
Modelisation Exemples introductifs
Exemple : echanges thermiques dans une enceinte
Considerons les variations de temperature a l’interieur d’une enceinte
La deperdition et le bilan energetique s’ecrivent
Ppertes = k(T − Text) et P − Ppertes = cdT
dt
Nous avons le modele :P = k(T − Text) + cT
Y. Ariba - Icam, Toulouse. Automatique Lineaire 8 / 131
![Page 17: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/17.jpg)
Modelisation Exemples introductifs
Exemple : echanges thermiques dans une enceinte
Considerons les variations de temperature a l’interieur d’une enceinte
La deperdition et le bilan energetique s’ecrivent
Ppertes = k(T − Text) et P − Ppertes = cdT
dt
Nous avons le modele :P = k(T − Text) + cT
Y. Ariba - Icam, Toulouse. Automatique Lineaire 8 / 131
![Page 18: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/18.jpg)
Modelisation Exemples introductifs
Exemple : echanges thermiques dans une enceinte
Considerons les variations de temperature a l’interieur d’une enceinte
La deperdition et le bilan energetique s’ecrivent
Ppertes = k(T − Text) et P − Ppertes = cdT
dt
Nous avons le modele :P = k(T − Text) + cT
Y. Ariba - Icam, Toulouse. Automatique Lineaire 8 / 131
![Page 19: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/19.jpg)
Modelisation Exemples introductifs
Exemple : dynamique des populations
Considerons l’evolution de la taille N d’une population
Premier modele :
N(t) = αN(t)− βN(t)
avec
α : taux de naissances
β : taux de deces
Second modele :
N(t) = rN(t)
(1−
N(t)
K
)avec
r = α− βK : nombre max d’individus
0temps t
taill
e po
pula
tion
N(t
)
α > β
N0
N0
β > α
0temps t
taill
e po
pula
tion
N(t
)
K
Y. Ariba - Icam, Toulouse. Automatique Lineaire 9 / 131
![Page 20: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/20.jpg)
Modelisation Exemples introductifs
Exemple : dynamique des populations
Considerons l’evolution de la taille N d’une population
Premier modele :
N(t) = αN(t)− βN(t)
avec
α : taux de naissances
β : taux de deces
Second modele :
N(t) = rN(t)
(1−
N(t)
K
)avec
r = α− βK : nombre max d’individus
0temps t
taill
e po
pula
tion
N(t
)
α > β
N0
N0
β > α
0temps t
taill
e po
pula
tion
N(t
)
K
Y. Ariba - Icam, Toulouse. Automatique Lineaire 9 / 131
![Page 21: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/21.jpg)
Modelisation Exemples introductifs
Exemple : dynamique des populations
Considerons l’evolution de la taille N d’une population
Premier modele :
N(t) = αN(t)− βN(t)
avec
α : taux de naissances
β : taux de deces
Second modele :
N(t) = rN(t)
(1−
N(t)
K
)avec
r = α− βK : nombre max d’individus
0temps t
taill
e po
pula
tion
N(t
)
α > β
N0
N0
β > α
0temps t
taill
e po
pula
tion
N(t
)
K
Y. Ariba - Icam, Toulouse. Automatique Lineaire 9 / 131
![Page 22: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/22.jpg)
Modelisation Exemples introductifs
Exemple : dynamique des populations
Considerons l’evolution de la taille N d’une population
Premier modele :
N(t) = αN(t)− βN(t)
avec
α : taux de naissances
β : taux de deces
Second modele :
N(t) = rN(t)
(1−
N(t)
K
)avec
r = α− βK : nombre max d’individus
0temps t
taill
e po
pula
tion
N(t
)
α > β
N0
N0
β > α
0temps t
taill
e po
pula
tion
N(t
)
K
Y. Ariba - Icam, Toulouse. Automatique Lineaire 9 / 131
![Page 23: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/23.jpg)
Modelisation Exemples introductifs
Exemple : dynamique des populations
Considerons l’evolution de la taille N d’une population
Premier modele :
N(t) = αN(t)− βN(t)
avec
α : taux de naissances
β : taux de deces
Second modele :
N(t) = rN(t)
(1−
N(t)
K
)avec
r = α− βK : nombre max d’individus
0temps t
taill
e po
pula
tion
N(t
)
α > β
N0
N0
β > α
0temps t
taill
e po
pula
tion
N(t
)
K
Y. Ariba - Icam, Toulouse. Automatique Lineaire 9 / 131
![Page 24: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/24.jpg)
Modelisation Notion de systeme
Modeles mathematiques
Ici, les modeles sont des equations differentielles
RCv(t) + v(t) = e(t)
⇒ lineaire
mx(t) + cx(t) + kx(t) = f(t)
⇒ lineaire
x = v cos θy = v sin θ
θ = ω
P = k(T − Text) + cT
⇒ lineaire
N(t) = rN(t)
(1− N(t)
K
)
? meme formalisme mathematique.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 10 / 131
![Page 25: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/25.jpg)
Modelisation Notion de systeme
Modeles mathematiques
Ici, les modeles sont des equations differentielles
RCv(t) + v(t) = e(t) ⇒ lineaire
mx(t) + cx(t) + kx(t) = f(t) ⇒ lineaire
x = v cos θy = v sin θ
θ = ω
P = k(T − Text) + cT ⇒ lineaire
N(t) = rN(t)
(1− N(t)
K
)
? meme formalisme mathematique.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 10 / 131
![Page 26: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/26.jpg)
Modelisation Notion de systeme
Modeles mathematiques
Ici, les modeles sont des equations differentielles
RCv(t) + v(t) = e(t) ⇒ lineaire
mx(t) + cx(t) + kx(t) = f(t) ⇒ lineaire
x = v cos θy = v sin θ
θ = ω
P = k(T − Text) + cT ⇒ lineaire
N(t) = rN(t)
(1− N(t)
K
)
? meme formalisme mathematique.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 10 / 131
![Page 27: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/27.jpg)
Modelisation Notion de systeme
Approche systeme
L’automaticien adopte une approche generique
l’entree - u(t) - represente le signal de commande qui va agir sur le systeme.
la sortie - y(t) - caracterise le comportement du systeme. Il est la reponse dusysteme a une excitation exterieure.
la perturbation - d(t) - represente les signaux d’entree incontroles affectant lesysteme. Il specifie les interactions avec son environnement.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 11 / 131
![Page 28: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/28.jpg)
Modelisation Notion de systeme
Approche systeme
L’automaticien adopte une approche generique
l’entree - u(t) - represente le signal de commande qui va agir sur le systeme.
la sortie - y(t) - caracterise le comportement du systeme. Il est la reponse dusysteme a une excitation exterieure.
la perturbation - d(t) - represente les signaux d’entree incontroles affectant lesysteme. Il specifie les interactions avec son environnement.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 11 / 131
![Page 29: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/29.jpg)
Modelisation Notion de systeme
Approche systeme
Meme formalisme pour chaque cas :
Y. Ariba - Icam, Toulouse. Automatique Lineaire 12 / 131
![Page 30: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/30.jpg)
Modelisation Notion de systeme
Les systemes linaires invariants
Dans de ce cours, on s’interesse exclusivement aux SLI
Relation entree-sortie : u(t)→ y(t)
⇒ equations differentielles lineairesa coefficients constants
Y. Ariba - Icam, Toulouse. Automatique Lineaire 13 / 131
![Page 31: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/31.jpg)
Modelisation Notion de systeme
Les systemes linaires invariants
Dans de ce cours, on s’interesse exclusivement aux SLI
Relation entree-sortie : u(t)→ y(t)
⇒ equations differentielles lineairesa coefficients constants
Y. Ariba - Icam, Toulouse. Automatique Lineaire 13 / 131
![Page 32: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/32.jpg)
Modelisation Fonction de transfert
La transformee de Laplace
La transformee de Laplace (TL) d’une fonction f(t) est la fonction du nombrecomplexe s :
f(s) =
∫ ∞0
e−stf(t)dt
On note : f(s) = Lf(t)
et f(t) = L−1
f(s)
.
Proprietes :
Linearite : Laf(t) + bg(t)
= af(s) + bg(s), a et b constants.
Derivation : Lddtf(t)
= sf(s)− f(0).
Integration : L∫ t
0f(θ)dθ
=
1
sf(s).
Valeur finale : limt→∞
f(t) = lims→0
sf(s).
Y. Ariba - Icam, Toulouse. Automatique Lineaire 14 / 131
![Page 33: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/33.jpg)
Modelisation Fonction de transfert
La transformee de Laplace
La transformee de Laplace (TL) d’une fonction f(t) est la fonction du nombrecomplexe s :
f(s) =
∫ ∞0
e−stf(t)dt
On note : f(s) = Lf(t)
et f(t) = L−1
f(s)
.
Proprietes :
Linearite : Laf(t) + bg(t)
= af(s) + bg(s), a et b constants.
Derivation : Lddtf(t)
= sf(s)− f(0).
Integration : L∫ t
0f(θ)dθ
=
1
sf(s).
Valeur finale : limt→∞
f(t) = lims→0
sf(s).
Y. Ariba - Icam, Toulouse. Automatique Lineaire 14 / 131
![Page 34: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/34.jpg)
Modelisation Fonction de transfert
Table des transformees
Rappelons que les fonctions temporelles f(t) ci-dessous ne sont definies que pour t ≥ 0.
Fonction Dom. temporel f(t) Trans. de Laplacef(s)
echelon 11
s
rampe t1
s2
puissance n-iemetn
n!
1
sn+1
exponentielledecroissante
e−at 1
s+ a
Y. Ariba - Icam, Toulouse. Automatique Lineaire 15 / 131
![Page 35: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/35.jpg)
Modelisation Fonction de transfert
Table des transformees
Rappelons que les fonctions temporelles f(t) ci-dessous ne sont definies que pour t ≥ 0.
Fonction Dom. temporel f(t) Trans. de Laplacef(s)
sinus sinωtω
s2 + ω2
cosinus cosωts
s2 + ω2
decroissanceexponentielle d’un
sinus
e−at
sin(bt)b
(s+ a)2 + b2
decroissanceexponentielle d’un
cosinus
e−at
cos(bt)s+ a
(s+ a)2 + b2
decroissanceexponentielle d’unepuissance n-ieme
tn
n!e−at 1
(s+ a)n+1
Y. Ariba - Icam, Toulouse. Automatique Lineaire 16 / 131
![Page 36: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/36.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
La TL permet de determiner la fonction de transfert d’un systeme.
5y(3)(t) + 2y(t) + y(t) + 4y(t) = u(t) + 2u(t)
y L
5s3y(s) + 2s2y(s) + sy(s) + 4y(s) = su(s) + 2u(s)yG(s) =
y(s)
u(s)=
s+ 2
5s3 + 2s2 + s+ 4
Nouvelle relation entre la sortie et l’entree
y(s) = G(s)u(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 17 / 131
![Page 37: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/37.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
La TL permet de determiner la fonction de transfert d’un systeme.
5y(3)(t) + 2y(t) + y(t) + 4y(t) = u(t) + 2u(t)y L
5s3y(s) + 2s2y(s) + sy(s) + 4y(s) = su(s) + 2u(s)yG(s) =
y(s)
u(s)=
s+ 2
5s3 + 2s2 + s+ 4
Nouvelle relation entre la sortie et l’entree
y(s) = G(s)u(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 17 / 131
![Page 38: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/38.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
La TL permet de determiner la fonction de transfert d’un systeme.
5y(3)(t) + 2y(t) + y(t) + 4y(t) = u(t) + 2u(t)y L
5s3y(s) + 2s2y(s) + sy(s) + 4y(s) = su(s) + 2u(s)yG(s) =
y(s)
u(s)=
s+ 2
5s3 + 2s2 + s+ 4
Nouvelle relation entre la sortie et l’entree
y(s) = G(s)u(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 17 / 131
![Page 39: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/39.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
Cas general :
any(n)(t) + · · ·+ a1y(t) + a0y(t) = bmu(m)(t) + · · ·+ b1u(t) + b0u(t)
y L
ansny(s) + · · ·+ a1sy(s) + a0y(s) = bmsmu(s) + · · ·+ b1su(s) + b0u(s)
yFonction de transfert : G(s) =
y(s)
u(s)=bmsm + · · ·+ b1s+ b0
ansn + · · ·+ a1s+ a0
La fonction de transfert est definie comme le rapport entre la TL de la sortie surla TL de l’entree, pour des conditions initiales nulles.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 18 / 131
![Page 40: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/40.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
Cas general :
any(n)(t) + · · ·+ a1y(t) + a0y(t) = bmu(m)(t) + · · ·+ b1u(t) + b0u(t)
y L
ansny(s) + · · ·+ a1sy(s) + a0y(s) = bmsmu(s) + · · ·+ b1su(s) + b0u(s)
yFonction de transfert : G(s) =
y(s)
u(s)=bmsm + · · ·+ b1s+ b0
ansn + · · ·+ a1s+ a0
La fonction de transfert est definie comme le rapport entre la TL de la sortie surla TL de l’entree, pour des conditions initiales nulles.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 18 / 131
![Page 41: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/41.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
Cas general :
any(n)(t) + · · ·+ a1y(t) + a0y(t) = bmu(m)(t) + · · ·+ b1u(t) + b0u(t)
y L
ansny(s) + · · ·+ a1sy(s) + a0y(s) = bmsmu(s) + · · ·+ b1su(s) + b0u(s)
yFonction de transfert : G(s) =
y(s)
u(s)=bmsm + · · ·+ b1s+ b0
ansn + · · ·+ a1s+ a0
La fonction de transfert est definie comme le rapport entre la TL de la sortie surla TL de l’entree, pour des conditions initiales nulles.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 18 / 131
![Page 42: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/42.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
Cas general :
any(n)(t) + · · ·+ a1y(t) + a0y(t) = bmu(m)(t) + · · ·+ b1u(t) + b0u(t)
y L
ansny(s) + · · ·+ a1sy(s) + a0y(s) = bmsmu(s) + · · ·+ b1su(s) + b0u(s)
yFonction de transfert : G(s) =
y(s)
u(s)=bmsm + · · ·+ b1s+ b0
ansn + · · ·+ a1s+ a0
La fonction de transfert est definie comme le rapport entre la TL de la sortie surla TL de l’entree, pour des conditions initiales nulles.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 18 / 131
![Page 43: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/43.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
Autre exemple :
y(t) + 4y(t)− 2y(t) = 6u(t)
y L
s2y(s) + 4sy(s)− 2y(s) = 6u(s)yG(s) =
y(s)
u(s)=
6
s2 + 4s− 2
Y. Ariba - Icam, Toulouse. Automatique Lineaire 19 / 131
![Page 44: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/44.jpg)
Modelisation Fonction de transfert
Fonctions de transfert
Autre exemple :
y(t) + 4y(t)− 2y(t) = 6u(t)y L
s2y(s) + 4sy(s)− 2y(s) = 6u(s)yG(s) =
y(s)
u(s)=
6
s2 + 4s− 2
Y. Ariba - Icam, Toulouse. Automatique Lineaire 19 / 131
![Page 45: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/45.jpg)
Modelisation Fonction de transfert
Definitions :
L’ordre d’un systeme est le degre du polynome du denominateur (n).
Le systeme est dit strictement propre (propre), si n > m (n = m).
Les poles sont les racines du denominateur.
Les zeros sont les racines du numerateur.
Exemple : systeme masse-ressort
La relation entre la force et la position de la masse est decrite par l’equation
mx(t) + cx(t) + kx(t) = f(t).
Sa fonction de transfert s’ecrit
x(s)
f(s)=
1
ms2 + cs+ k.
C’est un systeme : d’ordre deux ; sans zero ; les 2 poles sont solutions dems2 + cs+ k = 0.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 20 / 131
![Page 46: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/46.jpg)
Modelisation Fonction de transfert
Definitions :
L’ordre d’un systeme est le degre du polynome du denominateur (n).
Le systeme est dit strictement propre (propre), si n > m (n = m).
Les poles sont les racines du denominateur.
Les zeros sont les racines du numerateur.
Exemple : systeme masse-ressort
La relation entre la force et la position de la masse est decrite par l’equation
mx(t) + cx(t) + kx(t) = f(t).
Sa fonction de transfert s’ecrit
x(s)
f(s)=
1
ms2 + cs+ k.
C’est un systeme : d’ordre deux ; sans zero ; les 2 poles sont solutions dems2 + cs+ k = 0.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 20 / 131
![Page 47: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/47.jpg)
Modelisation Schema fonctionnel
Schema fonctionnel
Representation graphique des systemes complexes
Le formalisme des fonctions de transfert est tres pratique pour ce type derepresentation.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 21 / 131
![Page 48: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/48.jpg)
Modelisation Schema fonctionnel
Schema fonctionnel
Mise en serie de deux systemes
y2(s) = G2(s)y1(s)
= G2(s)G1(s)u1(s)
Transfert equivalent :
F (s) = G2(s)×G1(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 22 / 131
![Page 49: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/49.jpg)
Modelisation Schema fonctionnel
Schema fonctionnel
Mise en serie de deux systemes
y2(s) = G2(s)y1(s)
= G2(s)G1(s)u1(s)
Transfert equivalent :
F (s) = G2(s)×G1(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 22 / 131
![Page 50: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/50.jpg)
Modelisation Schema fonctionnel
Schema fonctionnel
Mise en parallele de deux systemes
y(s) = y1(s) + y2(s)
=(G1(s) +G2(s)
)u(s)
Transfert equivalent :
F (s) = G1(s) +G2(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 23 / 131
![Page 51: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/51.jpg)
Modelisation Schema fonctionnel
Schema fonctionnel
Mise en parallele de deux systemes
y(s) = y1(s) + y2(s)
=(G1(s) +G2(s)
)u(s)
Transfert equivalent :
F (s) = G1(s) +G2(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 23 / 131
![Page 52: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/52.jpg)
Modelisation Schema fonctionnel
Schema fonctionnel
Interconnexion feedback
y(s) = G1(s)u(s)
u(s) = e(s)−G2(s)y(s)
Transfert equivalent :
F (s) =G1(s)
1 +G1(s)G2(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 24 / 131
![Page 53: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/53.jpg)
Modelisation Schema fonctionnel
Schema fonctionnel
Interconnexion feedback
y(s) = G1(s)u(s)
u(s) = e(s)−G2(s)y(s)
Transfert equivalent :
F (s) =G1(s)
1 +G1(s)G2(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 24 / 131
![Page 54: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/54.jpg)
Modelisation Schema fonctionnel
Exemple 1
+
+
Transfert equivalent :
y(s) = G3(s)(G1(s) +G2(s)
)︸ ︷︷ ︸
F (s)
u(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 25 / 131
![Page 55: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/55.jpg)
Modelisation Schema fonctionnel
Exemple 1
+
+
Transfert equivalent :
y(s) = G3(s)(G1(s) +G2(s)
)︸ ︷︷ ︸
F (s)
u(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 25 / 131
![Page 56: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/56.jpg)
Modelisation Schema fonctionnel
Exemple 2
Transfert equivalent :
y(s) = F1(s)u1(s) + F2(s)u2(s)
avecF1(s) = G1(s) et F2(s) =
(G2(s)−G1(s)
)G3(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 26 / 131
![Page 57: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/57.jpg)
Modelisation Schema fonctionnel
Exemple 2
Transfert equivalent :
y(s) = F1(s)u1(s) + F2(s)u2(s)
avecF1(s) = G1(s) et F2(s) =
(G2(s)−G1(s)
)G3(s)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 26 / 131
![Page 58: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/58.jpg)
Reponse d’un systeme
Sommaire
1 Modelisation
Exemples introductifs
Notion de systeme
Fonction de transfert
Schema fonctionnel2 Reponse d’un systeme
Reponse temporelle
Reponse frequentielle
Resume3 L’Automatique
Definition
Applications
Architecture de commande
Notion de stabilite
Objectif du cours
4 Stabilite et performances d’un
asservissement
Stabilite
Precision
Rapidite
Marges de stabilite5 Synthese de correcteurs
Introduction
Synthese directe
Action proportionnelle
Action integrale
Action derivee
Combinaisons d’actions
Proportionnel-Integral-Derive
(PID)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 27 / 131
![Page 59: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/59.jpg)
Reponse d’un systeme Reponse temporelle
Reponse temporelle
Calcul explicite de la sortie y(t) a une entree e(t) donnee.
Systèmey(t)e(t)
Resolution avec les fonctions de transfert :
1 Exprimer la sortie y(s) = G(s)e(s).
2 Effectuer une decomposition en elements simples de y(s).
3 Appliquer la transformee de Laplace inverse (a l’aide de la table)
y(s)L−1
−−−→ y(t)
pour obtenir le signal dans le domaine temporel.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 28 / 131
![Page 60: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/60.jpg)
Reponse d’un systeme Reponse temporelle
Reponse temporelle
Calcul explicite de la sortie y(t) a une entree e(t) donnee.
Systèmey(t)e(t)
Resolution avec les fonctions de transfert :
1 Exprimer la sortie y(s) = G(s)e(s).
2 Effectuer une decomposition en elements simples de y(s).
3 Appliquer la transformee de Laplace inverse (a l’aide de la table)
y(s)L−1
−−−→ y(t)
pour obtenir le signal dans le domaine temporel.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 28 / 131
![Page 61: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/61.jpg)
Reponse d’un systeme Reponse temporelle
Exemple : calculons la reponse a un echelon unite du systeme suivant
Exprimons le signal de sortie
y(s) =2s+ 1
(s+ 1)(s+ 5)e(s) =
2s+ 1
(s+ 1)(s+ 5)
1
s
La reponse se decompose en
y(s) =A
s+
B
s+ 1+
C
s+ 5
Le calcul des coefficients donne :A = sy(s)
∣∣s=0
= 15
B = (s+ 1)y(s)∣∣s=−1
= 14
C = (s+ 5)y(s)∣∣s=−5
= − 920
Y. Ariba - Icam, Toulouse. Automatique Lineaire 29 / 131
![Page 62: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/62.jpg)
Reponse d’un systeme Reponse temporelle
Exemple : calculons la reponse a un echelon unite du systeme suivant
Exprimons le signal de sortie
y(s) =2s+ 1
(s+ 1)(s+ 5)e(s) =
2s+ 1
(s+ 1)(s+ 5)
1
s
La reponse se decompose en
y(s) =A
s+
B
s+ 1+
C
s+ 5
Le calcul des coefficients donne :A = sy(s)
∣∣s=0
= 15
B = (s+ 1)y(s)∣∣s=−1
= 14
C = (s+ 5)y(s)∣∣s=−5
= − 920
Y. Ariba - Icam, Toulouse. Automatique Lineaire 29 / 131
![Page 63: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/63.jpg)
Reponse d’un systeme Reponse temporelle
Exemple : calculons la reponse a un echelon unite du systeme suivant
Exprimons le signal de sortie
y(s) =2s+ 1
(s+ 1)(s+ 5)e(s) =
2s+ 1
(s+ 1)(s+ 5)
1
s
La reponse se decompose en
y(s) =A
s+
B
s+ 1+
C
s+ 5
Le calcul des coefficients donne :A = sy(s)
∣∣s=0
= 15
B = (s+ 1)y(s)∣∣s=−1
= 14
C = (s+ 5)y(s)∣∣s=−5
= − 920
Y. Ariba - Icam, Toulouse. Automatique Lineaire 29 / 131
![Page 64: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/64.jpg)
Reponse d’un systeme Reponse temporelle
Exemple : calculons la reponse a un echelon unite du systeme suivant
Exprimons le signal de sortie
y(s) =2s+ 1
(s+ 1)(s+ 5)e(s) =
2s+ 1
(s+ 1)(s+ 5)
1
s
La reponse se decompose en
y(s) =A
s+
B
s+ 1+
C
s+ 5
Le calcul des coefficients donne :A = sy(s)
∣∣s=0
= 15
B = (s+ 1)y(s)∣∣s=−1
= 14
C = (s+ 5)y(s)∣∣s=−5
= − 920
Y. Ariba - Icam, Toulouse. Automatique Lineaire 29 / 131
![Page 65: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/65.jpg)
Reponse d’un systeme Reponse temporelle
A l’aide de la table, nous obtenons la reponse temporelle
y(t) =1
5+
1
4e−t −
9
20e−5t
−1 0 1 2 3 4 5 6
0
0.2
0.4
0.6
0.8
1
Time (s)
Am
plitu
de
e(t)
y(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 30 / 131
![Page 66: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/66.jpg)
Reponse d’un systeme Reponse temporelle
Caracteristiques temporelles
Certaines caracteristiques du systeme sont releves a partir de la reponse indicielle.(≡ reponse a un echelon)
t
y(t)
y(∞)
±5% y(∞)
90% y(∞)
10% y(∞)
tr5%tm tp
D1
y(∞)D1
tmtr5%
tp
: valeur finale: 1er dépassement: temps de montée: temps de réponse: temps du pic
F Ces caracteristiques donnent des indices de performances.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 31 / 131
![Page 67: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/67.jpg)
Reponse d’un systeme Reponse temporelle
Cas des systemes du 1er ordre
modele canonique :
y(s) =k
τs+ 1e(s)
pour lequel on definit τ la constante de temps et k le gain statique.
Reponse a un echelon d’amplitude e0 : y(t) = e0k(
1− e−1τt)
t
y(t)
y(∞)= k e0
95% y(∞)
tr5%=3ττ
63% y(∞)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 32 / 131
![Page 68: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/68.jpg)
Reponse d’un systeme Reponse temporelle
Cas des systemes du 1er ordre
modele canonique :
y(s) =k
τs+ 1e(s)
pour lequel on definit τ la constante de temps et k le gain statique.
Reponse a un echelon d’amplitude e0 : y(t) = e0k(
1− e−1τt)
t
y(t)
y(∞)= k e0
95% y(∞)
tr5%=3ττ
63% y(∞)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 32 / 131
![Page 69: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/69.jpg)
Reponse d’un systeme Reponse temporelle
Cas des systemes du 2nd ordre
modele canonique :
y(s) =Kω2
n
s2 + 2ζωns+ ω2n
e(s)
ou l’on definit ζ le coefficient d’amortissement, ωn la pulsation propre et K le gainstatique. 3 cas pour la reponse indicielle :
cas ζ > 1 : regime aperiodique
y(t) = Ke0
[1 + p2
p1−p2e−ζωnt + p1
p2−p1e−ζωnt
]
cas ζ = 1 : regime critique
y(t) = Ke0
[1−
(1− p1t
)e−ζωnt
]
cas 0 < ζ < 1 : regime pseudo-periodique (ωp = ωn√
1− ζ2)
y(t) = Ke0
[1− e−ζωnt
(cos(ωpt) +
ζ√1− ζ2
sin(ωpt))]
Y. Ariba - Icam, Toulouse. Automatique Lineaire 33 / 131
![Page 70: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/70.jpg)
Reponse d’un systeme Reponse temporelle
Cas des systemes du 2nd ordre
modele canonique :
y(s) =Kω2
n
s2 + 2ζωns+ ω2n
e(s)
ou l’on definit ζ le coefficient d’amortissement, ωn la pulsation propre et K le gainstatique. 3 cas pour la reponse indicielle :
cas ζ > 1 : regime aperiodique
y(t) = Ke0
[1 + p2
p1−p2e−ζωnt + p1
p2−p1e−ζωnt
]
cas ζ = 1 : regime critique
y(t) = Ke0
[1−
(1− p1t
)e−ζωnt
]
cas 0 < ζ < 1 : regime pseudo-periodique (ωp = ωn√
1− ζ2)
y(t) = Ke0
[1− e−ζωnt
(cos(ωpt) +
ζ√1− ζ2
sin(ωpt))]
Y. Ariba - Icam, Toulouse. Automatique Lineaire 33 / 131
![Page 71: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/71.jpg)
Reponse d’un systeme Reponse temporelle
0 5 10 15 20 25 30 35 40 450
0.2
0.4
0.6
0.8
1
1.2
ζ = 4
Ke0=1ζ = 1
Cas ζ ≥ 1
(ici : e0 = 1, K = 1 et ωn = 1)
pas de depassement
quand ζ ou ωn , les temps dereponse et de montee
0 5 10 15 20 25 300
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
ζ=0.1
ζ=0.8Ke
0=1
ωn=1
Cas 0 < ζ < 1
(ici : e0 = 1, K = 1 et ωn = 1)
depassement : D1 = 100e− ζπ√
1−ζ2
depassement a t =π
ωp
temps de reponse : tr5% '3
ζωn
Y. Ariba - Icam, Toulouse. Automatique Lineaire 34 / 131
![Page 72: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/72.jpg)
Reponse d’un systeme Reponse temporelle
0 5 10 15 20 25 30 35 40 450
0.2
0.4
0.6
0.8
1
1.2
ζ = 4
Ke0=1ζ = 1
Cas ζ ≥ 1
(ici : e0 = 1, K = 1 et ωn = 1)
pas de depassement
quand ζ ou ωn , les temps dereponse et de montee
0 5 10 15 20 25 300
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
ζ=0.1
ζ=0.8Ke
0=1
ωn=1
Cas 0 < ζ < 1
(ici : e0 = 1, K = 1 et ωn = 1)
depassement : D1 = 100e− ζπ√
1−ζ2
depassement a t =π
ωp
temps de reponse : tr5% '3
ζωn
Y. Ariba - Icam, Toulouse. Automatique Lineaire 34 / 131
![Page 73: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/73.jpg)
Reponse d’un systeme Reponse temporelle
Exemple : Soit le systeme d’ordre 2
y(s) =2
s2 + s+ 4u(s)
Quelle est l’allure de sa reponse indicielle ?
Identifions ses parametres caracteristiques :
ωn = 2 K = 0.5 ζ = 0.25
Nous pouvons donc conclure pour la reponse
regime pseudo-periodique,
valeur finale y(∞) = 0.5,
premier depassement de 44% a t = 1.62s,
temps de reponse tr5% ' 6s.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 35 / 131
![Page 74: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/74.jpg)
Reponse d’un systeme Reponse temporelle
Exemple : Soit le systeme d’ordre 2
y(s) =2
s2 + s+ 4u(s)
Quelle est l’allure de sa reponse indicielle ?
Identifions ses parametres caracteristiques :
ωn = 2 K = 0.5 ζ = 0.25
Nous pouvons donc conclure pour la reponse
regime pseudo-periodique,
valeur finale y(∞) = 0.5,
premier depassement de 44% a t = 1.62s,
temps de reponse tr5% ' 6s.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 35 / 131
![Page 75: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/75.jpg)
Reponse d’un systeme Reponse temporelle
Exemple : Soit le systeme d’ordre 2
y(s) =2
s2 + s+ 4u(s)
Quelle est l’allure de sa reponse indicielle ?
Identifions ses parametres caracteristiques :
ωn = 2 K = 0.5 ζ = 0.25
Nous pouvons donc conclure pour la reponse
regime pseudo-periodique,
valeur finale y(∞) = 0.5,
premier depassement de 44% a t = 1.62s,
temps de reponse tr5% ' 6s.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 35 / 131
![Page 76: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/76.jpg)
Reponse d’un systeme Reponse temporelle
Reponse indicielle :
0 2 4 6 8 10 120
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8Step Response
Time (seconds)
Am
plitu
de
Y. Ariba - Icam, Toulouse. Automatique Lineaire 36 / 131
![Page 77: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/77.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Reprenons le circuit RC
Regime sinusoıdale : e(t) = em cos(ωt+ φe)
v(t) = vm cos(ωt+ φv)⇒
e = emejφe
v = vmejφv
Y. Ariba - Icam, Toulouse. Automatique Lineaire 37 / 131
![Page 78: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/78.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Reprenons le circuit RC
Regime sinusoıdale : e(t) = em cos(ωt+ φe)
v(t) = vm cos(ωt+ φv)⇒
e = emejφe
v = vmejφv
Y. Ariba - Icam, Toulouse. Automatique Lineaire 37 / 131
![Page 79: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/79.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Pour R = 1kΩ et C = 200µF , appliquons une tension e(t) = cos(8t).
1 1.5 2 2.5 3 3.5 4 4.5 5
−1
−0.5
0
0.5
1
temps (s)
e(t)v(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 38 / 131
![Page 80: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/80.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Analysons le circuit. Loi d’Ohm : u = Zi
ZR = R et ZC =1
jωC
Appliquons le pont diviseur de tension
v =ZC
ZC + ZRe
La transmittance du circuit s’ecrit donc :
T =
1jωC
1jωC
+R=
1
jωRC + 1.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 39 / 131
![Page 81: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/81.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Analysons le circuit. Loi d’Ohm : u = Zi
ZR = R et ZC =1
jωC
Appliquons le pont diviseur de tension
v =ZC
ZC + ZRe
La transmittance du circuit s’ecrit donc :
T =
1jωC
1jωC
+R=
1
jωRC + 1.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 39 / 131
![Page 82: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/82.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Analysons le circuit. Loi d’Ohm : u = Zi
ZR = R et ZC =1
jωC
Appliquons le pont diviseur de tension
v =ZC
ZC + ZRe
La transmittance du circuit s’ecrit donc :
T =
1jωC
1jωC
+R=
1
jωRC + 1.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 39 / 131
![Page 83: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/83.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Tracons le diagramme de Bode de la transmittance
−30
−25
−20
−15
−10
−5
0
5M
agni
tude
(dB
)
10−1 100 101 102−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 40 / 131
![Page 84: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/84.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 1
Tracons le diagramme de Bode de la transmittance
−30
−25
−20
−15
−10
−5
0
5M
agni
tude
(dB
)
10−1 100 101 102−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
ω = 8
Y. Ariba - Icam, Toulouse. Automatique Lineaire 40 / 131
![Page 85: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/85.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses du circuit pour ω = 0.8, 4, 8, 40
10 15 20 25 30 35 40 45
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=0.8)v(t)
2 3 4 5 6 7 8 9 10
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=4)v(t)
3 3.5 4 4.5 5 5.5 6 6.5 7
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=8)v(t)
10 10.1 10.2 10.3 10.4 10.5 10.6 10.7 10.8
−1
−0.5
0
0.5
1
e(t) (ω=40)v(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 41 / 131
![Page 86: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/86.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses du circuit pour ω = 0.8, 4, 8, 40
10 15 20 25 30 35 40 45
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=0.8)v(t)
2 3 4 5 6 7 8 9 10
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=4)v(t)
3 3.5 4 4.5 5 5.5 6 6.5 7
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=8)v(t)
10 10.1 10.2 10.3 10.4 10.5 10.6 10.7 10.8
−1
−0.5
0
0.5
1
e(t) (ω=40)v(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 41 / 131
![Page 87: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/87.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses du circuit pour ω = 0.8, 4, 8, 40
10 15 20 25 30 35 40 45
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=0.8)v(t)
2 3 4 5 6 7 8 9 10
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=4)v(t)
3 3.5 4 4.5 5 5.5 6 6.5 7
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=8)v(t)
10 10.1 10.2 10.3 10.4 10.5 10.6 10.7 10.8
−1
−0.5
0
0.5
1
e(t) (ω=40)v(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 41 / 131
![Page 88: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/88.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses du circuit pour ω = 0.8, 4, 8, 40
10 15 20 25 30 35 40 45
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=0.8)v(t)
2 3 4 5 6 7 8 9 10
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=4)v(t)
3 3.5 4 4.5 5 5.5 6 6.5 7
−1
−0.5
0
0.5
1
temps (s)
e(t) (ω=8)v(t)
10 10.1 10.2 10.3 10.4 10.5 10.6 10.7 10.8
−1
−0.5
0
0.5
1
e(t) (ω=40)v(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 41 / 131
![Page 89: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/89.jpg)
Reponse d’un systeme Reponse frequentielle
−25
−20
−15
−10
−5
0
5M
agni
tude
(dB
)
10−1
100
101
102
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
ω=8ω=4ω=0.8
ω=40
Y. Ariba - Icam, Toulouse. Automatique Lineaire 42 / 131
![Page 90: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/90.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle : exemple introductif 2
Oscillations forcees d’un pendule
Y. Ariba - Icam, Toulouse. Automatique Lineaire 43 / 131
![Page 91: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/91.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses de la masse pour differentes oscillations du point Aω = 1, 6, 11, 50 (excitation moteur)
15 20 25 30 35 40
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=1)x(t)
10 11 12 13 14 15 16−4
−3
−2
−1
0
1
2
3
4
temps (s)
z(t) (ω=6)x(t)
20 20.5 21 21.5 22 22.5 23
−1
−0.5
0
0.5
1
z(t) (ω=11)x(t)
20 20.1 20.2 20.3 20.4 20.5 20.6
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=50)x(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 44 / 131
![Page 92: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/92.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses de la masse pour differentes oscillations du point Aω = 1, 6, 11, 50 (excitation moteur)
15 20 25 30 35 40
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=1)x(t)
10 11 12 13 14 15 16−4
−3
−2
−1
0
1
2
3
4
temps (s)
z(t) (ω=6)x(t)
20 20.5 21 21.5 22 22.5 23
−1
−0.5
0
0.5
1
z(t) (ω=11)x(t)
20 20.1 20.2 20.3 20.4 20.5 20.6
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=50)x(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 44 / 131
![Page 93: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/93.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses de la masse pour differentes oscillations du point Aω = 1, 6, 11, 50 (excitation moteur)
15 20 25 30 35 40
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=1)x(t)
10 11 12 13 14 15 16−4
−3
−2
−1
0
1
2
3
4
temps (s)
z(t) (ω=6)x(t)
20 20.5 21 21.5 22 22.5 23
−1
−0.5
0
0.5
1
z(t) (ω=11)x(t)
20 20.1 20.2 20.3 20.4 20.5 20.6
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=50)x(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 44 / 131
![Page 94: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/94.jpg)
Reponse d’un systeme Reponse frequentielle
Testons les reponses de la masse pour differentes oscillations du point Aω = 1, 6, 11, 50 (excitation moteur)
15 20 25 30 35 40
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=1)x(t)
10 11 12 13 14 15 16−4
−3
−2
−1
0
1
2
3
4
temps (s)
z(t) (ω=6)x(t)
20 20.5 21 21.5 22 22.5 23
−1
−0.5
0
0.5
1
z(t) (ω=11)x(t)
20 20.1 20.2 20.3 20.4 20.5 20.6
−1
−0.5
0
0.5
1
temps (s)
z(t) (ω=50)x(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 44 / 131
![Page 95: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/95.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle
L’analyse frequentielle consiste a etudier la reponse d’un systeme lineaire a des entreessinusoıdales.
F(s) Y(s)U(s)
u(t) = u0 sin(ωt) u(t) = u0 sin(ωt) y(t) = y0 sin(ωt+φ) y(t) = y0 sin(ωt+φ)
0 2 4 6 8 10 12 14 16 18 20
−1
−0.5
0
0.5
1
1.5
entrée
sortie
Le signal de sortie est sinusoıdale, de meme pulsation, d’amplitude differente etpresente un dephasage.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 45 / 131
![Page 96: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/96.jpg)
Reponse d’un systeme Reponse frequentielle
Reponse frequentielle
L’analyse frequentielle consiste a etudier la reponse d’un systeme lineaire a des entreessinusoıdales.
F(s) Y(s)U(s)
u(t) = u0 sin(ωt) u(t) = u0 sin(ωt) y(t) = y0 sin(ωt+φ) y(t) = y0 sin(ωt+φ)
0 2 4 6 8 10 12 14 16 18 20
−1
−0.5
0
0.5
1
1.5
entrée
sortie
Le signal de sortie est sinusoıdale, de meme pulsation, d’amplitude differente etpresente un dephasage.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 45 / 131
![Page 97: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/97.jpg)
Reponse d’un systeme Reponse frequentielle
0 2 4 6 8 10 12 14 16 18 20
−1
−0.5
0
0.5
1
1.5
u0
y0
y(t)
u(t)T
T
∆ t
La reponse frequentielle est caracterisee par
son amplification :y0
u0
son dephasage : ±360∆t
T
L’amplification et le dephasage induit par le systeme dependent de la pulsation ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 46 / 131
![Page 98: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/98.jpg)
Reponse d’un systeme Reponse frequentielle
Exemple : considerons le systeme
F (s) =1/2
s+ 1
Observons sa reponse aux 3 entrees :
u1 = sin(0.05 t)
u2 = sin(1.5 t)
u3 = sin(10 t)
0 50 100 150 200 250 300 350 400 450 500
−1
−0.5
0
0.5
1
u1(t)
y1(t)
0 2 4 6 8 10 12 14 16
−1
−0.5
0
0.5
1 u2(t)
y2(t)
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
−1
−0.5
0
0.5
1 u3(t)
y3(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 47 / 131
![Page 99: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/99.jpg)
Reponse d’un systeme Reponse frequentielle
Exemple : considerons le systeme
F (s) =1/2
s+ 1
Observons sa reponse aux 3 entrees :
u1 = sin(0.05 t)
u2 = sin(1.5 t)
u3 = sin(10 t)
0 50 100 150 200 250 300 350 400 450 500
−1
−0.5
0
0.5
1
u1(t)
y1(t)
0 2 4 6 8 10 12 14 16
−1
−0.5
0
0.5
1 u2(t)
y2(t)
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
−1
−0.5
0
0.5
1 u3(t)
y3(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 47 / 131
![Page 100: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/100.jpg)
Reponse d’un systeme Reponse frequentielle
On montre que :
l’amplification = |F (jω)|,
le dephasage = argF (jω).
F (jω) est obtenue en remplacant la variable de Laplace s par jω.
Reprenons l’exemple precedent : F (jω) =1/2
jω + 1
pour ω = 0.05 rad/s : |F (j0.05)| = 0.5 et argF (j0.05) = −2.86.
pour ω = 1.5 rad/s : |F (j1.5)| = 0.277 et argF (j1.5) = −56.3.
pour ω = 10 rad/s : |F (j10)| = 0.05 et argF (j10) = −84.3.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 48 / 131
![Page 101: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/101.jpg)
Reponse d’un systeme Reponse frequentielle
On montre que :
l’amplification = |F (jω)|,
le dephasage = argF (jω).
F (jω) est obtenue en remplacant la variable de Laplace s par jω.
Reprenons l’exemple precedent : F (jω) =1/2
jω + 1
pour ω = 0.05 rad/s : |F (j0.05)| = 0.5 et argF (j0.05) = −2.86.
pour ω = 1.5 rad/s : |F (j1.5)| = 0.277 et argF (j1.5) = −56.3.
pour ω = 10 rad/s : |F (j10)| = 0.05 et argF (j10) = −84.3.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 48 / 131
![Page 102: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/102.jpg)
Reponse d’un systeme Reponse frequentielle
F (jω) est une fonction complexe de la pulsation ω (soit la frequence f =ω
2π).
On distingue trois representations graphiques :
les diagrammes de Bode,
la representation de Nyquist,
la representation de Black-Nichols.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 49 / 131
![Page 103: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/103.jpg)
Reponse d’un systeme Reponse frequentielle
Representation dans le diagramme de Bode
Representation du transfert F (jω) en fonction de ω sur deux graphes :
le diagramme de gain qui represente le module de F (jω) en decibel⇒ |F |dB = 20log|F (jω)|.
le diagramme de phase qui represente l’argument de F (jω) en degre ou radian.
|F |dB
ω
arg(F)
0,1 1 10 100 1000
ω0,1 1 10 100 1000
Y. Ariba - Icam, Toulouse. Automatique Lineaire 50 / 131
![Page 104: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/104.jpg)
Reponse d’un systeme Reponse frequentielle
Representation dans le diagramme de Bode
Exemple : Reponse frequentielle de F (p) =1/2
p+ 1⇒ representation de
F (jω) =1/2
jω + 1
Bode Diagram
Frequency (rad/sec)
−50
−40
−30
−20
−10
0M
agni
tude
(dB
)
10−2 10−1 100 101 102−90
−45
0
Pha
se (
deg)
20log(0.277)
−56 deg
ω=1.5
Pour une sollicitation sinusoıdale de pulsation ω = 1.5 :
amplification de 0.277,
dephasage de −56 deg.Y. Ariba - Icam, Toulouse. Automatique Lineaire 51 / 131
![Page 105: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/105.jpg)
Reponse d’un systeme Reponse frequentielle
Representation dans le diagramme de Bode
Propriete : le diagramme d’une fonction donnee est obtenu a partir de la somme destraces elementaires.
Soit F (p) = F1(p)× F2(p)× F3(p),
gain → 20log|F (jω)| = 20log|F1(jω)|+ 20log|F2(jω)|+ 20log|F3(jω)|,
phase → arg(F (jω)
)= arg
(F1(jω)
)+ arg
(F2(jω)
)+ arg
(F3(jω)
).
Methodes :
trace asymptotique, analyse du transfert lorsque ω → 0 et +∞ ainsi qu’enquelques points particuliers,
trace complet, a l’aide de logiciel de calcul numerique tel que MatlabR©,
trace experimental, solliciter le procede avec des entrees sinusoıdales de differentespulsation ω. Mesurer en sortie l’amplification/attenuation et le dephasage dusignal a l’oscilloscope.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 52 / 131
![Page 106: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/106.jpg)
Reponse d’un systeme Reponse frequentielle
Representation dans le diagramme de Bode
Propriete : le diagramme d’une fonction donnee est obtenu a partir de la somme destraces elementaires.
Soit F (p) = F1(p)× F2(p)× F3(p),
gain → 20log|F (jω)| = 20log|F1(jω)|+ 20log|F2(jω)|+ 20log|F3(jω)|,
phase → arg(F (jω)
)= arg
(F1(jω)
)+ arg
(F2(jω)
)+ arg
(F3(jω)
).
Methodes :
trace asymptotique, analyse du transfert lorsque ω → 0 et +∞ ainsi qu’enquelques points particuliers,
trace complet, a l’aide de logiciel de calcul numerique tel que MatlabR©,
trace experimental, solliciter le procede avec des entrees sinusoıdales de differentespulsation ω. Mesurer en sortie l’amplification/attenuation et le dephasage dusignal a l’oscilloscope.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 52 / 131
![Page 107: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/107.jpg)
Reponse d’un systeme Reponse frequentielle
Traces elementaires
Rappel
Soit un nombre complexe z = a+ jb.
a0
Im
Re
b z
|z|θ
• z = |z|ejθ • |z| =√a2 + b2 • tan(θ) =
b
a
Y. Ariba - Icam, Toulouse. Automatique Lineaire 53 / 131
![Page 108: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/108.jpg)
Reponse d’un systeme Reponse frequentielle
Traces elementaires
Terme constant : F (p) = k, (k > 0)
Fonction de transfert complexe : F (jω) = k
|F (jω)|db = 20log(k)
φ = 0
|F(jω)|dB
ω0+90
-90
φ20log|k|
ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 54 / 131
![Page 109: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/109.jpg)
Reponse d’un systeme Reponse frequentielle
Traces elementaires
Terme constant : F (p) = k, (k > 0)
Fonction de transfert complexe : F (jω) = k
|F (jω)|db = 20log(k)
φ = 0
|F(jω)|dB
ω0+90
-90
φ20log|k|
ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 54 / 131
![Page 110: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/110.jpg)
Reponse d’un systeme Reponse frequentielle
Traces elementaires
Derivateur : F (p) = p
Fonction de transfert complexe : F (jω) = jω
|F (jω)|db = 20log(ω)
φ = +90
|F(jω)|dB
ω0+90
-90
φ+20
-201 10
0.1
ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 55 / 131
![Page 111: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/111.jpg)
Reponse d’un systeme Reponse frequentielle
Traces elementaires
Derivateur : F (p) = p
Fonction de transfert complexe : F (jω) = jω
|F (jω)|db = 20log(ω)
φ = +90
|F(jω)|dB
ω0+90
-90
φ+20
-201 10
0.1
ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 55 / 131
![Page 112: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/112.jpg)
Reponse d’un systeme Reponse frequentielle
Traces elementaires
Integrateur : F (p) =1
p
Fonction de transfert complexe : F (jω) =1
jω
|F (jω)|db = −20log(ω)
φ = arg(−j 1ω
) = −90
|F(jω)|dB
ω0+90
-90
φ+20
-201
10
0.1 ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 56 / 131
![Page 113: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/113.jpg)
Reponse d’un systeme Reponse frequentielle
Traces elementaires
Integrateur : F (p) =1
p
Fonction de transfert complexe : F (jω) =1
jω
|F (jω)|db = −20log(ω)
φ = arg(−j 1ω
) = −90
|F(jω)|dB
ω0+90
-90
φ+20
-201
10
0.1 ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 56 / 131
![Page 114: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/114.jpg)
Reponse d’un systeme Reponse frequentielle
Fonction du premier ordre
F (p) =1
1 + τp
Fonction de transfert complexe : F (jω) =1
1 + jτω
Module :
|F (jω)| =1
|1 + jτω|=
1√
1 + τ2ω2
Argument :
arg(F (jω)
)= arg
(1
1 + jτω
)= arg
(1)− arg
(1 + jτω
)= −arctan
τω
1
|F (jω)|db = −20log
(√1 + τ2ω2
)φ = −arctan
(τω)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 57 / 131
![Page 115: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/115.jpg)
Reponse d’un systeme Reponse frequentielle
Fonction du premier ordre
Courbe de gain : Courbe de phase :
• ω 1
τ|F |db → −20log(1) = 0
• ω 1
τ|F |db → −20log(τω)
(pente a −20dB/dec)
• ω =1
τ|F |db = −20log(
√2) w −3dB
• ω 1
τφ→ 0
• ω 1
τφ→ −
π
2
• ω =1
τφ = −arctan(1) = −
π
4
|F(jω)|dBω
0-45
-90
φ
-3
1/τ ω0
-20db/decade
1/τ
Y. Ariba - Icam, Toulouse. Automatique Lineaire 58 / 131
![Page 116: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/116.jpg)
Reponse d’un systeme Reponse frequentielle
Fonction du premier ordre
Courbe de gain : Courbe de phase :
• ω 1
τ|F |db → −20log(1) = 0
• ω 1
τ|F |db → −20log(τω)
(pente a −20dB/dec)
• ω =1
τ|F |db = −20log(
√2) w −3dB
• ω 1
τφ→ 0
• ω 1
τφ→ −
π
2
• ω =1
τφ = −arctan(1) = −
π
4
|F(jω)|dBω
0-45
-90
φ
-3
1/τ ω0
-20db/decade
1/τ
Y. Ariba - Icam, Toulouse. Automatique Lineaire 58 / 131
![Page 117: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/117.jpg)
Reponse d’un systeme Reponse frequentielle
Fonction du second ordre
F (p) =1
1 + 2ζωnp+ 1
ω2np2
Fonction de transfert complexe : F (jω) =1
(1− ω2
ω2n
) + 2ζ ωωnj
Module :
|F (jω)| =1∣∣∣(1− ω2
ω2n
) + 2ζ ωωnj∣∣∣ =
1√(1− ω2
ω2n
)2 + (2ζ ωωn
)2
Argument :
arg(F (jω)
)= −arg
((1−
ω2
ω2n
) + 2ζω
ωnj
)= −arctan
2ζ ωωn
1− ω2
ω2n
|F (jω)|db = −20log
√(1− ω2
ω2n
)2 + (2ζ ωωn
)2
φ = −arctan2ζ ωωn
1− ω2
ω2nY. Ariba - Icam, Toulouse. Automatique Lineaire 59 / 131
![Page 118: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/118.jpg)
Reponse d’un systeme Reponse frequentielle
Fonction du second ordre
Courbe de gain : Courbe de phase :
• ω ωn |F |db → −20log(1) = 0
• ω ωn |F |db → −40log( ωωn
)
(pente a −40dB/dec)
• ω = ωn |F |db = −20log(2ζ)
• ω ωn φ→ 0
• ω ωn φ→ −arctan 2ζωn−ω = −π
• ω = ωn φ→ −arctan 2ζ0 = −π2
|F(jω)|dB
ω0
-90
-180
φ
ω0
-40db/decade
ωn ωnωr
20log(Q)
ς>0.707
ς<0.707
-20log(2ς)
• Resonance si ζ < 0.707 : ωr = ωn√
1− 2ζ2 et Q =1
2ζ√
1− ζ2
Y. Ariba - Icam, Toulouse. Automatique Lineaire 60 / 131
![Page 119: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/119.jpg)
Reponse d’un systeme Reponse frequentielle
Fonction du second ordre
Courbe de gain : Courbe de phase :
• ω ωn |F |db → −20log(1) = 0
• ω ωn |F |db → −40log( ωωn
)
(pente a −40dB/dec)
• ω = ωn |F |db = −20log(2ζ)
• ω ωn φ→ 0
• ω ωn φ→ −arctan 2ζωn−ω = −π
• ω = ωn φ→ −arctan 2ζ0 = −π2
|F(jω)|dB
ω0
-90
-180
φ
ω0
-40db/decade
ωn ωnωr
20log(Q)
ς>0.707
ς<0.707
-20log(2ς)
• Resonance si ζ < 0.707 : ωr = ωn√
1− 2ζ2 et Q =1
2ζ√
1− ζ2
Y. Ariba - Icam, Toulouse. Automatique Lineaire 60 / 131
![Page 120: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/120.jpg)
Reponse d’un systeme Reponse frequentielle
Exemple
Diagramme de Bode de la fonction de transfert F (p) =20(10p+1)
(100p+1)(p+1)
Fonction de transfert : 20
25
25.5
26
26.5
27
27.5
Mag
nitu
de (
dB)
100 101−1
−0.5
0
0.5
1
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Fonction de transfert : 1100p+1
−40
−30
−20
−10
0
Mag
nitu
de (
dB)
10−4 10−3 10−2 10−1 100−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
ω=0.01
Fonction de transfert : 10p+ 1
0
10
20
30
40
Mag
nitu
de (
dB)
10−3 10−2 10−1 100 1010
45
90
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
ω=0.1
Fonction de transfert : 1p+1
−40
−30
−20
−10
0
Mag
nitu
de (
dB)
10−2 10−1 100 101 102−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
ω=1
Y. Ariba - Icam, Toulouse. Automatique Lineaire 61 / 131
![Page 121: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/121.jpg)
Reponse d’un systeme Reponse frequentielle
Exemple
Diagramme de Bode de la fonction de transfert F (p) =20(10p+1)
(100p+1)(p+1)
Somme des traces
−40
−30
−20
−10
0
10
20
30
Mag
nitu
de (
dB)
10−4 10−3 10−2 10−1 100 101 102−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 61 / 131
![Page 122: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/122.jpg)
Reponse d’un systeme Reponse frequentielle
Representation de Nyquist
Representation graphique du transfert F (jω) dans le plan complexe.
F (jω) = Re[F (jω)
]+ j Im
[F (jω)
]
La courbe est parametree par la pulsation ω et doit etre orientee suivant le sens des ωcroissants.
Im
Re
F(jω)
Re[F(jω)]
Im[F(jω)]
ω=0
ω=∞
Y. Ariba - Icam, Toulouse. Automatique Lineaire 62 / 131
![Page 123: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/123.jpg)
Reponse d’un systeme Reponse frequentielle
Representation de Black-Nichols
Representation graphique du gain par rapport au dephasage.
En abscisses, le dephasage en degre : arg[F (jω)
]En ordonnees, le gain en dB : 20 log
∣∣F (jω)∣∣
La courbe est parametree par la pulsation ω et doit etre orientee suivant le sens des ωcroissants.
20 log |F(jω)|
ω=0
ω=∞
φ(ω)
-180° -90°
Y. Ariba - Icam, Toulouse. Automatique Lineaire 63 / 131
![Page 124: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/124.jpg)
Reponse d’un systeme Resume
Resume
Y. Ariba - Icam, Toulouse. Automatique Lineaire 64 / 131
![Page 125: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/125.jpg)
L’Automatique
Sommaire
1 Modelisation
Exemples introductifs
Notion de systeme
Fonction de transfert
Schema fonctionnel2 Reponse d’un systeme
Reponse temporelle
Reponse frequentielle
Resume3 L’Automatique
Definition
Applications
Architecture de commande
Notion de stabilite
Objectif du cours
4 Stabilite et performances d’un
asservissement
Stabilite
Precision
Rapidite
Marges de stabilite5 Synthese de correcteurs
Introduction
Synthese directe
Action proportionnelle
Action integrale
Action derivee
Combinaisons d’actions
Proportionnel-Integral-Derive
(PID)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 65 / 131
![Page 126: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/126.jpg)
L’Automatique Definition
L’Automatique ou la Theorie de la Commande
Discipline traitant de la modelisation, l’analyse et la commande des systemesdynamiques. Elle a pour objectif l’etude et la conception de dispositifsfonctionnant sans intervention humaine.
Apres les phases :
de modelisation,
d’analyse (etude du comportementale : temporelle et frequentielle),
la question que se pose l’automaticien :
comment agir sur le signal d’entree u(t) pourcontroler le signal de sortie y(t) ?
Y. Ariba - Icam, Toulouse. Automatique Lineaire 66 / 131
![Page 127: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/127.jpg)
L’Automatique Definition
L’Automatique ou la Theorie de la Commande
Discipline traitant de la modelisation, l’analyse et la commande des systemesdynamiques. Elle a pour objectif l’etude et la conception de dispositifsfonctionnant sans intervention humaine.
Apres les phases :
de modelisation,
d’analyse (etude du comportementale : temporelle et frequentielle),
la question que se pose l’automaticien :
comment agir sur le signal d’entree u(t) pourcontroler le signal de sortie y(t) ?
Y. Ariba - Icam, Toulouse. Automatique Lineaire 66 / 131
![Page 128: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/128.jpg)
L’Automatique Definition
L’Automatique ou la Theorie de la Commande
Discipline traitant de la modelisation, l’analyse et la commande des systemesdynamiques. Elle a pour objectif l’etude et la conception de dispositifsfonctionnant sans intervention humaine.
Apres les phases :
de modelisation,
d’analyse (etude du comportementale : temporelle et frequentielle),
la question que se pose l’automaticien :
comment agir sur le signal d’entree u(t) pourcontroler le signal de sortie y(t) ?
Y. Ariba - Icam, Toulouse. Automatique Lineaire 66 / 131
![Page 129: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/129.jpg)
L’Automatique Applications
Application thermodynamique
Regulation de temperature
Entree : alimentation resistance u(t).
Sortie : temperature de l’enceinte θ(t).
Perturbations : echange calorifique avec le milieu ambiant θair(t),introduction d’objets.
u(t)
θ(t)
R Capteur de température
θair(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 67 / 131
![Page 130: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/130.jpg)
L’Automatique Applications
Application spatiale
Systeme de controle d’attitude et d’orbite
Controle de l’orientation d’un satellite pour le pointage.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 68 / 131
![Page 131: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/131.jpg)
L’Automatique Applications
Application spatiale
Systeme de controle d’attitude et d’orbite
Controle de l’orientation d’un satellite pour le pointage.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 68 / 131
![Page 132: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/132.jpg)
L’Automatique Applications
Application robotique
Asservissement en vitesse et en position d’un robot mobile
Commande des moteurs pour le suivi de trajectoire.
Base mobile iRobot R©
Y. Ariba - Icam, Toulouse. Automatique Lineaire 69 / 131
![Page 133: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/133.jpg)
L’Automatique Applications
Application robotique
Asservissement en vitesse et en position d’un robot mobile
Commande des moteurs pour le suivi de trajectoire.
Base mobile iRobot R©
Y. Ariba - Icam, Toulouse. Automatique Lineaire 69 / 131
![Page 134: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/134.jpg)
L’Automatique Applications
Application thermodynamique
Pilotage d’un turboreacteur
Commande de l’injection de carburant pour le controle de la poussee, la pression,la temperature...
ADMISSION ÉCHAPPEMENT
Entrée d'air Chambres de combustion
Section froide Section chaude
Y. Ariba - Icam, Toulouse. Automatique Lineaire 70 / 131
![Page 135: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/135.jpg)
L’Automatique Applications
Application thermodynamique
Pilotage d’un turboreacteur
Commande de l’injection de carburant pour le controle de la poussee, la pression,la temperature...
ADMISSION ÉCHAPPEMENT
Entrée d'air Chambres de combustion
Section froide Section chaude
Y. Ariba - Icam, Toulouse. Automatique Lineaire 70 / 131
![Page 136: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/136.jpg)
L’Automatique Applications
Exemple academique : le pendule.
L’angle d’inclinaison est regi par l’equation :(modele simpliste)
lθ − x cos θ = g sin θ
? Objectif : stabiliser la position verticale haute.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 71 / 131
![Page 137: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/137.jpg)
L’Automatique Applications
Exemple academique : le pendule.
L’angle d’inclinaison est regi par l’equation :(modele simpliste)
lθ − x cos θ = g sin θ
? Objectif : stabiliser la position verticale haute.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 71 / 131
![Page 138: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/138.jpg)
L’Automatique Applications
Exemple academique : le pendule.
L’angle d’inclinaison est regi par l’equation :(modele simpliste)
lθ − x cos θ = g sin θ
? Objectif : stabiliser la position verticale haute.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 71 / 131
![Page 139: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/139.jpg)
L’Automatique Applications
Exemple academique : le pendule.
L’angle d’inclinaison est regi par l’equation :(modele simpliste)
lθ − x cos θ = g sin θ
? Objectif : stabiliser la position verticale haute.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 71 / 131
![Page 140: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/140.jpg)
L’Automatique Applications
c© Quanser
Y. Ariba - Icam, Toulouse. Automatique Lineaire 72 / 131
![Page 141: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/141.jpg)
L’Automatique Applications
c© Quanser
Y. Ariba - Icam, Toulouse. Automatique Lineaire 73 / 131
![Page 142: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/142.jpg)
L’Automatique Architecture de commande
Architecture de commande
Comment agir sur u(t) pour controler y(t) ?
1iere approche naturelle : commande en boucle ouverte
yc(t) est la consigne pour y(t),
necessite une parfaite connaissance du modele,
le systeme et son environnement doivent etre previsibles,
ne permet pas de compenser des perturbations non-connues,
la loi de commande est independante de l’evolution du systeme.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 74 / 131
![Page 143: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/143.jpg)
L’Automatique Architecture de commande
Architecture de commande
Comment agir sur u(t) pour controler y(t) ?
1iere approche naturelle : commande en boucle ouverte
yc(t) est la consigne pour y(t),
necessite une parfaite connaissance du modele,
le systeme et son environnement doivent etre previsibles,
ne permet pas de compenser des perturbations non-connues,
la loi de commande est independante de l’evolution du systeme.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 74 / 131
![Page 144: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/144.jpg)
L’Automatique Architecture de commande
Architecture de commande
Comment agir sur u(t) pour controler y(t) ?
1iere approche naturelle : commande en boucle ouverte
yc(t) est la consigne pour y(t),
necessite une parfaite connaissance du modele,
le systeme et son environnement doivent etre previsibles,
ne permet pas de compenser des perturbations non-connues,
la loi de commande est independante de l’evolution du systeme.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 74 / 131
![Page 145: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/145.jpg)
L’Automatique Architecture de commande
Architecture de commande
Principe fondamental en Automatique : la commande en boucle fermee
⇒ principe de contre-reaction
yc(t) est la consigne pour y(t),
la sortie du systeme est mesuree,
necessite la mise en place d’un capteur,
la loi de commande reagit a l’evolution du systeme,
principe d’asservissement du systeme.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 75 / 131
![Page 146: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/146.jpg)
L’Automatique Architecture de commande
Architecture de commande
Principe fondamental en Automatique : la commande en boucle fermee
⇒ principe de contre-reaction
yc(t) est la consigne pour y(t),
la sortie du systeme est mesuree,
necessite la mise en place d’un capteur,
la loi de commande reagit a l’evolution du systeme,
principe d’asservissement du systeme.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 75 / 131
![Page 147: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/147.jpg)
L’Automatique Architecture de commande
Structure classique d’un asservissement
Système(physique)
Signal de CommandeConsigne Sortie réelle
CapteurSortie
mesurée
+
-
Erreurd'asservissement Contrôleur
loi de commande
Système de commande
Perturbations
En bleu : le systeme physique / dispositif / procede a commander.
Cadre rouge : systeme de commande, realise par l’ingenieur.
Le systeme de commande est mis en oeuvre par une carte electronique.
Attention a differencier le systeme de commande et le systeme commande.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 76 / 131
![Page 148: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/148.jpg)
L’Automatique Architecture de commande
Structure classique d’un asservissement
Exemple.
Position(angle)
Système
Commande du moteur
Soit un systeme a commander : un bras manipulateur.
Un capteur permet de mesurer la position du bras.
La position mesuree est comparee a une position de reference.
Un controleur genere la commande pour corriger l’erreur d’asservissement.
L’etage de puissance permet d’adapter le signal de commande a la puissancede fonctionnement du moteur.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 77 / 131
![Page 149: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/149.jpg)
L’Automatique Architecture de commande
Structure classique d’un asservissement
Exemple.
Position(angle)
Système
Commande du moteur
tension représentative de la position
mesurée capteurpotentiomètre
+-
u
Soit un systeme a commander : un bras manipulateur.
Un capteur permet de mesurer la position du bras.
La position mesuree est comparee a une position de reference.
Un controleur genere la commande pour corriger l’erreur d’asservissement.
L’etage de puissance permet d’adapter le signal de commande a la puissancede fonctionnement du moteur.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 77 / 131
![Page 150: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/150.jpg)
L’Automatique Architecture de commande
Structure classique d’un asservissement
Exemple.
tension représentative de la position
souhaitée Position(angle)+
-
Système
Commande du moteur
tension représentative de la position
mesurée capteurpotentiomètre
+-
u
Soit un systeme a commander : un bras manipulateur.
Un capteur permet de mesurer la position du bras.
La position mesuree est comparee a une position de reference.
Un controleur genere la commande pour corriger l’erreur d’asservissement.
L’etage de puissance permet d’adapter le signal de commande a la puissancede fonctionnement du moteur.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 77 / 131
![Page 151: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/151.jpg)
L’Automatique Architecture de commande
Structure classique d’un asservissement
Exemple.
tension représentative de la position
souhaitée Position(angle)+
-
Système
Commande du moteur
tension représentative de la position
mesurée
Contrôleur
capteurpotentiomètre
+-
u
Soit un systeme a commander : un bras manipulateur.
Un capteur permet de mesurer la position du bras.
La position mesuree est comparee a une position de reference.
Un controleur genere la commande pour corriger l’erreur d’asservissement.
L’etage de puissance permet d’adapter le signal de commande a la puissancede fonctionnement du moteur.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 77 / 131
![Page 152: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/152.jpg)
L’Automatique Architecture de commande
Structure classique d’un asservissement
Exemple.
tension représentative de la position
souhaitée Position(angle)+
-
Système
Commande du moteur
tension représentative de la position
mesurée
Contrôleur
Etage depuissance
capteurpotentiomètre
+-
u
Soit un systeme a commander : un bras manipulateur.
Un capteur permet de mesurer la position du bras.
La position mesuree est comparee a une position de reference.
Un controleur genere la commande pour corriger l’erreur d’asservissement.
L’etage de puissance permet d’adapter le signal de commande a la puissancede fonctionnement du moteur.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 77 / 131
![Page 153: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/153.jpg)
L’Automatique Notion de stabilite
Notion de stabilite
Considerons le systeme de fonction de transfert
y(s)
u(s)=
1
s+ a.
Sa reponse temporelle a un echelon unite est
y(t) =1
a
(1− e−at
).
cas a > 0
0 1 2 3 4 5 60
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Time (sec)
Am
plitu
de y
(t)
le terme e−at tend vers 0
cas a < 0
0 1 2 3 4 5 60
50
100
150
200
250
300
350
400
450
500
Time (sec)
Am
plitu
de y
(t)
le terme e−at tend vers +∞
⇒ Propriete de stabilite
Y. Ariba - Icam, Toulouse. Automatique Lineaire 78 / 131
![Page 154: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/154.jpg)
L’Automatique Notion de stabilite
Notion de stabilite
Considerons le systeme de fonction de transfert
y(s)
u(s)=
1
s+ a.
Sa reponse temporelle a un echelon unite est
y(t) =1
a
(1− e−at
).
cas a > 0
0 1 2 3 4 5 60
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Time (sec)
Am
plitu
de y
(t)
le terme e−at tend vers 0
cas a < 0
0 1 2 3 4 5 60
50
100
150
200
250
300
350
400
450
500
Time (sec)
Am
plitu
de y
(t)
le terme e−at tend vers +∞
⇒ Propriete de stabilite
Y. Ariba - Icam, Toulouse. Automatique Lineaire 78 / 131
![Page 155: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/155.jpg)
L’Automatique Notion de stabilite
Autre exemple, considerons le systeme de fonction de transfert
y(s)
u(s)=
1
s2 + as+ 1, |a| < 2.
Sa reponse temporelle a un echelon unite est
y(t) = 1− e−a2t(
cosω0t +a/2
ω0
sinω0t), ω0 =
√1−
a2
4.
cas 0 < a < 2
0 5 10 15 20 25 30 350
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
Time (sec)
Am
plitu
de y
(t)
le terme e−a2t tend vers 0
cas −2 < a < 0
0 5 10 15 20 25 30 35−800
−600
−400
−200
0
200
400
600
800
1000
1200
Time (sec)
Am
plitu
de y
(t)
le terme e−a2t tend vers +∞
⇒ Propriete de stabilite
Y. Ariba - Icam, Toulouse. Automatique Lineaire 79 / 131
![Page 156: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/156.jpg)
L’Automatique Notion de stabilite
Autre exemple, considerons le systeme de fonction de transfert
y(s)
u(s)=
1
s2 + as+ 1, |a| < 2.
Sa reponse temporelle a un echelon unite est
y(t) = 1− e−a2t(
cosω0t +a/2
ω0
sinω0t), ω0 =
√1−
a2
4.
cas 0 < a < 2
0 5 10 15 20 25 30 350
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
Time (sec)
Am
plitu
de y
(t)
le terme e−a2t tend vers 0
cas −2 < a < 0
0 5 10 15 20 25 30 35−800
−600
−400
−200
0
200
400
600
800
1000
1200
Time (sec)
Am
plitu
de y
(t)
le terme e−a2t tend vers +∞
⇒ Propriete de stabilite
Y. Ariba - Icam, Toulouse. Automatique Lineaire 79 / 131
![Page 157: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/157.jpg)
L’Automatique Notion de stabilite
Definition
Un systeme est dit stable si pour toute entree bornee la sortie est bornee.
Theoreme
Un systeme de fonction de transfert F (s) est stable si et seulement si tous les poles deF (s) sont a partie reelle negative, c’est-a-dire qu’ils appartiennent au demi-plan gauchedu plan complexe.
Exemples
1
s− 2
⇒ Instable
4
s+ 0.5
⇒ Stable
3
(s+ 1)(s+ 3)
⇒ Stable
1
s(5s+ 1)
⇒ Instable
10
s2 + 2s+ 2
⇒ Stable
2s− 1
(s+ 1)(s+ 2)(s− 6)
⇒ Instable
Y. Ariba - Icam, Toulouse. Automatique Lineaire 80 / 131
![Page 158: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/158.jpg)
L’Automatique Notion de stabilite
Definition
Un systeme est dit stable si pour toute entree bornee la sortie est bornee.
Theoreme
Un systeme de fonction de transfert F (s) est stable si et seulement si tous les poles deF (s) sont a partie reelle negative, c’est-a-dire qu’ils appartiennent au demi-plan gauchedu plan complexe.
Exemples
1
s− 2
⇒ Instable
4
s+ 0.5
⇒ Stable
3
(s+ 1)(s+ 3)
⇒ Stable
1
s(5s+ 1)
⇒ Instable
10
s2 + 2s+ 2
⇒ Stable
2s− 1
(s+ 1)(s+ 2)(s− 6)
⇒ Instable
Y. Ariba - Icam, Toulouse. Automatique Lineaire 80 / 131
![Page 159: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/159.jpg)
L’Automatique Notion de stabilite
Definition
Un systeme est dit stable si pour toute entree bornee la sortie est bornee.
Theoreme
Un systeme de fonction de transfert F (s) est stable si et seulement si tous les poles deF (s) sont a partie reelle negative, c’est-a-dire qu’ils appartiennent au demi-plan gauchedu plan complexe.
Exemples
1
s− 2⇒ Instable
4
s+ 0.5⇒ Stable
3
(s+ 1)(s+ 3)⇒ Stable
1
s(5s+ 1)⇒ Instable
10
s2 + 2s+ 2⇒ Stable
2s− 1
(s+ 1)(s+ 2)(s− 6)⇒ Instable
Y. Ariba - Icam, Toulouse. Automatique Lineaire 80 / 131
![Page 160: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/160.jpg)
L’Automatique Notion de stabilite
Analyse de la stabilite entree-sortie
Comment analyser la stabilite d’un systeme d’entree u(s) et de sortie y(s) ?
y(s)
u(s)= G(s) =
bmsm + · · ·+ b1s+ b0
ansn + · · ·+ a1s+ a0
Methodes :
Calcul direct des poles ⇒ trouver les racines de l’equation caracteristique
ansn
+ · · ·+ a1s+ a0 = 0.
Si elles sont toutes a partie reelle negative alors le systeme est stable.
Critere algebrique de Routh.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 81 / 131
![Page 161: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/161.jpg)
L’Automatique Notion de stabilite
Analyse de la stabilite entree-sortie
Comment analyser la stabilite d’un systeme d’entree u(s) et de sortie y(s) ?
y(s)
u(s)= G(s) =
bmsm + · · ·+ b1s+ b0
ansn + · · ·+ a1s+ a0
Methodes :
Calcul direct des poles ⇒ trouver les racines de l’equation caracteristique
ansn
+ · · ·+ a1s+ a0 = 0.
Si elles sont toutes a partie reelle negative alors le systeme est stable.
Critere algebrique de Routh.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 81 / 131
![Page 162: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/162.jpg)
L’Automatique Notion de stabilite
Critere de Routh
Construction d’un tableau a partir des coefficients du polynome caracteristique
D(s) = ansn
+ · · ·+ a1s+ a0
Procedure :
1 Condition necessaire : tous les coeff. ai doivent etre strictement de meme signe.
2 Construction du tableau
pn an an−2 an−4 . . .
pn−1 an−1 an−3 an−5 . . .
pn−2 b1 b2 b3 . . .
pn−3 c1 c2 . . .
.
.
....
.
.
....
p0 α
avec
b1 = −1
an−1
∣∣∣∣ an an−2an−1 an−3
∣∣∣∣ b2 = −1
an−1
∣∣∣∣ an an−4an−1 an−5
∣∣∣∣ b3 = −1
an−1
∣∣∣∣ an an−6an−1 an−7
∣∣∣∣c1 = −
1
b1
∣∣∣∣ an−1 an−3b1 b2
∣∣∣∣ c2 = −1
b1
∣∣∣∣ an−1 an−5b1 b3
∣∣∣∣
3 Le systeme est stable ssi tous les coefficients de la 1ere colonne sont de meme signe.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 82 / 131
![Page 163: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/163.jpg)
L’Automatique Notion de stabilite
Critere de Routh
Construction d’un tableau a partir des coefficients du polynome caracteristique
D(s) = ansn
+ · · ·+ a1s+ a0
Procedure :
1 Condition necessaire : tous les coeff. ai doivent etre strictement de meme signe.
2 Construction du tableau
pn an an−2 an−4 . . .
pn−1 an−1 an−3 an−5 . . .
pn−2 b1 b2 b3 . . .
pn−3 c1 c2 . . .
.
.
....
.
.
....
p0 α
avec
b1 = −1
an−1
∣∣∣∣ an an−2an−1 an−3
∣∣∣∣ b2 = −1
an−1
∣∣∣∣ an an−4an−1 an−5
∣∣∣∣ b3 = −1
an−1
∣∣∣∣ an an−6an−1 an−7
∣∣∣∣c1 = −
1
b1
∣∣∣∣ an−1 an−3b1 b2
∣∣∣∣ c2 = −1
b1
∣∣∣∣ an−1 an−5b1 b3
∣∣∣∣
3 Le systeme est stable ssi tous les coefficients de la 1ere colonne sont de meme signe.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 82 / 131
![Page 164: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/164.jpg)
L’Automatique Notion de stabilite
Critere de Routh
Construction d’un tableau a partir des coefficients du polynome caracteristique
D(s) = ansn
+ · · ·+ a1s+ a0
Procedure :
1 Condition necessaire : tous les coeff. ai doivent etre strictement de meme signe.
2 Construction du tableau
pn an an−2 an−4 . . .
pn−1 an−1 an−3 an−5 . . .
pn−2 b1 b2 b3 . . .
pn−3 c1 c2 . . .
.
.
....
.
.
....
p0 α
avec
b1 = −1
an−1
∣∣∣∣ an an−2an−1 an−3
∣∣∣∣ b2 = −1
an−1
∣∣∣∣ an an−4an−1 an−5
∣∣∣∣ b3 = −1
an−1
∣∣∣∣ an an−6an−1 an−7
∣∣∣∣c1 = −
1
b1
∣∣∣∣ an−1 an−3b1 b2
∣∣∣∣ c2 = −1
b1
∣∣∣∣ an−1 an−5b1 b3
∣∣∣∣
3 Le systeme est stable ssi tous les coefficients de la 1ere colonne sont de meme signe.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 82 / 131
![Page 165: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/165.jpg)
L’Automatique Notion de stabilite
Critere de Routh
Construction d’un tableau a partir des coefficients du polynome caracteristique
D(s) = ansn
+ · · ·+ a1s+ a0
Procedure :
1 Condition necessaire : tous les coeff. ai doivent etre strictement de meme signe.
2 Construction du tableau
pn an an−2 an−4 . . .
pn−1 an−1 an−3 an−5 . . .
pn−2 b1 b2 b3 . . .
pn−3 c1 c2 . . .
.
.
....
.
.
....
p0 α
avec
b1 = −1
an−1
∣∣∣∣ an an−2an−1 an−3
∣∣∣∣ b2 = −1
an−1
∣∣∣∣ an an−4an−1 an−5
∣∣∣∣ b3 = −1
an−1
∣∣∣∣ an an−6an−1 an−7
∣∣∣∣c1 = −
1
b1
∣∣∣∣ an−1 an−3b1 b2
∣∣∣∣ c2 = −1
b1
∣∣∣∣ an−1 an−5b1 b3
∣∣∣∣
3 Le systeme est stable ssi tous les coefficients de la 1ere colonne sont de meme signe.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 82 / 131
![Page 166: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/166.jpg)
L’Automatique Notion de stabilite
Critere de Routh
Construction d’un tableau a partir des coefficients du polynome caracteristique
D(s) = ansn
+ · · ·+ a1s+ a0
Procedure :
1 Condition necessaire : tous les coeff. ai doivent etre strictement de meme signe.
2 Construction du tableau
pn an an−2 an−4 . . .
pn−1 an−1 an−3 an−5 . . .
pn−2 b1 b2 b3 . . .
pn−3 c1 c2 . . .
.
.
....
.
.
....
p0 α
avec
b1 = −1
an−1
∣∣∣∣ an an−2an−1 an−3
∣∣∣∣ b2 = −1
an−1
∣∣∣∣ an an−4an−1 an−5
∣∣∣∣ b3 = −1
an−1
∣∣∣∣ an an−6an−1 an−7
∣∣∣∣c1 = −
1
b1
∣∣∣∣ an−1 an−3b1 b2
∣∣∣∣ c2 = −1
b1
∣∣∣∣ an−1 an−5b1 b3
∣∣∣∣
3 Le systeme est stable ssi tous les coefficients de la 1ere colonne sont de meme signe.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 82 / 131
![Page 167: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/167.jpg)
L’Automatique Notion de stabilite
Exemple 1 :
Soit la fonction de transfert F (s) =s+ 4
s4 + s3 + 4s2 + 2s+ 1, stable ?
s4 1 4 1 0
s3 1 2 0 0
s2 2 1 0
s1 32 0
s0 1
b1 = − 11 (2− 4) = 2
b2 = − 11 (0− 1) = 1
c1 = − 12 (1− 4) = 3
2
c2 = 0
α = − 23 (0− 3
2 ) = 1
⇒ Systeme stable.
Exemple 2 :
Soit la fonction de transfert F (s) =7
3s3 + s2 + 2s+ 4, stable ?
s3 3 2 0
s2 1 4 0
s1 −10 0
s0 4
b1 = − 11 (12− 2) = −10
α = − 1−10 (0 + 40) = 4
⇒ Systeme instable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 83 / 131
![Page 168: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/168.jpg)
L’Automatique Notion de stabilite
Exemple 1 :
Soit la fonction de transfert F (s) =s+ 4
s4 + s3 + 4s2 + 2s+ 1, stable ?
s4 1 4 1 0
s3 1 2 0 0
s2 2 1 0
s1 32 0
s0 1
b1 = − 11 (2− 4) = 2
b2 = − 11 (0− 1) = 1
c1 = − 12 (1− 4) = 3
2
c2 = 0
α = − 23 (0− 3
2 ) = 1
⇒ Systeme stable.
Exemple 2 :
Soit la fonction de transfert F (s) =7
3s3 + s2 + 2s+ 4, stable ?
s3 3 2 0
s2 1 4 0
s1 −10 0
s0 4
b1 = − 11 (12− 2) = −10
α = − 1−10 (0 + 40) = 4
⇒ Systeme instable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 83 / 131
![Page 169: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/169.jpg)
L’Automatique Objectif du cours
Cadre de travail
Ce cours est dedie aux systemes continus lineaires invariants.
Il s’agit d’une classe de modeles particuliere liantl’entree u(t) et la sortie y(t).
Tres utilisee car bien connue.
Système
(ou processus)
y(t)u(t)
d(t)
Continuite : les signaux evoluent continument (ou presque) avec le temps,
x(t) ∈ R avec t ∈ R+.
Linearite : le systeme verifie les principes de superposition et d’homogeneite :u(t) = u1(t) + u2(t)
Systeme−−−−−−→ y(t) = y1(t) + y2(t),
α ∈ R, αu(t)Systeme−−−−−−→ αy(t).
Invariance : le systeme est invariant lorsque
∀τ ∈ R, u(t)Systeme−−−−−−→ y(t) ⇒ u(t− τ)
Systeme−−−−−−→ y(t− τ).
Y. Ariba - Icam, Toulouse. Automatique Lineaire 84 / 131
![Page 170: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/170.jpg)
L’Automatique Objectif du cours
Cadre de travail
Ce cours est dedie aux systemes continus lineaires invariants.
Il s’agit d’une classe de modeles particuliere liantl’entree u(t) et la sortie y(t).
Tres utilisee car bien connue.
Système
(ou processus)
y(t)u(t)
d(t)
Continuite : les signaux evoluent continument (ou presque) avec le temps,
x(t) ∈ R avec t ∈ R+.
Linearite : le systeme verifie les principes de superposition et d’homogeneite :u(t) = u1(t) + u2(t)
Systeme−−−−−−→ y(t) = y1(t) + y2(t),
α ∈ R, αu(t)Systeme−−−−−−→ αy(t).
Invariance : le systeme est invariant lorsque
∀τ ∈ R, u(t)Systeme−−−−−−→ y(t) ⇒ u(t− τ)
Systeme−−−−−−→ y(t− τ).
Y. Ariba - Icam, Toulouse. Automatique Lineaire 84 / 131
![Page 171: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/171.jpg)
L’Automatique Objectif du cours
Quelques remarques :
Attention a ne pas confondre avec l’Automatisme ou l’informatiqueindustrielle.
Composante Mathematiques appliquees importante.
Pluridisciplinaire par ses applications.
L’Automatique est une science cachee ?
Y. Ariba - Icam, Toulouse. Automatique Lineaire 85 / 131
![Page 172: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/172.jpg)
Stabilite et performances d’un asservissement
Sommaire
1 Modelisation
Exemples introductifs
Notion de systeme
Fonction de transfert
Schema fonctionnel2 Reponse d’un systeme
Reponse temporelle
Reponse frequentielle
Resume3 L’Automatique
Definition
Applications
Architecture de commande
Notion de stabilite
Objectif du cours
4 Stabilite et performances d’un
asservissement
Stabilite
Precision
Rapidite
Marges de stabilite5 Synthese de correcteurs
Introduction
Synthese directe
Action proportionnelle
Action integrale
Action derivee
Combinaisons d’actions
Proportionnel-Integral-Derive
(PID)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 86 / 131
![Page 173: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/173.jpg)
Stabilite et performances d’un asservissement Stabilite
Stabilite d’un asservissement
Comment analyser la stabilite d’un systeme asservi ?
Methodes :
Ecrire la fonction de transfert globale equivalente T (s) =C(s)G(s)
1 + C(s)G(s),
⇒ Appliquer l’une des deux methodes vues precedemment.
Critere du revers (critere graphique).
Si le systeme en boucle ouverte est stable et a minimum dephase (poles et zeros a partie reelle strictement negative) alors lesysteme asservi est stable si et seulement si le point critique
(−1 , 0) est laisse a gauche quand on parcourt le lieu detransfert de la boucle ouverte dans le plan de Nyquist dans lesens des ω croissants.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 87 / 131
![Page 174: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/174.jpg)
Stabilite et performances d’un asservissement Stabilite
Stabilite d’un asservissement
Comment analyser la stabilite d’un systeme asservi ?
Methodes :
Ecrire la fonction de transfert globale equivalente T (s) =C(s)G(s)
1 + C(s)G(s),
⇒ Appliquer l’une des deux methodes vues precedemment.
Critere du revers (critere graphique).
Si le systeme en boucle ouverte est stable et a minimum dephase (poles et zeros a partie reelle strictement negative) alors lesysteme asservi est stable si et seulement si le point critique
(−1 , 0) est laisse a gauche quand on parcourt le lieu detransfert de la boucle ouverte dans le plan de Nyquist dans lesens des ω croissants.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 87 / 131
![Page 175: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/175.jpg)
Stabilite et performances d’un asservissement Stabilite
Exemple 1 :
Les fonctions de transfert du procede et du correcteur sont de la forme
G(s) =10
10s2 + 18s+ 8et C(s) =
4
s2 + 3s+ 2.
Tracons le lieu de Nyquist de la fonction de transfert en boucle ouverte C(s)G(s).
−1.5 −1 −0.5 0 0.5 1 1.5 2 2.5−2
−1.5
−1
−0.5
0
0.5Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ le systeme en boucle fermee est stable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 88 / 131
![Page 176: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/176.jpg)
Stabilite et performances d’un asservissement Stabilite
Exemple 1 :
Les fonctions de transfert du procede et du correcteur sont de la forme
G(s) =10
10s2 + 18s+ 8et C(s) =
4
s2 + 3s+ 2.
Tracons le lieu de Nyquist de la fonction de transfert en boucle ouverte C(s)G(s).
−1.5 −1 −0.5 0 0.5 1 1.5 2 2.5−2
−1.5
−1
−0.5
0
0.5Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ le systeme en boucle fermee est stable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 88 / 131
![Page 177: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/177.jpg)
Stabilite et performances d’un asservissement Stabilite
Exemple 1 :
Les fonctions de transfert du procede et du correcteur sont de la forme
G(s) =10
10s2 + 18s+ 8et C(s) =
4
s2 + 3s+ 2.
Tracons le lieu de Nyquist de la fonction de transfert en boucle ouverte C(s)G(s).
−1.5 −1 −0.5 0 0.5 1 1.5 2 2.5−2
−1.5
−1
−0.5
0
0.5Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ le systeme en boucle fermee est stable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 88 / 131
![Page 178: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/178.jpg)
Stabilite et performances d’un asservissement Stabilite
Exemple 2 :
Les fonctions de transfert du procede et du correcteur sont de la forme
G(s) =5
s2 + 0.4s+ 1et C(s) =
s+ 3
10s+ 1.
Tracons le lieu de Nyquist de la fonction de transfert en boucle ouverte C(s)G(s).
−4 −2 0 2 4 6 8 10 12 14 16−8
−7
−6
−5
−4
−3
−2
−1
0
1
2Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ le systeme en boucle fermee est instable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 89 / 131
![Page 179: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/179.jpg)
Stabilite et performances d’un asservissement Stabilite
Exemple 2 :
Les fonctions de transfert du procede et du correcteur sont de la forme
G(s) =5
s2 + 0.4s+ 1et C(s) =
s+ 3
10s+ 1.
Tracons le lieu de Nyquist de la fonction de transfert en boucle ouverte C(s)G(s).
−4 −2 0 2 4 6 8 10 12 14 16−8
−7
−6
−5
−4
−3
−2
−1
0
1
2Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ le systeme en boucle fermee est instable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 89 / 131
![Page 180: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/180.jpg)
Stabilite et performances d’un asservissement Stabilite
Exemple 2 :
Les fonctions de transfert du procede et du correcteur sont de la forme
G(s) =5
s2 + 0.4s+ 1et C(s) =
s+ 3
10s+ 1.
Tracons le lieu de Nyquist de la fonction de transfert en boucle ouverte C(s)G(s).
−4 −2 0 2 4 6 8 10 12 14 16−8
−7
−6
−5
−4
−3
−2
−1
0
1
2Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ le systeme en boucle fermee est instable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 89 / 131
![Page 181: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/181.jpg)
Stabilite et performances d’un asservissement Stabilite
Performance d’un asservissement
Il existe differents criteres pour caracteriser un asservissement.
En plus de la stabilite, d’autres proprietes peuvent etre interessantes :
la precision.
la rapidite.
la marge de stabilite.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 90 / 131
![Page 182: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/182.jpg)
Stabilite et performances d’un asservissement Precision
Precision
La precision est determinee par l’erreur d’asservissement en regime permanent :
ε(t) = e(t)− y(t) lorsque t→∞
On definit :
erreur statique εs
Lorsque l’entree est un echelon
e(t) = e0, ∀t ≥ 0
t
e(t)
0
e0
y(t)
εs
erreur de trainage εt
Lorsque l’entree est une rampe
e(t) = e0t, ∀t ≥ 0
t0
εte(t) y(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 91 / 131
![Page 183: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/183.jpg)
Stabilite et performances d’un asservissement Precision
Precision
La precision est determinee par l’erreur d’asservissement en regime permanent :
ε(t) = e(t)− y(t) lorsque t→∞On definit :
erreur statique εs
Lorsque l’entree est un echelon
e(t) = e0, ∀t ≥ 0
t
e(t)
0
e0
y(t)
εs
erreur de trainage εt
Lorsque l’entree est une rampe
e(t) = e0t, ∀t ≥ 0
t0
εte(t) y(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 91 / 131
![Page 184: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/184.jpg)
Stabilite et performances d’un asservissement Precision
Precision
La precision est determinee par l’erreur d’asservissement en regime permanent :
ε(t) = e(t)− y(t) lorsque t→∞On definit :
erreur statique εs
Lorsque l’entree est un echelon
e(t) = e0, ∀t ≥ 0
t
e(t)
0
e0
y(t)
εs
erreur de trainage εt
Lorsque l’entree est une rampe
e(t) = e0t, ∀t ≥ 0
t0
εte(t) y(t)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 91 / 131
![Page 185: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/185.jpg)
Stabilite et performances d’un asservissement Precision
L’erreur en regime permanent s’exprime par limt→∞
e(t)− y(t).
Selon le theoreme de la valeur finale limt→∞
ε(t) = lims→0
s ε(s).
En pratique, on utilise la transformee de Laplace :
ε(s) = e(s)− y(s) =
(1− F (s)
)e(s)
ou F (s) =G(s)
1 +G(s)est la fonction de transfert en boucle fermee.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 92 / 131
![Page 186: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/186.jpg)
Stabilite et performances d’un asservissement Precision
L’erreur en regime permanent s’exprime par limt→∞
e(t)− y(t).
Selon le theoreme de la valeur finale limt→∞
ε(t) = lims→0
s ε(s).
En pratique, on utilise la transformee de Laplace :
ε(s) = e(s)− y(s) =
(1− F (s)
)e(s)
ou F (s) =G(s)
1 +G(s)est la fonction de transfert en boucle fermee.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 92 / 131
![Page 187: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/187.jpg)
Stabilite et performances d’un asservissement Precision
L’erreur en regime permanent s’exprime par limt→∞
e(t)− y(t).
Selon le theoreme de la valeur finale limt→∞
ε(t) = lims→0
s ε(s).
En pratique, on utilise la transformee de Laplace :
ε(s) = e(s)− y(s) =
(1− F (s)
)e(s)
ou F (s) =G(s)
1 +G(s)est la fonction de transfert en boucle fermee.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 92 / 131
![Page 188: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/188.jpg)
Stabilite et performances d’un asservissement Precision
Exemple 1 : G(s) =10
s+ 20
Quelques calculs donnent : F (s) =10
s+ 30et ε(s) =
s+ 20
s+ 30e(s).
On en deduit :
erreur statique εs = lims→0
ss+ 20
s+ 30
e0
s=
2
3e0,
erreur de trainage εt = lims→0
ss+ 20
s+ 30
e0
s2= +∞.
Exemple 2 : G(s) =Kω2
n
s2 + 2ζωns+ ω2n
Quelques calculs donnent :
F (s) =Kω2
n
p2 + 2ζωnp+ ω2n(K + 1)
et ε(s) =s2 + 2ζωns+ ω2
n
s2 + 2ζωns+ ω2n(K + 1)
e(s).
On en deduit :
erreur statique εs = lims→0
s
(1− F (s)
)e0
s=
1
K + 1e0,
erreur de trainage εt = lims→0
s
(1− F (s)
)e0
s2= +∞.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 93 / 131
![Page 189: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/189.jpg)
Stabilite et performances d’un asservissement Precision
Exemple 1 : G(s) =10
s+ 20
Quelques calculs donnent : F (s) =10
s+ 30et ε(s) =
s+ 20
s+ 30e(s).
On en deduit :
erreur statique εs = lims→0
ss+ 20
s+ 30
e0
s=
2
3e0,
erreur de trainage εt = lims→0
ss+ 20
s+ 30
e0
s2= +∞.
Exemple 2 : G(s) =Kω2
n
s2 + 2ζωns+ ω2n
Quelques calculs donnent :
F (s) =Kω2
n
p2 + 2ζωnp+ ω2n(K + 1)
et ε(s) =s2 + 2ζωns+ ω2
n
s2 + 2ζωns+ ω2n(K + 1)
e(s).
On en deduit :
erreur statique εs = lims→0
s
(1− F (s)
)e0
s=
1
K + 1e0,
erreur de trainage εt = lims→0
s
(1− F (s)
)e0
s2= +∞.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 93 / 131
![Page 190: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/190.jpg)
Stabilite et performances d’un asservissement Precision
Exemple 1 : G(s) =10
s+ 20
Quelques calculs donnent : F (s) =10
s+ 30et ε(s) =
s+ 20
s+ 30e(s).
On en deduit :
erreur statique εs = lims→0
ss+ 20
s+ 30
e0
s=
2
3e0,
erreur de trainage εt = lims→0
ss+ 20
s+ 30
e0
s2= +∞.
Exemple 2 : G(s) =Kω2
n
s2 + 2ζωns+ ω2n
Quelques calculs donnent :
F (s) =Kω2
n
p2 + 2ζωnp+ ω2n(K + 1)
et ε(s) =s2 + 2ζωns+ ω2
n
s2 + 2ζωns+ ω2n(K + 1)
e(s).
On en deduit :
erreur statique εs = lims→0
s
(1− F (s)
)e0
s=
1
K + 1e0,
erreur de trainage εt = lims→0
s
(1− F (s)
)e0
s2= +∞.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 93 / 131
![Page 191: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/191.jpg)
Stabilite et performances d’un asservissement Precision
Exemple 1 : G(s) =10
s+ 20
Quelques calculs donnent : F (s) =10
s+ 30et ε(s) =
s+ 20
s+ 30e(s).
On en deduit :
erreur statique εs = lims→0
ss+ 20
s+ 30
e0
s=
2
3e0,
erreur de trainage εt = lims→0
ss+ 20
s+ 30
e0
s2= +∞.
Exemple 2 : G(s) =Kω2
n
s2 + 2ζωns+ ω2n
Quelques calculs donnent :
F (s) =Kω2
n
p2 + 2ζωnp+ ω2n(K + 1)
et ε(s) =s2 + 2ζωns+ ω2
n
s2 + 2ζωns+ ω2n(K + 1)
e(s).
On en deduit :
erreur statique εs = lims→0
s
(1− F (s)
)e0
s=
1
K + 1e0,
erreur de trainage εt = lims→0
s
(1− F (s)
)e0
s2= +∞.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 93 / 131
![Page 192: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/192.jpg)
Stabilite et performances d’un asservissement Precision
Regles generales :
Un gain statique eleve en boucle ouverte permet d’obtenir une erreur de position(en asservissement) plus faible.
Loi des integrateurs : l’erreur en regime permanent, pour une entree e(s), est nullesi la boucle ouverte comprend au moins autant d’integrateurs que le signal e(s).
- S’il n’y a pas d’integrateur pur dans la FTBO, l’erreur de position est finie etl’erreur de trainage est infinie.
- S’il y a un integrateur pur dans la FTBO, l’erreur de position est nulle et l’erreurde trainage est finie.
- S’il y a deux integrateur pur dans la FTBO, l’erreur de position est nulle etl’erreur de trainage est nulle.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 94 / 131
![Page 193: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/193.jpg)
Stabilite et performances d’un asservissement Precision
Exemple 3 :
Quelques calculs donnent : F (s) =10
10 + s(20 + s)et ε(s) =
s(20 + s)
10 + s(20 + s)e(s).
erreur statique εs = lims→0
ss(20 + s)
10 + s(20 + s)
e0
s= 0,
erreur de trainage εt = lims→0
ss(20 + s)
10 + s(20 + s)
e0
s2= 2e0.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 95 / 131
![Page 194: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/194.jpg)
Stabilite et performances d’un asservissement Rapidite
Rapidite
La performance en rapidite d’un asservissement est caracterisee par le temps de reponseet le temps de monte.
Le temps de reponse est letemps mis par le signal de sortiepour atteindre sa valeur finale an% pres (souvent 5%) sansressortir de cet intervalle.
t
y(t)
y(∞)
±5% y(∞)
tr5%
Le temps de montee corresponda l’intervalle de temps dans lequella sortie passe de 10% a 90% de lavaleur finale.
t
y(t)
y(∞)90% y(∞)
10% y(∞)
tmY. Ariba - Icam, Toulouse. Automatique Lineaire 96 / 131
![Page 195: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/195.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Marges de stabilite
Soit le systeme a commander : G(s) =k
s(s+ 1)(s+ 2).
Stabilite du systeme asservi F (s) =k
s(s+ 1)(s+ 2) + k?
s3
1 2 0
s2
3 k 0
s1 6− k
30
s0
k
⇒ systeme stable si 0 < k < 6.
−3 −2.5 −2 −1.5 −1 −0.5 0 0.5−10
−8
−6
−4
−2
0
2
Nyquist Diagram
Real Axis
Imag
inar
y A
xis
k=1
k=4
F Notion de marges de stabilite : quantifient “l’eloignement”par rapport au seuilcritique d’instabilite.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 97 / 131
![Page 196: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/196.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Marges de stabilite
Soit le systeme a commander : G(s) =k
s(s+ 1)(s+ 2).
Stabilite du systeme asservi F (s) =k
s(s+ 1)(s+ 2) + k?
s3
1 2 0
s2
3 k 0
s1 6− k
30
s0
k
⇒ systeme stable si 0 < k < 6.
−3 −2.5 −2 −1.5 −1 −0.5 0 0.5−10
−8
−6
−4
−2
0
2
Nyquist Diagram
Real Axis
Imag
inar
y A
xis
k=1
k=4
F Notion de marges de stabilite : quantifient “l’eloignement”par rapport au seuilcritique d’instabilite.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 97 / 131
![Page 197: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/197.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Marges de stabilite
Soit le systeme a commander : G(s) =k
s(s+ 1)(s+ 2).
Stabilite du systeme asservi F (s) =k
s(s+ 1)(s+ 2) + k?
s3
1 2 0
s2
3 k 0
s1 6− k
30
s0
k
⇒ systeme stable si 0 < k < 6.
−3 −2.5 −2 −1.5 −1 −0.5 0 0.5−10
−8
−6
−4
−2
0
2
Nyquist Diagram
Real Axis
Imag
inar
y A
xis
k=1
k=4
F Notion de marges de stabilite : quantifient “l’eloignement”par rapport au seuilcritique d’instabilite.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 97 / 131
![Page 198: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/198.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Marge de gain
La marge de gain se definit comme le gain qu’il faut apporter au systeme en boucleouverte pour destabiliser le systeme asservi (en boucle fermee).
Mesure :
La marge de gain d’un systeme asservi est donnee par la formule :
∆G = −20 log |G(jω−π)|
ou ω−π est la pulsation pour laquelle la FTBO G(jω) est dephase de −180.
−2 −1 0 1 2 3 4−3
−2.5
−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis
AC
Y. Ariba - Icam, Toulouse. Automatique Lineaire 98 / 131
![Page 199: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/199.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Marge de phase
La marge de phase se definit comme le dephasage qu’il faut apporter au systeme enboucle ouverte pour destabiliser le systeme asservi (en boucle fermee).
Mesure :
La marge de phase d’un systeme asservi est donnee par la formule :
∆φ = π + argG(jω0dB)
ou ω0dB est la pulsation pour laquelle la FTBO G(jω) a un gain unitaire (0 en decibel).
−2 −1 0 1 2 3 4−3
−2.5
−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis
C
A
Y. Ariba - Icam, Toulouse. Automatique Lineaire 99 / 131
![Page 200: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/200.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Exemple : Soit l’asservissement suivant
avec G(s) =5
s3 + 3.5s2 + 3.5s+ 1. L’asservissement est-il stable ?
−2 −1 0 1 2 3 4 5 6−4
−3.5
−3
−2.5
−2
−1.5
−1
−0.5
0
0.5
1
1.5Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ Par application du critere du revers : stable
Y. Ariba - Icam, Toulouse. Automatique Lineaire 100 / 131
![Page 201: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/201.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Exemple : Soit l’asservissement suivant
avec G(s) =5
s3 + 3.5s2 + 3.5s+ 1. L’asservissement est-il stable ?
−2 −1 0 1 2 3 4 5 6−4
−3.5
−3
−2.5
−2
−1.5
−1
−0.5
0
0.5
1
1.5Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ Par application du critere du revers : stable
Y. Ariba - Icam, Toulouse. Automatique Lineaire 100 / 131
![Page 202: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/202.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Exemple : Soit l’asservissement suivant
avec G(s) =5
s3 + 3.5s2 + 3.5s+ 1. L’asservissement est-il stable ?
−2 −1 0 1 2 3 4 5 6−4
−3.5
−3
−2.5
−2
−1.5
−1
−0.5
0
0.5
1
1.5Nyquist Diagram
Real Axis
Imag
inar
y A
xis
⇒ Par application du critere du revers : stable
Y. Ariba - Icam, Toulouse. Automatique Lineaire 100 / 131
![Page 203: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/203.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Quelle est sa marge de stabilite ?
Des mesures sur le lieu de G(s) dans Bode ou Nyquist donnent
|G(jω)| = 1 pour la pulsation ω = 1.24 rad/s,
arg(G(jω)
)= −180 pour la pulsation ω = 1.87 rad/s.
Calculs des marges
marge de gain :
∆G = −20 log5
|(jw)3 + 3.5(jw)2 + 3.5(jw) + 1|, avec ω = 1.87
= 7.04 dB
marge de phase :
∆φ = π + arg5
(jw)3 + 3.5(jw)2 + 3.5(jw) + 1, avec ω = 1.24
= 0.51 rad (29.2 deg)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 101 / 131
![Page 204: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/204.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Quelle est sa marge de stabilite ?
Des mesures sur le lieu de G(s) dans Bode ou Nyquist donnent
|G(jω)| = 1 pour la pulsation ω = 1.24 rad/s,
arg(G(jω)
)= −180 pour la pulsation ω = 1.87 rad/s.
Calculs des marges
marge de gain :
∆G = −20 log5
|(jw)3 + 3.5(jw)2 + 3.5(jw) + 1|, avec ω = 1.87
= 7.04 dB
marge de phase :
∆φ = π + arg5
(jw)3 + 3.5(jw)2 + 3.5(jw) + 1, avec ω = 1.24
= 0.51 rad (29.2 deg)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 101 / 131
![Page 205: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/205.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Quelle est sa marge de stabilite ?
Des mesures sur le lieu de G(s) dans Bode ou Nyquist donnent
|G(jω)| = 1 pour la pulsation ω = 1.24 rad/s,
arg(G(jω)
)= −180 pour la pulsation ω = 1.87 rad/s.
Calculs des marges
marge de gain :
∆G = −20 log5
|(jw)3 + 3.5(jw)2 + 3.5(jw) + 1|, avec ω = 1.87
= 7.04 dB
marge de phase :
∆φ = π + arg5
(jw)3 + 3.5(jw)2 + 3.5(jw) + 1, avec ω = 1.24
= 0.51 rad (29.2 deg)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 101 / 131
![Page 206: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/206.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Celles-ci peuvent aussi etre directement mesurees sur le lieu de Nyquist de G(s)
−1.5 −1 −0.5 0 0.5 1 1.5−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis
a
distance a ' 0.44 : marge de gain ∆G = 20 log(a) ' 7.05 dB.
angle b ' 29 : marge de gain ∆φ = b ' 29 deg.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 102 / 131
![Page 207: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/207.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Celles-ci peuvent aussi etre directement mesurees sur le lieu de Nyquist de G(s)
−1.5 −1 −0.5 0 0.5 1 1.5−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis
a
distance a ' 0.44 : marge de gain ∆G = 20 log(a) ' 7.05 dB.
angle b ' 29 : marge de gain ∆φ = b ' 29 deg.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 102 / 131
![Page 208: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/208.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Celles-ci peuvent aussi etre directement mesurees sur le lieu de Nyquist de G(s)
−1.5 −1 −0.5 0 0.5 1 1.5−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis
distance a ' 0.44 : marge de gain ∆G = 20 log(a) ' 7.05 dB.
angle b ' 29 : marge de gain ∆φ = b ' 29 deg.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 102 / 131
![Page 209: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/209.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Celles-ci peuvent aussi etre directement mesurees sur le lieu de Nyquist de G(s)
−1.5 −1 −0.5 0 0.5 1 1.5−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis b
distance a ' 0.44 : marge de gain ∆G = 20 log(a) ' 7.05 dB.
angle b ' 29 : marge de gain ∆φ = b ' 29 deg.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 102 / 131
![Page 210: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/210.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Celles-ci peuvent aussi etre directement mesurees sur le lieu de Nyquist de G(s)
−1.5 −1 −0.5 0 0.5 1 1.5−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis b
distance a ' 0.44 : marge de gain ∆G = 20 log(a) ' 7.05 dB.
angle b ' 29 : marge de gain ∆φ = b ' 29 deg.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 102 / 131
![Page 211: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/211.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Celles-ci peuvent aussi etre directement mesurees sur le lieu de Nyquist de G(s)
−1.5 −1 −0.5 0 0.5 1 1.5−2
−1.5
−1
−0.5
0
0.5
1Nyquist Diagram
Real Axis
Imag
inar
y A
xis
a
b
distance a ' 0.44 : marge de gain ∆G = 20 log(a) ' 7.05 dB.
angle b ' 29 : marge de gain ∆φ = b ' 29 deg.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 102 / 131
![Page 212: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/212.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Ou encore, sur le diagramme de Bode de G(s)
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
10−2 10−1 100 101−270
−180
−90
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 103 / 131
![Page 213: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/213.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Ou encore, sur le diagramme de Bode de G(s)
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
10−2 10−1 100 101−270
−180
−90
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 103 / 131
![Page 214: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/214.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Ou encore, sur le diagramme de Bode de G(s)
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
10−2 10−1 100 101−270
−180
−90
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
∆ G
Y. Ariba - Icam, Toulouse. Automatique Lineaire 103 / 131
![Page 215: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/215.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Ou encore, sur le diagramme de Bode de G(s)
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
10−2 10−1 100 101−270
−180
−90
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 103 / 131
![Page 216: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/216.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Ou encore, sur le diagramme de Bode de G(s)
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
10−2 10−1 100 101−270
−180
−90
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
∆ φ
Y. Ariba - Icam, Toulouse. Automatique Lineaire 103 / 131
![Page 217: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/217.jpg)
Stabilite et performances d’un asservissement Marges de stabilite
Ou encore, sur le diagramme de Bode de G(s)
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
10−2 10−1 100 101−270
−180
−90
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
∆ φ
∆ G
Y. Ariba - Icam, Toulouse. Automatique Lineaire 103 / 131
![Page 218: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/218.jpg)
Synthese de correcteurs
Sommaire
1 Modelisation
Exemples introductifs
Notion de systeme
Fonction de transfert
Schema fonctionnel2 Reponse d’un systeme
Reponse temporelle
Reponse frequentielle
Resume3 L’Automatique
Definition
Applications
Architecture de commande
Notion de stabilite
Objectif du cours
4 Stabilite et performances d’un
asservissement
Stabilite
Precision
Rapidite
Marges de stabilite5 Synthese de correcteurs
Introduction
Synthese directe
Action proportionnelle
Action integrale
Action derivee
Combinaisons d’actions
Proportionnel-Integral-Derive (PID)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 104 / 131
![Page 219: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/219.jpg)
Synthese de correcteurs Introduction
Introduction
Phase de modelisation :
Chercher un modele representant le comportement du systeme.
relation entree-sortie ?
equation differentielle, fonction de transfert, schema-bloc...
Phase d’analyse :
Analyse des proprietes du modele et ses performances (en bo ou bf).
reponses temporelles et frequentielles
stabilite
analyse d’un asservissement pour un correcteur donne
⇒ Ce que vous avez vu jusqu’a maintenant.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 105 / 131
![Page 220: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/220.jpg)
Synthese de correcteurs Introduction
Introduction
Phase de modelisation :
Chercher un modele representant le comportement du systeme.
relation entree-sortie ?
equation differentielle, fonction de transfert, schema-bloc...
Phase d’analyse :
Analyse des proprietes du modele et ses performances (en bo ou bf).
reponses temporelles et frequentielles
stabilite
analyse d’un asservissement pour un correcteur donne
⇒ Ce que vous avez vu jusqu’a maintenant.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 105 / 131
![Page 221: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/221.jpg)
Synthese de correcteurs Introduction
Introduction
Phase de modelisation :
Chercher un modele representant le comportement du systeme.
relation entree-sortie ?
equation differentielle, fonction de transfert, schema-bloc...
Phase d’analyse :
Analyse des proprietes du modele et ses performances (en bo ou bf).
reponses temporelles et frequentielles
stabilite
analyse d’un asservissement pour un correcteur donne
⇒ Ce que vous avez vu jusqu’a maintenant.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 105 / 131
![Page 222: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/222.jpg)
Synthese de correcteurs Introduction
Phase de synthese :
Conception d’un systeme de commande.
rechercher un correcteur
stabilisation
amelioration des performances
⇒ Le systeme en boucle fermee doit satisfaire un certain cdc.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 106 / 131
![Page 223: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/223.jpg)
Synthese de correcteurs Introduction
Un systeme de commande a pour objectif de doter le systeme asservi de certainesproprietes telles que
la stabilite du systeme asservi,
la rapidite de la reponse temporelle (regime transitoire),
la precision en regime permanent,
la robustesse (marges de stabilite, rejet de perturbation).
Dans ce cours, seule la synthese par methodes frequentielles classiques est consideree.Elles sont tres utilisees en entreprises du fait de
leur aspect pratique,
l’existence de techniques de synthese simples ou empiriques,
la possibilite de regler les gains intuitivement.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 107 / 131
![Page 224: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/224.jpg)
Synthese de correcteurs Introduction
Un systeme de commande a pour objectif de doter le systeme asservi de certainesproprietes telles que
la stabilite du systeme asservi,
la rapidite de la reponse temporelle (regime transitoire),
la precision en regime permanent,
la robustesse (marges de stabilite, rejet de perturbation).
Dans ce cours, seule la synthese par methodes frequentielles classiques est consideree.Elles sont tres utilisees en entreprises du fait de
leur aspect pratique,
l’existence de techniques de synthese simples ou empiriques,
la possibilite de regler les gains intuitivement.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 107 / 131
![Page 225: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/225.jpg)
Synthese de correcteurs Introduction
Pour la suite de l’expose, nous garderons les notations suivantes :
C(s) represente le correcteur.
G(s) represente le systeme a commander.
u(s) est le signal de commande.
F (s) represente la fonction de transfert en boucle fermee.
Dans cette configuration la FTBO correspond a la chaıne directe : C(s)G(s).
L’erreur d’asservissement est notee ε(s) = e(s)− y(s).
Y. Ariba - Icam, Toulouse. Automatique Lineaire 108 / 131
![Page 226: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/226.jpg)
Synthese de correcteurs Synthese directe
Synthese directe : modele impose
Une 1iere methode directe et pratique consiste a imposer un modele pour la FTBF.
Qu’est-ce qu’un modele “ideal” ?
pour un premier ordre : Fd(s) =1
τs+ 1
0
0.2
0.4
0.6
0.8
1
Step Response
Time (sec)
Am
plitu
de
τ
63%
la constante de temps τ definit la rapidite : tr5% = 3τ
Y. Ariba - Icam, Toulouse. Automatique Lineaire 109 / 131
![Page 227: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/227.jpg)
Synthese de correcteurs Synthese directe
Synthese directe : modele impose
Une 1iere methode directe et pratique consiste a imposer un modele pour la FTBF.
Qu’est-ce qu’un modele “ideal” ?
pour un premier ordre : Fd(s) =1
τs+ 1
0
0.2
0.4
0.6
0.8
1
Step Response
Time (sec)
Am
plitu
de
τ
63%
la constante de temps τ definit la rapidite : tr5% = 3τ
Y. Ariba - Icam, Toulouse. Automatique Lineaire 109 / 131
![Page 228: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/228.jpg)
Synthese de correcteurs Synthese directe
Synthese directe : modele impose
Une 1iere methode directe et pratique consiste a imposer un modele pour la FTBF.
Qu’est-ce qu’un modele “ideal” ?
pour un premier ordre : Fd(s) =1
τs+ 1
0
0.2
0.4
0.6
0.8
1
Step Response
Time (sec)
Am
plitu
de
τ
63%
la constante de temps τ definit la rapidite : tr5% = 3τ
Y. Ariba - Icam, Toulouse. Automatique Lineaire 109 / 131
![Page 229: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/229.jpg)
Synthese de correcteurs Synthese directe
pour un second ordre : Fd(s) =1
s2
ω2n
+ 2ζωns+ 1
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
Step Response
Time (sec)
Am
plitu
de ζ=0.8
ζ=0.1
ωn=1
l’amortissement ζ definit le depassement : D1 = 100 e− ζπ√
1−ζ2
la pulsation propre et l’amortissement definissent la rapidite tr5% ' 3ζωn
Y. Ariba - Icam, Toulouse. Automatique Lineaire 110 / 131
![Page 230: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/230.jpg)
Synthese de correcteurs Synthese directe
Methode :
1 Specifier une FTBF ayant une dynamique souhaitee “ideale”.
type premier ordre
Fd(s) =1
τs+ 1
type second ordre
Fd(s) =1
s2
ω2n
+ 2ζωns+ 1
pour les ordres superieurs on peut envisager : un pole multiples ou un secondordre dominant.
2 En identifiant la FTBF avec la fonction desiree, calculer le correcteur :
Fd(p) = F (p) ⇒ Fd(p) =C(p)G(p)
1 + C(p)G(p)
⇒ C(p) =Fd(p)
G(p)(
1− Fd(p))
Y. Ariba - Icam, Toulouse. Automatique Lineaire 111 / 131
![Page 231: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/231.jpg)
Synthese de correcteurs Synthese directe
Methode :
1 Specifier une FTBF ayant une dynamique souhaitee “ideale”.
type premier ordre
Fd(s) =1
τs+ 1
type second ordre
Fd(s) =1
s2
ω2n
+ 2ζωns+ 1
pour les ordres superieurs on peut envisager : un pole multiples ou un secondordre dominant.
2 En identifiant la FTBF avec la fonction desiree, calculer le correcteur :
Fd(p) = F (p) ⇒ Fd(p) =C(p)G(p)
1 + C(p)G(p)
⇒ C(p) =Fd(p)
G(p)(
1− Fd(p))
Y. Ariba - Icam, Toulouse. Automatique Lineaire 111 / 131
![Page 232: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/232.jpg)
Synthese de correcteurs Synthese directe
3 Le correcteur doit etre propre pour etre realisable.(ordre numerateur ≤ ordre denominateur)
4 Une phase de simulation est necessaire pour valider le fonctionnement(comportement attendu, pas de saturation...).
Remarques :
L’approche est valide seulement si le modele du systeme est bien connu.
Le systeme a commander doit posseder des poles et zeros stables.
L’approche est adaptee plutot pour des systemes d’ordre faible.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 112 / 131
![Page 233: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/233.jpg)
Synthese de correcteurs Synthese directe
3 Le correcteur doit etre propre pour etre realisable.(ordre numerateur ≤ ordre denominateur)
4 Une phase de simulation est necessaire pour valider le fonctionnement(comportement attendu, pas de saturation...).
Remarques :
L’approche est valide seulement si le modele du systeme est bien connu.
Le systeme a commander doit posseder des poles et zeros stables.
L’approche est adaptee plutot pour des systemes d’ordre faible.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 112 / 131
![Page 234: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/234.jpg)
Synthese de correcteurs Synthese directe
3 Le correcteur doit etre propre pour etre realisable.(ordre numerateur ≤ ordre denominateur)
4 Une phase de simulation est necessaire pour valider le fonctionnement(comportement attendu, pas de saturation...).
Remarques :
L’approche est valide seulement si le modele du systeme est bien connu.
Le systeme a commander doit posseder des poles et zeros stables.
L’approche est adaptee plutot pour des systemes d’ordre faible.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 112 / 131
![Page 235: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/235.jpg)
Synthese de correcteurs Synthese directe
ExempleSoit le systeme a commander de fonction de transfert :
G(s) =6
s2 + 2s+ 8.
Simulons sa reponse indicielle :
Step Response
Time (sec)
Am
plitu
de
0 1 2 3 4 5 60
0.2
0.4
0.6
0.8
1
1.2
Depassement = 30% ; temps de reponse a 5% = 2.78s ; erreur de position = 25%.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 113 / 131
![Page 236: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/236.jpg)
Synthese de correcteurs Synthese directe
ExempleSoit le systeme a commander de fonction de transfert :
G(s) =6
s2 + 2s+ 8.
Simulons sa reponse indicielle :
Step Response
Time (sec)
Am
plitu
de
0 1 2 3 4 5 60
0.2
0.4
0.6
0.8
1
1.2
Depassement = 30% ; temps de reponse a 5% = 2.78s ; erreur de position = 25%.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 113 / 131
![Page 237: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/237.jpg)
Synthese de correcteurs Synthese directe
Choisissons un modele de reference du second ordre :
Fd(s) =1
s2
ω2n
+ 2ζωns+ 1
.
L’amortissement ζ est fixe a 0.8 de sorte que le depassement soit inferieur a2%.
La pulsation propre ωn est fixee a 1.87 de sorte que le temps de reponse soitde l’ordre de 2s.(formule approximative : tr5% ' 3
ζωn)
La FTBF desiree s’ecrit donc :
Fd(s) =1
0.284s2 + 0.853s+ 1.
Enfin, nous pouvons en deduire le correcteur :
C(s) =s2 + 2s+ 86ω2ns2 + 12ζ
ωns
=s2 + 2s+ 8
1.707s2 + 5.12s
Y. Ariba - Icam, Toulouse. Automatique Lineaire 114 / 131
![Page 238: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/238.jpg)
Synthese de correcteurs Synthese directe
Choisissons un modele de reference du second ordre :
Fd(s) =1
s2
ω2n
+ 2ζωns+ 1
.
L’amortissement ζ est fixe a 0.8 de sorte que le depassement soit inferieur a2%.
La pulsation propre ωn est fixee a 1.87 de sorte que le temps de reponse soitde l’ordre de 2s.(formule approximative : tr5% ' 3
ζωn)
La FTBF desiree s’ecrit donc :
Fd(s) =1
0.284s2 + 0.853s+ 1.
Enfin, nous pouvons en deduire le correcteur :
C(s) =s2 + 2s+ 86ω2ns2 + 12ζ
ωns
=s2 + 2s+ 8
1.707s2 + 5.12s
Y. Ariba - Icam, Toulouse. Automatique Lineaire 114 / 131
![Page 239: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/239.jpg)
Synthese de correcteurs Synthese directe
Choisissons un modele de reference du second ordre :
Fd(s) =1
s2
ω2n
+ 2ζωns+ 1
.
L’amortissement ζ est fixe a 0.8 de sorte que le depassement soit inferieur a2%.
La pulsation propre ωn est fixee a 1.87 de sorte que le temps de reponse soitde l’ordre de 2s.(formule approximative : tr5% ' 3
ζωn)
La FTBF desiree s’ecrit donc :
Fd(s) =1
0.284s2 + 0.853s+ 1.
Enfin, nous pouvons en deduire le correcteur :
C(s) =s2 + 2s+ 86ω2ns2 + 12ζ
ωns
=s2 + 2s+ 8
1.707s2 + 5.12s
Y. Ariba - Icam, Toulouse. Automatique Lineaire 114 / 131
![Page 240: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/240.jpg)
Synthese de correcteurs Synthese directe
Simulons la reponse indicielle du systeme asservi par notre correcteur :
F (s) =C(s)G(s)
1 + C(s)G(s)
Step Response
Time (sec)
Am
plitu
de
0 1 2 3 4 5 60
0.2
0.4
0.6
0.8
1
1.2
Depassement = 1.51% ; temps de reponse a 5% = 1.81s ; erreur de position = 0%.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 115 / 131
![Page 241: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/241.jpg)
Synthese de correcteurs Synthese directe
Simulons la reponse indicielle du systeme asservi par notre correcteur :
F (s) =C(s)G(s)
1 + C(s)G(s)
Step Response
Time (sec)
Am
plitu
de
0 1 2 3 4 5 60
0.2
0.4
0.6
0.8
1
1.2
Depassement = 1.51% ; temps de reponse a 5% = 1.81s ; erreur de position = 0%.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 115 / 131
![Page 242: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/242.jpg)
Synthese de correcteurs Action proportionnelle
Action proportionnelle
Ce type de correcteur est un simple amplificateur de gain reglable :
u(t) = kpε(t) ⇒ C(s) = kp.
10−1
100
101
102
−270
−180
−90
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
−50
0
Mag
nitu
de (
dB)
Avantages et Inconvenients :
Si kp , le systeme repond plus rapidement, l’erreur de position diminue.
Si kp , les marges de stabilite augmentent.
Si kp , peut entrainer des oscillations et un depassement prejudiciable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 116 / 131
![Page 243: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/243.jpg)
Synthese de correcteurs Action proportionnelle
Action proportionnelle
Ce type de correcteur est un simple amplificateur de gain reglable :
u(t) = kpε(t) ⇒ C(s) = kp.
10−1
100
101
102
−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
−70
−60
−50
−40
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
Avantages et Inconvenients :
Si kp , le systeme repond plus rapidement, l’erreur de position diminue.
Si kp , les marges de stabilite augmentent.
Si kp , peut entrainer des oscillations et un depassement prejudiciable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 116 / 131
![Page 244: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/244.jpg)
Synthese de correcteurs Action proportionnelle
Action proportionnelle
Ce type de correcteur est un simple amplificateur de gain reglable :
u(t) = kpε(t) ⇒ C(s) = kp.
Bode Diagram
Frequency (rad/sec)
−70
−60
−50
−40
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
10−1
100
101
102
−270
−180
−90
0
Pha
se (
deg)
Avantages et Inconvenients :
Si kp , le systeme repond plus rapidement, l’erreur de position diminue.
Si kp , les marges de stabilite augmentent.
Si kp , peut entrainer des oscillations et un depassement prejudiciable.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 116 / 131
![Page 245: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/245.jpg)
Synthese de correcteurs Action proportionnelle
Exemple :Soit le systeme a commander G(s) = 1
s2+s+1et le systeme de commande C(s) = kp. Le
systeme en boucle fermee s’exprime par
F (s) =kp
s2 + s+ 1 + kp.
Cette nouvelle fonction de transfert est caracterisee par
Kstatique =kp
kp + 1, ωn =
√1 + kp, ζ =
1
2√
1 + k.
Donc si kp alors ζ , les marges et l’erreur statique εs = 11+kp
.
−80
−60
−40
−20
0
20
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10−2
10−1
100
101
102
−180
−135
−90
−45
0
Pha
se (
deg)
kp=0.5
0 2 4 6 8 10 120
0.2
0.4
0.6
0.8
1
1.2
Step Response
Time (sec)
Am
plitu
de
kp=0.5
Y. Ariba - Icam, Toulouse. Automatique Lineaire 117 / 131
![Page 246: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/246.jpg)
Synthese de correcteurs Action proportionnelle
Exemple :Soit le systeme a commander G(s) = 1
s2+s+1et le systeme de commande C(s) = kp. Le
systeme en boucle fermee s’exprime par
F (s) =kp
s2 + s+ 1 + kp.
Cette nouvelle fonction de transfert est caracterisee par
Kstatique =kp
kp + 1, ωn =
√1 + kp, ζ =
1
2√
1 + k.
Donc si kp alors ζ , les marges et l’erreur statique εs = 11+kp
.
−80
−60
−40
−20
0
20
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10−2
10−1
100
101
102
−180
−135
−90
−45
0
Pha
se (
deg)
kp=0.5
0 2 4 6 8 10 120
0.2
0.4
0.6
0.8
1
1.2
Step Response
Time (sec)
Am
plitu
de
kp=0.5
Y. Ariba - Icam, Toulouse. Automatique Lineaire 117 / 131
![Page 247: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/247.jpg)
Synthese de correcteurs Action proportionnelle
Exemple :Soit le systeme a commander G(s) = 1
s2+s+1et le systeme de commande C(s) = kp. Le
systeme en boucle fermee s’exprime par
F (s) =kp
s2 + s+ 1 + kp.
Cette nouvelle fonction de transfert est caracterisee par
Kstatique =kp
kp + 1, ωn =
√1 + kp, ζ =
1
2√
1 + k.
Donc si kp alors ζ , les marges et l’erreur statique εs = 11+kp
.
−80
−60
−40
−20
0
20
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10−2
10−1
100
101
102
−180
−135
−90
−45
0
Pha
se (
deg)
kp=0.5
0 2 4 6 8 10 120
0.2
0.4
0.6
0.8
1
1.2
Step Response
Time (sec)
Am
plitu
dekp=0.5
Y. Ariba - Icam, Toulouse. Automatique Lineaire 117 / 131
![Page 248: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/248.jpg)
Synthese de correcteurs Action proportionnelle
Exemple :Soit le systeme a commander G(s) = 1
s2+s+1et le systeme de commande C(s) = kp. Le
systeme en boucle fermee s’exprime par
F (s) =kp
s2 + s+ 1 + kp.
Cette nouvelle fonction de transfert est caracterisee par
Kstatique =kp
kp + 1, ωn =
√1 + kp, ζ =
1
2√
1 + k.
Donc si kp alors ζ , les marges et l’erreur statique εs = 11+kp
.
−80
−60
−40
−20
0
20
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
kp=1
kp=0.5
0 2 4 6 8 10 12 140
0.2
0.4
0.6
0.8
1
1.2
Step Response
Time (sec)
Am
plitu
dekp=0.5
kp=1
Y. Ariba - Icam, Toulouse. Automatique Lineaire 117 / 131
![Page 249: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/249.jpg)
Synthese de correcteurs Action proportionnelle
Exemple :Soit le systeme a commander G(s) = 1
s2+s+1et le systeme de commande C(s) = kp. Le
systeme en boucle fermee s’exprime par
F (s) =kp
s2 + s+ 1 + kp.
Cette nouvelle fonction de transfert est caracterisee par
Kstatique =kp
kp + 1, ωn =
√1 + kp, ζ =
1
2√
1 + k.
Donc si kp alors ζ , les marges et l’erreur statique εs = 11+kp
.
−80
−70
−60
−50
−40
−30
−20
−10
0
10
20
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−210
−110
010
110
2−180
−135
−90
−45
0
Pha
se (
deg)
kp=1
kp=0.5
kp=5
0 2 4 6 8 10 12 140
0.2
0.4
0.6
0.8
1
1.2
Step Response
Time (sec)
Am
plitu
de
kp=1
kp=0.5
kp=5
Y. Ariba - Icam, Toulouse. Automatique Lineaire 117 / 131
![Page 250: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/250.jpg)
Synthese de correcteurs Action integrale
Action integrale
Ce correcteur introduit un integrateur qui ajoute un pole nul a la fonction de transferten boucle ouverte :
u(t) =
∫ t
0
ε(t)dt ⇒ C(s) =1
s.
−60
−50
−40
−30
−20
−10
0
10
20
30
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−110
010
110
2−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Avantages et Inconvenients :
Erreur statique nulle en boucle fermee (BF),
Diagramme de phase decale de −90 : marges de stabilites degradees,
Peut provoquer des oscillations et du depassement.Y. Ariba - Icam, Toulouse. Automatique Lineaire 118 / 131
![Page 251: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/251.jpg)
Synthese de correcteurs Action integrale
Action integrale
Ce correcteur introduit un integrateur qui ajoute un pole nul a la fonction de transferten boucle ouverte :
u(t) =
∫ t
0
ε(t)dt ⇒ C(s) =1
s.
−60
−50
−40
−30
−20
−10
0
10
20
30
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−110
010
110
2−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Avantages et Inconvenients :
Erreur statique nulle en boucle fermee (BF),
Diagramme de phase decale de −90 : marges de stabilites degradees,
Peut provoquer des oscillations et du depassement.Y. Ariba - Icam, Toulouse. Automatique Lineaire 118 / 131
![Page 252: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/252.jpg)
Synthese de correcteurs Action integrale
Exemple :Soit le systeme a commander G(s) = 0.8
s2+s+1et le systeme de commande C(s) = 1/s.
Le systeme en boucle fermee s’exprime par
F (p) =0.8
s3 + s2 + s+ 0.8.
F est un ordre 3 : la marge de gain est finie,Courbe plus proche du point critique : marges degradees, oscillations et depassement,Pas d’erreur statique. schema de gauche : bode de la BO ; schema de droite : reponseindicielle de la BF.
−80
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)0 10 20 30 40 50 60 70
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6Step Response
Time (sec)
Am
plitu
de
Y. Ariba - Icam, Toulouse. Automatique Lineaire 119 / 131
![Page 253: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/253.jpg)
Synthese de correcteurs Action integrale
Exemple :Soit le systeme a commander G(s) = 0.8
s2+s+1et le systeme de commande C(s) = 1/s.
Le systeme en boucle fermee s’exprime par
F (p) =0.8
s3 + s2 + s+ 0.8.
F est un ordre 3 : la marge de gain est finie,Courbe plus proche du point critique : marges degradees, oscillations et depassement,Pas d’erreur statique. schema de gauche : bode de la BO ; schema de droite : reponseindicielle de la BF.
−80
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)0 10 20 30 40 50 60 70
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6Step Response
Time (sec)
Am
plitu
de
Y. Ariba - Icam, Toulouse. Automatique Lineaire 119 / 131
![Page 254: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/254.jpg)
Synthese de correcteurs Action integrale
Exemple :Soit le systeme a commander G(s) = 0.8
s2+s+1et le systeme de commande C(s) = 1/s.
Le systeme en boucle fermee s’exprime par
F (p) =0.8
s3 + s2 + s+ 0.8.
F est un ordre 3 : la marge de gain est finie,Courbe plus proche du point critique : marges degradees, oscillations et depassement,Pas d’erreur statique. schema de gauche : bode de la BO ; schema de droite : reponseindicielle de la BF.
−80
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)0 10 20 30 40 50 60 70
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6Step Response
Time (sec)
Am
plitu
de
Y. Ariba - Icam, Toulouse. Automatique Lineaire 119 / 131
![Page 255: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/255.jpg)
Synthese de correcteurs Action integrale
Exemple :Soit le systeme a commander G(s) = 0.8
s2+s+1et le systeme de commande C(s) = 1/s.
Le systeme en boucle fermee s’exprime par
F (p) =0.8
s3 + s2 + s+ 0.8.
F est un ordre 3 : la marge de gain est finie,Courbe plus proche du point critique : marges degradees, oscillations et depassement,Pas d’erreur statique. schema de gauche : bode de la BO ; schema de droite : reponseindicielle de la BF.
−80
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
avec correction intégrale
0 10 20 30 40 50 60 700
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6Step Response
Time (sec)
Am
plitu
de
avec correction intégrale
Y. Ariba - Icam, Toulouse. Automatique Lineaire 119 / 131
![Page 256: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/256.jpg)
Synthese de correcteurs Action derivee
Action derivee
Ce correcteur introduit un derivateur qui ajoute un zero nul a la fonction de transferten boucle ouverte :
u(t) =dε(t)
dt⇒ C(s) = s.
−50
−40
−30
−20
−10
0
10
20
30
40
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−180
−135
−90
−45
0
45
90
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Avantages et Inconvenients :
Ajoute de la phase : marges susceptibles d’etre augmentees,
A tendance a accelerer la reponse du systeme : temps de monte diminue,
Le gain statique est diminue : precision degradee.Y. Ariba - Icam, Toulouse. Automatique Lineaire 120 / 131
![Page 257: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/257.jpg)
Synthese de correcteurs Action derivee
Action derivee
Ce correcteur introduit un derivateur qui ajoute un zero nul a la fonction de transferten boucle ouverte :
u(t) =dε(t)
dt⇒ C(s) = s.
−50
−40
−30
−20
−10
0
10
20
30
40
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−180
−135
−90
−45
0
45
90
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
Avantages et Inconvenients :
Ajoute de la phase : marges susceptibles d’etre augmentees,
A tendance a accelerer la reponse du systeme : temps de monte diminue,
Le gain statique est diminue : precision degradee.Y. Ariba - Icam, Toulouse. Automatique Lineaire 120 / 131
![Page 258: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/258.jpg)
Synthese de correcteurs Combinaisons d’actions
Correcteur Proportionnel-Integral (PI)
Ce correcteur a pour objectif de tirer profit des avantages de l’effet de I sans sesinconvenients :
C(s) = kp1 + τis
τis.
1/ω
ω
φ(ω)
20 log(kp)
τ i
20 logC(ω)
−90
Idee du correcteur :
Utiliser l’avantage de l’integrateur en basses frequences : precision infinie,
L’action integrale ne doit plus avoir d’effet dans les frequences elevees, enparticulier dans la region du point critique,
Y. Ariba - Icam, Toulouse. Automatique Lineaire 121 / 131
![Page 259: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/259.jpg)
Synthese de correcteurs Combinaisons d’actions
Reglage intuitif du correcteur :
Ajuster le gain proportionnel en fonction de l’objectif et des caracteristiques dusysteme avant correction :
le diminuer pour augmenter les marges de stabilite,l’augmenter pour ameliorer la rapidite du systeme.
Ensuite, la partie I est ajoutee en reglant le zero 1/τi de facon a ce que lacorrection ne se fasse qu’en basses frequences.
−100
−80
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
10−1
100
101
102
−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
sans correction
Y. Ariba - Icam, Toulouse. Automatique Lineaire 122 / 131
![Page 260: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/260.jpg)
Synthese de correcteurs Combinaisons d’actions
Reglage intuitif du correcteur :
Ajuster le gain proportionnel en fonction de l’objectif et des caracteristiques dusysteme avant correction :
le diminuer pour augmenter les marges de stabilite,l’augmenter pour ameliorer la rapidite du systeme.
Ensuite, la partie I est ajoutee en reglant le zero 1/τi de facon a ce que lacorrection ne se fasse qu’en basses frequences.
−100
−80
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−110
010
110
2−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
sans correction
avec la partie"proportionnel" seul
Y. Ariba - Icam, Toulouse. Automatique Lineaire 122 / 131
![Page 261: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/261.jpg)
Synthese de correcteurs Combinaisons d’actions
Reglage intuitif du correcteur :
Ajuster le gain proportionnel en fonction de l’objectif et des caracteristiques dusysteme avant correction :
le diminuer pour augmenter les marges de stabilite,l’augmenter pour ameliorer la rapidite du systeme.
Ensuite, la partie I est ajoutee en reglant le zero 1/τi de facon a ce que lacorrection ne se fasse qu’en basses frequences.
−100
−80
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−310
−210
−110
010
110
2−270
−225
−180
−135
−90
−45
0
Pha
se (
deg)
avec le PI "complet"
sans correction
avec la partie"proportionnel" seul
Y. Ariba - Icam, Toulouse. Automatique Lineaire 122 / 131
![Page 262: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/262.jpg)
Synthese de correcteurs Combinaisons d’actions
Correcteur a Avance de Phase
Ce correcteur a pour objectif d’apporter de la phase autour du point critique afind’augmenter les marges de stabilite
C(s) = kp1 + aTs
1 + Ts, a > 1.
ω
φ(ω)
20 logC(ω)
1/aT 1/T
20 log(kp)
20 log(kp.a)
90
mφ
ωm
ω
Y. Ariba - Icam, Toulouse. Automatique Lineaire 123 / 131
![Page 263: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/263.jpg)
Synthese de correcteurs Combinaisons d’actions
Le principe repose sur le reglage de a et T tels que le correcteur apporte de la phase(plage [ 1
aT ,1T ]) autour du point critique.
Methode de reglage du correcteur :
Regler le gain kp pour ajuster la precision ou la rapidite ou les marges.
Mesurer la marge de phase (apres la correction prop. kp) et en deduire la quantitede phase necessaire
φm = ∆φdesiree −∆φ.
A partir de φm, a peut etre calcule
a =1 + sinφm
1− sinφm
Enfin, nous avons la relation
ωm =1
T√a,
il s’agit alors de calculer T tel que ωm coıncide avec ω0db (apres correction prop.) :
T =1
ω0db√a.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 124 / 131
![Page 264: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/264.jpg)
Synthese de correcteurs Combinaisons d’actions
Le principe repose sur le reglage de a et T tels que le correcteur apporte de la phase(plage [ 1
aT ,1T ]) autour du point critique.
Methode de reglage du correcteur :
Regler le gain kp pour ajuster la precision ou la rapidite ou les marges.
Mesurer la marge de phase (apres la correction prop. kp) et en deduire la quantitede phase necessaire
φm = ∆φdesiree −∆φ.
A partir de φm, a peut etre calcule
a =1 + sinφm
1− sinφm
Enfin, nous avons la relation
ωm =1
T√a,
il s’agit alors de calculer T tel que ωm coıncide avec ω0db (apres correction prop.) :
T =1
ω0db√a.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 124 / 131
![Page 265: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/265.jpg)
Synthese de correcteurs Combinaisons d’actions
Exemple :
Soit la fonction de transfert
G(s) =100
(1 + s)2
On souhaite que le systeme en boucle fermee ait une marge de phase Mφ = 45 :
Calculons la marge de phase avant correction
G(ω) = 100
(1+ω20db
)= 1 ⇒ ω0db = 9.95rad/s
Mφ = π − 2arctanω0db = 11
La marge est insuffisante, il faut donc remonter la phase de φm = 34 a lapulsation ω0db
a =1 + sinφm
1− sinφm=
1 + sin34
1− sin34= 3.54
Puis, on regle T afin d’ajouter la quantite φm au bon endroit
1
T√a
= ω0db ⇒ T = 0.053
Souvent, on choisit φm = 1.2φnecessaire pour compenser le decalage de ω0db aprescorrection.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 125 / 131
![Page 266: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/266.jpg)
Synthese de correcteurs Combinaisons d’actions
Exemple :
Soit la fonction de transfert
G(s) =100
(1 + s)2
On souhaite que le systeme en boucle fermee ait une marge de phase Mφ = 45 :
Calculons la marge de phase avant correction
G(ω) = 100
(1+ω20db
)= 1 ⇒ ω0db = 9.95rad/s
Mφ = π − 2arctanω0db = 11
La marge est insuffisante, il faut donc remonter la phase de φm = 34 a lapulsation ω0db
a =1 + sinφm
1− sinφm=
1 + sin34
1− sin34= 3.54
Puis, on regle T afin d’ajouter la quantite φm au bon endroit
1
T√a
= ω0db ⇒ T = 0.053
Souvent, on choisit φm = 1.2φnecessaire pour compenser le decalage de ω0db aprescorrection.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 125 / 131
![Page 267: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/267.jpg)
Synthese de correcteurs Combinaisons d’actions
Exemple :
Soit la fonction de transfert
G(s) =100
(1 + s)2
On souhaite que le systeme en boucle fermee ait une marge de phase Mφ = 45 :
Calculons la marge de phase avant correction
G(ω) = 100
(1+ω20db
)= 1 ⇒ ω0db = 9.95rad/s
Mφ = π − 2arctanω0db = 11
La marge est insuffisante, il faut donc remonter la phase de φm = 34 a lapulsation ω0db
a =1 + sinφm
1− sinφm=
1 + sin34
1− sin34= 3.54
Puis, on regle T afin d’ajouter la quantite φm au bon endroit
1
T√a
= ω0db ⇒ T = 0.053
Souvent, on choisit φm = 1.2φnecessaire pour compenser le decalage de ω0db aprescorrection.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 125 / 131
![Page 268: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/268.jpg)
Synthese de correcteurs Combinaisons d’actions
Exemple :
Soit la fonction de transfert
G(s) =100
(1 + s)2
On souhaite que le systeme en boucle fermee ait une marge de phase Mφ = 45 :
Calculons la marge de phase avant correction
G(ω) = 100
(1+ω20db
)= 1 ⇒ ω0db = 9.95rad/s
Mφ = π − 2arctanω0db = 11
La marge est insuffisante, il faut donc remonter la phase de φm = 34 a lapulsation ω0db
a =1 + sinφm
1− sinφm=
1 + sin34
1− sin34= 3.54
Puis, on regle T afin d’ajouter la quantite φm au bon endroit
1
T√a
= ω0db ⇒ T = 0.053
Souvent, on choisit φm = 1.2φnecessaire pour compenser le decalage de ω0db aprescorrection.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 125 / 131
![Page 269: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/269.jpg)
Synthese de correcteurs Combinaisons d’actions
Exemple :
Soit la fonction de transfert
G(s) =100
(1 + s)2
On souhaite que le systeme en boucle fermee ait une marge de phase Mφ = 45 :
Calculons la marge de phase avant correction
G(ω) = 100
(1+ω20db
)= 1 ⇒ ω0db = 9.95rad/s
Mφ = π − 2arctanω0db = 11
La marge est insuffisante, il faut donc remonter la phase de φm = 34 a lapulsation ω0db
a =1 + sinφm
1− sinφm=
1 + sin34
1− sin34= 3.54
Puis, on regle T afin d’ajouter la quantite φm au bon endroit
1
T√a
= ω0db ⇒ T = 0.053
Souvent, on choisit φm = 1.2φnecessaire pour compenser le decalage de ω0db aprescorrection.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 125 / 131
![Page 270: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/270.jpg)
Synthese de correcteurs Combinaisons d’actions
Schema de gauche : bode de la BO ; schema de droite : reponse indicielle de la BF
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
10−1
100
101
102
103
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
sans correction
0 1 2 3 4 5 60
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8Step Response
Time (sec)
Am
plitu
de
Y. Ariba - Icam, Toulouse. Automatique Lineaire 126 / 131
![Page 271: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/271.jpg)
Synthese de correcteurs Combinaisons d’actions
Schema de gauche : bode de la BO ; schema de droite : reponse indicielle de la BF
−60
−40
−20
0
20
40
Mag
nitu
de (
dB)
10−1
100
101
102
103
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
sans correction
avec correction
0 1 2 3 4 5 60
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8Step Response
Time (sec)
Am
plitu
de
Y. Ariba - Icam, Toulouse. Automatique Lineaire 126 / 131
![Page 272: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/272.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Proportionnel-Integral-Derive (PID)
Ce correcteur PID est une combinaison des actions PI et PD
C(s) = kp
(1 + Tis
Tis
)(1 + Tds).
ω
ω
20 log(kp)
20 logC(ω)
φ(ω)
−90
1/τi
90
1/τd
Idee du correcteur :
Ajouter du gain en bf (precision),
Ajouter de la phase pres du pointcritique (Mφ).
Le correcteur est propre,
Attenuer l’effet du bruit (moins degain en hf).
Y. Ariba - Icam, Toulouse. Automatique Lineaire 127 / 131
![Page 273: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/273.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Proportionnel-Integral-Derive (PID)
Ce correcteur PID est une combinaison des actions PI et PD + un filtre
C(s) = kp
(1 + Tis
Tis
)(1 + Tds)
1
(1 + Tfs).
ω20 log(kp)
20 logC(ω)
φ(ω)
−90
1/τi
90
1/τd
ω
1/τf
Idee du correcteur :
Ajouter du gain en bf (precision),
Ajouter de la phase pres du pointcritique (Mφ).
Le correcteur est propre,
Attenuer l’effet du bruit (moins degain en hf).
Y. Ariba - Icam, Toulouse. Automatique Lineaire 127 / 131
![Page 274: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/274.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Proportionnel-Integral-Derive (PID)
Methode de reglage du correcteur :
Etudier le systeme en BO et ses caracteristiques
(marges, precision, rapidite...).
Regler le gain proportionnel kp afin d’obtenir un premier asservissementsatisfaisant
(en termes de marges, depassement, oscillations, rapidite).
Regler la constante de temps de l’action integrale de sorte qu’elle n’agisse qu’en bf
(pour ne pas penaliser les marges).
Regler la constante de temps de l’action derivee de sorte qu’elle n’agisse qu’autourdu point critique
(pour ne pas penaliser la precision).
Regler la constante de temps du filtre de sorte qu’il n’enleve pas de phase pres dupoint critique,
L’ordre du filtre peut etre augmente afin de mieux attenuer le bruit en hf.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 128 / 131
![Page 275: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/275.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Proportionnel-Integral-Derive (PID)
Methode de reglage du correcteur :
Etudier le systeme en BO et ses caracteristiques
(marges, precision, rapidite...).
Regler le gain proportionnel kp afin d’obtenir un premier asservissementsatisfaisant
(en termes de marges, depassement, oscillations, rapidite).
Regler la constante de temps de l’action integrale de sorte qu’elle n’agisse qu’en bf
(pour ne pas penaliser les marges).
Regler la constante de temps de l’action derivee de sorte qu’elle n’agisse qu’autourdu point critique
(pour ne pas penaliser la precision).
Regler la constante de temps du filtre de sorte qu’il n’enleve pas de phase pres dupoint critique,
L’ordre du filtre peut etre augmente afin de mieux attenuer le bruit en hf.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 128 / 131
![Page 276: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/276.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Exemple : Soit le systeme G(s) =5
s2 + s+ 1.
Une etude de ce systeme nous montre qu’il est stable en BF avec une erreur deposition εs = 16.7%, un depassement de 51% et une marge de phase Mφ = 27.8 aω0db = 2.33rad/s.
On choisit un premier correcteur proportionnel afin de baisser l’erreur et accelererle systeme : on prend kp = 2 pour avoir une erreur de 9%. Il en resulte une nouvellemarge de phase Mφ = 18.9 a ω0db = 3.23rad/s et un depassement de 61.9%.
On place la constante de temps du I avant le point critique Ti = 3 1/3.23. Le
correcteur devient C(s) = 21 + 3s
3s.
L’erreur en regime permanent est maintenant nul mais la nouvelle marge de phaseest de Mφ = 12.9 a ω0db = 3.24rad/s. Il s’agit ensuite de placer la constante detemps de la partie D avant le point critique de facon a laisser le temps aucorrecteur de remonter la phase jusqu’a obtenir une margesatisfaisante :Ti Td = 0.4 > 1/3.24. On obtient Mφ = 70.4.
Pour finir, afin de synthetiser un correcteur realisable et pour mieux filtrer lebruit, on ajoute un filtre avec une constante de temps relativement faibleTf = 0.1 < 1/3.24. Il en resulte Mφ = 45.9, ce qui reste satisfaisant.
Y. Ariba - Icam, Toulouse. Automatique Lineaire 129 / 131
![Page 277: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/277.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Diagramme de Bode de la fonction de transfert en boucle ouverte (C(s)G(s))
−60
−40
−20
0
20
40
60
Mag
nitu
de (
dB)
10−2
10−1
100
101
102
−180
−135
−90
−45
0
Pha
se (
deg)
Bode Diagram
Frequency (rad/sec)
C(p)=2
Y. Ariba - Icam, Toulouse. Automatique Lineaire 130 / 131
![Page 278: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/278.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Diagramme de Bode de la fonction de transfert en boucle ouverte (C(s)G(s))
−60
−40
−20
0
20
40
60
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−210
−110
010
110
2−180
−135
−90
−45
0
Pha
se (
deg)
C(p)=2
C(p)=2(3p+1)/(3p)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 130 / 131
![Page 279: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/279.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Diagramme de Bode de la fonction de transfert en boucle ouverte (C(s)G(s))
−60
−40
−20
0
20
40
60
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−210
−110
010
110
2−180
−135
−90
−45
0
Pha
se (
deg)
C(p)=2
C(p)=2(3p+1)/(3p)
C(p)=2(3p+1)(0.4p+1)/(3p)
Y. Ariba - Icam, Toulouse. Automatique Lineaire 130 / 131
![Page 280: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/280.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Diagramme de Bode de la fonction de transfert en boucle ouverte (C(s)G(s))
−60
−40
−20
0
20
40
60
Mag
nitu
de (
dB)
Bode Diagram
Frequency (rad/sec)10
−210
−110
010
110
210
3−180
−135
−90
−45
0
Pha
se (
deg)
C(p)=2
C(p)=2(3p+1)/(3p)
C(p)=2(3p+1)(0.4p+1)/(3p)
PID complet
Y. Ariba - Icam, Toulouse. Automatique Lineaire 130 / 131
![Page 281: Yassine Ariba Dpt GEI - Icam, Toulouse.homepages.laas.fr/yariba/enseignement/slides-regul.pdf · Objectif du cours 4 Stabilit e et performances d’un asservissement Stabilit e Pr](https://reader031.fdocument.pub/reader031/viewer/2022022020/5b9d41d109d3f2ed218bd4d3/html5/thumbnails/281.jpg)
Synthese de correcteurs Proportionnel-Integral-Derive (PID)
Reponses temporelles de la BF a un echelon unite
0 5 10 150
0.5
1
1.5Step Response
Time (sec)
Am
plitu
de
C(p)=2
0 5 100
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8Step Response
Time (sec)
Am
plitu
de
C(p)=2(3p+1)/3p
0 5 100
0.5
1
1.5Step Response
Time (sec)
Am
plitu
de C(p)=2(3p+1)(0.4p+1)/3p
0 5 100
0.5
1
1.5Step Response
Time (sec)
Am
plitu
de PID complet
Y. Ariba - Icam, Toulouse. Automatique Lineaire 131 / 131