Chapitre

5 Étude temporelle
des systèmes linéaires
5.1 Systèmes du 1er ordre

5.1.1 Définition

Un système du premier ordre est décrit par une équation différentielle linéaire à coefficients constants
du premier ordre :

ds(t) τ : constante de temps

τ· + s(t) = K.e(t) avec
dt K : gain statique

5.1.2 Fonction de transfert

L L
On pose : e(t) −
→ E(p) et s(t) −
→ S(p). On se place dans les conditions de Heaviside.
En appliquant la transformation de Laplace à l’équation précédente, on obtient :

τ · p · S(p) + S(p) = K · E(p)

On en déduit la fonction de transfert sous sa forme canonique et le schéma-bloc d’un système du

premier ordre :
S(p) K E(p) K S(p)
H(p) = =
E(p) 1+τ·p 1+τ·p

5.1.3 Étude temporelle - Réponse à un échelon

On se propose d’étudier la réponse temporelle de la sortie s(t) pour une entrée de type échelon
e(t) = E0 · H(t) .
Pour déterminer la réponse temporelle de la sortie s(t), il est nécessaire de résoudre l’équation diffé-
rentielle. On peut pour cela, soit utiliser les outils mathématiques classiques, soit utiliser les propriétés
de la transformée de Laplace.
Pour résoudre complètement il est nécessaire de préciser les conditions initiales, nous les suppose-
rons nulles (conditions de Heaviside).
La résolution ne pose pas de problème, on connaît la forme de l’équation différentielle à résoudre,

ds(t)
τ· + s(t) = K.e(t)
dt
et la solution est, pour un système partant du repos, de la forme :
 
−t
s(t) = K · E0 · 1 − e τ 

1
 5 Étude temporelle des systèmes linéaires

Nous allons retrouver ce résultat à partir de la transformation de Laplace.

On détermine dans un premier temps la sortie dans le domaine symbolique :

S(p) = H(p) · E(p)

K E0 L E0
S(p) = · avec e(t) − →
1+τ·p p p
A B
S(p) = + (décomposition en fractions simples)
p 1+τ·p
K · E0 K · E0 · τ
S(p) = − (par identification)
p 1+τ·p

À partir du tableau des transformées page 67 ;

K · E0 L−1
−−→ K · E0 · H(t)
p
−t
K · E0 · τ L−1
−−→ K · E0 · e τ · H(t)
1+τ·p

On retrouve bien la réponse temporelle (figure 5.1).

 
−t
s(t) = K · E0 · 1 − e τ  H(t)

s(t)
s∞ = K · E0
0,95 · K · E0

0,63 · K · E0

t
τ
3τ

FIGURE 5.1 – Caractéristiques graphiques d’un 1er ordre

a ) Étude et propriétés de s(t)

L’étude de cette fonction ne pose aucun problème particulier mais nous allons utiliser les outils de
la transformation de Laplace pour retrouver quelques propriétés.

Asymptote finale : à partir du théorème de la valeur finale 1 , nous allons vérifier que la réponse

1. Le théorème de la valeur finale n’est utilisable que si la fonction est convergente.

2
5.1 Systèmes du 1er ordre

possède une asymptote horizontale. Le théorème de la valeur finale s’écrit :

lim s(t) = lim p · S(p)

t→∞ p→0
( )
K E0
lim s(t) = lim p · ·
t→∞ p→0 1+τ·p p
lim s(t) = K · E0
t→∞

La sortie tend asymptotiquement vers une valeur finie K · E0 lorsque t → ∞.

Valeur initiale : le théorème de la valeur initiale appliqué à s(t) permet d’obtenir la valeur à l’instant
t = 0 de la sortie.

lim s(t) = lim p · S(p)

t→0 p→∞
( )
K E0
lim s(t) = lim p · ·
t→0 p→∞ 1+τ·p p
lim s(t) = 0
t→0

Tangente à l’origine : pour déterminer la tangente à l’origine, il faut déterminer la dérivée de s(t)
pour t = 0. On peut utiliser le théorème de la valeur initiale appliqué à la dérivée pour obtenir
cette valeur.
Nous savons que pour une fonction temporelle f(t) dont la transformée de Laplace est F(p), on
peut écrire :
( )
d f(t)
L = p · F(p) − f(0+ )
dt
( )
d f(t)
L = p · F(p) si les conditions initiales sont nulles.
dt

