Discussion :

[conanran]

Bonjour,
En effet, j'ai un petit code qui consiste à:

données d'entrée: un vecteur r et un vecteur T, en effet, pour chaque élément r(i) on associe un élément T(i)
on a une matrice A,
L'objectif consiste à chercher la valeur proche de r pour chaque élément de A, A(i,j), et à calculer la valeur correspond de vecteur T.

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
20
clear all; clc;
 
A=[20 4 5 5 7
   2 5 0 6 8
   1 15 5 4 10];
 
r=[1 2 3 4 5 6 7 8];
T=[5 4 8 7 1 2 4 5];
%%
i=1;
while i<size(A,1)+1
        j=1;
        while j<size(A,2)+1
            X=A(i,j);
             [pasbesoin,idx]=min(abs(r(1,:)-X));            
             B(i,j)=T(1,idx); 
             j=j+1;
        end
        i=i+1;
end

Avez vous une idée pour éviter les deux boucles While.
Merciiiii

[Jerome Briot]

Comme ceci :

Code :

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

1
2
3
4
5
6
7
8
A = [20 4 5 5 7
    2 5 0 6 8
    1 15 5 4 10];
 
r = [1 2 3 4 5 6 7 8];
T = [5 4 8 7 1 2 4 5];
 
B = interp1(r, T, A, 'nearest', 'extrap');

[conanran]

Merci
Exactement ça

[conanran]

Bonjour,

En effet, je veux modifier un peu mon code, maintenant, j'ai une matrice A en trois dimensions A(:,:,m) et T une matrice de T(m,..).
De la même façon, l'objectif consiste à chercher la valeur proche de r pour chaque élément de A(:,:,1,m), et à calculer la valeur correspond de vecteur T(1:m,..) .
J'ai essayé avec le code suivant mais il n'est pas marché.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
A=zeros(3,5,2)
A(:,:,1) = [20 4 5 5 7
    2 5 0 6 8
    1 15 5 4 10];
A(:,:,2) = [20 4 3 5 3
    2 5 0 7 8
    1 15 5 4 5];
 
r = [1 2 3 4 5 6 7 8];
T = [5 4 8 7 1 2 4 5
       4 5 7 8 4 6 1 4];
 
B(:,:,1:m)= interp1(r, T(1:m,:), A(:,:,1:m), 'nearest', 'extrap');

il donne comme erreur:

Error using griddedInterpolant
The grid vectors do not define a grid of points that match the given values.

Error in interp1 (line 158)
F = griddedInterpolant(Xext,V,method);

Merci pour votre aide