Discussion :

[nilsmatlab]

Bonjour j'ai ce message d'erreur lorsque je tente de faire une interpolation.

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

Code :

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

1
2
3
4
for i=1:z 
 d2(i)=interp1(time,d1,time2(i));
 P2(i)=interp1(d1,P,d2(i));
end

time, d1 et P sont des vecteurs issus de mes données, comportant 7580 valeurs.
time2 comporte les mêmes valeurs que time mais avec un pas plus petit pour plus de précision.
d2 est le vecteur d1 interpolé pour qu'il fasse la même longueur que time 2 (dimension > 40 000 valeurs).
Lorsque j'essaye de déterminer P2 par interpolation, j'obtiens le message d'erreur ci dessus.

Je pense que c'est dû au vecteur d1. En effet celui ci comporte plusieurs fois la même valeur à des pas de temps différents.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
294
294
294
295
295
295
295
295
295
296
296

par exemple.

Est ce que l'erreur vient d'ici (à noter que les valeurs de d1 sont toujours croissantes).
Je ne peux pas faire

Code :

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

unique

, car dans ce cas d1 et P n'aurait pas la même longueur. Je ne peux pas non plus recalculer P avec une boucle

Code :

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

for i=1:length(d1bis)

car il me manquerait des valeurs.

Comment pourrais-je faire pour réussir mon interpolation ?

Je vous remercie

[shaiHulud]

Bonjour,

Tout d'abord une première remarque, il me semble que tu peux éviter la boucle sur i en faisant directement

Code :

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

1
2
3
 
d2=interp1(time,d1,time2);
P2=interp1(d1,P,d2);

Ensuite le fond du problème, c'est effectivement que d1 ne prend pas de valeurs uniques et qu'il ne sait donc pas comment interpoler.

Cela dit vu la construction que tu fais, tu devrais pouvoir contourner le problème:
- Assure toi que pour deux valeurs identiques de d1, P prends la même valeur. (Si ce n'est pas le cas, ton problème semble mal posé et il faut que tu nous expliques ce que représente d1 et P et comment se fait il que tu ne connaisses pas les dates auxquelles correspondent ces valeur)
- Filtre par valeurs uniques avant la 2nde interpolation

Code :

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

1
2
3
4
 
d2=interp1(time,d1,time2);
[b,m,n] = unique(d1);
P2=interp1(d1(m),P(m),d2);

Une seconde solution consiste à rendre les dates distinctes (arbitrairement, selon leur ordre dans les matrices, voici un exemple):

Code :

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

1
2
3
4
5
6
7
8
 
d2=interp1(time,d1,time2);
dt= diff(d1);
dt_step= min(dt(dt>0));
N= length(d1);
tie_order= (dt_step/N):(dt_step/N):dt_step;
d1_untied= d1+ order;
P2=interp1(d1_untied,P(m),d2);

[nilsmatlab]

Merci VV33D, j'ai essayé la 1ere méthode, j'obtiens des résultats mais il ne sont pas cohérents.
Pour les valeurs issus des données j'obtiens. (d1 et d2 représentent des distances, et P et P2 représentent l'effort subit par mon véhicule à chaque distance). P est différent pour une même distance dans mes données car il dépend d'autres paramètres qui eux évolue pour une même distance (c'est juste que la précision de la distance n'est pas assez précise, les valeurs sont arrondies à l'unité)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
d1	P
0	1219.544105
0	942.9848732
1	597.9365095
1	463.5609502
1	439.0669896
1	441.6812616
2	443.8598216
2	236.6611456
2	102.7405776
2	-33.79426237

Quand pour ma simulation j'obtiens :

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
19
d2	P2
0	942.9848732
0	942.9848732
0	942.9848732
0.2	842.7241509
0.6	642.2027063
1	441.6812616
1	441.6812616
1	441.6812616
1	441.6812616
1	441.6812616
1	441.6812616
1	441.6812616
1	441.6812616
1.2	346.5861568
1.6	156.3959472
2	-33.79426237
2	-33.79426237
2	-33.79426237

Il associe tout le temps la dernière valeur de chaque distance à P2. De plus mon interpolation de d2 semble fausse également

[Jean Dumoncel]

Bonjour,

P est différent pour une même distance dans mes données car il dépend d'autres paramètres qui eux évolue pour une même distance (c'est juste que la précision de la distance n'est pas assez précise, les valeurs sont arrondies à l'unité)

et tu n'as pas accès aux valeurs non arrondies? Car tout tes problèmes viennent de là...

Par exemple, tu dis que l'interpolation de d2 est fausse :

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
19
d2
0	
0
0	
0.2
0.6
1	
1	
1	
1	
1
1	
1	
1	
1.2	
1.6	
2	
2
2

Quelles valeurs aurais-tu souhaité obtenir?

[shaiHulud]

Je ne comprends pas pourquoi tu ne fais pas juste

Code :

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

1
2
 
P2=interp1(d1,P,time2)

Vu que tu as l'air de dire que les valeurs P correspondent aux "temps" d1

[nilsmatlab]

Bonjour,

et tu n'as pas accès aux valeurs non arrondies? Car tout tes problèmes viennent de là...

Quelles valeurs aurais-tu souhaité obtenir?

Non, je récupère les valeurs comme ça et je ne peux pas les modifier.
Avoir des valeurs qui augmentent progressivement. Du coup j'ai recréer le vecteur d2 avec un pas que j'ai fixé

Code :

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

d2 = [d1(1):0.15:d1(end)];

, mais ça fait vraiment "bidouille", parce que le 0.15 ne provient de nulle part

Je ne comprends pas pourquoi tu ne fais pas juste
Code :Sélectionner tout - Visualiser dans une fenêtre à part
1
2
P2=interp1(d1,P,time2)

C'est ce que j'ai fait en tout premier, et c'est là que Matlab m'a sorti l'erreur qui est dû au fait que d1 comporte plusieurs valeurs identiques...

[Jean Dumoncel]

Du coup j'ai recréer le vecteur d2 avec un pas que j'ai fixé

Code :

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

d2 = [d1(1):0.15:d1(end)];

, mais ça fait vraiment "bidouille", parce que le 0.15 ne provient de nulle part

En même temps tu fais une interpolation (approximation) sur des données très approximatives... je vois pas bien ce que tu peux faire avec ces données

[nilsmatlab]

Mais il n'existe pas de solutions pour faire une interpolation lorsqu'il y a plusieurs x identiques ?

[Jean Dumoncel]

Un exemple indépendant de tes données : voici ce que interp1 "voit" :

des points qui ont la même abscisse mais des ordonnées différentes... comment veux-tu qu'il interpole quoi que ce soit? Quel règle il dois appliquer? Comment ferais-tu toi à la main pour interpoler? Il n'existe pas de méthodes générales dans ce cas, donc avec les infos que tu nous as donné, non l'interpolation n'est pas possible.

[nilsmatlab]

Je vous remercie, je vais essayer de contourner le problème, car je ne dispose pas de plus de données également.