Tracé de diagrammes de Bode

**syl55** · 22/05/2008, 15h33

Bonjour,

j'ai écrit un programme (.m) pour tracé les diagrammes de Bode d'une fonction de transfert mais rien ne se passe quand je clique sur run.
j'ai le message d'erreur suivant : "Error: Missing variable or function."

Voici le programme :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
% définition des différentes fonctions de transfert de la PLL. (en phase)
 
clear all;
close all;
 
w=logspace(-1,5,500);
 
Kc=1;                          % gain du comparateur de phase
 
Kf=1;                          % fonction de transfert du filtre passe bas de gain Kf, constante de temps Tau                              
Tau=1;
F=Kf/(1+Tau*i*w);        
 
Ko=1;                          % gain du VCO
 
G=1/i*w;                       % intégrateur
 
% tracé des diagrammes de Bode en boucle ouverte et en boucle fermée
 
figure(1);
 
A=Kc*F*Ko*G;                   % fonction de transfert en boucle ouverte
 
subplot(2,1,1);
semilogx(w,20*log10(abs(A)));
grid; xlabel('log(w)'); ylabel('G'); title('Diagramme de Bode de la FTBO en gain');
subplot(2,1,2);
semilogx(w,unwrap(angle(A))*180/pi);
grid; xlabel('log(w)'); ylabel('Phi'); title('Diagramme de Bode de la FTBO en phase');
 
 
figure(2);
 
B=(Kc*F*Ko*G)/(1+Kc*F*Ko*G);    % fonction de transfert en boucle fermée (retour unitaire)
 
subplot(2,1,1);
semilogx(w,20*log10(abs(B)));
grid; xlabel('log(w)'); ylabel('G'); title('Diagramme de Bode de la FTBF en gain');
subplot(2,1,2);
semilogx(w,unwrap(angle(B))*180/pi);
grid; xlabel('log(w)'); ylabel('Phi'); title('Diagramme de Bode de la FTBF en phase');

Si quelqu'un pourrait m'aider se serait vraiment gentil car c'est mon premier programme et je ne gère pas trop !
merci beaucoup.

**Jerome Briot** · 22/05/2008, 17h39

Envoyé par syl55

j'ai le message d'erreur suivant : "Error: Missing variable or function."

Es-tu sûr ? Je n'ai pas le même message sur ma machine

Quel est le message d'erreur complet ?

**syl55** · 22/05/2008, 20h57

oui je suis sur, pour preuve : (stabilité est le nom du fichier .m)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
??? stabilité
            |
Error: Missing variable or function.
 
>>

j'ai la version 7.0.0 (R14) de Matlab. Qu'est ce que vous avez vous comme message ?

merci encore !

**Jerome Briot** · 23/05/2008, 00h30

Je pense que MATLAB n'apprécie pas les caractères accentués ni dans les noms de fonction, ni dans les noms de variables

**syl55** · 23/05/2008, 10h25

j'ai changé le nom de mon fichier ( je l'ai appelé stabilite2.m mais je n'ai rien changé dans le corps du programme) et là je n'ai plus le même message d'erreur !

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
??? Error using ==> mrdivide
Matrix dimensions must agree.
 
Error in ==> stabilite2 at 12
F=Kf/(1+Tau*i*w);

Pourriez vous m'indiquer ce que cela signifie ? quel est le problème avec mon expression de F ?

merci beaucoup de votre aide !

**syl55** · 23/05/2008, 11h06

Je me suis inspiré de ce programme :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
%=======================================================================================
% b) plotting of frequency response as a Bode plot using standard complex number
% and plotting functions
%=======================================================================================
 
% open figure 2 for second frequency response plots
figure(2);
 
% create vector of 500 logarithmically spaced (i.e., same number of points per decade)
% frequencies from 10^-1 = 0.1rad/sec to 10^5 = 100000rad/sec
w = logspace(-1,5,500);
 
% define transfer function
H = 5000*j*w./(500 + j*w);
 
% divide figure window into two rows, one column, and plot magnitude response in top graph,
% phase response in bottom graph
subplot(2,1,1);
semilogx(w,20*log10(abs(H)));
grid; ylabel('|H(j\omega)|'); title('Bode Plot: Magnitude Response');
subplot(2,1,2);
semilogx(w,unwrap(angle(H))*180/pi);
grid; xlabel('\omega (rad/sec)'); ylabel('\angleH(j\omega) (\circ)'); title('Bode Plot: Phase Response');

