Discussion :

[payetonbob]

Salut,

Je travaille sur le traitement audio.
Je cherche à déterminer la ou les fréquences fondamentales d'un signal vocal. Cela me permettrait d'enrichir le signal en amplifiant les harmoniques ou en ajoutant d'autres harmoniques.

Est-ce que vous avez des idées ?

Merci

[progfou]

Les harmoniques se calculent avec la transformée de Fourier.
Sous Matlab, c'est la commande fft qui permet de l'avoir.

Après, connaître la fondamentale se fait très vite .

Si ça ne correspond pas à ce que tu cherches, précises un peu .

[payetonbob]

Les harmoniques se calculent avec la transformée de Fourier.
Sous Matlab, c'est la commande fft qui permet de l'avoir.

Après, connaître la fondamentale se fait très vite .

Si ça ne correspond pas à ce que tu cherches, précises un peu .

En fait, ce que je veux faire, c'est enrichir le signal vocal par des harmoniques.
Donc, il faut que je détecte la fondamentale. Seulement le problème, c'est que dans un son vocal, tu n'as pas vraiment de fondamentale parce que c'est pas un son pur. Tu as plusieurs fréquences prépondérantes.
Pour l'instant j'ai découpé mon signal en plusieurs parti à l'aide de la fenetre de hamming, et je cherche la fondamentale sur chaque fenetre. Graphiquement, on la voit, mais comment faire pour récupérer directement sa valeur sans lire le graphique ?

Et pour le fenêtrage, je n'arrive à appliquer la fenetre que sur le début du signal, pas sur la suite.

Je mets le code, peut être que ce sera plus compréhensible.

Merci

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
20
21
22
23
24
25
26
27
28
29
--
clear all
[y,FS,NBITS]=WAVREAD('modulator22');
t=1024;
N=floor(length(y)/t);  %nombre de fenêtres distinctes!! de 30ms
M=3*N;                 %nombre de fenêtres de 80ms
 
subplot(2,3,1)
plot(y)
X=abs(fftshift(fft(y,FS)));
subplot(2,3,2)
plot(X)
 
y1=y(1:N*t);
 
for i=1:(M-2);
    temp = hamming(t).*y1(((i-1)*80+1):(((i-1)*80)+t)); %signal fenetré
    xham=(fftshift(abs(fft(temp,FS))));                 %fft du signal fenetré
end
 
subplot(2,3,3)
plot(temp)
 
subplot(2,3,4)
plot(xham)
 
subplot(2,3,5)
psd(temp)
--

Sinon, je ne vois pas non plus comment utiliser la fonction psd, quels sont les arguments à spécifier ?

Merci d'avance