. d s(t)
En appliquant le théorème de la valeur initiale à s(t) =
dt
( )
. d s(t) K E0 K
lim s(t) = lim = lim p · S(p) = lim p ·
2 2
· = · E0
t→0 t→0 dt p→∞ p→∞ 1+τ·p p τ

K · E0 s∞
La pente à l’origine est non nulle et égale à = . On constate que la tangente à l’origine
τ τ
coupe l’asymptote horizontale (K · E0 ) au bout du temps τ (voir figure 5.2a).
Temps de réponse : on appelle temps de réponse à 5% (T5% ), le temps mis pour que la sortie atteigne
la valeur finale à ±5% près soit :

s(∞) − s(T5% )
= 0,05
s(∞)
( −T5%
)
K · E0 − K · E0 · 1 − e τ
= 0,05
K · E0
−T5%
e τ = 0,05

On en déduit

T5% ≈ 3 · τ

On note que le temps de réponse ne dépend pas de la consigne d’entrée ni du gain mais uni-
quement de la constante de temps τ, plus la constante de temps est importante, plus le temps

3
 5 Étude temporelle des systèmes linéaires

de réponse est important.

On retrouve sur les figures 5.1 et 5.2 les différentes caractéristiques de la réponse temporelle à un
échelon d’un système du 1erordre.

s(t) s(t)
K=2

K=1.6

K=1
K=0.8

t t
τ τ3 τ1 τ2
T5% = 3τ 3τ3 3τ1 3τ2
(a) Réponse à un échelon unitaire d’un 1er ordre pour diffé- (b) Réponse à un échelon d’un 1er ordre pour diffé-
rentes valeurs du gain K rentes valeurs de la constante de temps τ

FIGURE 5.2 – Réponse d’un système du 1er ordre

5.1.4 Réponse à une impulsion de Dirac

Pour une impulsion de Dirac (figure 5.3) :

e(t) = δ(t) alors K s(t)
τ K
Tangente en 0 :
t τ2
K K −
S(p) = H1 (p) = ⇒ s(t) = · e τ
1+τ·p τ

Remarque : compte tenu de la disconti-

nuité de l’impulsion de Dirac, le théorème de t
la valeur initiale n’est pas applicable. τ

FIGURE 5.3 – Réponse à une impulsion de Dirac

5.1.5 Réponse à une rampe

Réponse à une rampe e(t) = a · t · H(t) (fi-

gure 5.4) alors s(t)
K a
S(p) = H1 (p) · E(p) = · 2
1+τ·p p
  
t
−
s(t) = K · a · t − τ · 1 − e τ  Asymptote de
pente K · a t
τ
L’asymptote de pente K · a coupe l’axe des
abscisses en (τ, 0). FIGURE 5.4 – Réponse à une rampe de pente a

4
5.2 Systèmes du 2d ordre

5.2 Systèmes du 2d ordre

5.2.1 Définition

Un système du second ordre est décrit par une équation différentielle linéaire à coefficients constants
du second ordre :

d2 s(t) d s(t)
a2 · 2
+ a1 · + a0 · s(t) = b0 · e(t)
dt dt

Lorsque tous les coefficients sont de même signe, on préfère l’écrire sous sa forme canonique :

ωn > 0 : pulsation propre, en rad · s−1

1 d2 s(t) 2 · ξ d s(t)
· + · + s(t) = K · e(t) avec ξ > 0 : : facteur d’amortissement, sans dimension
ω2n dt2 ωn dt
K : gain statique

Remarque : en physique, vous utiliserez une autre représentation canonique :

ωn > 0 : pulsation propre

1 d2 s(t) 1 d s(t) 1
· + · + s(t) = K · e(t) avec Q = : facteur de qualité, sans dimension
ω2n dt 2 Q · ωn dt 2·ξ
K : gain statique

5.2.2 Fonction de transfert

L L
On pose : e(t) −
→ E(p) et s(t) −
→ S(p). On se place dans les conditions de Heaviside.
En appliquant la transformation de Laplace à l’équation précédente, on obtient :

p2 2·ξ
2
· S(p) + · p · S(p) + S(p) = K · E(p)
ωn ωn
( 2 )
p 2·ξ
+ · p + 1 · S(p) = K · E(p)
ω2n ωn

On en déduit la fonction de transfert sous sa forme canonique et le schéma-bloc d’un système du

