Discussion :

[neptune2024]

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 :

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

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

2 ème règle

Code :

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

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

[Jerome Briot]

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 ?

[neptune2024]

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 :

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

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

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

Code :

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

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 :

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

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.

[neptune2024]

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

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
 
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