Discussion :

[chris1495]

Bonjour, j'ai plusieurs problèmes de programmation donc je fais appel à vous pour venir à mon soutien :

- j'aurais tout d'abord voulu savoir comment détecter le temps lorsqu'une courbe (ici, un vecteur) franchit un certain seuil (qui dans mon cas est fixé à 12) que ce soit en front montant ou descendant, ne garder que les valeurs en dessous mais étant les plus proches de ce seuil et insérer à ce vecteur que l'on vient de créer les valeurs de début et de fin du temps (qui correspondent au début et à la fin de l'expérimentation)

- ensuite j'aurais aimé savoir comment former des couples de points à partir d'un vecteur

-enfin je voudrais savoir comment "découper" un vecteur par plusieurs couples de points

Merci d'avance pour votre aide

[soft001]

Montres nous ce que tu as fait

[chris1495]

Malheureusement, ceci constitue le début de mon programme Matlab et je n'ai pas pu faire de capture d'écran pour vous éclaircir mais je peux apporter quelques précisions si vous le souhaitez.

Mon vecteur possède donc plusieurs pics et ces pics passent par un seuil que j'ai fixé et j'aimerais en fait récupérer les valeurs inférieures à 12 mais étant le plus proche de 12 lorsque ces pics sont en front montant ou descendant (ce qui par la suite créera un nouveau vecteur qui me permettra de créer (mais je ne sais pas actuellement comment) des couples de points).

[Jean Dumoncel]

Bonjour,

Malheureusement, ceci constitue le début de mon programme Matlab et je n'ai pas pu faire de capture d'écran pour vous éclaircir

Tu n'as pas besoin de faire une capture d'écran, il suffit de faire un copiez-coller de ton code dans un message sur le forum, sans oublier d'utiliser les balises code :


Pour ton problème, tu peux regarder la fonction find.

[chris1495]

Pour l'instant, je n'ai pas de code puisque cette étape doit constituer le début de mon programme (si ce n'est simplement un vecteur en entrée).

J'ai déjà pensé à la fonction find mais le problème, c'est comment ne sélectionner que les valeurs les plus proches de 12 et non pas toutes celles qui se situent en dessous

[Jean Dumoncel]

Il y a plusieurs façons de procéder :

• soit tu parcours chaque valeur de ton vecteur, et tu détermines les passages au-dessus et en-dessous de 12, et en fonction, tu sauvegardes les valeurs correspondantes.
• soit tu travailles sur le vecteur complet et tu utilises des fonctions comme diff et find pour déterminer les changements de croissance/décroissance de ton vecteur

[chris1495]

Ok, merci pour l'indication. Il y a peut-être moyen de bricoler quelque chose avec les indications que tu viens de me fournir. En fait, pour mon premier "problème", je voudrais effectuer un seuillage (à 12) et ne garder que les valeurs les plus proches mais inférieures à ce seuil (exceptées les valeurs de début et fin d'expérience).

Sinon, concernant les autres questions, quelqu'un aurait-il une idée ?

[Jean Dumoncel]

Tes autres questions ne sont pas très claires...

Si tu pouvais essayer de reexpliquer en donnant un exemple de vecteur et le résultat que tu souhaites?

[chris1495]

Ok, prenons comme exemple le vecteur temps T :

T=[0 5 10 15 20 25 30 35]

Comment à partir de ce vecteur pourrais-je former les couples de points suivants : 0-5 ; 10-15 ; 20-25 ; 30-35 ?

Comment ensuite à partir de ces couples de points pourrais-je redécouper un autre vecteur ? (pour ensuite reformer plusieurs nouveaux vecteurs qui ne seraient que des parties du vecteur "découpé" par les intervalles (ou couples de points) ci-dessus).

[soft001]

Tu veux dire comme ça :

Code :

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

1
2
3
T=[0 5 10 15 20 25 30 35]
[m,n]=size(T);
T1=reshape(T',2,n/2)

[chris1495]

Pas vraiment. En fait, il faudrait fractionner le vecteur temps T pour former des couples de points comme ci-dessus, mais comment ? ^^,

[soft001]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
>> T1(:,1)
 
ans =
 
     0
     5
 
>> T1(:,2)
 
ans =
 
    10
    15

[Invité]

Si j'ai bien compris, le code de soft001 peut tout de même t'aider :

Code :

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

1
2
3
4
5
6
7
T=[0 5 10 15 20 25 30 35]
[m,n]=size(T);
T1=reshape(T',2,n/2)
pas = ...
for i = 1:n/2
   T1(1,i):pas:T1(2,i) 
end

[chris1495]

Merci à tous ! Grâce à vous petit à petit j'avance. Comment fait-on pour "fragmenter" un vecteur à partir de ces couples de points ?

[Invité]

Pourrais-tu définir plus précisément ce que tu souhaites faire ?

[chris1495]

Bien sûr ^^,. Je vais de nouveau prendre un exemple pour illustrer.
Nous avons un vecteur A :
A = [1 3 2 6 9 1 7 5]
Il y a aussi le vecteur temps T défini précédemment de la manière suivante :
T = [0 5 10 15 20 25 30 35].
Nous avons réussi à former les couples de points suivants : 0-5 ; 10-15 ; 20-25 ; 30-35
A partir de ces couples de points, comment fragmenter le vecteur A pour former :
- le vecteur B = [1 3] et qui correspond à la valeur du vecteur A entre le temps T compris entre 0 et 5
- le vecteur C = [2 6] et qui correspond à la valeur du vecteur A entre le temps T compris entre 10 et 15
- le vecteur D = [9 1] et qui correspond à la valeur du vecteur A entre le temps T compris entre 20 et 25
- le vecteur E = [7 5] et qui correspond à la valeur du vecteur A entre le temps T compris entre 30 et 35

Voilà ^^,

[Invité]

Rien ne change : c'est le même principe que précédemment.

[ol9245]

Bien sûr ^^,. Je vais de nouveau prendre un exemple pour illustrer.
Nous avons un vecteur A :
A = [1 3 2 6 9 1 7 5]
Il y a aussi le vecteur temps T défini précédemment de la manière suivante :
T = [0 5 10 15 20 25 30 35].
Nous avons réussi à former les couples de points suivants : 0-5 ; 10-15 ; 20-25 ; 30-35
A partir de ces couples de points, comment fragmenter le vecteur A pour former :
- le vecteur B = [1 3] et qui correspond à la valeur du vecteur A entre le temps T compris entre 0 et 5
- le vecteur C = [2 6] et qui correspond à la valeur du vecteur A entre le temps T compris entre 10 et 15
- le vecteur D = [9 1] et qui correspond à la valeur du vecteur A entre le temps T compris entre 20 et 25
- le vecteur E = [7 5] et qui correspond à la valeur du vecteur A entre le temps T compris entre 30 et 35

Voilà ^^,

Ca me parait tordu ton truc.
et ce qui se passe entre 5 et 10 secondes c'est pour les cochons ?

pour sélectionner un front montant qui franchit la valeur 12 :

Code :

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

k = find(M<12 & M([2:end 1])>12) ;

Ca te donne l'indice du point justa avant de franchir le front.

Si tu cherche la date de cet événement dans un vecteur temps :

Code :

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

t = Temps(k) ;

L'avantage de cette approche c'est que tu détectee d'un seul coup TOUS les fronts montants.
Si tu veux que le premier :

Code :

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

premier = k(1) ;

si tu veux pas que ton code t'insulte si il n'a rien trouver :

Code :

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

1
2
3
4
if isempty(k)
disp('pas de front montant')
return
end

[chris1495]

Ce que j'ai donné était un exemple, ce n'est pas ce que j'obtiens réellement.

Et pour détecter les fronts descendants aussi ?

[chris1495]

En fait, c'est bon, j'ai trouvé pour les fronts descendants ^^,.

Mais je reviens alors au problème suivant.
A partir de ces couples de points, comment fragmenter le vecteur A (A = [1 3 2 6 9 1 7 5]) pour former :
- le vecteur A1 = [1 3] et qui correspond à la valeur du vecteur A entre le temps T compris entre 0 et 5
- le vecteur A2 = [2 6] et qui correspond à la valeur du vecteur A entre le temps T compris entre 10 et 15
- le vecteur A3 = [9 1] et qui correspond à la valeur du vecteur A entre le temps T compris entre 20 et 25
- le vecteur A4 = [7 5] et qui correspond à la valeur du vecteur A entre le temps T compris entre 30 et 35

Attention, je me suis rendu compte que j'avais fait une légère erreur en prenant mon exemple car j'ai pris pour le temps un pas régulier alors qu'en réalité, il n'y a pas de concordance entre les valeurs qui composent le vecteur temps T.