Discussion :

[Webby1234]

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

[kmaniche]

Bonsoir,

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

Tu peux montrer ton code, on essayera de voir ...

[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

[phryte]

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

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
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])

[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

[yann.falevoz]

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