Discussion :

[Newenda]

Bonjour,

j'ai un signal auquel j'applique la FFT. J'aimerais calculer la fréquence centrale de ce spectre, c'est à dire ici, la fréquence pour laquelle la somme des amplitudes avant celle-ci soit identique à la somme après.



mais je ne vois pas comment faire, , si quelqu'un a une idée? Je pourrais faire la somme de tout, la diviser par deux et voir approximativement vers quelle indice j'atteins cette valeur mais peut être y'a t il plus propre?

Merci

[Jerome Briot]

Code :

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

1
2
idx = cumsum(amplitude) <= sum(amplitude)/2;
frequence_centrale = frequence(idx)

Non ?

[Newenda]

Code :

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

1
2
idx = cumsum(amplitude) <= sum(amplitude)/2;
frequence_centrale = frequence(idx)

Non ?

Tiens pratique CUMSUM
Euh disons que idx dans ce cas est le vecteur du 1 au dernier ! mais merci pour l'idée de base, je pense que comme ça c'est bon :

Code :

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

1
2
idx = find(cumsum(amplitude) <= sum(amplitude)/2 , 1 , 'last') ;
frequence_centrale = frequence(idx)

non?

[Vince2505]

Il faut calculer la moyenne des fréquences pondérées par leur valeurs d'amplitudes.
ex :
f1=10hz amplitude du pic 1000
f2=20hz amplitude du pic 200
f3=30hz amplitude du pic 500

frequence centrale : (1000*10+200*20+500*30)/(1000+200+500)=17.06hz

[patrick7413]

En gros tu veux calculer le centre de gravité spectral de ton signal, voici 1 code qui fait ce calcul, (je ne l'ai pas testé)

http://www.mathworks.com/matlabcentr...exchange/11345

[Newenda]

Il faut calculer la moyenne des fréquences pondérées par leur valeurs d'amplitudes.
ex :
f1=10hz amplitude du pic 1000
f2=20hz amplitude du pic 200
f3=30hz amplitude du pic 500

frequence centrale : (1000*10+200*20+500*30)/(1000+200+500)=17.06hz

en gros

Code :

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

freq_cent_pond = sum(spect.*freq)/sum(spect)

Mais quel sens physique explique de donner plus d'importance aux amplitudes des fréquences les plus hautes?

Merci

[Jerome Briot]

Mais quel sens physique explique de donner plus d'importance aux amplitudes des fréquences les plus hautes?

Si j'ai bien compris, ça revient à un simple calcul de barycentre