Salta al contenuto
Note per Studenti Esercizio 133 · composizione di rotazioni e cambio di riferimento

Esercizio 133composizione di rotazioni e cambio di riferimento

In questa pagina 5

Testo. Un albero viene ruotato di 90∘90^\circ attorno all'asse xx e poi di 90∘90^\circ attorno all'asse zz (rotazioni antiorarie, riferimento destrorso).

  1. Scrivere le matrici RxR_x e RzR_z delle due rotazioni e verificare che sono ortogonali.
  2. Calcolare P=RzRxP = R_zR_x (prima RxR_x, poi RzR_z) e RxRzR_xR_z. Le due rotazioni commutano?
  3. Dire di che rotazione si tratta PP: asse, angolo, verso. Calcolare P3P^3.
  4. Un vettore v=(1,2,3)v = (1, 2, 3) (per esempio una velocità) è dato nel riferimento originale. Il riferimento ruotato ha come versori degli assi Pe1,Pe2,Pe3Pe_1, Pe_2, Pe_3: calcolare le componenti di vv nel nuovo riferimento.

Teoria: Isometrie e matrici ortogonaliUna matrice Q è ortogonale se QᵀQ = I, cioè se ha colonne (e righe) ortonormali; equivale a conservare prodotto scalare, norme e angoli, ed è la matrice di una isometria lineare. Ha Q⁻¹ = Qᵀ e det Q = ±1. In R² le ortogonali sono le rotazioni (det 1) e le riflessioni rispetto a una retta (det −1); in R³ le rotazioni hanno un asse (autovalore 1) e un angolo con tr Q = 1 + 2 cos θ, le riflessioni rispetto a un piano hanno Q = I − 2nnᵀ/(n·n). Le isometrie affini sono x ↦ Qx + b.Isometrie e matrici ortogonali →.


1. Le due matrici

La rotazione di 90∘90^\circ attorno a xx lascia fisso e1e_1 e manda e2↦e3e_2 \mapsto e_3, e3↦−e2e_3 \mapsto -e_2 (guardando da xx positivo, yy va in zz). Le colonne sono le immagini dei versori: Rx=(10000−1010),Rz=(0−10100001)R_x = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix}, \qquad R_z = \begin{pmatrix} 0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{pmatrix} (RzR_z: e1↦e2e_1 \mapsto e_2, e2↦−e1e_2 \mapsto -e_1, e3e_3 fisso). Entrambe hanno colonne ortonormali (sono permutazioni di versori, con un segno), quindi sono ortogonali, e det⁡Rx=1⋅(0⋅0−(−1)⋅1)=1\det R_x = 1 \cdot (0\cdot0 - (-1)\cdot1) = 1, det⁡Rz=1\det R_z = 1: sono rotazioni. In ciascuna il blocco 2×22 \times 2 sul piano ortogonale all'asse è RθR_\theta con θ=π2\theta = \frac\pi2, cioè cos⁡θ=0\cos\theta = 0 e sin⁡θ=1\sin\theta = 1.

2. I prodotti

P=RzRx=(0−10100001)(10000−1010)=(001100010),P = R_zR_x = \begin{pmatrix} 0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{pmatrix}\begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix} = \begin{pmatrix} 0 & 0 & 1 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{pmatrix}, RxRz=(10000−1010)(0−10100001)=(0−1000−1100).R_xR_z = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix}\begin{pmatrix} 0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{pmatrix} = \begin{pmatrix} 0 & -1 & 0 \\ 0 & 0 & -1 \\ 1 & 0 & 0 \end{pmatrix} . (Il prodotto ABAB si legge "prima BB, poi AA": RzRxR_zR_x applica prima RxR_x.) Le due matrici sono diverse: le rotazioni nello spazio non commutano. Si vede anche sul versore e2e_2: con RzRxR_zR_x si ha Rxe2=e3R_x e_2 = e_3 e poi Rze3=e3R_z e_3 = e_3, quindi e2↦e3e_2 \mapsto e_3; con RxRzR_xR_z si ha Rze2=−e1R_z e_2 = -e_1 e poi Rx(−e1)=−e1R_x(-e_1) = -e_1, quindi e2↦−e1e_2 \mapsto -e_1. L'ordine delle rotazioni dell'albero conta. (Nel piano invece le rotazioni attorno allo stesso asse commutano.)

3. Natura di PP

PP è prodotto di rotazioni, quindi ortogonale con det⁡P=1\det P = 1: una rotazione. Le sue colonne dicono che PP permuta ciclicamente i versori: e1↦e2↦e3↦e1e_1 \mapsto e_2 \mapsto e_3 \mapsto e_1.

  • Asse. Ker⁡(P−I)\operatorname{Ker}(P - I): da P(x,y,z)=(z,x,y)=(x,y,z)P(x, y, z) = (z, x, y) = (x, y, z) segue x=y=zx = y = z. L'asse è ⟨(1,1,1)⟩\langle(1, 1, 1)\rangle (la diagonale del cubo).
  • Angolo. tr⁡P=0=1+2cos⁡θ\operatorname{tr} P = 0 = 1 + 2\cos\theta, quindi cos⁡θ=−12\cos\theta = -\frac12 e θ=±120∘\theta = \pm120^\circ.
  • Verso. Con u=(1,1,1)u = (1, 1, 1), w=(1,−1,0)w = (1, -1, 0) (ortogonale a uu) e Pw=(0,1,−1)Pw = (0, 1, -1): w×Pw=((−1)(−1)−0⋅1, 0⋅0−1⋅(−1), 1⋅1−(−1)⋅0)=(1,1,1)w \times Pw = \left((-1)(-1) - 0\cdot1,\ 0\cdot0 - 1\cdot(-1),\ 1\cdot1 - (-1)\cdot0\right) = (1, 1, 1) e u⋅(w×Pw)=3>0u \cdot (w \times Pw) = 3 > 0: antioraria. Quindi PP è la rotazione di +120∘+120^\circ attorno a (1,1,1)(1, 1, 1).