second ordre :
S(p) K K
H(p) = = E(p) S(p)
E(p) 2·ξ p2 2·ξ p2
1+ ·p+ 2 1+ ·p+ 2
ωn ωn ωn ωn
Le dénominateur de la fonction de transfert admet soit des racines réelles, soit une racine réelle
double, soit deux racines complexes conjuguées.

4 · ξ2 4 4
∆= 2
− 2 = 2 · (ξ2 − 1)
ωn ωn ωn

On constate que le signe de ∆ ne dépend que de ξ :

Cas ξ > 1 alors ∆ > 0 : le dénominateur admet alors deux racines réelles,
( √ ) ( √ )
r1 = ωn · −ξ − ξ2 − 1 et r2 = ωn · −ξ + ξ2 − 1

5
 5 Étude temporelle des systèmes linéaires

On pose :

1 1 1 1
τ1 = − =− ( √ ) et τ2 = − =− ( √ )
r1 ωn · −ξ − ξ2 − 1 r2 ωn · −ξ + ξ2 − 1

soit :
τ1 + τ2 1
ξ= √ et ωn = √
2 · τ1 · τ2 τ1 · τ2

Il est alors préférable d’écrire la fonction de transfert sous la forme d’un produit de deux fonc-
tions du premier ordre :

K
H(p) =
(1 + τ1 · p) · (1 + τ2 · p)

Cas ξ = 1 alors ∆ = 0 : le dénominateur admet alors une racine réelle double,

r = −ωn

La fonction de transfert peut alors se mettre sous la forme :

K K
H(p) = ( )2 =
p (1 + τ · p)2
1+
ωn

Cas 0 < ξ < 1 alors ∆ < 0 : le dénominateur admet alors deux racines complexes conjuguées.
( √ ) ( √ )
r1 = ωn · −ξ − j · 1 − ξ2 et r2 = ωn · −ξ + j · 1 − ξ2

Il est d’usage de laisser la fonction de transfert dans ce cas, sous la forme d’un polynôme du second
degré afin de ne pas faire apparaître les termes imaginaires :

K
H(p) =
2·ξ p2
1+ ·p+ 2
ωn ωn

5.2.3 Réponse temporelle à une entrée en échelon

L’allure de la réponse temporelle va donc dépendre du coefficient d’amortissement ξ.

Soit e(t) = E0 · H(t) avec H(t) la fonction de Heaviside alors :

K E0
S(t) = ·
2·ξ p2 p
1+ ·p+ 2
ωn ωn

On considère que ξ > 0 et ωn > 0.

Valeur finale

6
5.2 Systèmes du 2d ordre

La valeur finale ne dépend que du gain statique K et de l’amplitude de l’échelon.

 
 K E0 
lim s(t) = lim p · S(p) = lim  p· · 
t→∞ p→0 p→0  2·ξ p2 p
1+ ·p+ 2
ωn ωn
lim s(t) = K · E0
t→∞

Valeur initiale
 
 K E0 
lim s(t) = lim p · S(p) = lim  p· · 
t→0 p→∞ p→∞  2·ξ p2 p
1+ ·p+ 2
ωn ωn
lim s(t) = 0
t→0

Tangente à l’origine
La courbe de réponse pour une entrée en échelon débute en 0, avec une tangente horizontale.
 

.  2 K E0 
p ·
lim s(t) = lim p · p · S(p) = lim  · 
t→0 p→∞ p→∞ 2·ξ p2 p
1+ ·p+ 2
ωn ωn
.
lim s(t) = 0
t→0

Les valeurs initiale et finale ne dépendent pas de ξ ni de ωn si ξ > 0 et ωn > 0. L’allure de la

réponse temporelle dépend par contre de la nature des racines du dénominateur de la fonction de
transfert et de la nature des racines réelles ou complexes du dénominateur.

a ) Cas ξ > 1 : deux racines réelles

La fonction de transfert a donc la forme suivante :

K E0
S(p) = ·
(1 + τ1 · p) · (1 + τ2 · p) p

La décomposition en fraction simple s’écrit :

( )
1 τ21 1 τ22 1
S(p) = K · E0 · − · − ·
p τ1 − τ2 1 + τ1 · p τ2 − τ1 1 + τ2 · p

Du tableau des transformées inverses, on déduit :

 
t t
 τ1 − τ2 −

