Discussion :

[sabrina-e]

Bonjour tout le monde , j'aurais besoin d'aide,
je cherche a trouver la valeur d'une constante "epsilon" sur un ensemble de points d'une grille (xoy), je m'explique:
J'ai une grille en x,y à la quelle j'affecte pour chaque cellule une constante " epsilon" ; et je souhaite récupérer la valeur de cette constante en certains point de ma grille, pour ce la j'ai choisis de faire une interpolation comme suit:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
y=1:1:3;
x=1:1:3;
[X,Y]=meshgrid(x,y) % définition de ma grille;
epsilon=[1 3 2;3 5 6;11 23 26] % affectation des valeurs d'epsilon a chaque cellule de ma grille;
% Interpolation 
xo =1; yo =1;   % (x,y) origine 
xd = 3; yd = 3;   % (x,y) fin
nbpoints = 4;    % nombre 
xprofil = linspace(xo, xd, nbpoints)
yprofil = linspace(yo, yd, nbpoints)
epsilon_interpol = interp2(X,Y,epsilon, xprofil, yprofil)

le résultat de code me renvoit:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
xprofil =
 
    1.0000    1.6667    2.3333    3.0000
 
yprofil =
 
    1.0000    1.6667    2.3333    3.0000
 
epsilon_interpole =
 
    1.0000    3.6667   11.5556   26.0000

Or que sur ma grille epsilon est une constante et non pas une fonction continu , donc c'est pas ça ce que je souhaite avoir , plutôt je souhaite avoir ça :

Code :

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

1
2
3
epsilon_interpole =
 
    1 5 26 26

Merci !

[Jerome Briot]

Tu peux essayer comme ceci :

Code :

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

epsilon_interpol = interp2(X, Y, epsilon, xprofil, yprofil, 'nearest')

[sabrina-e]

Bonjour,
Merci beaucoup ça marche parfaitement !