Discussion :

[syl55]

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]

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]

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]

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]

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]

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.