s(t) = K · E0 1 − · e τ1 − · e τ2  · H(t)
τ1 − τ2 τ2 − τ1

On note sur la figure 5.5 :

τ1 + τ2
— l’influence de ξ > 1, plus ξ = √ est grand, plus le temps de réponse est important
2 · τ1 · τ2
pour ωn constant ;
— que la tangente à l’origine est horizontale et la réponse ne présente pas d’oscillation.

7
 5 Étude temporelle des systèmes linéaires

s(t)
K · E0
0,95 · K · E0

ξ>
1 cr
oiss
ant

FIGURE 5.5 – Réponse temporelle d’un second ordre avec ωn constant pour ξ > 1 (tracés continus) et
ξ = 1 (tracé en pointillé) à une entrée en échelon

b ) cas ξ = 1 : une racine double

K E0
S(t) = ·
(1 + τ · p) 2 p

Décomposition en fraction simple :

( )
1 τ τ
S(p) = K · E0 · − −
p (1 + τ · p) 2 (1 + τ · p)

D’où la réponse temporelle :

 
t
t+τ − 
s(t) = K · E0 · 1 − · e τ · H(t)
τ

en revenant à ξ et ωn , la réponse temporelle s’écrit :

( )
s(t) = K · 1 − (1 + ωn · t) · e−ωn ·t · H(t)

c ) Cas 0 < ξ < 1 : deux racines complexes conjuguées

La fonction de transfert s’écrit :
K E0
S(p) = ·
2·ξ p2 p
1+ ·p+ 2
ωn ωn

La réponse temporelle s’écrit alors :

( )
1 ( √ )
s(t) = K · E0 1 − √ e−ξ·ωn ·t · sin ωn 1 − ξ2 · t + φ · H(t) avec φ = arccos ξ
1 − ξ2

On pose en général ωp la pseudo pulsation avec

√
ωp = ωn · 1 − ξ2

8
5.2 Systèmes du 2d ordre

La réponse temporelle devient :

( )
1
s(t) = K · E0 1− √ e−ξωn t · sin (ωp · t + φ) · H(t)
1 − ξ2

que l’on peut aussi mettre sous la forme :

( ( ) )
ξ
s(t) = K · E0 1 − cos (ωp · t) + √ · sin (ωp · t) · e−ξ·ωn ·t · H(t)
1 − ξ2

s(t)

1,05 · K · E0
K · E0
0,95 · K · E0
ξ→
0a
ve
c 1>
ξ>
0

FIGURE 5.6 – Réponse temporelle d’un second ordre avec ωn constant pour 0 < ξ < 1 (tracés continus)
et ξ = 1 (tracé en pointillé) à une entrée en échelon

pseudo-périodeTp
s(t)

1er dépassement

1,05 · K · E0
K·E
0,95 · K · E00

t
Tp
Tpm =
2

FIGURE 5.7 – Réponse temporelle d’un second ordre pour 0 < ξ < 1 à une entrée en échelon

Les caractéristiques de cette réponse temporelle pseudo oscillatoire sont :

2·π
— La pseudo–période : Tp = √ .
ω n · 1 − ξ2

9
 5 Étude temporelle des systèmes linéaires

√
— La pseudo–pulsation : ωp = ωn · 1 − ξ2 .
Tp π
— L’instant du premier maximum : Tpm = = √ .
2 ω n · 1 − ξ2
−π · ξ
√
2
— L’amplitude relative du premier dépassement : D1% = e 1 − ξ .

Dépassements : L’amplitude des dépassements augmente d’autant plus fortement que le coeffi-
cient d’amortissement ξ diminue. Le nombre d’oscillation augmente lui aussi avec la diminution de
ξ. L’abaque page suivante (figure 5.8) permet de déterminer l’amplitude relative de chaque dépasse-
ment en fonction du coefficient d’amortissement.

D
1
0.8 D1
0.7
0.6 D2
0.5 D3
0.4
D4
0.3 D5
D6
0.2
D7
D8

0.1

0.05
D1
1
D3
D4
D2

z
0.01
05

5
6
7
01

1
0.

0.

0.
0.
0.
0.

0.
0.
0.

FIGURE 5.8 – Abaque des dépassements

