matlab find() indices de cellules pour particule, double 4D

Question

J'ai le problème suivant. Considérant que j'ai 4 particules dans une zone.

Les particules traversent cette zone avec une certaine vélocité. À un moment donné, cependant, je veux savoir dans quelle cellule se trouve une telle particule.

Mon exemple de code ressemble à ceci:

x = -10:1:10; % area 
y = -10:1:10; 

particle = zeros(1,2,1,4); % 4 Particle 
particle(1,1,1,1) = 9.5; % x-position of the first particle 
particle(1,2,1,1) = 9.5; % y-position of the first particle 
particle(1,1,1,2) = 7.25; 
particle(1,2,1,2) = 8.5; 
particle(1,1,1,3) = 6.25; 
particle(1,2,1,3) = 7.5; 
particle(1,1,1,4) = 5.25; 
particle(1,2,1,4) = 6.5; 

for xp=1:size(particle,3) 
    for yp=1:size(particle,4) 
     if particle(1,1,xp,yp) < x(1) 
      cell_x = 1; 
     elseif particle(1,1,xp,yp) > x(end) 
      cell_x = size(x,2); 
     else 
      cell_x = find(x <= particle(1,1,xp,yp),1,'last'); 
     end 
     if particle(1,2,xp,yp) < y(1) 
      cell_y = 1; 
     elseif particle(1,2,xp,yp) > y(end) 
      cell_y = size(y,2); 
     else 
      cell_y = find(y <= particle(1,2,xp,yp),1,'last'); 
     end 
    end 
end 

Mon code prend trop de temps pour un certain nombre de 52 x 10 particules. Je voulais omettre une boucle for en entrant yp comme un tableau dans particle(1,2,xp,yp), mais malheureusement la fonction find() ne peut pas gérer les tableaux. Ou ai-je tort? Avez-vous une idée de comment le faire fonctionner beaucoup plus vite pour que j'ai

cell_x = [20,18,17,16] 
cell_y = [20,19,18,17] 

Answer 1

Un exemple de code à mon commentaire

x = -10:1:10; % area x 

%input data 
Numberofdata=4; 
saved_x=[9.5 7.25 6.25 5.25]; 

%conversion from value to index 
idx_x1=floor(saved_x)+abs(min(x))+1; 
%search for index 
for k=1:Numberofdata 
idx_x2(k,1)=find(x <=saved_x(k),1,'last'); 
end 

% Output 
[idx_x1' idx_x2] 

20 20 
18 18 
17 17 
16 16 

Je sais qu'il est appliqué à un vecteur et non une matrice cependant, vous pouvez appliquer cette transformation à n'importe quelle dimension de votre matrice.

Dans le cas où vous ne changez pas x par pas de 1 mais quelque chose comme ça:

x=min_x:dx:max_x 

alors vous devez faire ce qui suit:

idx_x1=floor(saved_x/dx)+abs(min_x/dx)+1;