Il marche très bien sur ma machine et le mien est quasiment la copie conforme sauf que moi il ne marche pas ! j'ai remarqué cependant une différence entre ce programme et le mien c'est la présence d'un "." dans la fonction de transfert H, que signifie ce point ? mon message d'erreur vient t-il de là ?

merci.

**mr_samurai** · 23/05/2008, 11h14

Salut,

Le point signifie que l'opération sera faites terme à terme :

(x.*y) (i) est égale à x(i)*y(i)
Alors que :
x*y est le produit matriciel de x et y.

Ceci explique l'erreur que tu as obtenu.

++

**syl55** · 23/05/2008, 11h16

j'ai essayé de mettre des points un peu partout dans mon programme

et maintenant ça marche mais je ne sais pas si les diagrammes de bode qui s'affichent sont les bons car j'ai peut etre mis les points au mauvaise endroit.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
% definition des differentes fonctions de transfert de la PLL. (en phase)
 
clear all;
close all;
 
w=logspace(-1,5,500);
 
Kc=1;                          % gain du comparateur de phase
 
Kf=1;                          % fonction de transfert du filtre passe bas de gain Kf, constante de temps Tau                              
Tau=1;
F=Kf./(1+Tau*i*w);        
 
Ko=1;                          % gain du VCO
 
G=2*pi/i*w;                    % integrateur
 
% trace des diagrammes de Bode en boucle ouverte et en boucle fermee
 
figure(1);
 
A=Kc.*F*Ko.*G;                   % fonction de transfert en boucle ouverte
 
subplot(2,1,1);
semilogx(w,20*log10(abs(A)));
grid; xlabel('log(w)'); ylabel('G'); title('Diagramme de Bode de la FTBO en gain');
 
subplot(2,1,2);
semilogx(w,unwrap(angle(A))*180/pi);
grid; xlabel('log(w)'); ylabel('Phi'); title('Diagramme de Bode de la FTBO en phase');
 
 
figure(2);
 
B=(Kc.*F*Ko.*G)./(1+Kc.*F*Ko.*G);    % fonction de transfert en boucle fermee (retour unitaire)
 
subplot(2,1,1);
semilogx(w,20*log10(abs(B)));
grid; xlabel('log(w)'); ylabel('G'); title('Diagramme de Bode de la FTBF en gain');
 
subplot(2,1,2);
semilogx(w,unwrap(angle(B))*180/pi);
grid; xlabel('log(w)'); ylabel('Phi'); title('Diagramme de Bode de la FTBF en phase');

est ce que vous pourriez me dire si je n'ai pas fait de catastrophe avec mes points ?

