Discussion :

[hollowdeadoss]

Bonjour tout le monde,

J'ai un signal et je cherches a savoir combien de fois la constante "c" est depasse par ce signal.
J'ai donc rentre dans un vecteur "X" toutes les valeurs de mon signal et j'ai essaye de faire cela :

Code :

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

1
2
3
4
 
X1=X([1:end-1]);
X2=X([2:end]);
Tmp=find (X1 <=c & X2>c);

J'ai un resultat mais quand j'essaye avec un vecteur "test" du genre

Code :

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

test=[1 2 3 4 5 6];

Ca marche pas tres bien.
Si je prends c=4.2, ou n'importe quelle autre valeur, j'ai l'impression que c'est bon mais un peu par hasard (vraiment dependant de mon vecteur on va dire).

De plusla constante "c" est en realite un vecteur donc je pensais faire une boucle et pour chaque valeur de "c" compter le nombre de fois ou la valeur est depassee.

Avez vous des conseils ?

Hollow

[mr_samurai]

Le code est correcte, pourquoi "Ca marche pas tres bien" ?

++

[hollowdeadoss]

En fait c'est dans la conceptualisation de mon programme que j'ai "foire".

En effet sur ce signal je veux savoir combien de fois mon processus croise la constante "c".

La je regarde seulement la valeur ou il croise ... enfin avec le code que je t'ai montre, j'obtiens que un seul indice, or je vois bien sur mon signal qu'il croise la constante plusieurs fois ...

[Caro-Line]

Est-ce que tu pourrais essayer de définir mieux "dépasser", "croiser" ?
J'avoue que ce n'est pas très clair

Donne nous un exemple (avec une courbe non monotone ce sera peut-être plus parlant), avec le résultat que tu espères avoir.

[hollowdeadoss]

Ouai c'est vrai que c'est pas super clair

Alors imaginons que j'ai une sinusoide quelconque sur [-4pi 4pi] pour faire simple
Si la constante c=0.3 on va dire.
J'essaye de regarder le nombre de fois ou la sinusoide coupe 0.3 en montant ... je sais pas si c'est plus clair

Pour expliquer le coup du "montant" si on a la fonction carre, la courbe croise deux fois la valeur 0.3 mais une seule fois quand la fonctione st croissante.

Merci de votre aide

Hollow

[Caro-Line]

Par rapport à ce que tu écrit ton code semble toujours correct.
Donc il y a un problème entre ce que tu veux et ce que nous on comprend.

Peux-tu donner un exemple concret pour lequel ton code serait faux ?

[mr_samurai]

essayes comme ceci :

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
 
 
clc
close all
clear all
 
n = 40
X = 10 * rand(1,n);
IX = 1:n;
c = 8;
 
S = find((X(1:end-1) - c) .* (X(2:end) - c) < 0);
 
plot(IX,X)
hold on
scatter(IX(S),X(S))
plot(IX, c * ones(size(IX)))

++ bonne chance

[hollowdeadoss]

Voici mon code avec ma boucle

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
 
 
u=[-0.06:0.01:0.06]; %mes differents niveaux
ZZ=Z1(:,2); %les valeurs de mon signal
 
A1=ZZ([1:end-1]);
A2=ZZ([2:end]);
 
 
N=length(u);
lc=zeros(N,2); %de facon a avoir le nombre de depassement et le niveau %correspondant
 
for i=1:N
    lev=u(i);
    tmp=0;
    j=find(A1<=lev & A2>lev);
    TF=isempty(j);
    if TF==0
        NN=length(j);
        tmp=tmp+NN;
    end
    lc(i,1)=tmp;
    lc(i,2)=lev;
end