Coerenza: tre rotazioni di 120∘120^\circ fanno un giro completo, cioè P3=IP^3 = I (e infatti ciclicamente e1↦e2↦e3↦e1e_1 \mapsto e_2 \mapsto e_3 \mapsto e_1 in tre passi).

Per confronto, RxRzR_xR_z ha ancora traccia 00 e determinante 11: è una rotazione di 120∘120^\circ, ma l'asse è un altro, Ker⁡(RxRz−I)=⟨(1,−1,1)⟩\operatorname{Ker}(R_xR_z - I) = \langle(1, -1, 1)\rangle (da (−y, −z, x)=(x,y,z)(-y,\ -z,\ x) = (x, y, z) segue y=−xy = -x e z=xz = x, cioè (x,y,z)=x (1,−1,1)(x, y, z) = x\,(1, -1, 1)). Le due composizioni sono rotazioni dello stesso angolo attorno ad assi diversi.

4. Componenti nel riferimento ruotato

I versori del nuovo riferimento sono e1′=Pe1=e2e_1' = Pe_1 = e_2, e2′=Pe2=e3e_2' = Pe_2 = e_3, e3′=Pe3=e1e_3' = Pe_3 = e_1 (colonne di PP). Poiché sono ortonormali, le componenti di vv sono i prodotti scalari con essi: v′=(v⋅e1′, v⋅e2′, v⋅e3′)=(v⋅e2, v⋅e3, v⋅e1)=(2,3,1),v' = (v\cdot e_1',\ v\cdot e_2',\ v\cdot e_3') = (v\cdot e_2,\ v\cdot e_3,\ v\cdot e_1) = (2, 3, 1), cioè v′=PTv=(010001100)(123)=(2,3,1)v' = P^Tv = \begin{pmatrix} 0&1&0 \\ 0&0&1 \\ 1&0&0 \end{pmatrix}\begin{pmatrix} 1 \\ 2 \\ 3 \end{pmatrix} = (2, 3, 1). Controllo: ∥v′∥2=4+9+1=14=∥v∥2\|v'\|^2 = 4 + 9 + 1 = 14 = \|v\|^2: il modulo non dipende dal riferimento.


Errori comuni

  • Leggere RzRxR_zR_x da sinistra a destra ("prima RzR_z"): si applica prima la matrice a destra, che è quella adiacente al vettore.
  • Credere che le rotazioni nello spazio commutino.
  • Per le componenti nel riferimento ruotato usare PvPv invece di PTv=P−1vP^Tv = P^{-1}v.

Versione ripasso

Testo. Rotazioni di 90∘90^\circ attorno a xx e a zz: matrici, composizione, asse e angolo del prodotto, componenti di v=(1,2,3)v = (1, 2, 3) nel riferimento ruotato.

  • Rx=(10000−1010)R_x = \begin{pmatrix} 1&0&0 \\ 0&0&-1 \\ 0&1&0 \end{pmatrix}, Rz=(0−10100001)R_z = \begin{pmatrix} 0&-1&0 \\ 1&0&0 \\ 0&0&1 \end{pmatrix}: ortogonali, det⁡1\det 1.
  • P=RzRx=(001100010)P = R_zR_x = \begin{pmatrix} 0&0&1 \\ 1&0&0 \\ 0&1&0 \end{pmatrix} (prima RxR_x); RxRz≠PR_xR_z \ne P: non commutano (es. e2↦e3e_2 \mapsto e_3 contro e2↦−e1e_2 \mapsto -e_1).
  • PP: asse ⟨(1,1,1)⟩\langle(1, 1, 1)\rangle, tr⁡P=0⇒cos⁡θ=−12\operatorname{tr} P = 0 \Rightarrow \cos\theta = -\frac12, verso positivo (u⋅(w×Pw)=3u \cdot (w \times Pw) = 3): rotazione di +120∘+120^\circ; P3=IP^3 = I. RxRzR_xR_z: ancora 120∘120^\circ, asse ⟨(1,−1,1)⟩\langle(1, -1, 1)\rangle.
  • Componenti nel nuovo riferimento: v′=PTv=(2,3,1)v' = P^Tv = (2, 3, 1), ∥v′∥=∥v∥=14\|v'\| = \|v\| = \sqrt{14}.

Errori tipici: ordine dei fattori; credere che commutino; PvPv al posto di PTvP^Tv. Passaggi. RzRxR_z R_x: la colonna jj è Rz(Rxej)R_z(R_x e_j): e1↦e1↦e2e_1 \mapsto e_1 \mapsto e_2, e2↦e3↦e3e_2 \mapsto e_3 \mapsto e_3, e3↦−e2↦e1e_3 \mapsto -e_2 \mapsto e_1. Verso di PP: w=(1,−1,0)w = (1, -1, 0), Pw=(0,1,−1)Pw = (0, 1, -1), w×Pw=(1,1,1)w \times Pw = (1, 1, 1), u⋅(w×Pw)=3>0u \cdot (w \times Pw) = 3 > 0. Componenti: v′=PTv=(v⋅e2, v⋅e3, v⋅e1)=(2,3,1)v' = P^Tv = (v\cdot e_2,\ v \cdot e_3,\ v \cdot e_1) = (2, 3, 1).

Teoria collegata