Discussion :

[jems87]

Bonjour,

Je souhaite réaliser un filtre passe-bande sur son .wav, écouter la musique avant et après filtrage pour voir la différence tout en agrémentant le programme de courbes de réponses fréquentielle et autres...

J'ai du mal à réaliser le filtre proprement dit. J'ai essayé avec filter, puis fir1 mais je ne sais pas trop comment faire au final.

Voila ou en est rendu mon code pour l'instant :

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
42
43
44
45
clear
clc
[x,Fe,Nbits]=wavread('Light_My_Fire.wav');  %On importe le fichier audio
Fe      %Fréquence d'échantillonnage 
Nbits   %Nombre de bits
g = x(:,1); %voie gauche = toutes les lignes, colonne 1
d = x(:,2); %voie droite = toutes les lignes, colonne 2
col=['r','g','b'];
subplot(3,1,1)
plot(g), hold on %On affiche l'enveloppe du signal audio de la voie gauche
title('Signal audio voie gauche','FontSize',16);
grid on;
subplot(3,1,2)
plot(d), hold on %On affiche l'enveloppe du signal audio de la voie gauche
title('Signal audio voie droite','FontSize',16);
grid on;
subplot(3,1,3)
plot(x), hold on %On affiche l'enveloppe du signal audio de la voie gauche
title('Signal audio en entier','FontSize',16);
grid on;
 
 
 
%sound(g,Fe) %restituer le son de la voie gauche
%sound(d,Fe) %restituer le son de la voie droite
%sound(x,Fe) %restituer tout le son
 
figure
z=fft(x);
title('Transformé de Fourier du signal','FontSize',16);
grid on;
 
 
%Appliquons lui maintenant un filtre passe-bande
 
wn=2*2*[1e3 2e3]/Fe
y=fir1(2,wn);   %Filtre d'ordre 2
figure
freqz(y,1);
title('Réponse fréquentielle du signal filtré','FontSize',16);
grid on;
figure
impz(1,y);
title('Réponse impulsionnelle du signal filtré','FontSize',16);
grid on;

Avez-vous des remarques, est-ce que mon filtre est bon?

[phryte]

Bonjour.
Je te propose ces modifications :

Code :

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

1
2
3
4
5
figure(2)
z=fft(x);
title('Transformé de Fourier du signal','FontSize',16);
plot(abs(z(1:length(z)/2)))
grid on;

Code :

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

y=fir1(48,wn);   %Filtre d'ordre 2

Les filtres à réponse impulsionnelle finie nécessitent un ordre élevé.
Après si tu veux filtrer le signal par un passe-bande tu fais (avec n pair) :

Code :

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

1
2
[num,den]=butter(n,[f1 f2]/(Fe/2))
x_filt=filter(num,den,x)

[jems87]

Super!

ça marche très bien avec le filtre butterworth.
Par contre avec fir1 ça ne fonctionne pas, en même temps c'est logique parce que mon signal d'entrée x n'intervient pas du tout dans la fonction fir1

Code :

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

1
2
3
wn=2*2*[1e3 2e3]/Fe;
y=fir1(48,wn);   %Filtre d'ordre 48
sound(y,Fe)

Comment faire?

Merci

[phryte]

Bonjour.

Par contre avec fir1 ça ne fonctionne pas

La primitive "fir1" ne filtre pas mais permet la synthèse du filtre RIF (ou RII).
Exemple (en fréquences normalisées):

Code :

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

1
2
3
clear
num=fir1(48,[0.3 0.7]);
freqz(num,1,512);

C'est-à-dire calcule Gain et Phase du filtre (mais pas la fonction de transfert.
C'est pourquoi je t'ai conseillé Butterworth.
Tu peux aussi utiliser "filter":

Code :

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

y=filter(num,den,x)

[jems87]

D'accord je comprend mieux. Merci

Je m'éloigne un peu du sujet maintenant mais comme l'amplitude et la phase donne des courbes assez moches (à cause du fait que l'ordre fir1 soit à 48), j'essaye d'ajouter une courbe de tendance avec le Basic fitting mais j'ai vraiment du mal et l'aide Mathwork ne m'avance pas trop.

Encore un coup de main?

[phryte]

une courbe de tendance

Qu' as-tu programmé ?

[jems87]

ben c'est sur la courbe de la réponse fréquentielle c'est pas joli c'est la dessus que je veux mettre une courbe de tendance mais je vois pas comment fonctionne le Basic Fitting (quel paramètre sélectionner(cubic,4éme ordre...) par exemple)