Temps de réponse : Le temps de réponse à 5% évolue lui aussi de manière notable en fonction du
coefficient d’amortissement ξ.
L’abaque de la figure 5.10 donne le temps de réponse à
5% (en fait T5% · ωn ) pour un système du second ordre. s
ξ ≈ 0, 7
On constate sur cette abaque deux parties : +5%
— pour ξ > 0,7, le temps de réponse augmente lorsque −5%
ξ augmente ;
— pour ξ < 0,7, le temps de réponse augmente lorsque
ξ diminue.
Le temps de réponse est minimal pour ξ ≈ 0,7. Le temps ξ=1
de réponse sans dépassement lui est minimal pour ξ = 1.
−π · ξ
√
2
ξ = 0,7 est la solution de D1% = e 1 − ξ = 0,05. t
La courbe réponse est alors tangente à la limite supé- T r0,7 T r1
rieure de détermination de T5% .
FIGURE 5.9 – Temps de réponse minimal

10
5.3 Intégrateur

Tr · ωn
80

50
K
H(p) = p2
30 2·z
1+ ωn p + ω2n

10

z
0.05

0.1

0.2

0.3

0.4
0.5

0.7

1.5

4
5

10
FIGURE 5.10 – Abaque des temps de réponse d’un second ordre

5.3 Intégrateur

On appelle système intégrateur, un système

décrit par l’équation différentielle du premier s s(t) = K · E0 · t · H(t)
ordre :

d s(t) e(t) = E0 · H(t)

= K · e(t)
dt
∫ +∞
s(t) = K · e(u)du
0 t
La sortie est l’intégrale de l’entrée. FIGURE 5.11 – Réponse temporelle d’un intégrateur

Dans les conditions de Heaviside :

S(p) K E(p) K S(p)
H(p) = =
E(p) p p

5.4 Dérivateur

On appelle système dérivateur, un système décrit par l’équation différentielle du premier ordre :

d e(t)
s(t) = K ·
dt
Dans les conditions de Heaviside :
S(p) E(p) S(p)
H(p) = = K·p K·p
E(p)

11
 5 Étude temporelle des systèmes linéaires

5.5 Retard

On appelle système retard, un système décrit par l’équation différentielle :

s(t) = e(t − τ)

Dans les conditions de Heaviside :

S(p) E(p) S(p)
H(p) = = e−τ·p e−τ·p
E(p)

5.6 Influence de la boucle fermée

On se propose dans cette partie d’évaluer l’influence du bouclage sur les caractéristiques tempo-
relles des systèmes du premier et du second ordre.

Soit un système bouclé à retour unitaire défini E(p) ε(p) S(p)

par le schéma-bloc ci-contre : + T (p)
−

La formule de Black, nous permet de déterminer rapidement la fonction de transfert en boucle

fermée :
S(p) T (p)
BF(p) = =
E(p) 1 + T (p)

5.6.1 1er ordre

Ko
T (p) =
1 + τo · p

avec Ko le gain en boucle ouverte et τo la constante de temps en boucle ouverte.

Alors
Ko
S(p) 1 + τo · p Ko
BF(p) = = =
E(p) Ko 1 + Ko + τo · p
1+
1 + τo · p
Ko
1 + Ko
BF(p) = τo
1+ ·p
1 + Ko

On constate qu’un système du premier ordre bouclé est aussi un système du premier ordre avec les
caractéristiques suivantes :
Ko
— gain en boucle fermée : Kf =
1 + Ko
τo
— constante de temps en boucle fermée : τf =
1 + Ko
Le système bouclé est toujours plus rapide que le système en boucle ouverte, et d’autant plus rapide
que le gain en boucle ouverte est important et le gain tend vers 1 si Ko → ∞

12
5.7 Systèmes complexes

Pour une entrée en échelon e(t) = E0 · H(t), la sortie tend vers

( )
Kf E0
lim s(t) = lim p · S(p) = lim p · · = Kf · E0
t→∞ p→0 p→0 1 + τf · p p
Ko
lim s(t) = · E0
t→∞ 1 + Ko
et l’erreur indicielle εi

εi = lim ε(t) = lim p · ε(p) = lim (p · (E(p) − S(p)))

