Discussion :

[Nico76131]

Bonjour à tous,

J'ai l'évolution de la vitesse d'un camion en fonction du temps et je cherche à la linéariser. C'est à dire, je voudrais supprimer les pics d'accélération ou de décélération et toutefois respecter l'allure globale de la vitesse.

Comme sur la courbe bleue où l'on voit ces nombreuses variations de vitesse.
Voici le premier programme que j'ai tapé : (où G3(:,1) est ma colonne vitesse)

(G3(:,2) est la colonne temps)

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
G4(1,1)=0; 
G4(1,2)=G3(1,2);
 
for i=1:length(G3)-20
 
    if G3(i,1)<G3(i+1,1) && G3(i+1,1)>G3(i+20,1) || G3(i+1,1)<G3(i+20,1);
        G4(i+1,1)=G3(i,1);
        G4(i+1,2)=G3(i+1,2);
 
    elseif G3(i,1)>G3(i+1,1) && G3(i+1,1)<G3(i+20,1) || G3(i+1,1)>G3(i+20,1);
        G4(i+1,1)=G3(i,1);
        G4(i+1,2)=G3(i+1,2);
 
    else G4(i+1,1)=G3(i,1);
        G4(i+1,2)=G3(i,2);
 
    end
end

En deuxième pièce jointe la nouvelle courbe de vitesse après avoir appliqué le programme.
J'ai vu que ce programme me supprimait quelque valeurs mais malheureusement pas toutes celles espérées. J'ai donc répété ce programme en changeant le pas de comparaison (ici : 20).

Quel autre programme ou quelle amélioration pourriez vous me conseiller?

Merci.

[Jerome Briot]

Pourrais-tu nous fournir un jeu de données sous forme d'un fichier mat (mis dans une archive rar ou zip) ?

Sinon, dans un premier temps, tu peux aussi t'inspirer des différents codes de cette discussion : Détection de pics

[Nico76131]

Voici une banque de données
Nouveau dossier compressé.zip

[Nico76131]

Merci pour le lien.

[Nico76131]

J'ai un petit problème, j'ai essayé de m'inspirer de la discussion sur les pics d'une courbe mais je n'ai pas réussi à faire un programme qui convient pour ce que je veux faire.

Quelqu'un aurait il une idée? (elle serait la très bien venue)

Merci

Nico76131

[Nico76131]

J'ai trouvé un programme en "langue française" qui pourrait je pense convenir à ce que je recherche pour mes données, seul petit bémol je ne sais pas comment l'écrire dans Matlab.

Je voudrais écrire :

V est mon vecteur original et V1 est mon nouveau vecteur.

pour tous les éléments de mon vecteur (

Code :

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

for i=1:length(V)

)

si V(i) est différente de V(i+1) et V(i+1)différente de V(i+2)
alors V(i+1)=V(i).
(

Code :

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

if V(i,1)~=V(i+1,1) && V(i+1,1)~=V(i+2,1);    V1(i,1)=V(i,1); V1(i+1,1)=V(i,1);

)

si V(i) est différente de V(i+1) et V(i+1)égale à V(i+2)
alors il faut regarder si V(i+2)=V(i+3).
si V(i+3)est différent de V(i+2)
alors V1(i+1)=V1(i+2)=V1(i+3)=V(i).

par contre si V(i+2)=V(i+3)
alors il faut regarder si V(i+3)=V(i+4)
si V(i+3)est différent de V(i+4)
alors V1(i+1)=V1(i+2)=V1(i+3)=V1(i+4)=V(i).

par contre si V(i+3)=V(i+4)
alors il faut regarder si V(i+4)=V(i+5)
si V(i+4)est différent de V(i+5)
alors V1(i+1)=V1(i+2)=V1(i+3)=V1(i+4)=V(i)

par contre si V(i+2)=V(i+3)
alors il faut regarder si V(i+3)=V(i+4)
si V(i+3)est différent de V(i+4)
alors V1(i+1)=V1(i+2)=V1(i+3)=V1(i+4)=V1(i+5)=V(i)
etc..... jusqu'à V(i+10)


ensuite il recommencera la boucle à partir de V(i+10) qu'il considèrera comme V(i).



J'ai regardé sur le forum et je ne sais pas du tout quelles types de boucles je pourrais combiner pour créer ce programme (select case, if else, do)???

Comment insérer plein de conditions dépendantes dans une boucle?
Est il possible de faire des if de if??

Merci pour vos propositions.

[membreComplexe12]

peux etre avec une fonction "filtre" qu iest propre à matlab (je ne me rappel plsu du mot clef, je crois que c'est tout simplement: filtre)

sinon tu fait une boucle sur toutes tes valeurs et tu supprimer les valeurs qui dépassent un certain seuil...

[Jerome Briot]

Voici une autre idée :

Tout d'abord, tu comptes le nombre de valeurs consécutives pour chaque palier :

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
load Exp_03_VGeschwindigkeitneu.mat
 
x = U1(:,2);
y = U1(:,1);
 
k = 1;
bin(k) = y(1);
num(k) = 1;
 
for n = 2:numel(y)
    if y(n)==y(n-1)
        num(k) = num(k)+1;
    else
        k = k+1;
        bin(k) = y(n);
        num(k) = 1;
    end
end

Ensuite, tu définis un seuil et tu modifies la valeur des paliers dont le nombre de points est inférieur au seuil en mettant cette valeur à la valeur du palier précédent (par exemple)

Code :

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

1
2
3
4
5
6
7
8
9
10
seuil = 4;
 
if num(1)<=seuil
    bin(1) = 0;
end
for n = 2:numel(bin)
    if num(n)<=seuil
       bin(n) = bin(n-1);
    end
end

Il ne te reste plus ensuite qu'à construire le nouveau signal :

Code :

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

1
2
3
4
5
6
newy = zeros(size(y,1),size(y,2));
k = 1;
for n = 1:numel(bin)
    newy(k:k+num(n)-1) = bin(n);
    k = k+num(n);
end

ou encore :

Code :

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

1
2
3
4
5
newy = [];
 
for n = 1:numel(bin)
    newy = [newy ; repmat(bin(n),num(n),1)];
end

On peut encore pousser un peu plus l'optimisation du code mais je ne pense pas que cela soit réellement nécessaire...

[membreComplexe12]

ouahou !
pas mal du tout cela!

[Nico76131]

Programme tout à fait adapté et très performant.
Merci