MONTIGNY Eric

Aide mmoire sur le diagramme de Bode

Bode Diagram
40

35

30
Magnitude (dB)

25

20

15

10

0
90
Phase (deg)

45

0
-2 -1 0 1 2
10 10 10 10 10
Frequency (rad/sec)

MONTIGNY Eric
Fvrier 2006
 MONTIGNY Eric

A. Rappels et gnralits sur le diagramme de Bode

Le diagramme de Bode est constitu de deux diagrammes :
- le module
- largument

a) Le module
Module
Frquences Echelle
logarithmique

Echelle
linaire

b) La phase
Phase

Frquences Echelle
logarithmique

Echelle
linaire
 MONTIGNY Eric

B. Lchelle logarithmique
Lorsque lon trace lvolution dun signal sur un grand domaine de frquence, il est prfrable que lvolution de
la frquence soit reprsente sur une chelle logarithmique.
Dans une chelle logarithmique une octave est lcart entre X et 2X, et une dcade est lcart entre X et 10X.
Pour ce trac on utilise soit un papier millimtr chelle logarithmique, soit on construit lchelle
10
logarithmique partir de lchelle linaire. On sait quune dcade reprsente doctave. Donc si on divise une
3
dcade en dix units, alors trois units vont reprsentes une octave :

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14

Linaire

Logarithmique

X 2X 4X 5X 8X 10X

Octave Octave Octave

Octave

Dcade

C. Rappel sur les systmes du premier ordre

N ( jw)
Un systme de transmittance H ( jw) = est du premier ordre, si le dnominateur D(jw) est un polynme
D( jw)
en jw de degrs 1, et le numrateur N(jw) est un polynme en jw de degrs 1 ou 0.

D. Diagramme de Bode
Le diagramme de Bode dune fonction complexe H(jw) est constitu de deux courbes en reprsentation semi-
logarithmique : H ( jw) et ARG(H(jw)).

Selon que la fonction se prsente sous la forme dun produit, ou dun quotient, on retiendra que :
Si H = A.B alors Bode(H) = Bode(A) + Bode(B)
A
Si H = alors Bode(H) = Bode(A) Bode(B)
B

E. Trac du module
Le trac de la courbe se fait point par point, notons quil convient de prendre des points caractristique
permettant de procder un trac des asymptotes, et ensuite on peut approximer la courbe relle.
1
Etape 1 : Considrons la fonction complexe tracer : H ( jw) =
1 + jw
1
Etape 2 : Intressons-nous au module : H ( jw) =
1 + w

Etape 3 : Gain : G DB = 20. log H ( jw) = 20. log
1
( ) (
= 20. log(1) 20. log 1 + w = 20. log 1 + w )
1 + w
F 0
Etape 4 : Basse frquence : G DB F 0
20. log 1 + 0 0dB( )
Etape 5 : Haute frquence : F +
G DB F +
( )
20. log 1 + w 20. log(w) (Pente de -20dB/dcade).
 MONTIGNY Eric

Il est ensuite possible de tracer le module, en remarquant que le changement de pente va se faire pour une
frquence donne (que lon pourra appeler la pulsation caractristique w0).
H MAX
Par dfinition, on sait que pour w = w0 on a H = , soit :
2
1 1
H ( jw) = = , ce qui implique que w0 = 1.
1 + w 2

GDB w0 10.w0

0dB f

-20dB/dcade

-20dB

Notons que le trac ci-dessus est uniquement asymptotique, mais rien ne vous empche de tracer la courbe relle approche
en vous aidant des asymptotes !

F. Trac de la phase (ou de largument)

Commenons par des rappels sur les nombres complexes :
z = a + jb
IMAGINAIRE
M
b

REEL
0
a
Comme lillustre cette figure, on note quil y a un triangle rectangle, ce qui permet dutiliser la trigonomtrie
0a est le cot adjacent langle
0M est lhypotnuse
0b (soit aM dans le triangle) est le cot oppos langle )
Etant donn que lon connat 0a et 0b et que nous recherchons , nous allons utiliser la tangente, car :
0b
tan =
0a

Si on formalise cette formule, on aura :

b
= tan 1
a

Nous avons donc tous les outils ncessaires pour le trac de largument, dans le diagramme de Bode.
 MONTIGNY Eric

Il est difficile de tracer la courbe point par point, nous allons donc user des asymptotes, en suivant les tapes
suivantes :
1
Etape 1 : Considrons la fonction complexe tracer : H ( jw) =
1 + jw
Etape 2 : Intressons-nous largument : Arg {H ( jw)}
0
Etape 4 : Basse frquence : Arg {H ( jw)}F 0
F 0
tan 1 0
1
Etape 5 : Haute frquence : Arg {H ( jw)}F +
F +

Traons cela dans un repre :

GDB 0,1.w0 w0 10.w0

0 f

- 45

- 90

On notera que le trac en trait fort reprsente les asymptotes brutes. Le trac gris illustre une asymptote plus prcise

G. Allure des fonctions importantes

Etudions certaines fonctions particulires, qui nous permettrons de voir lapplication de la thorie que nous
venons de voir :

a) La fonction H 1 ( jw) = jw
Etape 1 : Considrons la fonction complexe tracer : H 1 ( jw) = jw
Etape 2 : Intressons-nous au module : H 1 ( jw) = w = w
Etape 3 : Gain : G DB = 20. log H 1 ( jw) = 20. log{w}
F +
Etape 5 : Haute frquence : G DB 20. log(w) 20. log(w) (Pente de +20dB/dcade).
F +

Module fonction de la frequence

100

90

80

70

60

50

40

30

20

