Signal audio et fft

Version imprimable

25/01/2010, 18h50
Webby1234

Signal audio et fft

Bonjour à tous !!

J'ai un fichier audio au format wav.
Sous matlab, je le récupère avec un wavread et je stocke alors les échantillons dans y (variable utilisée dans l'aide)

Ensuite, pour avoir le spectre d'amplitude de ce signal je dois faire un abs(fft(y)) pour avoir les réels et abs(ifft(y)) pour les imaginaires.

Jusque là on est certainement d'accord.

Ce que je ne comprend pas c'est ce que fais la fonction fft.

Car je dois implémenter ce principe décris au dessus, sans utiliser matlab.
J'ai une fonction fft mais qui requiert à la fois les imaginaires et les réels, comment faire pour les prendre sachant que je n'ai que les y.

Merci
25/01/2010, 22h26
kmaniche

Bonsoir,

Est-ce que tu veux reprogrammer la fonction fft sous un autre langage.

Tu peux montrer ton code, on essayera de voir ...
25/01/2010, 22h50
Webby1234

En fait en Java, j'ai cette méthode :

public static Complex[] fft(Complex[] x)

avec ce constructeur :
public Complex(double real, double imag)

ET je ne sais pas quoi donner comme paramètres sachant que je n'ai que mon tableau de y

Bonjour.
Tu peux faire le calcul "à la main" un peu long !. Résultat dans sum :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clear
Fs=4410;
T=1/Fs;
t = 0:T:0.5;
N=length(t);
y = 2*sin(2*pi*50*t)+5*sin(2*pi*120*t)+ 0.5*randn(length(t));
f = linspace(0,Fs/2,N/2);
for k=1:N
    sum=0;
  for n=1:N
    sum =sum+ y(n)*exp(-j*2*pi*(k-1)*(n-1)/N);
end
 X(k) =sum/N;
end
h=semilogx(f,2*abs(X(1:floor(N/2))));
xlabel('Frequence [Hz]');
ylabel('Amplitude');grid
set(h,'LineWidth',1,'Color','r')
axis([1 1E4 0 6])

26/01/2010, 10h13
Webby1234

Merci, là dessus je suis d'accord.
Mais ma question reste la même quelle est la partie imaginaire et la aprtie réelle que je donne au calcul de la fft ?

Tkx
26/01/2010, 16h29
yann.falevoz

Citation:

Envoyé par Webby1234

pour avoir le spectre d'amplitude de ce signal je dois faire un abs(fft(y)) pour avoir les réels et abs(ifft(y)) pour les imaginaires.

Bonjour,

Ben, non, je pense pas que ça marche... pour la partie réelle : real(fft(y))
et pour la partie imaginaire : imag(fft(y))

le abs(fft(y)) c'est le module de ton spectre. Pour avoir la phase, il doit y avoir une fonction arg(fft(y)), ou sinon tu fais atan2(imag(fft(y))/real(fft(y))).

Le signal que tu lis dans le fichier wav est réel. Sa partie imaginaire est donc nulle...

Yann