Discussion :

[membreComplexe12]

salut tous,

j'ai une question surement simple pour vous mais pas pour moi...

données entrées
j'ai deux vecteurs qui contiennent des données expérimentales X (abscisses) et Y (coordonnées), du genre:

Code :

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

1
2
Xexp=1 2 3 4 5 6 7 8 9
Yexp=5 8 6 4 7 1 5 2 8

calcul
je fais un calcul avec un modele et voici les résultats que j'obtiens:

Code :

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

1
2
X=1 2.2 2.9 5.2 6.3 9
Y=5 8.2 6.3 7.2 1.3 8

resultats
ce que je voudrais comme résultats c'est un vecteur erreur qui est:

Code :

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

erreur(i)=Y(i)-Yexp(i)

probleme
le probleme est que pour faire ceci il faut deux choses:

1°) que les vecteurs Yexp et Y ont la meme dimension
2°) que les Y(i) et Yexp(i) soient pris au meme abscisses

le soucis et que je ne sais pas trop comment transformer mon vecteur Y pour que les calculs correspondent aux données expérimentales...

je pense que matlab doit avoir une fonction préprogrammée pour cela mais je ne la connais pas...

j'espere que vous pourrez m'aider

[Jean Dumoncel]

Bonjour,

pour l'interpolation, il existe la fonction interp1, mais tu ne peux pas directement calculer ton modèle pour les valeurs de Xexp? Cela aurait plus de sens...

[membreComplexe12]

Bonjour,

salut et merci de m'aider

pour l'interpolation, il existe la fonction interp1

ok super !

mais tu ne peux pas directement calculer ton modèle pour les valeurs de Xexp? Cela aurait plus de sens...

au fait c'est le schema numérique qui fait les calculs. (j'utilise ode23) je lui donne un temps initial un temps final et il calcul les increments lui meme et il n'y a pas de raison pour que ça tombe pile avec ce que j'ai comme donnée experimentale..


EDIT:
c'est bon j'ai trouvé comment faire. Merrci en tout cas pour ton aide

[membreComplexe12]

en fait je bloque... j'aurais besoin de votre aide svp:

Données de départ
experimental

Code :

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

1
2
3
4
5
6
 
Xexp   Yexp
1           10
2           20
3           30
4           40

calculé

Code :

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

1
2
3
4
5
 
Xcalcul   Ycalcul
1.1           15
2.1           25
3.1           35

ma premiere idée
mettre les abscisses expérimental dans ces calculée et vis versa

Code :

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

1
2
3
donneesExp1=zeros(length(Xexp{1})+length(Xcalcul{1}),2)
donneesExp1(:,1)=[Xexp{1};Xcalcul{1}];
donneesExp1(:,2)=[Yexp{1};zeros(length(Ycalcul{1}),1)];

je dois avant l'interpolation classer ces données par ordre croissant par rapport à la première colonne, pour ceci je fais:

Code :

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

1
2
 
donneesExp1=sortrows(donneesExp1)

par contre comment dois je faire pour faire une interpolation que les données rajoutées, comme ceci ? :

Code :

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

1
2
 
donneesExp1=interp1(donneesExp1(:,1),donneesExp1(:,2),Xcalcul{1})

==> en fait je n'ai pas compris comment faire pour interpoler qu'aux endroits où il y a des zeros (c'est à dire où j'ai inséré des abscisses)

ma deuxieme idée
je pense que la premiere idée est pas mal mais le probleem est que les données expérimentales sont trop nombreuses et je préférais plutot conserver les données expérimentales qui correspondent aux pas de temps calculés car ça me fait des vecteurs beaucoup moins grand et largement suffisant.

le probleme est que je ne vois pas trop comment faire...

[Jean Dumoncel]

au fait c'est le schema numérique qui fait les calculs. (j'utilise ode23) je lui donne un temps initial un temps final et il calcul les increments lui meme et il n'y a pas de raison pour que ça tombe pile avec ce que j'ai comme donnée experimentale..

Ah bon? pourtant dans la doc de ode23 :

To obtain solutions at the specific times t0, t1,...,tf (all increasing or all decreasing), use tspan = [t0,t1,...,tf].

Ensuite je pense que c'est une très mauvaise idée de mélanger tes données expérimentales et calculées, calculer l'erreur la-dessus n'aurait pas vraiment de sens.

[membreComplexe12]

Ah bon? pourtant dans la doc de ode23 :

as oui en effet !

Ensuite je pense que c'est une très mauvaise idée de mélanger tes données expérimentales et calculées, calculer l'erreur la-dessus n'aurait pas vraiment de sens.

en fait je voulais juste garder les abscisses experimentales pour pouvoir faire les interpolations pour ces abscisses

[Jean Dumoncel]

en fait je voulais juste garder les abscisses experimentales pour pouvoir faire les interpolations pour ces abscisses

Ok, dans ce cas tu peux faire simplement :

Code :

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

donneesInterpole=interp1(Xcalcul, Ycalcul, Xexp)

[membreComplexe12]

sinon pour l'interpolation que j'ai fais avec interp1 qu'en penses tu c'est bien comme ça qu'il faut faire ?

[Jean Dumoncel]

Oui, mais je persiste en pensant que tu vas dégrader ton modèle en procédant ainsi : tu crée un modèle que j'imagine assez complexe puis tu vas l'interpoler avec une fonction linéaire ou polynomiale. Donc l'erreur que tu calculeras ne refléteras pas exactement l'erreur entre ton modèle et tes valeurs expérimentales.

Pourquoi ne pas directement calculer ton modèle avec les abscisses qui t'intéressent comme suggérer plus haut?

[membreComplexe12]

tu as raison, c'est ce que je vais faire.
merci pour ton aide