Bonjour à tous,
J'ai ouvert il y a peu un sujet relatif à la densité spectrale de puissance DSP et au filtrage : ici ; cependant je n'ai reçu de réponse et je ne sais pas moi-même que faire des résultats obtenus. Je repars donc sur une autre piste.
Le problème est le suivant : je dispose de 6 pistes d'impédances qui, en résumé, sont 6 vecteurs constitués chacun d'une colonne. Je souhaiterais, à l'aide d'une analyse spectrale, pouvoir mesurer les similitudes éventuelles existantes entre ces vecteurs. Malheureusement je n'ai pas une idée claire de comment procéder. Malgré tout, j'ai débuté le travail par le calcul d'une fft et du fftshift pour chacun des 6 vecteurs. En me calquant ensuite sur ce que j'ai pu faire par le passé je représente le module du logarithme de la partie réelle des 6 vecteurs spectraux résultant. Pour rendre plus clair, mon propos voici le code utilisé et le résultat :
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
30
31
32
33 %-------------------------------------------------------------------------- %ANALYSE EN UTILISANT LA FFT. %-------------------------------------------------------------------------- valeurs_Imp1_fft=fftshift(fft2(valeurs_Imp1)); valeurs_Imp2_fft=fftshift(fft2(valeurs_Imp2)); valeurs_Imp3_fft=fftshift(fft2(valeurs_Imp3)); valeurs_Imp4_fft=fftshift(fft2(valeurs_Imp4)); valeurs_Imp5_fft=fftshift(fft2(valeurs_Imp5)); valeurs_Imp6_fft=fftshift(fft2(valeurs_Imp6)); %Visualisation des résultats issus du filtrage fft en prenant les modules %des logarithmes des parties réelles. figure(4) plot(abs(log(real(valeurs_Imp1_fft))),'c') hold on plot(abs(log(real(valeurs_Imp2_fft))),'r') hold on plot(abs(log(real(valeurs_Imp3_fft))),'g') hold on plot(abs(log(real(valeurs_Imp4_fft))),'b') hold on plot(abs(log(real(valeurs_Imp5_fft))),'k') hold on plot(abs(log(real(valeurs_Imp6_fft))),'m') hold on title('log du module de la partie réelle de la fft des impédances pour un patient donné') legend('Imp1','Imp2','Imp3','Imp4','Imp5','Imp6','Location','NorthEast') xlabel('fréquences') ylabel('impédances')
Je compte faire de même pour la partie imaginaire, module et phase. Mais au regard de cette première figure, qu'en pensez-vous ? Est-ce que quelque chose saute directement au yeux ?
Je sollicite plus généralement votre aide afin de savoir comment je pourrais concrètement mesurer les ressemblances entre 6 vecteurs. J'ai, en outre, penser à des calculs de corrélations mais je ne sais pas tout à fait comment procéder.
Je donc serais ravi si vous me donniez ne serait-ce que des idées d'études que je pourrais mener sur ces 6 vecteurs ou des calculs, et particulièrement des idées de fonctions matlab à utiliser et paramétrer.
Partager