10
0 1 2 3 4
10 10 10 10 10
 MONTIGNY Eric

Etape 1 : Considrons la fonction complexe tracer : H 1 ( jw) = jw

Etape 2 : Intressons-nous largument : Arg {H 1 ( jw)}
IMAGINAIRE

=/2

REEL
0

Quelque soit la frquence, la phase sera gale / 2, soit :

Phase fonction de la frequence
180

160

140

120

100

80

60

40

20

0
0 1 2 3 4
10 10 10 10 10

1
b) La fonction H 2 ( jw) =
jw
1
Etape 1 : Considrons la fonction complexe tracer : H 2 ( jw) =
jw
1 1
Etape 2 : Intressons-nous au module : H 2 ( jw) = =
w w
1
Etape 3 : Gain : G DB = 20. log H 2 ( jw) = 20. log = 20. log(w)
w
F +
Etape 5 : Haute frquence : G DB 20. log(w) 20. log(w) (Pente de -20dB/dcade).
F +

1
Etape 1 : Considrons la fonction complexe tracer : H 2 ( jw) =
jw
Etape 2 : Intressons-nous largument : Arg {H 2 ( jw)}
 MONTIGNY Eric

1 1
Remarquons que H 2 ( jw) = = j. , ce qui sera reprsente comme suit :
jw w
IMAGINAIRE

REEL
0 =-/2

-b

Donc la phase sera gale - / 2 :

Bode Diagram
5

0
Magnitude (dB)

-5

-10

-15

-20
-89

-89.5
Phase (deg)

-90

-90.5

-91
0 1
10 10
Frequency (rad/sec)
 MONTIGNY Eric

c) La fonction H 3 ( jw) = 1 + jw
Pour le module :
Etape 1 : Considrons la fonction complexe tracer : H 3 ( jw) = 1 + jw
Etape 2 : Intressons-nous au module : H 3 ( jw) = 1 + w
Etape 3 : Gain : G DB = 20. log H ( jw) = 20. log 1 + w{ }
Etape 4 : Basse frquence : F 0
G DB 20. log(1) 0dB
F 0

F +
Etape 5 : Haute frquence : G DB ( )
20. log w 20. log(w) (Pente de +20dB/dcade).
F +

Bode Diagram
40

35

30
Magnitude (dB)

25

20

15

10

0
90
Phase (deg)

45

0
-2 -1 0 1 2
10 10 10 10 10
Frequency (rad/sec)

Pour la phase :
Etape 1 : Considrons la fonction complexe tracer : H 3 ( jw) = 1 + jw
Etape 2 : Intressons-nous largument : Arg {H 3 ( jw)}
0
Etape 4 : Basse frquence : Arg {H ( jw)}
F 0
tan 1 0
1
w +
Etape 5 : Haute frquence : Arg {H ( jw)}
F +
tan 1
1 2
 MONTIGNY Eric

1
d) La fonction H 4 ( jw) =
1 + jw
Pour le module :
1
Etape 1 : Considrons la fonction complexe tracer : H 4 ( jw) =
1 + jw
1
Etape 2 : Intressons-nous au module : H 4 ( jw) =
1 + w
{
Etape 3 : Gain : G DB = 20. log H 4 ( jw) = 20. log 1 + w }
F 0
Etape 4 : Basse frquence : G DB 20. log(1) 0dB
F 0

F +
Etape 5 : Haute frquence : G DB ( )
20. log w 20. log(w) (Pente de -20dB/dcade).
F +

Pour la phase :
1
Etape 1 : Considrons la fonction complexe tracer : H 4 ( jw) =
1 + jw
Etape 2 : Intressons-nous largument : Arg {H 4 ( jw)}
0
Etape 4 : Basse frquence : Arg {H ( jw)}
F 0
tan 1 0
1
w
Etape 5 : Haute frquence : Arg {H ( jw)}
F +
tan 1
1 2

Bode Diagram
0

-5

-10
Magnitude (dB)

-15

-20

-25

-30

-35

-40
0
Phase (deg)

-45

-90
-2 -1 0 1 2
10 10 10 10 10
Frequency (rad/sec)
 MONTIGNY Eric

H. Trac de diagramme de Bode avec MATLAB

Il convient de se rappeler que MATLAB utilise les vecteurs et les matrices Il nest donc pas possible de
donner directement la fonction et dobtenir le trac du diagramme de Bode (chose que ferait MAPLE, mais il ne
dispose pas dune fonction spcialement ddie au diagramme de Bode). Il faut donc entrer les coefficients du
polynme. Voyons cela en tapes :

1
Prenons par exemple la fonction H 4 ( jw) = : on constate quau numrateur il y a un polynme (pas trs
1 + jw
toff certes !), et au numrateur il y a aussi un polynme :
1 a x 4 + a3 .x 3 + ... + a 0
H 4 ( jw) = 4 4
1 + jw b4 x + b3 .x 3 + ... + b0
En procdant par identification, on a :
a 0 = 1

b1 = 1
b = 1
0
Il faut saisir ces coefficients dans lordre dcroissant (du plus grand au plus petit donc !) :
> numerateur=[1];
> denominateur=[1 1];

Ensuite il faut dfinir que la fonction de transfert est le rapport dun polynme se trouvant au numrateur, sur un
polynme prsent au dnominateur (dont nous venons pralablement de dfinir les coefficients) :
> h = tf(numerateur,denominateur);

La dernire tape va consister utiliser linstruction bode() qui se charge delle-mme de tracer le module et
largument :
> bode(h);

La magie va ensuite oprer :

Je vous recommande de saisir les quatre instructions dans un M-FILE, ainsi il est trs facile de procder des
modifications.