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.