Interpolation de vecteurs de tailles différentes

Bonjour à tous,

Je vous explique mon problème : je voudrais avoir des valeurs pour toutes les heures de la journée, or, je n'ai que les valeurs toutes les trois heures. Je voudrais donc interpoler.

Notons le vecteur colonne X qui a 8 valeurs. La valeur d'indice 1 correspond à 00h, l'indice 2 à 3h etc...
Si je veux avoir non plus 8 mais 24 valeurs pour X, comment faire ?

J'ai créé un vecteur colonne Y de 24 valeurs, mais comment interpoler par pas de deux mon vecteur X ? A chaque fois cela ne fonctionne pas parce que X et Y ne sont pas de la même taille, ce que je comprends mais du coup je ne vois vraiment pas comment faire...

Pouvez-vous m'aider s'il vous plaît ?
Merci.

Cleindorie

soit Y ton vecteur de 8 valeurs correspondant à X (8 valeur aussi)
alors

Code:

---

1 2	X24 = 0:23; Y24 = interp1(X,Y,X24);

---

et ouala

Alors j'ai écrit :

Code:

---

1 2 3 4

X=valeurs(1:8); %ce sont mes 8 valeurs Y=Y; % un vecteur colonne qui va de 1 à 8 X24 = 0:23; Y24 = interp1(X,Y,X24);

---

J'ai le message d'erreur suivant :

??? Error using ==> interp1 at 261
The values of X should be distinct.

Error in ==> interpol at 4
Y24 = interp1(X,Y,X24);

Mes valeurs sont pourtant bien distinctes puisque

Code:

---

1 2 3 4 5 6 7 8

X = 8,71 3,67 17,01 83,61 83,61 157,60 99,25 46,69

---

et

Code:

---

Y=[1 2 3 4 5 6 7 8 ]

---

(en colonne)

Citation:

---

Envoyé par cleindorie

Mes valeurs sont pourtant bien distinctes puisque

---

Es-tu bien sûr ?

Citation:

---

Envoyé par cleindorie

Code:

---

1 2 3 4 5 6 7 8

X = 8,71 3,67 17,01 83,61 83,61 157,60 99,25 46,69

---

---

Et je présume que la virgule n'est pas utilisée comme séparateur décimal dans MATLAB.

Bonjour DUT, merci pour ta réponse, mais en ce qui concerne la virgule, c'est juste quand j'ai fait le copier-coller, il y a eu un petit bugg sinon dans matlab j'ai bien des "."
J'ai bien deux valeurs de X identiques, mais ces résultats viennent d'un calcul expérimental donc je ne peux pas les changer. En plus, là je fais un échantillon sur 8 valeurs pour en avoir 24, mais en tout normalement j'ai 1440 valeurs à interpoler, donc il est possible qu'il y ait des doublons...
J'ai essayé d'écrire un algorithme pour interpoler différemment mais là je suis complètement paumée, je n'arrive absolument pas à faire cette interpolation...

Ton problème concerne donc la prise en compte des doublons.

Reviens au phénomène physique que tu mesures. Pourquoi y a-t-il deux valeurs distinctes de Y pour une même valeur de X ?
Ces deux valeurs sont-elles "justes" ? L'une est-elle meilleure que l'autre ? Ou bien serait-il possible de calculer une valeur moyenne ?

Je ne peux pas enlever les doublons, ils sont au contraire importants pour moi. Si j'ai deux valeurs de Y pour une même valeur de X c'est parce que Y représente le temps, donc si X est constant dans le temps, c'est normal. Mes valeurs de X sont déjà des moyennes, je ne peux vraiment pas enlever les doublons. Mais c'est juste pour cela que je ne peux pas faire d'interpolation ? dans ce cas, est-il possible d'utiliser une lookup table ou quelque chose de similaire ?

Si Y représente le temps, il faut inverser X et Y, non ?

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

X = [8.71 3.67 17.01 83.61 83.61 157.60 99.25 46.69];   Y = [1 2 3 4 5 6 7 8 ];   Yi = linspace(1, 8, 24);   Xi = interp1(Y, X, Yi);   figure plot(Y, X, 'b-o'); hold on plot(Yi, Xi, 'r*');

---

Je suis désolée de t'embêter encore un peu DUT, c'est presque ça, mais pas tout à fait.
Je t'explique :

Ce que tu m'as proposé fonctionne mais j'ai quelque peu modifié pour avoir ce que je veux, seulement je ne sais pas si c'est correct.
J'ai mis ton code, ce qui me donnait 24 valeurs de Y allant de 1 à 8. Mais moi je voulais des valeurs allant de 1 à 24, donc si je remplaçais Yi = linspace(1, 8, 24) par Yi = linspace(1, 24, 24) j'avais tout bon pour Yi mais pas pour X forcément, puisqu'il ne pouvait interpoler après la huitième valeur.
J'ai donc remplacé Yi par Yi=[1:1:24]', après avoir lancé le programme pour garder les valeurs trouvées de Xi.
Cela me donne une figure correcte aussi, alors c'est un peu du bidouillage mais ça fonctionne.
Qu'en penses-tu ?

J'ai l'impression que tu t'embrouilles.

Le vecteur de temps Y n'est pas :

Code:

---

Y = [1 2 3 4 5 6 7 8 ];

---

mais, une fois exprimé en heures toutes les 3 heures :

Code:

---

Y = 0:3:21;

---

Et tu souhaites obtenir les valeurs toutes les heures de la journée :

Code:

---

Yi = 0:1:23;

---

Donc :

Code:

---

Xi = interp1(Y, X, Yi, 'linear', 'extrap');

---

Bien entendu, si tu ajoute au vecteur Y la valeur à minuit du jour suivant, l'extrapolation devient inutile.

Merci beaucoup DUT. Tu as raison, je m'étais embrouillée, merci encore pour ton aide.

A bientôt sur le forum.

Cleindorie