Discussion :

[sone47]

Boujour,
j'aurais besoin de renseignement sur le fonctionnement de freqz, en effet j'ai déja regardé les aides matlab mais je ne comprend pas pourquoi desfois on créé une variable a laquelle on affecte freqz desfois on en créé deux (dans leur exemple h et f) et d'autres fois on utilise freqz directement sans l'affecter a une variable pour tracer la réponse.
Quels sont les différences de ces utilisations?
Merci

[rostomus]

Bonjour,

Code :

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

1
2
3
4
5
[h,w] = freqz(b,a,l) returns the frequency response vector h and the
corresponding angular frequency vector w for the digital filter whose transfer
function is determined by the (real or complex) numerator and denominator
polynomials represented in the vectors b and a, respectively. The vectors h
and w are both of length l.

Code :

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

1
2
3
h = freqz(b,a,w) returns the frequency response vector h calculated
at the frequencies (in radians per sample) supplied by the vector w. The
vector w can have any length.

Code :

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

1
2
freqz(b,a,...) plots the magnitude and unwrapped phase of the
frequency response of the filter. The plot is displayed in the current figure window.

[sone47]

Merci bien mais la premiere chose que je suis allé voir sont les aides Matlab et c'est justement pour cela que je demandé de l'aide car je ne comprend pas á quel moment utiliser quel forme. Car mes résultats selon les cas different, ou ne correspondent pas avec ma théorie.
Pour faire simple, j'utilise un filtre IIR passe bande.
Ma reponse est donc:

Code :

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

yf=filter(a,b,xf);

avec a,et b coefficients et xf vecteur d'entrée.
Maintenant pour vérifier mes résultats je voudrait tracer la réponse en fréquence et la phase pour m'assurer que ce que je trouve correspond á la théorie c'est á ce moment la que je bloque sur freqz.
Je vous remercie davance.

[rostomus]

Bonjour,
je ne comprend pas comment vous voulez comparer entre "y=filter(b,a,x)" et "freqz".
"freqz" dessine la reponse fréquencielle de votre filtre en amplitude (dB) et en phase. mais si vous voulez récuperer cette réponse pour un traitement, vous utilisez "[h,w]=freqz(b,a,l)", et si vous voulez récuperer h pour des pulsations w données, vous utilisez "h = freqz(b,a,w)".

est-ce que ca repond a la question?

[sone47]

Bonjour
merci je testerai donc la solution [h,w]=freqz(b,a,l)
j'affiche ensuite h en fonction de w?
Merci

[Amokrane]

Bonjour
merci je testerai donc la solution [h,w]=freqz(b,a,l)
j'affiche ensuite h en fonction de w?
Merci

plot(w,abs(h));