roty

Rotation matrix for rotations around y-axis

Syntax

R = roty(ang)

Description

R = roty(ang) creates a 3-by-3 matrix used to rotated a 3-by-1 vector or 3-by-N matrix of vectors around the y-axis by ang degrees. When acting on a matrix, each column of the matrix represents a different vector. For the rotation matrix R and vector v, the rotated vector is given by R*v.

Examples

Rotation Matrix for 45° Rotation

Construct the matrix for a rotation of a vector around the y-axis by 45°. Then let the matrix operate on a vector.

R = roty(45)

R = 3×3

    0.7071         0    0.7071
         0    1.0000         0
   -0.7071         0    0.7071

v = [1;-2;4];
y = R*v

y = 3×1

    3.5355
   -2.0000
    2.1213

Under a rotation around the y-axis, the y-component of a vector is invariant.

Input Arguments

ang — Rotation angle

real-valued scalar

Rotation angle specified as a real-valued scalar. The rotation angle is positive if the rotation is in the counter-clockwise direction when viewed by an observer looking along the y-axis towards the origin. Angle units are in degrees.

Example: 30.0

Data Types: double

Output Arguments

R — Rotation matrix

real-valued orthogonal matrix

3-by-3 rotation matrix returned as

$$R_y(\beta )=\left[\begin{array}{ccc}\cos \beta & 0& \sin \beta \\ 0& 1& 0\\ -\sin \beta & 0& \cos \beta \end{array}\right]$$

for a rotation angle β.

More About

Rotation Matrices

Rotation matrices are used to rotate a vector into a new direction.

In transforming vectors in three-dimensional space, rotation matrices are often encountered. Rotation matrices are used in two senses: they can be used to rotate a vector into a new position or they can be used to rotate a coordinate basis (or coordinate system) into a new one. In this case, the vector is left alone but its components in the new basis will be different from those in the original basis. In Euclidean space, there are three basic rotations: one each around the x, y and z axes. Each rotation is specified by an angle of rotation. The rotation angle is defined to be positive for a rotation that is counterclockwise when viewed by an observer looking along the rotation axis towards the origin. Any arbitrary rotation can be composed of a combination of these three (Euler’s rotation theorem). For example, you can rotate a vector in any direction using a sequence of three rotations: $$v^{\prime }=Av=R_z(\gamma )R_y(\beta )R_x(\alpha )v$$.

The rotation matrices that rotate a vector around the x, y, and z-axes are given by:

• Counterclockwise rotation around x-axis

  $$R_x(\alpha )=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \alpha & -\sin \alpha \\ 0& \sin \alpha & \cos \alpha \end{array}\right]$$

• Counterclockwise rotation around y-axis

  $$R_y(\beta )=\left[\begin{array}{ccc}\cos \beta & 0& \sin \beta \\ 0& 1& 0\\ -\sin \beta & 0& \cos \beta \end{array}\right]$$

• Counterclockwise rotation around z-axis

  $$R_z(\gamma )=\left[\begin{array}{ccc}\cos \gamma & -\sin \gamma & 0\\ \sin \gamma & \cos \gamma & 0\\ 0& 0& 1\end{array}\right]$$

The following three figures show what positive rotations look like for each rotation axis:

For any rotation, there is an inverse rotation satisfying $$A^{-1}A=1$$. For example, the inverse of the x-axis rotation matrix is obtained by changing the sign of the angle:

$$R_x^{-1}(\alpha )=R_x(-\alpha )=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos \alpha & \sin \alpha \\ 0& -\sin \alpha & \cos \alpha \end{array}\right]=R_x^{\prime }(\alpha )$$

This example illustrates a basic property: the inverse rotation matrix is the transpose of the original. Rotation matrices satisfy A’A = 1, and consequently det(A) = 1. Under rotations, vector lengths are preserved as well as the angles between vectors.

We can think of rotations in another way. Consider the original set of basis vectors, $$i,j,k$$, and rotate them all using the rotation matrix A. This produces a new set of basis vectors $$i^{\prime },j{,}^{\prime }{k}^{\prime }$$ related to the original by:

$$\begin{array}{ll}{i}^{\prime }\hfill & =Ai\hfill \\ j^{\prime }\hfill & =Aj\hfill \\ k^{\prime }\hfill & =Ak\hfill \end{array}$$

Using the transpose, you can write the new basis vectors as a linear combinations of the old basis vectors:

$$\left[\begin{array}{c}{i}^{\prime }\\ j^{\prime }\\ k^{\prime }\end{array}\right]=A^{\prime }\left[\begin{array}{c}i\\ j\\ k\end{array}\right]$$

Now any vector can be written as a linear combination of either set of basis vectors:

$$v=v_{x}i+v_{y}j+v_{z}k={v^{\prime }}_x{i}^{\prime }+{v^{\prime }}_y{j}^{\prime }+{v^{\prime }}_z{k}^{\prime }$$

Using algebraic manipulation, you can derive the transformation of components for a fixed vector when the basis (or coordinate system) rotates. This transformation uses the transpose of the rotation matrix.

$$\left[\begin{array}{c}{v^{\prime }}_x\\ {v^{\prime }}_y\\ {v^{\prime }}_z\end{array}\right]=A^{-1}\left[\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right]=A^{\prime }\left[\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right]$$

The next figure illustrates how a vector is transformed as the coordinate system rotates around the x-axis. The figure after shows how this transformation can be interpreted as a rotation of the vector in the opposite direction.

References

[1] Goldstein, H., C. Poole and J. Safko, Classical Mechanics, 3rd Edition, San Francisco: Addison Wesley, 2002, pp. 142–144.

Extended Capabilities

C/C++ Code Generation

Generate C and C++ code using MATLAB® Coder™.

Usage notes and limitations:

Does not support variable-size inputs.