Discussion :

[med.doc]

Bonjour,

je dispose de l'expression d'une fonction de transfert :

L(f)=omega/(omega-2*pi*f)

et je souhaite tracer le réponse fréquentielle en amplitude et en phase sachant que ma Fc=0.843 et omega=2*pi*Fc*(1+i/175);

voilà ce que je fais sous MATLAB/Octave:

Code :

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

1
2
3
4
5
f=0.6:0.05:1.2;
L=omega_c./(omega_c-2*pi*f);
fe =0.843;
plot(f,abs(L))
plot(f,arg(L))

Malheureusement je ne trouve pas le résultat espéré.
Pouvez-m'aidez SVP

[Invité]

Bonjour,

La fonction arg n'existe pas, c'est angle.
Tel quel, tu ne visualiseras que la seconde courbe.

Code :

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

1
2
subplot(211), plot(f,abs(L))
subplot(212), plot(f,angle(L))

[med.doc]

bonjour, je vous remercie de votre réponse.
sous octave, le résultat de la fonction arg et angle est le même.

le problème c'est que les graphes que j'obtiens ne correspondent pas à ce que je dois trouver.

ce que je suis sensé trouver est explicité dans les slide 7 et 8 de ce lien :
http://syrte.obspm.fr/journees2010/powerpoint/zotov.pdf

merci de votre aide

[med.doc]

je dispose de l'expression d'une fonction de transfert :
L(f)=omega/(omega-2*pi*f)
et je souhaite tracer le réponse fréquentielle en amplitude et en phase
normalement ceci est le résultat, mais quand je le fais sous octave, ce n'est pas ce que je trouve
voici mon cade MATLAB :

Code :

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

1
2
3
4
5
6
7
8
9
10
Fc=0.843;
f=0.0:0.05:1.2;
omega_c=2*pi*Fc*(1+i/2*Q);
for i=1:length(f)
	L(i)=omega_c/(omega_c-2*pi*f(i));
end
figure(1)
plot(f,abs(L))
figure(2)
plot(f,angle(L))

les graphes ne donnent pas la même chose
pouvez m'aider svp ?
merci

[FLB]

Salut,
tu as oublié une parenthèse (2*Q) :

Code :

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

1
2
3
4
5
6
7
8
9
Fc=1/433;
Q = 175;
f=(0.6:0.01:1.2)./365;
omega_c=2*pi*Fc*(1+sqrt(-1)/(2*Q));
L=omega_c./(omega_c-2*pi*f);
figure(1)
plot(f*365,abs(L))
figure(2)
plot(f*365,angle(L))

[med.doc]

effectivement !
merci beaucoup