( j'ai mis des points dans les fonctions de transfert F, A et B)
merci

**mr_samurai** · 23/05/2008, 11h29

A priori, ton code devrait bien tourner. On ne fait pas de calcul matriciel pour les fonctions de transferts

.

De plus tu connais les expressions de tes fonctions de transfert. Et tu peux donc vérifier les pentes (-20 dB, -40dB, ...).

++ bon courage

**syl55** · 23/05/2008, 11h31

merci pour votre aide,

après avoir faire quelques tests je me suis rendu compte que les diagrammes de bode ne correspondaient pas aux fonctions de transfert ! cependant j'ai trouvé mon erreur, j'ai oublié de mettre des parenthèses à certains endroit. d'où le programme final pour ceux que ça intéresse (on ne sait jamais !)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
% definition des differentes fonctions de transfert de la PLL. (en phase)
 
clear all;
close all;
 
w=logspace(-3,3,500);
 
Kc=(1);                              % gain du comparateur de phase
 
Kf=(1);                              % fonction de transfert du filtre passe bas de gain Kf, constante de temps Tau                              
Tau=(1);
F=(Kf./(1+Tau.*j.*w));        
 
Ko=(1);                              % gain du VCO
 
G=(2*pi./(j.*w));                    % integrateur
 
% trace des diagrammes de Bode en boucle ouverte et en boucle fermee
 
figure(1);
 
A=Kc.*F.*Ko.*G;                       % fonction de transfert en boucle ouverte
 
subplot(2,1,1);
semilogx(w,20*log10(abs(A)));
grid; ylabel('G en dB'); title('Diagramme de Bode de la FTBO en gain');
 
subplot(2,1,2);
semilogx(w,unwrap(angle(A))*180/pi);
grid; xlabel('w en rad/s'); ylabel('Phi en degre'); title('Diagramme de Bode de la FTBO en phase');
 
 
figure(2);
 
B=(Kc.*F.*Ko.*G)./(1+Kc.*F.*Ko.*G);  % fonction de transfert en boucle fermee (retour unitaire)
 
subplot(2,1,1);
semilogx(w,20*log10(abs(B)));
grid; ylabel('G en dB'); title('Diagramme de Bode de la FTBF en gain');
 
subplot(2,1,2);
semilogx(w,unwrap(angle(B))*180/pi);
grid; xlabel('w en rad/s'); ylabel('Phi en degre'); title('Diagramme de Bode de la FTBF en phase');

Personne ne sait comment afficher les marges de phases et de gains sur les diagrammes ??

**phryte** · 23/05/2008, 16h30

L'instruction "margin" est faite pour ça.

**syl55** · 23/05/2008, 16h34

oui j'ai vu cette instruction dans l'aide de Matlab, mais je n'arrive pas à l'utiliser, pourriez m'indiquer la syntaxe dans le cas de mon programme svp ?

**phryte** · 23/05/2008, 17h04

Tu définis la fonction de transfert en BO avec l'instruction tf : BO = tf(...)
Puis tu exécutes margin(BO)

**syl55** · 23/05/2008, 17h12

voici mon programme :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
% définition des différentes fonctions de transfert de la PLL. (en phase)
 
clear all;
close all;
 
w=logspace(-3,3,500);
 
Kc=(1);                              % gain du comparateur de phase
 
Kf=(1);                              % fonction de transfert du filtre passe bas de gain Kf, constante de temps Tau                              
Tau=(1);
F=(Kf./(1+Tau.*j.*w));        
 
Ko=(1);                              % gain du VCO
 
G=(2*pi./(j.*w));                    % intégrateur
 
% tracé des diagrammes de Bode en boucle ouverte et en boucle fermée
 
figure(1);
 
A=Kc.*F.*Ko.*G;                      % fonction de transfert en boucle ouverte
 
subplot(2,1,1);
semilogx(w,20*log10(abs(A)));
grid; ylabel('G en dB'); title('Diagramme de Bode de la FTBO en gain');
 
subplot(2,1,2);
semilogx(w,unwrap(angle(A))*180/pi);
grid; xlabel('w en rad/s'); ylabel('Phi en degré'); title('Diagramme de Bode de la FTBO en phase');
 
 
figure(2);
 
B=(Kc.*F.*Ko.*G)./(1+Kc.*F.*Ko.*G);  % fonction de transfert en boucle fermée (retour unitaire)
 
subplot(2,1,1);
semilogx(w,20*log10(abs(B)));
grid; ylabel('G en dB'); title('Diagramme de Bode de la FTBF en gain');
 
subplot(2,1,2);
semilogx(w,unwrap(angle(B))*180/pi);
grid; xlabel('w en rad/s'); ylabel('Phi en degré'); title('Diagramme de Bode de la FTBF en phase');

je voudrais avoir les marges de phases et de gains des fonctions de transfert A et B, mais quand je tape margin(A); ou margin(B); cela ne marche pas.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
??? Error using ==> margin
Not enough input arguments.
 
Error in ==> stabilite2 at 31
margin(A);

savez vous comment corriger cette erreur ?
merci

**phryte** · 23/05/2008, 17h23

Margin utilise la FTBO en symbolique !

Tu peux la définir aussi avec le numérateur et le dénominateur sous formes matricielles :
BO=tf(N D);

**syl55** · 23/05/2008, 17h27

ça veut dire que je ne peux pas demander les marges avec mon programme ? car mes fonctions de transfert ne sont pas définies par FTBO=tf... ?

**phryte** · 23/05/2008, 17h35

Tu peux avoir les marges avec : [Gm,Pm,Wcg,Wcp] = margin(BO).
Ou le tracé avec : margin(BO).
Dans ton cas,si je ne me trompe pas :
Code :
BO=tf(Kf*Kc*Ko*2*pi,[Tau 1 0]);

**syl55** · 23/05/2008, 17h38

super ça marche.

c'est bon , il faut juste rajouter grid pour avoir la grille sur le graph;

merci à vous.

Tracé de diagrammes de Bode [Débutant]

MATLAB

Discussions similaires

Partager

Partager