Discussion :

[amarion]

Bonjour à tous,
Je reviens sur la discussion que j'avais lancé en début de semaine (10/03) pour remplir une matrice à parti d'un vecteur.
J'ai un vecteur pos de Nd (nombre de diffuseurs) lignes et 2 colonnes (coordonnées x et y). Les diffuseurs ont des coordonnées sur une grille entre [0 Tx] et [0 Ty].
L'idée était de créer une matrice posit contenant comme valeur le nombre de diffuseurs pour chaque noeud.

La solution suivante faisant bien ce qui est souhaité.

Code :

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

1
2
3
4
5
6
7
8
 
pos = round(2*rand(8,2)+1)
[upos,A,B] = unique(pos,'rows');
nB = histc(B,unique(B));
[x,y] = meshgrid(1:max(pos(:,1)),1:max(pos(:,2)));
posit = zeros(size(x)).';
idx = sub2ind(size(posit),upos(:,1),upos(:,2));
posit(idx) = nB

Néanmoins, je souhaite à présent que chaque diffuseur à une position (x,y) donnée ait un poids distribué gaussien avec randn. La matrice de sortie posit devra prendre en compte ces poids.

Exemple :

Code :

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

1
2
3
4
pos = 
[1 1
 3 2
 2 3]

le diffuseur de coordonnées (1,1) aura comme poids 0.3, (3,2) 0.5 et (2,3) 0.7.

Posit ressemblera donc à :

Code :

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

1
2
3
[0.3 0.0 0.0
 0.0 0.0 0.7
 0.0 0.5 0.0]

au lieu de:

Code :

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

1
2
3
[1 0 0
 0 0 1
 0 1 0]

comme avant.

S'il y a plusieurs diffuseurs par noeud, on additionnera les poids.

J'ai bien pensé à faire un find sur chaque coordonnées puis sommer les poids d'indice correspondant mais je ne suis pas sûr que ce soit optimum.

Merci d'avance,
Adrien

[Jerome Briot]

Essaie quelque chose comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
pos = [1 1
       3 2
       2 3]

poids = [0.3
         0.5
         0.7]
   
[upos,A,B] = unique(pos,'rows');
nB = histc(B,unique(B));
[x,y] = meshgrid(1:max(pos(:,1)),1:max(pos(:,2)));
posit = zeros(size(x)).';
idx = sub2ind(size(posit),upos(:,1),upos(:,2));
posit(idx) = poids(A)

[amarion]

On ne peut pas respecter la contrainte d'addition des poids de tous les diffuseurs placés sur le noeud en faisant ainsi.

[amarion]

Quelqu'un a-t-il mieux que ça?

[l,c]=(find(posit~=0));
for p=1:length(l),
posit(l(p),c(p))=sum(Amp(find((posbis(:,1)==c(p))&(posbis(:,2)==l(p)))));
end

la boucle est mauvaise ...

[Jerome Briot]

Comme ceci alors ?

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
21
22
23
24
25
26
27
pos = [1 1
       3 2
       1 1
       1 1
       3 2
       2 3]
 
poids = [0.1
         0.2
         0.1
         0.1
         0.3
         0.7]
 
[upos,A,B] = unique(pos,'rows');
nB = histc(B,unique(B));
[x,y] = meshgrid(1:max(pos(:,1)),1:max(pos(:,2)));
posit = zeros(size(x)).';
idx = sub2ind(size(posit),upos(:,1),upos(:,2));
 
nA = numel(A);
%p=zeros(nA,1);
for n=1:nA
   p(n)=sum(poids(B==n)); 
end
 
posit(idx) = p

[amarion]

Merci.
Je voulais absolument éviter une boucle codée en Matlab.
Tant pis, j'arrive tout de même à mes fins.
Adrien