Discussion :

[appolin]

Salut,
j'ai un signal audio transformé en fichier .mat dont je connais bien les caractéristiques: (1000 Hz et 94dB) ainsi que la fréquence d'échantillonnage : 65536 Hz.
on sait en plus que: dB=20*log(P/Po) avec Po=2*10exp(-5) Pascal

J'ai envie de retrouver la fréquence de 1000Hz et la magnitude de 94dB à partir du signal temporel que je détient.
Voici mon code en Matlab:

Code :

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

1
2
3
4
5
6
7
8
9
10
x = abs(fft(Channel_1_Data));
Fs=65536;
index = find(x == 0);
x(index) = 1e-17;
mag = 20*log10(x);
m=length(Channel_1_Data);
n=pow2(nextpow2(m));
f = (0:length(mag)-1)*Fs;
f = f(:);
plot(f,mag)

Voilà mes 2 signaux en pièce jointe( 1er signal=temporel et 2eme signal=frequentiel)
Je ne sais pas si cela est bien codé et en plus je ne sais pas comment retrouver mes 1000Hz et 94 dB

Merci d'avance pour vos réponses.

[acx01b]

salut

Code :

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

1
2
3
4
5
6
7
8
 
x = abs(fft(Channel_1_Data));
mag = 20*log10(x+1e-17);
 
Fs=65536;
f = (0:length(Channel_1_Data)-1)*Fs;
 
plot(f,mag)

m'a l'air bon (équivalent à ce que tu as mis)
donc c'est ton Channel_1_Data qui ne l'est pas
c'est bien un vecteur ligne ?