Discussion :

[jean_pierr]

Bonjour,

j ai un vecteur ligne X de taille 1*4000. qui contient 4000 angles.
j ai besoin de construire une matrice 3D qui contient 4000 matrices de rotations(pour chaque angle on a une matrice) sans utiliser les boucles car ils sont lourdes dans ce cas, merci de m'aider si vous avez une idee.

[soft001]

Quelle est la matrice associée à chaque angle ???

[jean_pierr]

merci a votre reponse.

j'ai 4000 angles enregistrer dans un vecteur x de taille 1*4000.

pour chaque angle je vais calculer une matrice de rotation R=[cos t -sin t; sin t cos t].

supposons que t=X(1)=90
R=[0 -1;1 0].
j 'ai besoin de calculer 4000 matrices comme ci dessus en un seul coup sans utiliser les boucles et on peut utiliser une matrice 3 D W(2,2,4000) pour enregistrer les resultat. merci de m'aider.

[FR119492]

Salut!

une matrice 3D

Une matrice est nécessairement à deux dimensions. Dans ton cas, il s'agit peut-être plutôt d'un tableau.
Jean-Marc Blanc

[jean_pierr]

par abus de language on dit matrice,

je sais bien qu'on dit un tensor pour une superposition des plusieurs matrices 2D.

comme dans mon cas je vais superposer 4000 matrices de taille 2*2.
merci

[soft001]

Tu regardes ce code :

Code :

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

1
2
3
t=[0 pi/4 pi]
R=[cos(t);-sin(t);sin(t);cos(t)]
R=reshape(R,2,2,3)

dans ton cas tu remplaces le 3 par numel(t) (égale 4000)

[jean_pierr]

Grand merci,

tu as une petite faute : R=[cos(t) ; sin(t) ; -sin(t) ; cos(t)]. grand merci c'est une bonne solution. Je crois que n'importe quelle boucle peut être changée par un calcul compact. Grand merci

[soft001]

Grand merci,
tu as une petite faute : R=[cos(t) ; sin(t) ; -sin(t) ; cos(t)]

Oui j'ai pas fait attention

Je crois que n'importe quelle boucle peut être changée par un calcul compact. Grand merci

Je ne pense pas, si non tu aurais dû pas trouver la syntaxe

Code :

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

1
2
for 
end

[shaiHulud]

Je crois que n'importe quelle boucle peut être changée par un calcul compact.

Oui, ca se démontre. Cela dit attention:

- C'est de moins en moins lisible pour les calculs compliqués (dès que le corps d'une itération dépend des itérations précédentes)

- Si par exemple l'itération dépend de toutes les itérations passées, la taille mémoire de l'opération vectorielle va exploser (j'ai souvent du transformer un code vectorisé en boucle pour éviter les 'out of memory')

- Les boucles ont été optimisées dans les dernières versions de matlab et rivalisent avec les cellfun/arrayfun

- les boucles peuvent être dispatchées sur plusieurs processeurs via parfor, ce qui démultiplie les performances . Donc dans le cas de boucles dont les itérations sont indépendantes, l'argument performance n'est plus du tout convaincant.

[jean_pierr]

j'ai pas compris parfaitement tes remarques!