Discussion :

[ARHAB]

Bonsoir,

J'ai un petit soucis lorsque j'essaye de faire une manipulation sous matlab, voici mon problème :

J'ai deux vecteurs x et z de même dimension, ainsi qu'un troisième vecteur y de même dimension et que je défini comme fonction des deux premiers. C'est à dire qu'au couple (xi,zi) est associé la valeur yi. Maintenant je construis les vecteurs X et Z de plus grande dimension que les précédents, obtenus en répétant plus d'une fois les éléments des vecteurs x et z, exemple : si x=(x1 x2 x3) alors X=( x1 x1 x2 x1 x3 x2 x3)
z=(z2 z1 z3) alors Z=(z2 z2 z1 z2 z3 z1 z3)

dans ce cas y=(y12 y21 y33) alors je cherche à avoir

Y=(y12 y12 y21 y12 y33 y21 y33)

Moi j'ai les vecteur x, z, y, X et Z et je cherche à obtenir le vecteur Y par une fonction matlab, sans utiliser la commande for. Merci

[Invité]

Bonjour,

un troisième vecteur y de même dimension et que je défini comme fonction des deux premiers.

Quelle est l'expression de cette fonction ?

Quel est le problème ? Tu obtiens des erreurs ? Lesquelles ?

[ARHAB]

Salut,

Voici la partie de mon code que je veux accélérer :

Code :

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

1
2
3
4
5
6
7
Gr=zeros(N,N) ;
for s=1:M*M
    TATA=abs(mnvx-XN(s))+abs(mnvz-ZM(s)) ;
    [r,c]=find(TATA==0) ;
    idx=sub2ind(size(Gr),r,c) ;
    Gr(idx)=y(s) ;
end

Ici: XN ZM et y sont de dimensions (M x M). Alors que mnvx, mnvz et Gr sont de dimensions (N x N). Avec M < N. Est-il possible de le réécrire en remplaçant la boucle for par une fonction MATLAB ? Merci

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
Gr = zeros(N,N) ;
 
for s = 1:M*M
    idx = (abs(mnvx-XN(s))+abs(mnvz-ZM(s)))==0; 
    Gr(idx) = y(s) ;
end

Voir la FAQ : Qu'est-ce que l'indexage logique ?

[ARHAB]

ça marche également mais c'est juste un peu moins long! Pour des valeurs de M importantes, le programme devient excessivement lent. y a t-il un moyen de se passer de la boucle for ?

Mon idée que je n'arrive pas à concrétiser est de faire ceci :

Code :

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

1
2
3
4
5
VXN=reshape(XN,M^2,1);
VZM=reshape(ZM,M^2,1);
Vy=reshape(y,M^2,1) ;
Vmnvx=reshape(mnvx,N^2,1);
Vmnvz=reshape(mnvz,N^2,1);

Ensuite à partir de VXN,VZM,Vy et Vmnvx, Vmnvz déduire par une fonction MATLAB le vecteur VGr de taille (N^2,1) et faire ensuite :

Code :

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

Gr=reshape(VGr,N,N) ;

Est ce que vous voyez ce que je veux dire ? En tous cas merci pour l'aide fournie jusque là.

[Jerome Briot]

Tu peux essayer avec bsxfun (si tu as la quantité de mémoire suffisante) :

Code :

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

1
2
3
4
A1 = bsxfun(@minus,repmat(mnvx(:),1,M*M),XN(:).');
A2 = bsxfun(@minus,repmat(mnvz(:),1,M*M),ZM(:).');
 
idx = (abs(A1)+abs(A2))==0;

idx est un tableau de dimension (N*N)x(M*M) où chaque colonne correspond à une itération de ton exemple.

Par contre, pour la suite, je ne suis pas sûr de bien comprendre ce que fait ton exemple. Certaines valeurs de Gr prennent différentes valeurs de y au cours des itérations