Discussion :

[Contact2012]

Bonjour;
J'ai un vecteur de nombre réels (données ECG) que je veux analyser et détecter les cycles qu'il contient. Mais, il n'est pas parfaitement régulier . Voir sa représentation sous Matlab jointe au message.
Pouvez vous m'aider à trouver un algorithme pouvant extraire les cycles de mon signal?
Merci pour vos aides.

[ToTo13]

Bonjour,

est ce que le problème est bien posé ?
Il me semble que dans un ECG, il faut regarder l'onde T. Elle suit les pics. Donc est ce que retrouver tous les pics positifs et/ou négatifs puis le "mini" pic suivant ne suffirait pas ?

[Contact2012]

Merci ToTo pour ta réponse.
En fait, lorsque j'ai un signal bien régulier, je procède par détection du complex QRS : un pic maximum très apparant par rapport aux autres ondes. Puis, une fois le signal ECG est décomposé en morceaux commenceant pas onde R, je fais translater pour que ca commence par l'onde P.
Mon programme de décomposition en cycle est le suivant :

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;
[xn,fs]=wavread('300-0-598-0--sig0.wav');
m=max(xn);
tres=0.7*m;
peaks=[];
% detection des pics dans le signal ecg
for i=1:length(xn)
    if xn(i)>tres
        peaks=[peaks i];
    end
end
%elimination des valeurs répétées dans peaks et calcul du cycle length
distance=[];
peaks2=[];
for j=1:length(peaks)-1
    % la longueur d'un cycle est sup à 200
    if (peaks(j+1)-peaks(j))>200
        distance = [distance peaks(j+1)-peaks(j)];
        peaks2=[peaks2 peaks(j)];
    end
end
cycle_len=fix(mean(distance));
% reconstruction du signal
xs=[];
for i=1:length(peaks2)-1
        cyc = xn(peaks2(i):peaks2(i+1));
        xs=[xs ; cyc];
end
plot(xs);

Mais avec le signal joint (qui n'est pas régulier), cet algorithme ne me donne pas de bon résultats.
Merci pour toute aide.

[FR119492]

Salut!
Le nom scientifique de ton type de problème est série temporelle ou série chronologique, en anglais time serie. Il y a une masse de publications sur ce sujet.
Jean-Marc Blanc

[Contact2012]

Salut;
Oui c'est vrai ca s'appelle "times series data" en anglais.
stp, tu peux me fournir des liens ou je pourrais trouver de l'aide?
Merci.

Bonjour;
SVP, est ce que vous pouvez me conseiller d'autres forums et/ou sites ou je pourrais trouver une solution pour mon pb?
Merci à tous.

[FR119492]

Salut!
Tu peux commencer par regarder dans Wikipedia en français et en anglais. Ce n'est en effet pas du tout mon domaine d'activité, mais tu dois apprendre à chercher l'information.
Jean-Marc Blanc

[cite olympique]

Merci ToTo pour ta réponse.
En fait, lorsque j'ai un signal bien régulier, je procède par détection du complex QRS : un pic maximum très apparant par rapport aux autres ondes. Puis, une fois le signal ECG est décomposé en morceaux commenceant pas onde R, je fais translater pour que ca commence par l'onde P.
Mon programme de décomposition en cycle est le suivant :

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;
[xn,fs]=wavread('300-0-598-0--sig0.wav');
m=max(xn);
tres=0.7*m;
peaks=[];
% detection des pics dans le signal ecg
for i=1:length(xn)
    if xn(i)>tres
        peaks=[peaks i];
    end
end
%elimination des valeurs répétées dans peaks et calcul du cycle length
distance=[];
peaks2=[];
for j=1:length(peaks)-1
    % la longueur d'un cycle est sup à 200
    if (peaks(j+1)-peaks(j))>200
        distance = [distance peaks(j+1)-peaks(j)];
        peaks2=[peaks2 peaks(j)];
    end
end
cycle_len=fix(mean(distance));
% reconstruction du signal
xs=[];
for i=1:length(peaks2)-1
        cyc = xn(peaks2(i):peaks2(i+1));
        xs=[xs ; cyc];
end
plot(xs);

Mais avec le signal joint (qui n'est pas régulier), cet algorithme ne me donne pas de bon résultats.
Merci pour toute aide.

[Contact2012]

le programme est déjà donné dans un post précédent (voir la discussion)