Spectre audio Scilab

Bonjour,
Je fais actuellement un projet tuteuré sur le traitement d'un signal sonore. Je dois traiter le signal avec scilab, représenter la FFT et le filtrer afin de garder uniquement le fondamental. Jusque là tout va bien.
Seulement j'ai un problème avec le spectre qui (à mon avis) n'est pas ce lui que je devrais trouver (je le considère à l'envers ...)! Si quelqu'un pouvait me dire si mon spectre est bon, ou si il y a une phrase à modifier dans le programme, ça serait vraiment sympa ! ça me bloque depuis maintenant 2 mois et je peux difficilement avoir de l'aide avec les professeurs ... Si jamais je n'ai pas été clair (ce que je peux comprendre ! =] ) posez vos questions ! Merci d'avance !!


Voici mon programme :


clear;


[y,fe,bits]=wavread("C:\Users\Thibault\Documents\Projet Tuteuré 2009\BANJO_HI.wav");fe,bits


playsnd(y,fe)

//Caractéristiques de départ
sig=y';
[L,C]=size(sig);
N=length(sig);

//Transformée de Fourrier du signal
spect=fft(sig,-1);


f=fe*(0:N-1)/N;
n=size(f,'*');
Te=1/fe
t=Te*(0:N-1);

Z=[spect(1:N)];
//définition du filtre du filtre passe bas, butteworth d'ordre 10

G=analpf(10,'butt',[0 0],2*%pi*1000);
GAIN=freq(G(2),G(3),%i*2*%pi*f);
G=[GAIN(1:N)];
//
//filtrage
SIGPB=G'.*spect;
sigpb=fft(SIGPB,1);

//Affichage du spectre
xset("window",2); xbasc(0); xset("font size",4);
plot2d(f,abs(Z/N))

//affichage des graphes
//N°1
xset("window",0); xbasc(0) ; xset("font size",4);
plot2d(f,abs(spect(1:n))/N)
plot2d(f,abs(G(1:n)))
xtitle("Représentation du filtre utilisé","fréquence(Hz)","amplitude (V)")

//N°2
xset("window",1); xbasc(1) ; xset("font size",4);
plot2d(f,abs(SIGPB(1:n))/n)
xtitle("Représentation du spectre filtré","fréquence (Hz)", "ampltide (V)")

//N°3
xset("window",3); xbasc(3) ; xset("font size",4);
plot2d(t,sig(1:N))
plot2d(t,sigpb,2)
xtitle("Signal + signal filtré en fonction du temps","temps (s)", "amplitude (V)")

playsnd(sigpb',fe)