t→∞ p→0 p→0
( )
Ko 1
εi = ( 1 − · E0 = · E0
1 + Ko 1 + Ko

La précision du système s’améliore donc lorsque Ko → ∞.

5.6.2 2d ordre

Ko
T (p) =
2 · ξo p2
1+ ·p+ 2
ωn0 ωn0

avec Ko le gain en boucle ouverte, ξo le coefficient d’amortissement en boucle ouverte et ωno la pul-
sation propre.
La fonction de transfert en boucle fermée s’écrit sous sa forme canonique :

Ko
1 + Ko
BF(p) =
2ξ0 p2
1+ ·p+
(1 + Ko ) · ωno (1 + Ko ) · ω2no

Le bouclage ne change pas l’ordre du système. Par identification, on obtient les caractéristiques du
second ordre en boucle fermée :
Ko
— gain statique en boucle fermée : Kf =
1 + Ko
√
— pulsation propre en boucle fermée : ωnf = ωno · 1 + Ko > ωno
ξo
— coefficient d’amortissement en boucle fermée : ξf = √ < ξo
1 + Ko
Comme pour le système du premier ordre, le gain statique tend vers 1 lorsque Ko augmente.
On constate aussi que le facteur d’amortissement diminue et que la pulsation propre augmente.
Cela implique qu’un système qui n’était pas oscillant en boucle ouverte peut le devenir en boucle
fermée et s’il était oscillant, l’amplitude des oscillations et leur fréquence augmentent.
En conclusion, le bouclage d’un système modifie ses caractéristiques, il permet d’améliorer la rapi-
dité et la précision mais peut rendre le système instable.

5.7 Systèmes complexes

5.7.1 Systèmes d’ordre supérieur à 2

Soit
K
H(p) =
1 + a1 · p + a2 · p2 + · · · + an · pn

13
 5 Étude temporelle des systèmes linéaires

Un système d’ordre élevé comporte un grand nombre de paramètres et étudier ses réponses tem-
porelles peut se révéler fastidieux.
L’objectif est de pouvoir déterminer le comportement d’un système en fonction de ses réponses
temporelles tout en considérant un nombre restreint de paramètres en le ramenant à un système re-
lativement simple.
Un système d’ordre élevé comporte un grand nombre de pôles (réels ou complexes conjugués) et
tous ne possèdent pas la même influence sur les réponses temporelles.
Soit une fonction de transfert d’ordre supérieur à 2.

S(p) 1
H(p) = = ( p) ( p)
E(p) (1 + p) · 1 + · 1+
5 20
La réponse temporelle s’écrit

25 −t 1 1
s(t) = 1 − ·e − · e−20·t + · e−5·t
19 57 3

Dans cette réponse temporelle, tous les éléments n’ont pas le même poids, le terme e20·t tend beau-
coup plus rapidement vers zéro que e−t et e−5·t , il peut donc être négligé par rapport aux deux autres.
Sur la figure 5.12 sont tracées les réponses temporelle pour un échelon unitaire des fonctions de
transfert
1 1
H2 (p) = ( p ) et H1 (p) = (1 + p)
(1 + p) · 1 +
5
On remarque qu’il n’y a peu de différence en régime permanent. La première simplification est
proche du comportement de H(p), par contre l’erreur pour la seconde simplification est notable au
début.

s s

t t
(a) de H(p) et H2 (p) (b) de H(p) et H1 (p)

FIGURE 5.12 – Comparaison des réponses temporelles

On peut déduire de cet exemple.

— Un système d’ordre élevé possède, la plupart du temps, 1 ou 2 pôles dominants et se « comporte »
donc comme un système du premier ou du deuxième ordre.
— On peut simplifier la fonction de transfert d’un système d’ordre élevé en ne conservant que
le (ou les) pôle(s) dominant(s) (approximation par un système du premier ou du deuxième
ordre).
— En pratique, un pôle peut être négligé dès que sa constante de temps est 4 à 5 fois inférieure à
la constante de temps du pôle dominant.
Cet exemple ne traite que de la simplification d’un système non oscillatoire mais les conclusions
sont généralisables aux systèmes avec un comportement oscillatoire amorti.

14
5.7 Systèmes complexes

5.7.2 Influence d’un zéro

Soit la fonction de transfert possédant un zéro 2 réel :

1+a·p
H(p) = K ·
1 + b1 · p + b2 · p2

Cette fonction peut s’écrire :

( )
1 p
H(p) = K · +a·
1 + b1 · p + b2 · p2 1 + b1 · p + b2 · p2

Elle se construit donc comme la somme de deux fonctions de transfert. La fonction sans le numéra-
teur à laquelle on ajoute la dérivée de cette fonction affectée du coefficient a.
La réponse temporelle à un échelon unitaire est donc :
( )
1 1 p 1
S(p) = K · · +a· ·
1 + b1 · p + b2 · p2 p 1 + b1 · p + b2 · p2 p
( )
1 1
On pose s1 (t) = L−1 · , alors
1 + b1 · p + b2 · p2 p
( )
d s1 (t)
s(t) = K · s1 (t) + a ·
dt

À partir du théorème de la valeur finale, on constate que le terme dérivé ne modifie pas la valeur
finale.

lim s(t) = lim p · S(p)

t→∞ p→0
( ( ))
1 1 p 1
= lim p · K · · +a· · =K
p→0 1 + b1 · p + b2 · p p
2 1 + b1 · p + b2 · p p
2

Le zéro a surtout un effet pendant le régime transitoire.

a ) Cas d’un zéro négatif

