An $n\times n$ orthogonal real matrix $A$ is a set ${A_{ij}}$ of $n^2$ real numbers that satisfy the constraints:
- $$\sum_k A_{ik} A_{kj} = \delta_{ij} $$
for all $1\leq i,j\leq n$. The equations (1.) represent
$$ n + \binom{n}{2} = \frac{n(n+1)}{2} $$
constraints on the $n^2$ components of $A$. Therefore an orthogonal matrix $A$ has
$$ n^2 - \frac{n(n+1)}{2} = \frac{n(n-1)}{2} $$
degrees of freedom. Is there a simple representation for these $n(n-1)/2$ parameters? In other words, is there a coordinate system where the constraints (1.) take a simple form? For example, in two dimensions, one can use polar coordinates:
$$ x = r \cos \phi $$ $$ y = r \sin \phi $$
Then the constraint of orthonormality of two vectors $(x_1,y_1)$, $(x_2,y_2)$ takes the simple form: $r_1=r_2=1$, $\phi_1 = \phi_2 \pm \pi/2 $, and a real orthogonal $2\times 2$ matrix has the general parametric form:
$$A = \pm\left(\begin{array}{cc} \cos\phi & \sin\phi\\ -\sin\phi & \cos\phi \end{array}\right)$$
Is there a relatively simple parametric representation like this for higher dimensions?