Aide mémoire 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)
A. Rappels et généralités sur le diagramme de Bode
Le diagramme de Bode est constitué de deux diagrammes :
- le module
- l’argument
a) Le module
Module
Fréquences Echelle
logarithmique
Echelle
linéaire
b) La phase
Phase
Fréquences Echelle
logarithmique
Echelle
linéaire
B. L’échelle logarithmique
Lorsque l’on trace l’évolution d’un signal sur un grand domaine de fréquence, il est préférable que l’évolution de
la fréquence soit représentée sur une échelle logarithmique.
Dans une échelle logarithmique une octave est l’écart entre X et 2X, et une décade est l’écart entre X et 10X.
Pour ce tracé on utilise soit un papier millimétré à échelle logarithmique, soit on construit l’échelle
10
logarithmique à partir de l’échelle linéaire. On sait qu’une décade représente d’octave. Donc si on divise une
3
décade en dix unités, alors trois unités vont représentées une octave :
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
Linéaire
Logarithmique
X 2X 4X 5X 8X 10X
Octave Octave Octave
Octave
Décade
C. Rappel sur les systèmes du premier ordre
N ( jw)
Un système de transmittance H ( jw) = est du premier ordre, si le dénominateur D(jw) est un polynôme
D( jw)
en jw de degrés 1, et le numérateur N(jw) est un polynôme en jw de degrés 1 ou 0.
D. Diagramme de Bode
Le diagramme de Bode d’une fonction complexe H(jw) est constitué de deux courbes en représentation semi-
logarithmique : H ( jw) et ARG(H(jw)).
Selon que la fonction se présente sous la forme d’un produit, ou d’un 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 qu’il convient de prendre des points caractéristique
permettant de procéder à un tracé des asymptotes, et ensuite on peut approximer la courbe réelle.
1
Etape 1 : Considérons la fonction complexe à tracer : H ( jw) =
1 + jw
1
Etape 2 : Intéressons-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 fréquence : G DB F ≈0
→ −20. log 1 + 0 ≈ 0dB( )
Etape 5 : Haute fréquence : F → +∞
G DB F ≈ +∞
( )
→ −20. log 1 + w² ≈ −20. log(w) (Pente de -20dB/décade).
Il est ensuite possible de tracer le module, en remarquant que le changement de pente va se faire pour une
fréquence donnée (que l’on pourra appeler la pulsation caractéristique w0).
H MAX
Par définition, 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/décade
-20dB
Notons que le tracé ci-dessus est uniquement asymptotique, mais rien ne vous empêche de tracer la courbée réelle approchée
en vous aidant des asymptotes !
F. Tracé de la phase (ou de l’argument)
Commençons par des rappels sur les nombres complexes :
♦ z = a + jb
IMAGINAIRE
M
b
φ
REEL
0
a
Comme l’illustre cette figure, on note qu’il y a un triangle rectangle, ce qui permet d’utiliser la trigonométrie…
♦ 0a est le coté adjacent à l’angle φ
♦ 0M est l’hypoténuse
♦ 0b (soit aM dans le triangle) est le coté opposé à l’angle φ)
Etant donné que l’on connaît 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 nécessaires pour le tracé de l’argument, dans le diagramme de Bode.
Il est difficile de tracer la courbe point par point, nous allons donc user des asymptotes, en suivant les étapes
suivantes :
1
Etape 1 : Considérons la fonction complexe à tracer : H ( jw) =
1 + jw
Etape 2 : Intéressons-nous à l’argument : Arg {H ( jw)}
0
Etape 4 : Basse fréquence : Arg {H ( jw)}F →0
F ≈0
→ tan −1 ≈ 0
1
Etape 5 : Haute fréquence : Arg {H ( jw)}F →+∞
F ≈ +∞
→
Traçons cela dans un repère :
GDB 0,1.w0 w0 10.w0
0° f
- 45°
- 90°
On notera que le tracé en trait fort représente les asymptotes ‘brutes’. Le tracé gris illustre une asymptote plus précise…
G. Allure des fonctions importantes
Etudions certaines fonctions particulières, qui nous permettrons de voir l’application de la théorie que nous
venons de voir :
a) La fonction H 1 ( jw) = jw
Etape 1 : Considérons la fonction complexe à tracer : H 1 ( jw) = jw
Etape 2 : Intéressons-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 fréquence : G DB → 20. log(w) ≈ 20. log(w) (Pente de +20dB/décade).
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
Etape 1 : Considérons la fonction complexe à tracer : H 1 ( jw) = jw
Etape 2 : Intéressons-nous à l’argument : Arg {H 1 ( jw)}
IMAGINAIRE
φ=π/2
REEL
0
Quelque soit la fréquence, 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 : Considérons la fonction complexe à tracer : H 2 ( jw) =
jw
1 1
Etape 2 : Intéressons-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 fréquence : G DB → −20. log(w) ≈ −20. log(w) (Pente de -20dB/décade).
F ≈ +∞
1
Etape 1 : Considérons la fonction complexe à tracer : H 2 ( jw) =
jw
Etape 2 : Intéressons-nous à l’argument : Arg {H 2 ( jw)}
1 1
Remarquons que H 2 ( jw) = = − j. , ce qui sera représentée 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)
c) La fonction H 3 ( jw) = 1 + jw
Pour le module :
Etape 1 : Considérons la fonction complexe à tracer : H 3 ( jw) = 1 + jw
Etape 2 : Intéressons-nous au module : H 3 ( jw) = 1 + w²
Etape 3 : Gain : G DB = 20. log H ( jw) = 20. log 1 + w²{ }
Etape 4 : Basse fréquence : F →0
G DB → 20. log(1) ≈ 0dB
F ≈0
F → +∞
Etape 5 : Haute fréquence : G DB ( )
→ 20. log w² ≈ 20. log(w) (Pente de +20dB/décade).
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 : Considérons la fonction complexe à tracer : H 3 ( jw) = 1 + jw
Etape 2 : Intéressons-nous à l’argument : Arg {H 3 ( jw)}
0
Etape 4 : Basse fréquence : Arg {H ( jw)}
F ≈0
→ tan −1 ≈ 0
1
w +π
Etape 5 : Haute fréquence : Arg {H ( jw)}
F ≈ +∞
→ tan −1 ≈
1 2
1
d) La fonction H 4 ( jw) =
1 + jw
Pour le module :
1
Etape 1 : Considérons la fonction complexe à tracer : H 4 ( jw) =
1 + jw
1
Etape 2 : Intéressons-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 fréquence : G DB → −20. log(1) ≈ 0dB
F ≈0
F → +∞
Etape 5 : Haute fréquence : G DB ( )
→ −20. log w² ≈ −20. log(w) (Pente de -20dB/décade).
F ≈ +∞
Pour la phase :
1
Etape 1 : Considérons la fonction complexe à tracer : H 4 ( jw) =
1 + jw
Etape 2 : Intéressons-nous à l’argument : Arg {H 4 ( jw)}
0
Etape 4 : Basse fréquence : Arg {H ( jw)}
F ≈0
→ tan −1 ≈ 0
1
− w −π
Etape 5 : Haute fréquence : 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)
H. Tracé de diagramme de Bode avec MATLAB
Il convient de se rappeler que MATLAB utilise les vecteurs et les matrices… Il n’est donc pas possible de
donner directement la fonction et d’obtenir le tracé du diagramme de Bode (chose que ferait MAPLE, mais il ne
dispose pas d’une fonction spécialement dédiée au diagramme de Bode…). Il faut donc entrer les coefficients du
polynôme. Voyons cela en étapes :
1
Prenons par exemple la fonction H 4 ( jw) = : on constate qu’au numérateur il y a un polynôme (pas très
1 + jw
étoffé certes !), et au numérateur il y a aussi un polynôme :
1 a x 4 + a3 .x 3 + ... + a 0
H 4 ( jw) = ⇔ 4 4
1 + jw b4 x + b3 .x 3 + ... + b0
En procédant par identification, on a :
a 0 = 1
b1 = 1
b = 1
0
Il faut saisir ces coefficients dans l’ordre décroissant (du plus grand au plus petit donc !) :
> numerateur=[1];
> denominateur=[1 1];
Ensuite il faut définir que la fonction de transfert est le rapport d’un polynôme se trouvant au numérateur, sur un
polynôme présent au dénominateur (dont nous venons préalablement de définir les coefficients) :
> h = tf(numerateur,denominateur);
La dernière étape va consister à utiliser l’instruction bode() qui se charge d’elle-même de tracer le module et
l’argument :
> bode(h);
La magie va ensuite opérer :
Je vous recommande de saisir les quatre instructions dans un M-FILE, ainsi il est très facile de procéder à des
modifications.