s
s(t)
s1 (t)

.
0.5 · s1 (t)
t

FIGURE 5.13 – Influence d’un zéro négatif

Sur la figure 5.13 on a tracé la réponse temporelle à un échelon de la fonction de transfert :

1 + 0.5 · p 1 0.5 · p
H(p) = 2
= 2
+
1+p+p 1+p+p 1 + p + p2
.
On retrouve la réponse temporelle s(t) = s1 (t) + 0.5 · s1 (t).

2. Un zéro est une racine du numérateur.

15
 5 Étude temporelle des systèmes linéaires

On constate sur ce graphe que l’effet du numérateur avec un zéro négatif est d’accélérer la réponse
temporelle mais aussi d’amplifier les dépassements.

b ) Cas d’un zéro positif

Sur la figure 5.14 on a tracé la fonction de transfert :

1 − 0.5 · p 1 0.5 · p
H(p) = 2
= 2
+
1+p+p 1+p+p 1 + p + p2

présentant un pôle positif.

On constate que la réponse temporelle s(t) présente un « undershot » c’est-à-dire un dépassement
vers le bas au démarrage.

s1 (t)
s(t)

.
0.5 · s1 (t)

FIGURE 5.14 – Influence d’un zéro positif

5.7.3 Pôles positifs

Jusqu’ici, nous n’avons traité que le cas de pôles réels négatifs ou de pôles complexes à partie réelle
négative. Quel est l’effet d’un pôle réel positif sur le comportement du système ?
Soit
s t
K
H(p) = avec τ > 0
(1 + τ · p) · (1 − τ · p)

Étudions la réponse temporelle à un échelon :

e(t) = E0 · H(t).

K E0
S(p) = ·
(1 + τ · p) · (1 − τ · p) p FIGURE 5.15 – Influence d’un pôle
réel positif
d’où la décomposition en fractions simples :
( )
1 τ τ
S(p) = K · E0 · − +
p 2 · (1 + τ · p) 2 · (1 − τ · p)

et la réponse temporelle
  
t t
1 −
s(t) = K · E0 · 1 − e τ + e τ  H(t)
2

Cette fonction temporelle diverge (figure 5.15), elle tend −∞.

16
5.7 Systèmes complexes

De la même manière, un système possédant deux pôles

s
complexes conjugués à partie réelle positive diverge aussi
mais en oscillant. Le graphe 5.16 présente l’allure de la
réponse temporelle à un échelon unitaire de :

1 t
H(p) = .
1 − p + 10 · p2

FIGURE 5.16 – Influence de pôles

complexes conjugués à
partie réelle positives]

5.7.4 Réponse temporelle et position des pôles

On peut regrouper les différents résultats sur la figure 5.17 sur laquelle est représentée la réponse
à une impulsion de Dirac (e(t) = δ(t)) en fonction de la nature des pôles de la fonction de transfert.

CONVERGENCE DIVERGENCE

Racine nulle Racine nulle

simple ℑ multiple
Racines com- Racine réelle
plexes à partie 1 positive
réelle négative
0 1

0 1 0
0 1

1
1

0
0 1

0
0 1

0
0 1

0
0 1

Racine réelle 1

0
0 1
1

0
0 1 Racines com-
négative plexes à partie
réelle positive
Racines imagi- Racines imagi-
naires simples naires doubles

FIGURE 5.17 – Réponse temporelle en fonction de la position des pôles

Un critère important de réglage de fonction de transfert d’un système asservi est de s’assurer que
tous les pôles de la fonction de transfert en boucle fermée sont à partie réelle négative.

17
 5 Étude temporelle des systèmes linéaires

