Discussion :

[helaaa]

bonjour,

j'ai un vecteur de valeurs contenant des NaN

X=[NaN NaN 5 4 NaN 9 NaN 19 NaN 9 NaN 8 9 5 4 NaN 5 6 3 8 7 5 NaN NaN 8 NaN 4 5 NaN NaN 3 8 NaN NaN NaN]

j'ai essayé d'utiliser interp1 mais ça ne marche pas, je pense que le problème est que mon vecteur peut commencer et finir par des NaN

comment je peut procéder dans ce cas?

[mihaispr]

Salut!

A mon avis je dois qu'il s'agit des valeurs NaN.

Etape1: Essaye de remplacer les valeurs NaN par des valeurs =0 comme ca:

Code :

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

1
2
 
x(find(isnan(x)==1))=0 % remplacer les valeurs NaN avec 0

Etape 2: Utiliser puis la fonction interp 1

Code :

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

1
2
 
>>doc interp1

J'espere que ca t'aides!

Amicalment,

Michel

[helaaa]

merci pour la proposition mais le problème c'est juste qu'aux bornes il n'interpole pas, j'ai fais ceci, ça marche bien au milieu mais pas aux bornes

Code :

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

1
2
3
X=[NaN NaN 5 4 NaN 9 NaN 19 NaN 9 NaN 8 9 5 4 NaN 5 6 3 8 7 5 NaN NaN 8 NaN 4 5 NaN NaN 3 8 NaN NaN NaN]
idx=find(~isnan(X))
Yi=interp1(idx,X(idx),1:length(X))

[Jean Dumoncel]

Salut,

par défaut, interp1 effectue une interpolation linéaire, donc il ne peut interpoler que entre 2 points qui sont définis, les valeurs de début et de fin ne sont pas encadrer, donc interp1 ne peut pas les calculer.

edit : dans la doc de interp1 :

For the 'nearest', 'linear', and 'v5cubic' methods, interp1(x,Y,xi,method) returns NaN for any element of xi that is outside the interval spanned by x. For all other methods, interp1 performs extrapolation for out of range values.

yi = interp1(x,Y,xi,method,'extrap') uses the specified method to perform extrapolation for out of range values.

Comme tu veux extrapoler, il faut utiliser l'argument 'extrap' dans interp1.

[helaaa]

j'ai fais la même chose en ajoutant

method,'extrap'

mais j'ai ce message d'erreur

??? Undefined function or variable 'method'.

[Jerome Briot]

yi = interp1(x,Y,xi,method) interpolates using alternative methods:

• 'nearest' Nearest neighbor interpolation
• 'linear' Linear interpolation (default)
• 'spline' Cubic spline interpolation
• 'pchip' Piecewise cubic Hermite interpolation
• 'cubic' (Same as 'pchip')
• 'v5cubic' Cubic interpolation used in MATLAB® 5. This method does not extrapolate. Also, if x is not equally spaced, 'spline' is used/

[helaaa]

et ceci ne marche pas?? je l'ai trouvé dans le help de interp1

yi = interp1(x,Y,xi,method,'extrap')

[Caro-Line]

Il faut que tu initialises la variable method avec une des valeurs données dans la documentation.

[helaaa]

oui je suis bête
j'avais la tête ailleurs j'essayais de chercher la méthode d'interpolation linéaire utilisée par matlab, j'essaye de lire la fonction interp1 mais j'arrive pas à comprendre quelle formule on utilise
d'après ce lienhttp://fr.wikipedia.org/wiki/Interpo..._lin%C3%A9aire on a trois formules possibles mais c'est très compliqué de déterminer laquelle matlab utilise pour l'interpolation 'linear'

[Jerome Briot]

d'après ce lienhttp://fr.wikipedia.org/wiki/Interpo..._lin%C3%A9aire on a trois formules possibles mais c'est très compliqué de déterminer laquelle matlab utilise pour l'interpolation 'linear'

trois formulations équivalentes

Prends une feuille de papier, un crayon et développe (ou factorise) les équations

Sinon, pour répondre peu être plus précisément à ta question, en regardant le code de INTERP1, il semble que la formulation utilisée soit la seconde (formule de Taylor-Young au premier ordre)

[helaaa]

oui c'est ce que j'ai fais, c'est Taylor young au premier ordre