Discussion :

[Hatorii]

Bonjour,

Je ne comprends pas un point du code suivant. Comment le filtrage peut fonctionner alors que la fréquence d’échantillonnage Fs n'est pas donnée en entrée de fir1() ou/et de filter() ?
On a la valeur 0.48 (la deuxième entrée de fir1) qui est la fréquence de coupure normalisée et équivaut à 0.48*Fs/2 . Comment ceci fonctionne puisque Fs peut changer et donc la fréquence de coupure 0.48*Fs/2 aussi ?

Il faut que Fs soit connue pour connaitre la fréquence de coupure et ainsi connaitre la largeur de la réponse impulsionnel du filtre !? Ou bien je me trompe ?

Code :

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

1
2
3
4
5
6
7
 
t = (0:length(y)-1)/Fs;
 
bhi = fir1(34,0.48,'high',chebwin(35,30));
freqz(bhi,1)
 
outhi = filter(bhi,1,y);

[Hatorii]

Personne pour m'aider?

[Cheetor]

Déjà, je crois qu'il manque la commande suivante :

Code :

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

load Chirp

Chirp est un fichier ".mat" contenant un signal pulsé, caractérisé par y et une fréquence d'échantillonnage Fs. Ce qui permet au code de fonctionner, au moins pour récupérer les informations de temps.
Après, fir1 a l'air de fonctionner avec des rapports avec la fréquence d'échantillonnage, ce qui explique qu'il n'y ait pas besoin d'avoir Fs. (Je peux me tromper sur ce point).

[Hatorii]

oui il manque "load chirp". Mais fir1 ne va tout de même pas chercher l'info dedans si ?

Fs n'est mise en entrée ni de fir1(), ni de filter() et si j'utilise mon propre signal, la valeur Fs contenue dans "chirp.mat" ne sera plus correcte. Donc je saisie pas comment il récupère l'info.

[Cheetor]

Fir1 utilise un rapport, donc si tu fais 0.48*Fs/Fs, tu obtiens 0.48 par exemple.

Je me demande d'ailleurs, mais là aussi, ça serait à confirmer si le passage par fir1 et filter n'utilise pas le concept de synthèse de filtre (du moins normalisation du filtre). Plus d'infos ici => http://xcotton.pagesperso-orange.fr/..._de_filtre.pdf

(Edit : en fait, non, là, seule la normalisation par rapport à la fréquence d'échantillonnage est présente dans FiR1, le lien n'a donc aucun rapport avec le schmilblick).

[Hatorii]

excuse moi mais je comprends pas où tu veux en venir. Car dans tout les cas, normalisation ou pas, Fs doit être connue pour permettre un rapport et retrouver la fréquence de coupure. Mais je ne saisie pas comment fir1 et/ou filter s'en sortent sans demander Fs en entrée.

[Cheetor]

Fir1 n'utilise pas la fréquence d'échantillonnage, mais il accepte en entrée des valeurs dite-normalisées (dans l'exemple que tu utilises, 0.48 est la fréquence de coupure normalisée par rapport à ta fréquence de sampling, le rapport est déjà fait par rapport à ta fréquence d'échantillonnage).

D'ailleurs, si tu regardes les figures issus de l'exemple, on constate que l'axe des x n'est pas en fréquence, mais en fréquence normalisée. Si tu veux obtenir le signal avec un axe en fréquence, il te faut alors multiplier ton axe des x, par rapport à ta fréquence Fs.

Si tu veux utiliser ton propre signal, il te faut connaitre Fs, mais il te faut aussi calculer les valeurs normalisées, pour pouvoir utiliser le filtre.

[Hatorii]

Mon point est là justement. J'ai bien compris ce que veux dire le 0.48, mais pour que cette valeur corresponde à la fréquence de coupure souhaité, il doit y avoir un calcul du type Fc = 0.48*x quelque part. Autrement le 0.48 ne veut rien dire.
ors ici le x correspond à Fe/2 (la fréquence de coupure est calculé à partir de le fréquence de Nyquist ), mais je ne saisie pas ou cette infos est trouvé par fir1 ou filter.

Ps : En réalité, l'exemple de Matlab montre l'axe des abscisses en (pi*radians)/sample allant de 0 à 1. Ce n'est donc pas la fréquence mais la pulsation qui est normalisé. Et ceci de 0 a pi (et pas de 0 à 2pi puis ce que l'information est redondante après pi), mais ça ne change pas le problème.

[Hatorii]

J'ai résolu mon problème je crois. Je ne sais pas si c'est ce que tu voulais dire, mais enfaite il n'y a jamais besoin de Fs, même si j'utilise mon propre signal.

fir1 prend en entrée le 0.48 qui correspond (comme on l'a déjà dit) en version dénormalisé à Fc = 0.48*Fs/2=0.48*Fn avec respectivement les fréquences de coupure, d'échantillonnage et de Nyquist.
Or, on peut calculer Fn/Fc = 1/0.48 = 2.08. Ce qui correspond à un nombre d'échantillon. Donc quelque soit la fréquence d'échantillonnage (et donc de coupure), ce nombre ne change pas.

Ainsi, le modèle à moyenne mobile qui est utilisé pour le filtrage avec Fc = 0.48*Fn, sera toujours une porte contenant deux échantillons.