1
6. Stabilit des systmes discrets et chantillonns
6.1.
Dfinition
6.2.
Critres algbriques de stabilit
6.2.1. Critre simplifi de Jury (1961)
6.2.2. Critre de Routh modifi
6.2.3. Cas des systmes boucls
6.2.4. Application
6.3.
Mthode du lieu des racines (lieu de EVANS)
6.3.1. Principe de la mthode
6.3.2. Rgles de construction
6.3.3. Application
6.1.
DEFINITION
Comme pour les systmes continus, on dfinit la stabilit dun systme chantillonn par sa
capacit revenir de lui-mme sa position dquilibre aprs avoir t excit par une
perturbation. Il convient donc que sa rponse impulsionnelle vrifie la relation1 :
lim h(kT ) = 0
Or H ( z ) =
{hk } =
hk .z k = i =d0
k =0
i =0
ai z i
=
bi z i
N ( z)
D( z )
En supposant que les ples z Pi de la transmittance, soient simples, rels ou complexes 2 2
conjugus, la rponse impulsionnelle est de la forme :
Cd
C1
C2
H ( z)
=
+
+ ... +
z
z z P1 z z P 2
z z Pd
soit hk =
d
i =0
k
Ci z Pi
Comme hk tend vers zro si i z Pi < 1 on noncera le thorme suivant :
La transmittance H(z) dun systme numrique est stable si et seulement si tous ses ples
sont situs l'intrieur du cercle de rayon unit. Il est dautant plus stable que ses ples
sont prs de lorigine.
En d'
autres termes la stabilit du systme numrique est acquise ds lors que le module de
chacun de ses ples est infrieur 1.
Le systme est juste oscillant si ses ples sont de module 1 (i.e. situs sur le cercle unit)2.
On parle aussi de BIBO stabilit (Bounded Input Bounded Output = entre borne sortie borne) si C tel que
la rponse impulsionnelle soit C.
2 Rappel : lintgrateur admet le ple z = 1.
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
m (z)
Plan z
1
intgration
(z)
1
juste oscillant
Figure 1 : Zone de stabilit des fonctions de transfert numriques
Sachant que z = eTp posons p = + j . La condition de stabilit scrit z = eT < 1 ce qui
quivaut < 0 . Ainsi on retrouve la condition de stabilit nonce pour les systmes
continus.
6.2.
CRITERES ALGEBRIQUES DE STABILITE
En gnral le calcul des racines de l'
quation caractristique du systme (racines de D(z))
n'
est pas facile en l'
absence de moyen de calcul appropri (le recours au logiciel MATLAB
est recommand). Aussi utilise-t-on des critres algbriques qui, sans avoir rsoudre
l'
quation caractristique, permettent de savoir si ses racines sont contenues dans le cercle
unit du plan complexe z.
6.2.1. CRITERE SIMPLIFIE DE JURY (1961)
Soit la transmittance physiquement ralisable H(z) :
n
H ( z) =
n 1
N ( z ) an z + an 1z + ... + a1z + a0
=
=
D( z ) bd z d + bd 1z d 1 + ... + b1z + b0
i =0
d
i =0
ai z i
bi z i
Considrons le dnominateur D(z) de degr d.
D( z ) = bd z d + bd 1 z d 1 + ... + bi z i + ... + b1 z + b0 avec bi
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
Dressons le tableau suivant :
b0
b1
b2
...
bd k
...
bd 2
bd 1
bd
2
3
bd
c0
bd 1
c1
bd 2
c2
...
...
bk
cd k
...
...
b2
b0
cd 2
b1
cd 1
cd 2
d1
cd 3
d2
...
ck 1
...
c1
c0
cd 1
d0
...
...
...
dd 2
d d 3
dd 4
...
...
...
dd 2
d0
...
...
p0
p3
q0
...
...
p1
p2
q1
...
...
p2
p1
q2
...
...
...
...
...
...
...
...
2d-5
2d-4
2d-3
Avec :
ck =
b0
bd
bd k
bk
dk =
q0 =
p0
p3
p3
p0
q1 =
c0
cd 1
p0
p3
cd 1 k
ck
p2
p1
ek =
q2 =
d0
dd 2
p0
p3
dd 2 k
dk
p1
p2
Enonc du critre de JURY :
Pour que toutes les racines de D(z) = 0 soient situes l'
intrieur du cercle unit il faut et il
suffit que les (d + 1) conditions suivantes soient satisfaites :
b0 < bd
c0 > cd 1
D(1) > 0
D(1)
> 0 pour d pair
< 0 pour d impair
et (d 1) contraintes
d0 > dd 2
e0 > ed 3
................
q0 > q 2
Cas particuliers :
Systme de second ordre
Systme du troisime ordre
D( z ) = z 2 + b1z + b0 (d = 2)
D(1) = 1 + b1 + b0 > 0
D(1) = 1 b1 + b0 > 0
D( z ) = z 3 + b2 z 2 + b1z + b0
D(1) = 1 + b2 + b1 + b0 > 0
D(1) = 1 + b2 b1 + b0 < 0
b0 < 1
b0 < 1
(d = 3)
b02 1 > b0b2 b1
6.2.2. CRITERE DE ROUTH MODIFIE
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
Comme nous l'
avons voqu loccasion de ltude des systmes de commande continus, le
critre de ROUTH indique si un polynme contient des racines partie relle positive. Ce
critre ne s'
applique donc pas aux transmittances en z puisque dans ce cas il sagit de vrifier
que le module des racines de D(z) est <1. Cependant il s'
applique aux polynmes D(w)
obtenus par substitution de la variable w (transformation de MBIUS) la variable z avec :
z 1
1+ w
z=
z +1
1 w
En effet cette transformation fait correspondre au cercle de rayon unit du plan complexe z
(zone de stabilit dans ce plan) le demi-plan gauche de la variable complexe w.
w=
z = e Tp
Plan p
Plan z
z 1
z + 1
w =
Plan w
m
1
+2/
2/
Figure 2 : Plan complexe et zone de stabilit
Le critre s'
nonce comme suit :
D( w) = bd w d + bd 1w d 1 + ... + bi wi + ... + b1w + b0 a toutes ses racines partie relle ngative
si : i tous les bi 0 et de mme signe
tous les termes de la premire colonne du tableau de ROUTH sont de mme signe.
Formons le tableau de ROUTH :
wd
w d 1
w d 2
wd 3
bd
bd 1
w d 4
..........
e1
w
w0
Comnum_6_txt.doc
c1
d1
..........
..........
..........
bd 2
bd 3
c2
d2
bd 4
bd 5
c3
..........
bd 6
bd 7
c4
..........
..........
..........
..........
..........
..........
..........
..........
..........
..........
..........
..........
..........
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
Avec :
c1 =
bd 1bd 2 bd bd 3
bd 1
d1 =
c1bd 3 c2bd 1
c1
c2 =
bd 1bd 4 bd bd 5
bd 1
d2 =
c1bd 5 c3bd 1
c1
6.2.3. CAS DES SYSTEMES BOUCLES
On calcule le dnominateur de la FTBF(z) gal :
Dbf ( z ) = Dbo( z ) + Nbo( z )
et lon applique les critres algbriques.
6.2.4. APPLICATION
On considre l'
asservissement de position reprsent ci aprs. La priode d'
chantillonnage T
est gale une seconde. Le correcteur numrique est un simple gain K.
e(t)
+
T=1s
C(z) = K
s(t)
B0 (p)
p (1+
1+ p)
p
Figure 3 : Asservissement lmentaire de position commande numrique
La fonction de transfert du processus associe au BOZ est donne par lexpression suivante.
G (z) = (1 z 1 )
1
2
p (1 + p)
0,37 z + 0,26
( z 1)( z 0,37)
et
FTBO( z ) =
Nbo( z )
= K.G( z)
Dbo( z )
Calculons la FTBF(z).
FTBF ( z ) =
K (0,37 z + 0,26)
z 2 + (0,37 K 1,37) z + 0,37 + 0,26 K
Nbf ( z )
Nbo( z )
=
Dbf ( z ) Dbo( z ) + Nbo( z )
Appliquons le critre de JURY D ( z ) = z 2 + (0,37 K 1,37) z + 0,37 + 0,26 K
D (1) = 1 + 0,37 K 1,37 + 0,37 + 0,26 K = 0,63K > 0
D (1) = 1 0,37 K + 1,37 + 0,37 + 0,26 K = 2,74 0,11K > 0
0 < K < 2,42
0,37 + 0,26 K < 1
Pour K = 2,42 il vient Dbf ( z ) = z 2 0,4 z + 1 = ( z e1,37 j )( z e 1,37 j )
z P1 = e1,37 j = eTp avec p = j os = 1,37 j car T = 1 seconde. Ainsi : os = 1,37 rd/s
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
Appliquons le critre de ROUTH
G ( w) =
2w
.
1+ w
1
2
p (1 + p)
= [G(z)]z = 1+ w =
1 w
1 w 1 + 0,164w
2 w 1 + 2,164 w
D ( w) = (4,328 0,164 K ) w 2 + (2 0,836 K ) w + K
Formons le tableau de ROUTH.
w2
4,328 0,164 K
K < 26,39
w1
2 0,836 K
K < 2,4
w0
0 < K < 2,4
K>0
A partir du tableau de ROUTH on peut encore calculer la pulsation des auto-oscillations
obtenues pour K = 2,392. La dmarche est la suivante :
(4,328 0,164 K ) w 2 + K = 0 avec K = 2,4. Cette quation donne une valeur de w.
w = tanh
Tp
avec p = j os .
2
Dans le cas de notre problme wos = j. tan
6.3.
T os
2
soit
os = 1,33 rad/s
METHODE DU LIEU DES RACINES (LIEU DE EVANS)
6.3.1. PRINCIPE DE LA METHODE
Considrons un systme boucl retour unitaire :
(t)
e(t)
C(z)
mn
T
B0(p)
m(t)
s(t)
P(p)
_
G(z) =
(z)
E(z)
+
B0(p). P(p)
S(z)
C(z)
G(z)
Figure 4 : Transposition de la boucle de commande dans le domaine z
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
Nbo( z )
du systme de commande, il s'
agit d'
apprcier sa stabilit
Dbo( z )
en boucle ferme. Elle dpend de la distribution des ples de la FTBF (z ) dans le plan
complexe z . Lorsquon parle de lieu de racines ou des ples il faut entendre lieu des ples
de la FTBF (z ) lorsquun paramtre (en gnral le facteur de gain) varie de 0 linfini.
Connaissant la FTBO ( z ) =
Ces ples sont les racines du dnominateur de la FTBF (z ) . Calculons cette FTBF(z).
FTBF ( z ) =
C ( z )G ( z )
FTBO ( z )
Nbo( z )
Nbf ( z )
=
=
=
1 + C ( z )G ( z ) 1 + FTBO( z ) Dbo( z ) + Nbo( z ) Dbf ( z )
La stabilit de la FTBF(z) dpend des racines de l'
quation :
[1 + FTBO( z )] = Dbo( z ) + Nbo( z ) = 0
Ecrivons la FTBO(z ) sous la forme canonique permettant de mettre en vidence ses ples, ses
zros et son facteur de gain k :
n
FTBO ( z ) = k
(z z Zi )
i =1
d
(z z Pi )
Nbo( z )
; attention le facteur de gain k K = gain statique.
Dbo( z )
i =1
Plaons-nous dans le plan complexe z. Soit un point M d'
affixe z dans ce plan. Ce point
reprsente un ple de la FTBF (z ) la condition que, pour cette valeur de z,
[1 + FTBO( z)] = 0 soit vrifie. A chaque valeur de k on peut ainsi associer les solutions de :
[1 + FTBO( z)] = 0
i =1
i =1
(z z Pi ) + k (z z Zi ) = Dbo( z ) + Nbo( z ) = 0
infini les points reprsentatifs des ples de la FTBF (z )
Si le facteur de gain k varie de zro l'
dcrivent une courbe appele lieu de EVANS ou lieu des racines (roots locus).
Pour tudier la stabilit d'
un systme de commande par la mthode de EVANS il convient :
de dterminer la FTBO (z ) du systme concern ;
dcrire la FTBO (z ) selon la forme canonique indique ci dessus ;
de construire le lieu de EVANS la main, selon les rgles qui seront explicites par
la suite, ou en utilisant un logiciel adapt (e.g. MATLAB) ;
de paramtrer le lieu
en fonction du facteur de gain k ;
d'
analyser le trac de
en fonction de k.
ainsi ralis afin de prciser les zones de stabilit de FTBF(z)
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
6.3.2. REGLES DE CONSTRUCTION
Rappelons quelques lments concernant la reprsentation dune fonction de transfert dans le
plan complexe.
Soit titre d'
exemple FTBO( z ) =
k ( z z Z1 )
( z z P1 )( z z P 2 )( z z P 3 )
Traons les ples et zros de cette transmittance dans le plan complexe z. Pour tout point M
d'
affixe z M il est possible de mesurer le module et l'
argument de chaque ple et zro et ainsi
de dterminer la transmittance FTBO ( z M ) au facteur de gain k prs.
m
Z P2
X
X
ZP3
X ZM
d2
0 X
d1
d3
[FTBO( z )]z = z M
=k
1
d1d 2 d 3
.e j (1 1 2 3 )
Zz1
ZP1
a. Conditions des angles et des modules
n
(z z Zi )
i =1
d
1 + FTBO ( z ) = 0
(z z Pi )
1
k
k >0
avec
1 1 j( 2 +1)
= .e
k
k
i =1
(z M
Gnralisons :
i =1
d
i
i =1
d
(z M z Pi ) d i
i =1
Si z M
z Zi )
i =1
i =1
1 j ( 2 +1)
e
k
i =1
(lieu de EVANS) alors les relations suivantes sont vrifies :
condition des angles :
i i = (2 + 1)
i =1
i =1
condition des modules :
i
i =1
d
di
1
k
i =1
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
b. Axe de symtrie
Les racines de [1 + FTBO( z )] = 0 sont relles ou 2 2 conjugues. L'axe rel du plan
complexe z est donc axe de symtrie.
c. Points de dpart et d'arrive
Le facteur de gain k est le paramtre de rglage. Il varie de zro l'
infini.
k =0
z = z Pi
les d ples de la FTBO(z) sont les points de dpart de
k=
z = z Zi
les n zros de la FTBO(z) sont les points d'arrive de
d. Directions asymptotiques
Il existe d points de dpart et n points d'
arrive. Comme d n il existe ( d n) directions
asymptotiques.
Avec :
d
Direction asymptotique
a =
i =1
a =
z Pi
n
i =1
z Zi
d n
(2 + 1)
d n
e. Parties de l'axe rel au lieu
Tout point M de laxe rel
appartient
si on compte un
nombre impair de ples et de
zros droite de M.
Cette appartenance se dmontre
facilement par la condition des
angles.
f. Points de branchement
Ces points correspondent des racines doubles du lieu. L'
exemple ci contre matrialise le
point de branchement dont on dsire calculer l'
abscisse. On appliquera la mthode de calcul
suivante :
n
d
1
1
=
i =1 b z Zi
i =1 b z Pi
d
ou encore :
dk
= 0 avec k =
dz
( z z Pj )
j =1
n
( z z Zi )
i =1
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
10
6. Stabilit des systmes discrets et chantillonns
m
X
zZ1
zP1
zP2
X
g. Trac du lieu d'EVANS de quelques transmittances
On consultera l'
annexe A.6.1.
6.3.3. APPLICATION
Considrons l'
asservissement de position tudi en 6.2.4.
e(t)
+
G (z) = (1 z 1 )
T=1s
1
2
p (1 + p )
C(z) = K
Bo(p)
s(t)
p (1+ p)
0,37 z + 0,26
( z 1)( z 0,37)
a. Mthode graphique et analytique
Identifions la FTBO(z)
FTBO( z ) =
0,37 K ( z + 0,7)
k ( z + 0,7)
=
( z 1)( z 0,37) ( z 1)( z 0,37)
avec le facteur de gain k = 0,37. K
Elments pour tracer le lieu de
n=1
1 point d'
arrive ; d = 2
2 points de dpart ;
( d n) = 1
1 direction asymptotique
Equation du lieu (hors axe rel)
Dbf ( z ) = dnominateur de la FTBF(z)
Dbf ( z ) = z 2 + ( k 1,37) z + (0,37 + 0,7 k )
( x + 0,7) 2 + y 2 = (1,35) 2
k = 1,37 2 x
La premire quation indique que le lieu est un cercle.
La seconde quation permet de paramtrer le lieu.
dk
Points de branchement :
=0
b1 = 2,05 et b2 = 0,65
dz
Soit z = x + jy
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
11
6. Stabilit des systmes discrets et chantillonns
Intersection avec le cercle unit
( x + 0,7) 2 + y 2 = (1,35) 2
x2 + y2 = 1
Trac du lieu de EVANS
x = 0,2375
m
kos= 0,898
Sachant que k = 1,37 2 x
Soit Kos = 2,42
zZ1= 0,7
kos = 0,895
Cercle unit
Analyse succincte des rsultats
Le systme est stable pour :
0 < k < 0,898 soit 0 < K < 2,42
La pulsation des auto-oscillations est :
zos = 0,2375 + j 0,9714 = e j os
os = 1,33 rd / s
X
zP2= 1
b1= 2,05
b1= 0,65
zP1= 0,37
Figure 5 : Trac du lieu de Evans
b. Utilisation du logiciel MATLAB:
Soit lasservissement de position tudi en 6.2.4. La priode dchantillonnage est de 1s.
>> Gdp=tf([1],[1 1 0]);
>> Gdz=c2d(Gdp,1)
Transfer function:
0.3679 z + 0.2642
--------------------------z^2 - 1.368 z + 0.3679
Sampling time: 1
>> Ndz=[0.3679 0.2642];
>> Ddz=[1 1.368 0.3679];
>> Gdz=tf(Ndz,Ddz,-1)
On dfinit G ( p) = 1 / p(1 + p)
On convertit le systme continu associ un
bloqueur dordre 0 en un systme discret par la
commande (c2d) avec T = 1 s.
0n obtient directement Gdz cest dire G (z ) .
Dans ce cas MATLAB sait que Gdz est la fonction
de transfert dun systme discrtis au pas de
une seconde.
On aurait pu ne connatre que les polynmes Ndz
et Ddz. Dans ce cas on obtient la fonction de
transfert Gdz par linstruction ci-contre. Le
paramtre 1 indique que le systme est numrique de priode non prcise.
On obtient la forme de EVANS :
>> zpkGdz=zpk(Gdz) % on obtient la forme zpk de Gdz
Zero/pole/gain:
0.3679 (z+0.7181)
----------------------(z-1) (z-0.3678)
Sampling time: unspecified
Traons le lieu de EVANS :
>> rlocus(Gdz),zgrid % on obtient la figure la figure ci-dessous
Attention le gain affich par le LTI viewer est gal au gain statique.
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
12
6. Stabilit des systmes discrets et chantillonns
On explicitera au chapitre 8 lintrt de la grille obtenue avec la commande zgrid .
Avec la souris on peut dplacer le curseur de mesure sur le lieu et lon apprcie les
performances de la boucle ferme. Ainsi pour K = 2,42 le systme nest plus amorti et il
oscille la pulsation de 1,33 rad/s. Avec le zoom on peut travailler finement sur cette figure.
On apprcie finement les performances de lasservissement en boucle ferme. Ainsi pour K =
0,34 on obtient un amortissement de 0,7 correspondant un dpassement est de 4,6 %. De
mme pour K = 1 on vrifie que le dpassement est suprieur 44 %.
On peut encore slectionner un point sur la figure et connatre, dans la fentre MATLAB, les
donnes attaches ce point (ples de la FTBF et facteur de gain k).
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
6. Stabilit des systmes discrets et chantillonns
13
>> [k,poles]=rlocfind(Gdz)
Select a point in the graphics window
selected_point =
0.6164 + 0.2633i
k = 0.3368
poles = 0.6221 + 0.2644i
0.6221 - 0.2644i
A titre dexercice calculons la FTBF du systme continu (Ftbfdp) et celle du systme
chantillonn (Ftbfdz) et comparons leur rponse indicielle grce au LTI viewer .
>> Ftbfdz=feedback(Gdz,1)
Transfer function:
0.3679 z + 0.2642
---------------------z^2 - z + 0.6321
Sampling time: 1
>> Ftbfdp=feedback(Gdp,1)
Transfer function:
1
-------------s^2 + s + 1
Attention la rponse indicielle du systme chantillonn en boucle ferme est donne aux
seuls instants dchantillonnage. Ce qui se passe entre ces instants nous l'
ignorons.
Pour le connatre on simulera le processus avec SIMULINK.
>> step(Ftbfdz,Ftbfdp)
On remarquera leffet dstabilisant de lchantillonnage.
Comnum_6_txt.doc
Mis jour le 11/02/2006
Cours de M. Cougnon JL
14
6. Stabilit des systmes discrets et chantillonns
ANNEXE A.6.1.
LIEU D'EVANS
Transmittance
Lieu d'EVANS
Im
1
z 1
k=
k=
z
z 1
Im
k =0
k=
z
(z 1) 2
Comnum_6_txt.doc
k =0
Im
Mis jour le 11/02/2006
k =0
X
X
Cours de M. Cougnon JL
15
6. Stabilit des systmes discrets et chantillonns
k=
z
(z a) 2
Im
Im
a,b
XXa
Im
k=
et
Mis jour le 11/02/2006
k =0
X
1
k =0
Xa
o
o
X
Comnum_6_txt.doc
z( z b)
( z a)( z a)
Xa
Im
k =0
o Xb
k=
z
( z a)( z a)
k =0
k=
z
( z a)( z b)
Cours de M. Cougnon JL
