matrice et rotation

bonjur tous le monde , je commence enfin a me depatouiller avec la 3D , mais je bloque sur les rotation ,
j ai une matrice (4*4) de rotation:
je mutliplis un point 3D :

(x (1 0 0 0
y * 0 cos(a) -sin(a) 0
z 0 sin(a) cos(a) 0
t) 0 0 0 1)

je dois transcrir sa en pascal
j ai cree un type identitter et t3Dmatrice
rotation := identite3D;
rotation[2,2]:=cos(a);
rotation[3,2]:=sin(a);
rotation[2,3]:=sin(a);
rotation[3,2]:=cos(a);

et sa me donne des resultats legerement delirant , quelqu un peux t il m expliquer ou je me suis tromper ?
merci beaucoup
@+ charly :)

il faudrait préciser tes notations : rotation par rapport à quoi et selon quel axe.

Et sinon rotation[2,2] ça veut dire quoi ?

Merci.

excuser moi :
rotation est une matrice , un tableau de reel , en pascal on note sa
matrice : array[0..4,0..4];
pour donner une egaliter ( pour l affecter un nombre a une partie du tableau ) , on peux faire directemment en dessinant le tableau ce qui me donenrais logiquemment :

const ((1,0,0,0),
(0,cos(a),-sin(a),0),
(0,sin(a),cos(a),0),
(0,0,0,1));

mais turbo pascal n accepte pas cette forme la ( a cause du cos et sin ) du coup on est obliger d affecter directemment la variable a une partie du tableau
rotation[2,2]:=cos(a); signifie que cos(a) est affecter a la 2 eme colone 2eme ligne ( ou l inverse ) !

normalement sa me fais une rotation autour des x(axe des absisse) , mais je ne comprends pas pkoi , sa me donne des resultat extravagant !!!
@+charly

Euh, je vais petetre paraitre pour un lourd mais as-tu penser à vérifier si les angles sont en radians ou en degrés ( on les oublis vite ceux là )

ton code en exemple ( pas celui en pascal ) fait un peu brouillon. ( désolé de dire cela ).

Bon en tout cas voici une rotation d'un objet sur une surface plane ( sufit de mettre une 3 eme dimension pour la 3D )

Code:

---

1 2 3 4

P = (x,y);   A = ( cos(a), -sin(a) sin(a), cos(a) )

---

et tu fais A.P = P' en sachant que a est l'angle de rotation, P est une matrice colonne d'ordre 2 et A une matrice carré d'ordre 2.


cela donnerait comme résultat :

Code:

---

1 2	x' = x * cos(a) + x * -sin(a) y' = y * sin(a) + y * cos(a)

---

Ceci donne une rotation par rapport à l'origine. Bon si tu veux autour d'un axe en 3D, ca revient au meme ( sauf si le plan de rotation n'est pas perpendiculaire par rapport à l'axe de rotation, mais là un calcul en plus et c'est fini. )