Discussion :

[rollo]

Bonjour à tous,

Je cherche à analyser un signal physiologique humain, la respiration.
Sur le signal brut (fréq échantillonnage 10 Hz, voir image), on peut voir que le signal est composé de 2 signaux :
- une onde "porteuse", quand le sujet inspire et expire (fréquence rapide, ~0.35Hz),
- une onde "modulante", quand le sujet a des périodes où il respire plus ou moins fort, à une fréquence + longue (~0.03 Hz).

Quand j'applique une FFT, l'analyse spectrale révèle bien le pic à la fréquence respiratoire (0.35 Hz) mais rien sur l'oscillation plus longue (0.03 Hz).
J'ai également appliqué plusieurs filtres passe-bas et bande passante (simple et butterworth) : rien non plus, pas de pic autour de 0.03 Hz.
Est-ce parce que l'amplitude de cette longue oscillation est trop faible (j'en doute, elle est quand même nettement visible) ou bien y a-t-il un réglages FFT ou filtre à affiner ?

Merci d'avance pour votre aide !
Eric

PS : données brutes en pièce jointes

[Ludelphi7]

Je ne comprends pas ton pb : Le spectre que tu présentes montre une raie proche du zéro hz, précisément la où tu cherches 0,01hz ...
ludelphi7

[tachmou]

Salut

Tu as vu la durée du signal que tu analyses ? Si tu cherches des fréquences à 0.01 Hz plus visibles, rajoutes moi des enregistrements de plus longue durée, ton 0.01 Hz sera beaucoup plus visible !

Et utilises peut-être une échelle log ?

Cordialement,

[Shagocode]

Bonjour,

j'ai déjà eu le même soucis de la fft trop "bête" pour voir la "modulante".
Du coup j'ai utilisé la transformée de hilbert pour chopper l'enveloppe :

Code :

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

1
2
3
4
5
6
 
mo = sin(2*pi*1e5*t); %signal qui module
testo = mo.*sin(2*pi*5e6*t); %signal modulé
Ht = hilbert(testo); % transformée de hilbert du signal
aHt = abs(Ht); % on prend le module
plot(abs(fft(aHt))) % -> ici on voit bien le beau pic à 1e5 Hz du signal qui module

voilà c'est à l'arrache mais ça marche bien