Comment verifier que les fréquences aux délà de 20Hz et 20KHz sont non audibles

**Contact2012** · 09/01/2009, 01h25

Bonjour;
Est ce que quelqu'un peut m'aider à mettre en oeuvre le principe que l'oreil humain ne peut pas s'apercevoir que des fréquences entre 20Hz et 20KHz.
Alors, j'ai pris un signal, j'ai calculé son FFT. puis j'ai supprimé les fréquences au déla de l'intervalle d'audibilité [20H,20KHz].
Le morceau de code est donné comme suit:

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
clear all;
[xn,fs]=wavread('rooster44100.wav');
sound(xn,fs);
% Analyse spectrale
Xn = fft(xn);
Xn_shift = fftshift(Xn);
A=abs(Xn);
P=angle(Xn);
N=length(Xn);
f = (-N/2:N/2-1).*fs/N;
plot(f,A);
%% retrouver le son original
s_orig=real(ifft(A.*exp(i*P)));
sound(s_orig,fs);
%supp les freq inf 20 Hz
a=find(f>-20 & f<20);
A2=A;
A2(a)=0;
figure(1);plot(f,A2);
s_mod=real(ifft(A2.*exp(i*P)));
sound(s_mod,fs);
%supp les freq sup 20 KHz
A1=A;
b = find(f>22000 | f<-22000);
A1(b)=0;
figure(2);plot(f,A1);
s_mod=real(ifft(A1.*exp(i*P)));
sound(s_mod,fs);

Mais le pb que je sens la modification faite malgré qu'elle a affecté l'intervalle non audible

.
Merci pour toute aide

**FR119492** · 09/01/2009, 14h11

Salut!

Comment vérifier que les fréquences aux delà de 20Hz et 20KHz sont non audibles

C'est un problème de physiologie et non d'informatique; ça n'a donc rien à voir avec Matlab et, par conséquent, rien à faire dans ce forum.
Pour résoudre ton problème, tu dois prendre un échantillon représentatif d'individus et déterminer expérimentalement leur audiogramme en les soumettant

à des sons purement sinusoïdaux de différentes fréquences et différentes intensités;
aux mêmes sons avec superposition d'un bruit.

Jean-Marc Blanc

**phryte** · 09/01/2009, 15h11

Bonjour.

sound(xn,fs);

b = find(f>22000 | f<-22000);

AS-tu bien fs >= 44 000 Hz ?

**Contact2012** · 09/01/2009, 15h43

Bonjour;
en fait, fs = c'est la fréquence d'échantillonage
et f=les fréquence dans la fft
C'est ca la question si j'ai bien compris?

SVP FR119492, vous voulez dire quoi par 'audiogramme'?
et aussi j'ai pas compris :

1. à des sons purement sinusoïdaux de différentes fréquences et différentes intensités;
2. aux mêmes sons avec superposition d'un bruit.

Merci d'avance.

**phryte** · 09/01/2009, 15h47

en fait, fs = c'est la fréquence d'échantillonage

Si tu testes une fréquence de 22000Hz il faut que fs soit > à 2 fois cette valeur ! SHANNON

**Contact2012** · 09/01/2009, 16h04

Ah oui, j'entends parler du théorème de shanon. Mais j'ai pas une grande idée sur son application.

oui fs=44100Hz (la fréq d'echatillonage du fichier audio en entrée)>=2fmax
fmax = étant la fréq maximale 22050Hz
C'est bien ça?
Mais ce que je veux verifier c'est le "domaine d'audibilité"
je veux supprimer les fréquences hors 20-20000Hz dans la fft et vérifier que c'est pas perceptible en audio (malgré que le sifgnal a été modifié mais c'est pas perceptible).
Merci pour vos conseils et aides.

**phryte** · 09/01/2009, 17h04

J'ai donné un exemple de filtrage par fft - ifft là :
http://www.developpez.net/forums/d54...es-signal-wav/

**Contact2012** · 09/01/2009, 17h07

Merci Phryte. Je vais voir le lien et vous tenir courant de mon avancement.
Bon travail

**Contact2012** · 09/01/2009, 22h40

Oui c'est le meme principe presque (malgrés que j'ai pas compris le gabarit

)
Moi je veux aussi supprimer des intervalles de fréq dans la fft. Puis reconstruire le signal pour voir l'effet de cette modification. En fait il parait que le signal n'a pas de fréq dans l'intervalle 0-20Hz et >20KHz

Je vais voir avec d'autres signaux.
Merci pour tous vos aides Phryte

**acx01b** · 10/01/2009, 13h49

salut

je propose

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
clear all;
close all;
 
fmin = 20;
fmax = 20000;
 
[xn,fs]=wavread('son.wav');
Xn = fftshift(fft(xn));
N=length(Xn);
f = (-N/2:N/2-1).*fs/N;
Xn(find(f >= -fmin & f <= fmin))=0;
Xn(find(f <= -fmax | f >= fmax))=0;
 
xn2 = real(ifft(ifftshift(Xn)));
sound(xn2,fs);
plot(f,abs(Xn));

**Contact2012** · 10/01/2009, 13h56

Salut;
Ben, je ne vois aucune différence avec le code que j'ai mis dans mon premier post.
peut tu m'expliquer stp c'est quoi l'apport.
Merci pour ton aide

**acx01b** · 10/01/2009, 14h04

moi j'en vois

en fait dans ton code tu confonds la version shift et la version non shift

Code :

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

1
2
f = (-N/2:N/2-1).*fs/N; % version shift
f = [(0:N/2), (-N/2+1:-1)] .*fs/N; % version non shift

**Contact2012** · 10/01/2009, 14h06

Oui, c'est plus simplifié en fait, et plus clair donc facile à comprendre mais ca donne le meme résultat que mon premier code. Je viens de le tester

**Contact2012** · 10/01/2009, 14h11

Ah oui, c'est pourquoi j'étais obligé à utiliser les fonctions angle et abs pour m'en sortir. Je connais pas la fonction ifftshift.
Merci bcp, j'ai remarqué comment t'as utilisé ifft et ifftshift pour retrouver le signal original

Si tu peux m'aider aussi, je dois tester un signal avec des fréquences dans l'intervalle non audible. comment je pourrais avoir un tel signal

Merci de toute aide.