5.8 Feuille de travaux dirigés n°5

Exercice 1 - Analyse réponse temporelle

Corrigé page 21

Les courbes figure 5.18, présentent la réponse de 3 systèmes à une entrée e(t) = 10 · u(t) avec u(t)
la fonction de Heaviside.

s(t)

13

12
3
11

10

9
2
8

7 1

t
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8

FIGURE 5.18 – Réponses temporelles

Q1. Déterminer : le temps de réponse à 5%, le gain statique, le dépassement relatif s’il existe pour
chacune des courbes réponses.
Q2. Proposer une fonction de transfert pour chacune des représentations, déterminer les coefficients
caractéristiques des réponses temporelles n°2 et n°3.
Q3. Que proposez-vous pour la réponse temporelle n°1 ?

Exercice 2 - Identification BF-BO

Corrigé page 21

Un système asservi est décrit par le schéma-bloc :

E(p) ϵ(p) U(p) S(p)
+ K O(p)
−

— Le gain K est un gain proportionnel réglé à K = 10.

— O(p) est la fonction de transfert de l’actionneur du système. Cette fonction n’est pas connue.

18
5.8 Feuille de travaux dirigés n°5

S(p)
Q1. Déterminer la fonction de transfert en boucle fermée F(p) = en fonction de K et O(p).
E(p)
On cherche à déterminer expérimentalement la fonction de transfert de l’actionneur. Pour cela on
soumet le système asservi à un échelon d’amplitude 4 volt : e(t) = 4 · H(t) avec H(t) la fonction de
Heaviside. La figure 5.19 présente la réponse temporelle relevée expérimentalement de s(t).

s en V
5

e(t)
4

t en s
0,1 0,2 0,3 0,4

FIGURE 5.19 – Réponse temporelle du système

Q2. À partir de cette représentation temporelle, justifier que la fonction de F(p) peut être modélisée par
une fonction de transfert du second ordre. Rappeler la forme canonique. Déterminer graphiquement
les paramètres caractéristiques sur la figure 5.19.
Q3. En déduire O(p).
Le cahier des charges du système précise que :
— l’erreur indicielle ϵi doit être inférieure à 8% ;
— le premier dépassement D1% doit être inférieur à 50%.
Q4. Déterminer K qui permette de respecter ces deux critères.
Q5. Déterminer la fonction de transfert F(p) avec pour K, la valeur entière qui permet de minimiser
l’erreur indicielle.
Q6. Tracer l’allure de la réponse temporelle sur le graphe précédent (figure 5.19). Préciser le temps
de réponse et l’erreur indicielle du système corrigé.

Exercice 3 - temps de réponse d’un second ordre

Corrigé page 23

Soit le système décrit par les équations :

d s(t)
+ 5 · s(t) = 2 · v(t)
dt
d v(t)
+ 2 · v(t) = u(t)
dt
u(t) = K · (e(t) − s(t))

Q1. Tracer le schéma-bloc correspondant à ces équations.

Q2. Déterminer la fonction de transfert en boucle fermée puis la mettre sous sa forme canonique.
Q3. Déterminer K pour que le temps de réponse soit minimal

19
 5 Étude temporelle des systèmes linéaires

Q3a. sans dépassement,

Q3b. avec dépassement.
Q4. Déterminer la valeur finale pour une entrée en échelon unitaire.
Q5. Tracer l’allure de la réponse temporelle à partir des abaques.
Exercice 4 - Identification temporelle
Corrigé page 23

On a sollicité différents systèmes avec un échelon unitaire, les différentes réponses temporelles sont
présentées dans la figure 5.20.
Q1. Déterminer pour chacun des graphes la fonction de transfert du système.

s
s

1 1

t t
0 0
0 1 0 1

(a) (b)

s s

1 1

t t
0 0
0 1 0 1

(c) (d)

FIGURE 5.20 – Réponses temporelles

Exercice 5 - 1er ordre et intégration

Corrigé page 23

Soit un système décrit par le schéma-bloc

E(p) ε(p) U(p) 1 S(p)
+ K
− p · (1 + 10 · p)

Q1. Déterminer la fonction de transfert en boucle fermée. Mettre sous forme canonique.
Q2. Déterminer l’erreur indicielle. Conclure sur l’effet d’une intégration sur la précision.
Q3. Déterminer K pour que le temps de réponse soit minimal.

20