Traitement du signal

Bonjour,

J'aimerai récupérer les parties d'un signal passé en paramètre (dont je ne connais ni la période d'échantillonnage, ni le nombre de points, la seule chose connu est son intervalle [0 1]) qui vérifie les caractéristiques suivants :
0< signal_X <0,3
et signal_X ne varie pas (0,02 ie 2%) pendant un temps supérieur à 30s.

1ère règle

Code:

---

a=find(0<Vxx_acel_pdl_rat & Vxx_acel_pdl_rat<0.3);

---

2 ème règle

Code:

---

1 2	T= (abs(diff(Vxx_acel_pdl_rat))<=0.02); toto=find(T==1);

---

après je ne vois pas trop comment faire apparaitre la notion de temps et si c'est possible.

Un exemple de signal_X (appuie pédale (0 à 100%), je souhaite retrouver les parties entourées en rouge.

Merci pour votre aide :D:D:D

Je ne suis pas sur de bien comprendre...

Si tu connais la durée totale de ton signal et le nombre de points mesurés, tu peux en déduire la fréquence d'échantillonnage, non ?

Oui effectivement, j'ai facilement récupérer la fréquence d'échantillonnage.
J'ai pris le raisonnement autrement :

- Récupérer les points du signal vérifiant : 0 < Signal_X < 0,3

Code:

---

a=find(0<Signal_X & Signal_X<0.3);

---

- Trouver les intervalles supérieur ou égal à 30s sur ce points.

Code:

---

1 2	b=find(diff(a)~=1); % Recherche les points consécutifs c=find((diff(b)*mean(diff(PériodeEchan)))=>30) ; % Recherche les intervalles supérieurs à 30 secondes.

---

- Sur les intervalles supérieurs à 30s, vérifier abs(min(Signal_X)-max(Signal_X)) =< 0,02.

Code:

---

Signal_X(a(b(c(1):a(b(c(2))) % 1 er intervalle

---

.

Glisser sur une période de 30 s sur le premier intervalle comme expliquer dans le jpeg associé, décaler d'un point re-vérifier et ainsi de suite jusqu'au dernier point l'intervalle.
Recommencer le traitement sur le second intervalle.

Avez-vous des idées sur la réalisation, je bloque au niveau du code réalisant le "glisser" de 30 secondes sur l'intervalle.

Merci de votre aide.

Bon cela ne fait exactement ce que je voulais faire mais je partage le code si des personnes souhaitent effectuer des choses similaires.

Code:

---

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

a=find(0<Signal_X & Signal_X <0.3); b=find(diff(a)~=1); Sample_time=mean(diff(TimeVector)); if length(b)==1 Inter_temp{1}=a(1):a(b(1)); else Inter_temp{1}=a(1):a(b(1)); for i=2:length(b)-2 Inter_temp{i}=a(b(i-1)+1):a(b(i)); end Inter_temp{length(b)}=a(b(length(b)-1)+1):a(b(length(b))); end k=1; Inter=''; for j=1:length(Inter_temp) if (length(Inter_temp{j})*Sample_time>=30) Inter{k}=Inter_temp{j}; k=k+1; end end if isempty(Inter) InterF=0; else InterF=cell2mat(Inter); end

---