Discussion :

[bakaratoun]

Bonjour,
J'ai une variable booleenne, vect définie comme suit:

Code :

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

1
2
Time=[0 1 2 5];
Vect=[0 1 0 1];

Connaissez vous une fonction f sous Matlab qui permet d'interpoler des valeurs booleenne (interpolation constante)

je précise: je cherche f telle que, quand je fais:

Code :

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

1
2
 
f(Time,Vect,1.5)

me donne: 1.

la seule fonction d'interpolation que je connais est interp1, mais elle fait une interpolation linéaire, . moi je voudrais bien une interpolation constante.

Merci

[duf42]

Bonjour,

Et si tu utilises l'option 'nearest' dans INTERP1, ca ne convient pas à ton usage?

Bonne apm,
Duf

[bakaratoun]

bonjour,

le problème avec nearest c'est que, quand je fais:

Code :

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

YI = interp1(Time,Vect,1.7,'nearest')

ca me donne 0 alors que c'est sensé me donner 1.

en effet, Time=[0 1 2 5] et Vect=[0 1 0 1].. donc la valeur entre 1 et 2 doit etre une constante égale à 1.

voilà,

[duf42]

Ah ok, j'avais compris un peu de travers

Je sais pas si c'est la réponse à ton problème mais tu peux peut-être jeter un coup d'oeil à la fonction STAIRS

[bakaratoun]

d'accord.

Le seul problème c'est que stairs permet de tracer la courbe et ne fais pas le calcul.

dommage que:

Code :

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

a=stairs(Time,Vect,1.5)

ne marche pas.

[duf42]

Normalement tu peux récupérer les arguments de sortie de STAIRS (sans tracer de courbe), après je pense que c'est possible de l'adapter à ton besoin:

Code :

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

1
2
 
[x,y]=stairs(Time,Vect);

Après tu peux peut-être utiliser INTERP1 plus facilement...

[bakaratoun]

Ok, ca devrait marcher.
merci,

[bakaratoun]

bonjour,
en fait, je viens de m'en rendre compte, la focntion stairs renvoie un´vecteur temps qu'on ne peut pas interpoler vu qu'il contient deux fois le meme instant. . Donc, c'est un peu long à mettre en place cette interpolation "constante":

1)d'abord on créer les bon vecteurs avec stairs.
2)on corrige les instants redondants.
3)on interpole avec interp1.

Connaissez vous pas une fonction qui fait ca directement ?

emrci

[bakaratoun]

wé, em fait, c'est pas tr`s compliqué:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
function yi= interpConstant(time,value,xi)
 
[time_,value_]=stairs(time,value);
 
   for k=1:length(time_)-1
       if time_(k)==time_(k+1)
           time_(k)=time_(k)-(1e-10);
       end
   end
 
   yi=interp1(time_,value_,xi);
 
end