Euler angle velocity of 321 Sequence

The method of deriving the Euler angle velocity for a given sequence is to transform each of the derivatives into the reference frame. Remember that an Euler angle sequence is made up of three successive rotations. In other words, the angular velocity ψ needs one rotation, 𝜃 needs two and ϕ needs three.

⌊ ⌋ ⌊ ⌋ ⌊ ⌋ 0 0 ˙ψ ⃗ω = R1 (ψ )R2(𝜃)R3 (ϕ)⌈ 0 ⌉ + R1 (ψ)R2 (𝜃)⌈ ˙𝜃 ⌉ + R1 (ψ )⌈ 0 ⌉ ˙ϕ 0 0

Carrying out the matrix multiplication with R₁(ψ)R₂(𝜃)R₃(ϕ) being the Euler 321 sequence

⌊ ⌋ c𝜃 0 − s𝜃 R (ψ )R (𝜃) = ⌈ s s c s c ⌉ 1 2 ψ 𝜃 ψ ψ 𝜃 cψs𝜃 − sψ cψc𝜃

and

⌊ ⌋ 1 0 0 R (ψ) = ⌈ 0 c s ⌉ 1 ψ ψ 0 − sψ cψ

gives us

⌊ ⌋ ⌊ ˙ ⌋ ⌊ ⌋ ⌊ ⌋ ωx − s𝜃ϕ 0 ψ˙ ⌈ ωy ⌉ = ⌈ sψc𝜃ϕ˙⌉ + ⌈ cψ𝜃˙ ⌉ + ⌈ 0 ⌉ ωz cψc𝜃ϕ˙ − sψ𝜃˙ 0

Adding the vectors together yields

⌊ ⌋ ⌊ ⌋ ωx − s𝜃 ˙ϕ + ψ˙ ⌈ ωy ⌉ = ⌈ sψc𝜃ϕ˙+ cψ𝜃˙ ⌉ ω ˙ ˙ z cψc𝜃ϕ − sψ𝜃

Of course, we also wish to have the Euler angle velocities in terms of the angular velocities which requires us to solve the linear equations for them. Using a program like Matlab makes it easy for us to get

⌊ ⌋ ⌊ ⌋ ˙ϕ (ωys ψ + ωzcψ)sec(𝜃) ⌈ ˙𝜃 ⌉ = ⌈ ω c − ω s ⌉ ˙ y ψ z ψ ψ ωx + ωys ψt𝜃 + ωzc ψt𝜃

In matlab solving for the Euler Angle velocites can be done with the following commands. Using the notation Ax = b, we want to solve for x, such that x = A⁻¹b. For our problem then

syms wx wy wz phd thd psd SPS cth cps b x A;

b = [wx wy wz]’;

x = [phd thd psd]’;

A = [ -sth 0 1; sps*cth cps 0; cps*cth -sps 0];

and solve for the angle velocites with the command

x = inv(A)*b

Note that matlab spits out extra sine and cosine terms that just equal 1 through

2 2 sψ + cψ = 1

The shorthand notation used in this article is

sψ = sin(ψ )

cψ = cos(ψ )

tψ = tan(ψ )

"Euler angle velocity of 321 Sequence" is owned by bloftin.
(view preamble)