Discussion :

[ammoula1]

Bonjour ,
je suis débutante en Matlab ,
je fais un projet de traitement de la parole ,

je veux faire l'analyse spectrale d'un signal wav , j'ai écrit un code mais ce dernier marche avec les signaus de musique , mais en entrant un signal de parole cette exception apparaît :

Index exceeds matrix dimensions.
Error in ==>(E:...\...\ )
on line ...==> plot(f,20*log10(abs(Y(length(f)))+eps));

voici le code :

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
 %=====================================
 
[x,fs]=wavread('digits1.wav');
Y=fft(x.*hamming(length(x)));
 
% plot of spectrum of buttom 15000Hz
 
hz=15000*length(Y)/fs;
f=(0:hz)*fs/length(Y);
plot(f,20*log10(abs(Y(length(f)))+eps));
legend('Spectrum');
xlabel('frequency(Hz)');
ylabel('Magnitude(dB)');
 
%===========================================

Toute aide est la bien venue ,
je suis totalement bloquée ,
Merci ,

[Caro-Line]

Dans la
Index exceeds matrix dimensions.

Comme a priori Y est de la même taille que f (enfin j'en ai bien l'impression), tu devrais vérifier si tu n'as pas défini quelque part une variable abs ou log10...

[ammoula1]

Bonjour ,

merci pour ton aide , mais j'ai pas défini quelque part une variable abs ou log10...
quelle est la solution lors ????

[Caro-Line]

Et bien il faut que tu comprennes pourquoi tu as ce message d'erreur en débuguant ton programme.
Cherche dans l'expression :

Code :

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

20*log10(abs(Y(length(f)))+eps)

où est le problème en la décomposant.

[Jerome Briot]

cette exception apparaît :

Il nous faut le message d'erreur exacte et complet

[fatenov]

Salut,

Tout d'abord je n'ai pas bien compris le but de cette ligne

Code :

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

hz=15000*length(Y)/fs;

Comment tu fixes ta fréquence max.

Mais bon, en essayant sur exemple, je pourrais te dire que t'as un problème de taille de vecteurs. ton vecteur de fréquence f (déterminé par hz) est plus long que le vecteur Y qui a la même taille de ton signal.wav
Lorsque tu fais appel à

Code :

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

20*log10(abs(Y(length(f)))+eps)

le programme ne trouve aucune valeur dans cet index.

Mais je ne comprends toujours pas ton choix de l'axe des fréquences
Pourquoi ne pas faire

Code :

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

1
2
3
4
5
6
Y=abs(fft(x.*hamming(length(x))));
 
% plot of spectrum of buttom 15000Hz
N=length(x);
f=(-N/2:N/2-1)./N.*fs;
plot(f,20*log10(Y));

ça te donne un spectre classique de période =fs et à fréquences centrées à zéros. Tout dépond de ce que tu cherche finalement

[ammoula1]

Bonjour ,

Merci pour ton aide ,

Code :

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

hz=15000*length(Y)/fs;

en fait , je fais la multiplicaion par 15000 car je veux m'assurer que la représentation du signal s'étend jusqu'a 15000 Hz ,

Merci encore ,

[phryte]

Essaye avec :

Code :

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

 plot(f,20*log10(abs(Y(1:length(f)))+eps));