Discussion :

[membreComplexe12]

Salut tous,

je cherche à écrit un bout de code mais de manière optimisé sur matlab car j'aurais à utiliser ce bout de code très souvant.

Contexte :

1°) j'ai deux vecteurs qui sont stocké dans mon code : A1 et A2, il ont la même taille car A1 représente des abscisses et A2 les ordonnées correspondantes.

2°) j'ai un scalaire (par exemple = 37) qui est forcement compris entre les valeurs qui se trouvent dans le vecteur A2 et je veux connaitre le A1 correspondant.

Exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
A1=[1 2 3 4 5];
A2=[10 20 30 40 50];
 
for i=1:1:5
scalaire(i)=37;
%je fais donc une interpolation linéaire entre 30 et 40 pour connaitre le A1 %correspondant à la valeur de 3.5 et prendre le plus près des données %donnues
index=??
if (abs(A1(index))-A1reel)
...
end
end

[Gakusei]

Salut tous,


2°) j'ai un scalaire (par exemple = 3.7) qui est forcement compris entre les valeurs qui se trouvent dans le vecteur A2


A2=[10 20 30 40 50];

Tu veux dire A1?

Si c'est le cas, il me faudrai savoir est-ce que A1 est forcément composé d'entier?
Si oui:

Code :

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

1
2
3
4
5
 
scalaire = 3.7
int_scalaire =  floor(scalaire)
low_index = find(A1 == int_scalaire ) #indice inferieur 
interp = (scalaire-int_scalaire)*(A2(low_index+1)-A2(low_index)) + A2(low_index) #interpolation lineaire

Je pense que la fonction find est suffisamment optimisé.

[Merel]

Salut 21did21,

Ca n'est pas très clair ce que tu exposes...

Tu parles d'un scalaire (3.7) qui est forcement compris entre les valeurs qui se trouvent dans le vecteur A2. Or sur ton exemple ton vecteur A2 va de 10 à 50....Ca ne serait plutôt un scalaire de A1 et tu veux connaitre ton A2 correspondant?

De manière plus générale, je ne vois pas pourquoi tu veux faire une interpolation.
D'après ce que tu dis, tu chercherais juste le rapport qui existe entre tes deux axes.

Code :

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

3.7*A2(1)/A1(1)

suffit à connaitre ton A2 hypothétique.

Ou alors tu cherches le A2 correspondant à l'ordonnée d'une courbe qui a pour abcisse A1 ?

[membreComplexe12]

MERCI DE VOTRE AIDE

je suis vraiment désolé, je viens d'éditer mon message :

=> ce n'est pas 3.7 mais 37 que je voulais marquer.


comme je n'ai pas été clair je viens de faire un petit schema.

=> je veux calculer pour tout point de la courbe noir la distance Y2-A2 et savoir combien elle vaut.

=> je connais les ordonnées A2 et les abscisses A1 de la courbe noire et je connais les abscisses et les ordonnées de la courbe verte Y1 et Y2.

=> ce que je pensais faire c'est faire une boucle pour toutes les ordonnées de la courbe noir et voir a quelle ordonnées ça correspond sur la courbe verte et ensuite faire la différence des deux abscisses associées.

ce qui pose probleme est que les composantes que j'ai à disposition pour la courbe verte ne tombe pas forcement sur les ordonnées de la courbe noire.
-> d'où l'interpolation linéaire...

[tanguy4724]

regarde du côté de la fonction spline . Tu y trouveras peut être ton bonheur

[membreComplexe12]

regarde du côté de la fonction spline . Tu y trouveras peut être ton bonheur

en fait je connais la fonction splines mais comment l'utiliser ici ? je ré interpole la courbe verte à chaque itération de ma boucle pour avoir une ordonnée qui tombe "pile" sur l'ordonnées que je souhaite ?

=> ça fait un peu lourd non ?

[Jean Dumoncel]

Bonjour,

Qu'entends-tu par un peu lourd?
Il suffit d'interpoler pour les valeurs qu'il te manque, de toute manière, je ne vois pas trop quel autre choix tu pourrais avoir...

[membreComplexe12]

merci, j'y suis arivé