Fact-checked by Grok 2 weeks ago

Orthogonal matrix

In linear algebra, an orthogonal matrix is a real square matrix Q whose columns and rows form an orthonormal set of vectors, equivalently satisfying Q^T Q = I_n, where Q^T denotes the transpose of Q and I_n is the n \times n identity matrix.^[1] This condition implies that Q is invertible with inverse equal to its transpose, Q^{-1} = Q^T.^[2] Orthogonal matrices preserve the Euclidean inner product, ensuring that \| Q \mathbf{v} \| = \| \mathbf{v} \| for any vector \mathbf{v} and that angles between vectors are maintained up to sign, thus representing isometries such as rotations and reflections in Euclidean space.^[3] The determinant of an orthogonal matrix is either +1 or -1, with those having determinant +1 forming the special orthogonal group SO(n), which consists of proper rotations.^[2]^[1] The collection of all n \times n orthogonal matrices constitutes the orthogonal group O(n), a Lie group that is compact and plays a central role in the study of symmetries and transformations.^[1] Eigenvalues of orthogonal matrices lie on the unit circle in the complex plane, with absolute value 1.^[2] Orthogonal matrices are essential in numerical linear algebra for algorithms like QR decomposition and Householder reflections, which stabilize computations by preserving norms and orthogonality.^[4] They also arise in applications such as analytic geometry for axis rotations and quadratic form reductions, as well as in computer graphics for modeling rigid body motions.^[5]^[6]

Definition and Fundamentals

Definition

An orthogonal matrix is a square matrix Q over the real numbers such that its transpose Q^T, obtained by interchanging the rows and columns of Q, satisfies Q^T Q = I, where I is the identity matrix of the same dimension.^[7]^[8] This condition is equivalent to the columns (and similarly the rows) of Q forming an orthonormal basis for the Euclidean space, meaning each column has Euclidean norm 1 and the inner product between distinct columns is zero, with the inner product defined as \langle x, y \rangle = x^T y for real vectors x and y.^[7]^[9] A direct consequence of this definition is that orthogonal matrices preserve the Euclidean norm of vectors, so \| Q x \| = \| x \| for all vectors x, where the Euclidean norm is \| x \| = \sqrt{\langle x, x \rangle}.^[10]^[8] In the complex case, the analogous concept is a unitary matrix, which satisfies Q^* Q = I where Q^* is the conjugate transpose, serving as the complex extension of orthogonal matrices while preserving the Hermitian inner product.^[11]

Basic Properties

A square matrix Q is orthogonal if it satisfies Q^T Q = I, where I is the identity matrix and Q^T denotes the transpose of Q. This condition implies that the columns of Q are orthonormal. From this defining relation, it follows that Q is invertible, with the inverse given by Q^{-1} = Q^T. The proof proceeds by noting that Q^T Q = I establishes Q^T as a left inverse of Q; since Q is square, this left inverse is also the right inverse, confirming invertibility.^[12] A symmetric consequence of the definition is the relation Q Q^T = I, which similarly shows that the rows of Q are orthonormal and reinforces the invertibility property.^[12] The orthogonality condition leads to the preservation of the Euclidean inner product: for any vectors x and y, \langle Qx, Qy \rangle = (Qx)^T (Qy) = x^T Q^T Q y = x^T y = \langle x, y \rangle.^[12] Setting y = x yields norm preservation: \| Qx \| = \| x \|, since \| Qx \|^2 = \langle Qx, Qx \rangle = \langle x, x \rangle = \| x \|^2. This inner product preservation also implies that angles between vectors are maintained, as the cosine of the angle \theta between x and y is \cos \theta = \frac{\langle x, y \rangle}{\| x \| \| y \|}, and each component remains unchanged under the transformation by Q.^[12] The determinant of an orthogonal matrix satisfies |\det(Q)| = 1, so \det(Q) = \pm 1. This follows from \det(Q^T Q) = \det(I) = 1 and the property that \det(Q^T) = \det(Q), yielding [\det(Q)]^2 = 1.^[12] Regarding eigenvalues, if Q v = \lambda v for some nonzero vector v and scalar \lambda, then norm preservation gives \| Q v \| = \| v \|, so |\lambda| \| v \| = \| v \|, implying |\lambda| = 1. Thus, the eigenvalues of Q lie on the unit circle in the complex plane.^[12]^[13]

Examples and Interpretations

Concrete Examples

The identity matrix I of any dimension is the simplest orthogonal matrix, satisfying I^T I = I by definition, and it preserves the Euclidean norm of vectors as \|I \mathbf{v}\| = \|\mathbf{v}\| for any vector \mathbf{v}.^[14] A fundamental example in two dimensions is the rotation matrix Q = \begin{pmatrix} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{pmatrix}, which is orthogonal because its transpose is Q^T = \begin{pmatrix} \cos \theta & \sin \theta \\ -\sin \theta & \cos \theta \end{pmatrix}, and their product yields

Q^T Q = \begin{pmatrix} \cos^2 \theta + \sin^2 \theta & 0 \\ 0 & \sin^2 \theta + \cos^2 \theta \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} = I,

using the Pythagorean identity \cos^2 \theta + \sin^2 \theta = 1.^[15] For reflections, consider the matrix Q = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}, which flips the sign of the second coordinate; it is orthogonal since Q^T = Q and Q^2 = I, confirming Q^T Q = I. This arises as a Householder reflection with unit vector \mathbf{a} = (0, 1)^T, via Q = I - 2 \mathbf{a} \mathbf{a}^T.^[14] Permutation matrices form another class of orthogonal matrices with determinant \pm 1; for instance, the 2×2 matrix Q = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} swaps coordinates and satisfies orthogonality because each row and column has a single 1 with zeros elsewhere, ensuring Q^T Q = I.^[14]

Geometric Interpretations

Orthogonal matrices represent linear isometries of Euclidean space \mathbb{R}^n, which are transformations that preserve distances between points, thereby maintaining the geometric structure of the space. Specifically, multiplication by an orthogonal matrix Q defines a linear map T: \mathbb{R}^n \to \mathbb{R}^n given by T(\mathbf{x}) = Q\mathbf{x} that satisfies \|T(\mathbf{x}) - T(\mathbf{y})\| = \|\mathbf{x} - \mathbf{y}\| for all \mathbf{x}, \mathbf{y} \in \mathbb{R}^n. These isometries encompass both rotations and reflections: the special orthogonal group SO(n), consisting of orthogonal matrices with determinant 1, corresponds to proper rotations that rigidly turn the space without flipping it, while the full orthogonal group O(n), including matrices with determinant -1, incorporates reflections that reverse handedness.^[16]^[4]^[17] The determinant of an orthogonal matrix distinguishes between orientation-preserving and orientation-reversing transformations. Matrices in SO(n) with \det Q = 1 preserve the orientation of \mathbb{R}^n, meaning they map positively oriented bases to positively oriented bases, akin to rotations in 3D space that do not mirror objects. In contrast, those with \det Q = -1 reverse orientation, effectively including a reflection component that swaps the direction of right-handed coordinate systems. This distinction arises because the determinant measures the signed volume scaling factor, which remains \pm 1 for orthogonal matrices but indicates whether the transformation inverts spatial handedness.^[18]^[19] Orthogonal matrices act on the unit sphere S^{n-1} = \{\mathbf{x} \in \mathbb{R}^n : \|\mathbf{x}\| = 1\} by mapping it to itself, as they preserve the Euclidean norm \|\mathbf{x}\| = 1 for all vectors on the sphere. This isometry of the sphere highlights their role in maintaining angular relationships and surface geometry without distortion. For instance, rotation matrices in SO(n) interpret as pure rotations of vectors around the origin, altering directions but not lengths or angles between them, which geometrically corresponds to spinning objects in space while keeping their shape intact.^[16]^[15]

Constructions

Elementary Operations

One fundamental way to construct orthogonal matrices is through Householder reflections, which are reflections across hyperplanes orthogonal to a given vector. The orthogonal projector onto the span of a nonzero vector \mathbf{v} \in \mathbb{R}^n is given by P = \frac{\mathbf{v} \mathbf{v}^T}{\mathbf{v}^T \mathbf{v}}, which projects vectors onto the line spanned by \mathbf{v}.^[20] The corresponding Householder reflection matrix is then H = I - 2P, where I is the n \times n identity matrix, resulting in an orthogonal matrix that reflects vectors across the hyperplane perpendicular to \mathbf{v}.^[20] For a unit vector \mathbf{u} (so \mathbf{u}^T \mathbf{u} = 1), this simplifies to the general form

H = I - 2 \mathbf{u} \mathbf{u}^T.

Such reflections are symmetric (H^T = H) and involutory (H^2 = I), confirming their orthogonality.^[21] Another elementary construction involves Givens rotations, which generalize 2D rotation matrices to higher dimensions by rotating within specific planes. A Givens rotation matrix G is an orthogonal matrix obtained by modifying the identity matrix in the (i,j) and (j,i) positions with cosine c and sine s values satisfying c^2 + s^2 = 1, leaving all other entries unchanged; for example, in the i-th and j-th rows/columns,

G = \begin{pmatrix} \cdots & & c & & s & \\ \cdots & -s & & c & & \\ \cdots & & & \cdots & & \cdots \end{pmatrix}.

^[22] These are constructed to zero out specific entries in vectors or matrices, such as setting c = a / \sqrt{a^2 + b^2} and s = b / \sqrt{a^2 + b^2} for entries a and b.^[22] In higher dimensions, a Givens rotation acts as the identity outside the chosen 2D plane spanned by the i-th and j-th standard basis vectors, preserving orthogonality since the columns remain orthonormal.^[22] The orthogonal group O(n) can be generated by products of these elementary orthogonal matrices, specifically Householder reflections and Givens rotations. By the Cartan–Dieudonné theorem, every matrix in O(n) is a product of at most n reflections (Householder matrices), and rotations (Givens matrices) can be expressed as products of two reflections, allowing full generation of O(n) from these building blocks.^[17] This decomposition is particularly useful for algorithmic constructions of orthogonal matrices.^[17]

Low-Dimensional Cases

In two dimensions, orthogonal matrices with determinant 1, known as elements of the special orthogonal group SO(2), can be parameterized by a single angle θ, representing rotations around the origin. The rotation matrix takes the form

\begin{pmatrix} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{pmatrix},

where θ ranges from 0 to 2π, and this parameterization identifies SO(2) with the circle group, as the group operation corresponds to addition of angles modulo 2π.^[23]^[24] In three dimensions, rotations belonging to SO(3) require three parameters for parameterization, with common methods including Euler angles and unit quaternions. Euler angles (φ, θ, ψ) describe a sequence of rotations: first by φ about the z-axis, then by θ about the new x-axis, and finally by ψ about the new z-axis, yielding the overall rotation matrix as the product of the individual axis rotation matrices.^[25] Alternatively, unit quaternions q = (cos(θ/2), n sin(θ/2)), where n is a unit vector along the rotation axis and θ is the rotation angle, represent rotations via conjugation: a vector p is mapped to q p q^, with q^ the quaternion conjugate, offering a compact and singularity-free parameterization for SO(3).^[26] A direct formula for 3D rotation matrices is provided by Rodrigues' rotation formula, which constructs the matrix R for a rotation by angle θ about a unit axis vector ω = (ω_x, ω_y, ω_z) as

R = I + \sin \theta \, K + (1 - \cos \theta) K^2,

where I is the 3×3 identity matrix and K is the skew-symmetric cross-product matrix

K = \begin{pmatrix} 0 & -\omega_z & \omega_y \\ \omega_z & 0 & -\omega_x \\ -\omega_y & \omega_x & 0 \end{pmatrix}.

This formula efficiently computes the rotation while preserving orthogonality and determinant 1.^[27] Orthogonal matrices with determinant -1 in 3D include reflections across planes, which can be constructed using the Householder reflection formula. For a plane with unit normal vector v, the reflection matrix H is

H = I - 2 \frac{v v^T}{v^T v},

where v v^T is the outer product; this matrix reflects vectors across the plane perpendicular to v, fixing the plane pointwise while reversing the component along v.^[28]

Higher-Dimensional Generalizations

One approach to constructing orthogonal matrices in higher dimensions involves block-diagonal forms, where the matrix is assembled as a direct sum of smaller orthogonal matrices. Specifically, if Q_1, Q_2, \dots, Q_k are orthogonal matrices of sizes n_1 \times n_1, n_2 \times n_2, \dots, n_k \times n_k with n_1 + \dots + n_k = n, then the block-diagonal matrix Q = \operatorname{diag}(Q_1, Q_2, \dots, Q_k) satisfies Q Q^T = I_n, making it orthogonal. This method leverages known lower-dimensional constructions to build scalable higher-dimensional examples, though it preserves the overall structure without introducing novel interactions between blocks. Hadamard matrices provide another explicit construction for certain dimensions, particularly powers of 2. A Hadamard matrix H_n of order n is a square matrix with entries \pm 1 such that H_n H_n^T = n I_n, implying that \frac{1}{\sqrt{n}} H_n is an orthogonal matrix.^[29] The Sylvester construction recursively builds these matrices: starting from

H_1 = {{grok:render&&&type=render_inline_citation&&&citation_id=1&&&citation_type=wikipedia}}

, H_{2m} = \begin{pmatrix} H_m & H_m \\ H_m & -H_m \end{pmatrix}, yielding orthogonal matrices (up to scaling) for n = 2^k.^[29] Existence for other orders remains conjectural, with the Hadamard conjecture positing availability for all multiples of 4, but verified only up to moderate sizes. Clifford algebras offer a algebraic framework for generating orthogonal matrices, particularly through representations of the orthogonal group. The Clifford algebra Cl(n) over \mathbb{R}^n with the standard quadratic form is generated by an orthonormal basis \{e_i\} satisfying e_i e_j + e_j e_i = -2 \delta_{ij}, and its even subalgebra yields the spin group \operatorname{Spin}(n), a double cover of SO(n).^[30] Elements of \operatorname{Spin}(n) act on \mathbb{R}^n via the twisted adjoint representation, producing rotations as products of even numbers of reflections, in line with the Cartan-Dieudonné theorem; exponentiating bivectors like \exp(\theta/2 \, e_i e_j) generates plane rotations that combine into full orthogonal matrices.^[30] This approach is systematic for arbitrary n, facilitating computational generation in specific signatures. Parameterizing general orthogonal matrices in dimensions n \geq 4 presents significant challenges compared to lower dimensions. While SO(n) has dimension n(n-1)/2, allowing parameterization by that many angles via successive plane rotations (generalized Euler angles), no simple, singularity-free analogue to the three Euler angles of SO(3) exists; instead, decompositions like products of Givens rotations introduce complexities in uniformity and computational stability. Rational parameterizations are possible but often require solving Diophantine equations, limiting explicit constructions.^[31]

Core Properties

Matrix and Vector Properties

Orthogonal matrices play a fundamental role in preserving key structural properties of matrices and vectors under linear transformations. Specifically, if Q is an orthogonal matrix, then for any square matrix A, the transformation Q^T A Q is a similarity transformation, which preserves the eigenvalues of A. This follows from the general property of similarity transformations with invertible matrices, and since orthogonal matrices are invertible with Q^{-1} = Q^T, they maintain the spectrum of A.^[32]^[33] For symmetric matrices, this preservation is particularly significant. Consider a symmetric matrix A, so A = A^T. The eigenvalues of Q^T A Q are identical to those of A, as the similarity ensures that the characteristic polynomial remains unchanged:

\det(\lambda I - Q^T A Q) = \det(Q^T (\lambda I - A) Q) = \det(\lambda I - A).

This property allows orthogonal transformations to diagonalize symmetric matrices while retaining their spectral information, a cornerstone in applications like principal component analysis.^[34]^[35] Beyond eigenvalues, orthogonal matrices exhibit invariance properties that extend to vector norms and matrix traces. For any vector \mathbf{x}, the Euclidean norm is preserved under multiplication by an orthogonal matrix Q, satisfying \| Q \mathbf{x} \|_2 = \| \mathbf{x} \|_2. This isometry arises because Q^T Q = I, ensuring that inner products and lengths remain unchanged, which geometrically corresponds to rotations or reflections in Euclidean space.^[36]^[4] Similarly, the trace of a matrix is invariant under orthogonal conjugation: \operatorname{tr}(Q^T A Q) = \operatorname{tr}(A). This cyclic property of the trace, combined with the orthogonality condition, implies that such transformations do not alter the sum of the diagonal elements, providing a robust invariant for analyzing matrix similarities.^[37]^[35] Orthogonal matrices are integral to the singular value decomposition (SVD) of any real matrix M, expressed as M = U \Sigma V^T, where U and V are orthogonal matrices and \Sigma is a diagonal matrix containing the singular values. The orthogonal factors U and V capture the rotational components of the transformation, separating the scaling effects encoded in \Sigma, which enables efficient computations in data analysis and compression.^[38]^[39]

Group and Algebraic Structure

The orthogonal group O(n) is the set of all n \times n real matrices Q such that Q^T Q = I_n, forming a group under matrix multiplication since the product of two such matrices remains orthogonal: for Q_1, Q_2 \in O(n), (Q_1 Q_2)^T (Q_1 Q_2) = Q_2^T Q_1^T Q_1 Q_2 = Q_2^T Q_2 = I_n.^[40] The inverse of any Q \in O(n) is its transpose Q^T, which is also orthogonal, ensuring the group axioms of closure under inversion and the identity element I_n.^[40] A key subgroup is the special orthogonal group SO(n) = \{ Q \in O(n) \mid \det Q = 1 \}, which is normal in O(n) because the determinant is a continuous group homomorphism \det: O(n) \to \{ \pm 1 \} satisfying \det(Q_1 Q_2) = \det Q_1 \cdot \det Q_2, thereby preserving the multiplicative structure of \{ \pm 1 \}.^[40] As a Lie group, O(n) is compact with two connected components distinguished by the sign of the determinant, where SO(n) forms the connected component of the identity for n \geq 2.^[40] The irreducible representations of O(n) encompass the standard action on \mathbb{R}^n and its exterior powers, featuring weights like \pm \omega_1, \dots, \pm \omega_n in the defining representation, with the adjoint representation having highest weight \omega_1 + \omega_2 for types B_n and D_n.^[40] The associated Weyl group, which governs the structure of these representations through its action on weights via permutations and sign changes, is the hyperoctahedral group S_n \ltimes (\mathbb{Z}_2)^n for SO(2n+1) and S_n \ltimes (\mathbb{Z}_2)^{n-1} for SO(2n).^[40]

Canonical Forms

Over the real numbers, every orthogonal matrix Q \in O(n) is orthogonally similar to a block-diagonal matrix consisting of 1×1 blocks with entries \pm 1 and 2×2 rotation blocks of the form \begin{pmatrix} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{pmatrix} for \theta \in [0, \pi).^[41] This form arises from the real Schur decomposition adapted to the normal structure of orthogonal matrices, where the eigenvalues are on the unit circle, leading to real eigenvalues \pm 1 on the diagonal for 1×1 blocks and complex conjugate pairs corresponding to the 2×2 blocks.^[42] The number of 1×1 blocks with +1 and -1, along with the angles \theta in the 2×2 blocks, classify the equivalence class under orthogonal similarity.^[41] The Cartan-Dieudonné theorem provides another canonical decomposition, stating that every orthogonal transformation in O(n, \mathbb{R}) (except the identity) is the composition of at most n reflections about hyperplanes.^[17] Reflections are involutory orthogonal matrices represented by Householder transformations, and the theorem follows by induction on the dimension, reducing to cases where the transformation fixes a hyperplane or maps no nonzero vector to itself.^[17] Thus, any Q \in O(n) can be expressed as Q = s_1 s_2 \cdots s_k where each s_i is a reflection matrix and k \leq n, with the parity of k determining the determinant of Q.^[17] Involutory orthogonal matrices, satisfying Q^2 = I, form a special subclass where the canonical form simplifies further. Such matrices are symmetric (Q = Q^T) and diagonalizable by an orthogonal matrix to a diagonal form \operatorname{diag}(I_r, -I_{n-r}) for some r \in \{0, 1, \dots, n\}, consisting solely of 1×1 blocks with \pm 1 and no 2×2 rotation blocks.^[43] The integer r (the number of +1 eigenvalues) classifies these matrices up to orthogonal similarity, and they correspond to reflections or products of disjoint reflections.^[43]

Advanced Mathematical Aspects

Lie Algebra Connections

The Lie algebra of the special orthogonal group SO(n), denoted \mathfrak{so}(n), consists of all n \times n real skew-symmetric matrices, i.e., matrices A satisfying A^T = -A.^[44] This vector space, of dimension \frac{n(n-1)}{2}, forms the tangent space at the identity element of SO(n) and captures the infinitesimal generators of rotations.^[44] The Lie bracket on \mathfrak{so}(n) is given by the matrix commutator [A, B] = AB - BA, which endows it with the structure of a Lie algebra and encodes the non-commutativity of the group operation near the identity.^[44] The exponential map provides a fundamental connection between \mathfrak{so}(n) and SO(n), defined by \exp: \mathfrak{so}(n) \to \mathrm{SO}(n), where for a skew-symmetric matrix A, \exp(A) = \sum_{k=0}^\infty \frac{A^k}{k!}. This map is a surjective Lie group homomorphism, ensuring that every rotation matrix in SO(n) arises as the exponential of some element in its Lie algebra.^[44] In particular, for a skew-symmetric matrix K and scalar \theta, the one-parameter subgroup Q(\theta) = \exp(\theta K) generates a continuous family of rotations, preserving lengths and angles infinitesimally.^[44] The adjoint representation of SO(n) acts on \mathfrak{so}(n) by conjugation: for Q \in \mathrm{SO}(n) and X \in \mathfrak{so}(n), \mathrm{Ad}(Q)(X) = Q X Q^{-1}, which preserves the skew-symmetry and thus maps \mathfrak{so}(n) to itself.^[45] At the Lie algebra level, the infinitesimal adjoint action is \mathrm{ad}(X)(Y) = [X, Y] for X, Y \in \mathfrak{so}(n), reflecting how infinitesimal rotations transform the generators of the algebra via commutators.^[45] This representation, of dimension \frac{n(n-1)}{2}, is irreducible for n \geq 3 except for n=4, where it decomposes as the direct sum of two 3-dimensional irreducible representations, and plays a key role in understanding symmetries in representation theory.^[46]

Spin and Pin Extensions

The spin group Spin(n) serves as the universal double cover of the special orthogonal group SO(n), capturing the simply connected topology required for representing rotations in higher dimensions where SO(n) is not simply connected for n ≥ 3.^[30] It is constructed within the even subalgebra of the Clifford algebra Cl(n) associated to the Euclidean space ℝⁿ, consisting of those invertible even-grade elements that preserve the underlying vector space under the twisted adjoint action and have norm ±1.^[30] This embedding reflects the double-covering homomorphism from Spin(n) to SO(n), where rotations by angle θ in SO(n) correspond to elements in Spin(n) that trace paths closing after 4π rather than 2π, essential for consistent spin representations in physics.^[47] The pin group Pin(n) extends this construction to a double cover of the full orthogonal group O(n), incorporating both rotations and reflections.^[47] Defined as a subgroup of the multiplicative group of the Clifford algebra Cl(n), Pin(n) includes elements of both even and odd grades that normalize the vector space via the twisted adjoint action, with the projection map to O(n) having kernel {±1}.^[47] Reflections in O(n), generated by householder transformations orthogonal to unit vectors, lift to odd-grade elements like the basis vectors in the Clifford algebra, enabling Pin(n) to represent improper orthogonal transformations such as parity inversions.^[47] Spin(n) arises as the even-grade subgroup of Pin(n), distinguishing pure rotations from those combined with reflections.^[47] Spinor representations arise as the faithful irreducible representations of Spin(n) that do not factor through SO(n), transforming the Clifford module spaces under the group action.^[48] These representations act on spinors—vectors in the spinor space S, typically of dimension 2^{⌊n/2⌋}—via left multiplication in the Clifford algebra, where an element g ∈ Spin(n) maps a spinor ψ to gψ.^[48] Under the covering map to SO(n), spinors acquire a projective character, changing sign under 2π rotations, which resolves the half-integer spin of fermions in quantum mechanics.^[48] For even n = 2m, the spinor space decomposes into two chiral half-spinor representations, each of dimension 2^{m-1}, while for odd n = 2m+1, there is a single irreducible spinor representation of dimension 2^m.^[48] The development of spinors gained physical prominence in the late 1920s through Paul Dirac's formulation of the relativistic wave equation for electrons, which naturally incorporated spin-1/2 particles using four-component spinors transforming under the Lorentz group, a extension of orthogonal transformations in Minkowski space.^[49] This post-1928 work unified quantum mechanics with special relativity, predicting antimatter and laying the foundation for quantum field theory, where spinors describe fermion fields under orthogonal and Lorentz boosts.^[50]

Numerical and Computational Uses

Computational Advantages

Orthogonal matrices offer significant advantages in numerical computations primarily due to their inherent numerical stability. The condition number \kappa_2(Q) of an orthogonal matrix Q, defined as the ratio of the largest to smallest singular value, equals 1 in the Euclidean (2-)norm, which is the minimal possible value for any nonsingular matrix. This property ensures that perturbations in the input data lead to proportionally small changes in the output, preventing the amplification of rounding errors that is common in computations with ill-conditioned matrices. As a result, algorithms involving orthogonal transformations, such as those in linear algebra routines, maintain high accuracy even in finite-precision arithmetic. In floating-point arithmetic, multiplication by an orthogonal matrix preserves the relative error bound of the input vectors. Specifically, if x is a vector with relative error bounded by the machine epsilon \epsilon, then the computed Qx satisfies \|Qx - \mathrm{fl}(Qx)\| / \|Qx\| \approx \epsilon, because orthogonal transformations preserve the Euclidean norm \|Qx\| = \|x\|. Unlike multiplications with general matrices, where the condition number \kappa(A) > 1 can cause rounding errors to grow exponentially with repeated operations, orthogonal matrix multiplications exhibit no such growth, keeping errors bounded by O(n \epsilon) for an n \times n matrix in the worst case. This stability is particularly beneficial in iterative processes, such as optimization algorithms and physical simulations, where repeated applications of transformations accumulate errors minimally over many steps. The preservation of norms and inner products by orthogonal matrices further enhances their utility in computational settings, ensuring that geometric interpretations and distances remain accurate without distortion from numerical artifacts.

Decompositions and Algorithms

One key decomposition involving orthogonal matrices is the QR decomposition, which factors an m \times n matrix A (with m \geq n) into A = QR, where Q is an m \times m orthogonal matrix and R is an m \times n upper triangular matrix.^[51] This factorization can be obtained via the classical Gram-Schmidt process, which orthogonalizes the columns of A sequentially by subtracting projections onto the subspace spanned by prior columns to produce the columns of Q, with R storing the corresponding coefficients.^[51] The Gram-Schmidt approach, while conceptually straightforward, can suffer from numerical instability due to subtractive cancellation in finite precision.^[51] A more robust alternative employs Householder reflections, which triangularize A through a series of orthogonal transformations that zero out elements below the diagonal in each column.^[52] Each reflection is defined by a vector v such that the Householder matrix H = I - 2 \frac{v v^T}{v^T v} is orthogonal and symmetric, applied from the left to progressively eliminate subdiagonal entries.^[52] The product of these n Householder matrices yields Q, and the accumulated transformations produce R. This method, introduced by Householder in 1958, ensures better preservation of orthogonality in floating-point arithmetic.^[52] Another fundamental decomposition is the singular value decomposition (SVD), which expresses A as A = U \Sigma V^T, where U and V are orthogonal matrices whose columns are the left and right singular vectors, respectively, and \Sigma is a diagonal matrix with non-negative singular values on the diagonal in decreasing order.^[53] The SVD's theoretical foundation traces to the late 19th century, with independent discoveries by Beltrami in 1873 for the real case and Jordan in 1874 for the general case, revealing A's action as a rotation, scaling, and another rotation.^[53] Computationally, the SVD is often derived via iterative methods like the Golub-Reinsch algorithm, which builds on bidiagonalization followed by orthogonal transformations to diagonalize.^[54] Orthogonal matrices play a central role in eigenvalue algorithms through similarities that preserve eigenvalues while simplifying structure. The QR algorithm, a cornerstone for computing all eigenvalues and eigenvectors, relies on repeated QR decompositions to generate a sequence of orthogonal similarities converging to upper triangular (Schur) form. Developed by Francis in 1961, the basic unshifted iteration starts with a matrix A_0 = A, computes its QR factorization A_k = Q_k R_k, and sets A_{k+1} = R_k Q_k, ensuring A_{k+1} is similar to A_k via the orthogonal Q_k. For symmetric real matrices, efficiency is enhanced by first reducing A to symmetric tridiagonal form using Householder transformations, which eliminate off-tridiagonal elements while preserving symmetry and eigenvalues.^[52] The subsequent QR iterations on the tridiagonal matrix then deflate as eigenvalues appear on the diagonal. The iterative step on the tridiagonal T_k is

\begin{aligned} T_k &= Q_k R_k, \\ T_{k+1} &= R_k Q_k, \end{aligned}

where Q_k is orthogonal and bidiagonal (from Givens rotations or similar), promoting quadratic convergence near eigenvalues. Shifts, such as Wilkinson shifts, accelerate global convergence by targeting the bottom-right eigenvalue at each step.

Approximation Techniques

One common method for generating random orthogonal matrices distributed according to the Haar measure on the orthogonal group involves constructing a matrix with independent standard Gaussian entries and applying the QR decomposition, retaining only the orthogonal factor Q.^[55] This approach ensures uniformity because the QR factorization of a Gaussian random matrix yields an orthogonal Q that is Haar-distributed, a property rooted in the invariance of the Gaussian distribution under orthogonal transformations.^[56] The resulting Q satisfies Q^T Q = I and is widely used in Monte Carlo simulations and randomized algorithms for its computational efficiency and statistical properties.^[55] A key approximation problem is finding the nearest orthogonal matrix Q to a given square matrix A, which minimizes the Frobenius norm ||A - Q||_F subject to Q^T Q = I; this is known as the (unbalanced) orthogonal Procrustes problem.^[57] The solution is the unitary factor in the polar decomposition of A, providing the optimal orthogonal approximation in the Frobenius sense.^[58] Specifically, when A is invertible, this orthogonal factor is given by

Q = A (A^T A)^{-1/2},

where (A^T A)^{-1/2} denotes the unique positive definite square root of A^T A.^[58] This approximation is unique and arises in applications such as data alignment and attitude determination, where preserving orthogonality while minimizing deviation from A is essential.^[57] The Cayley transform provides another technique for parameterizing and approximating orthogonal matrices by mapping skew-symmetric matrices to the special orthogonal group, excluding matrices with -1 as an eigenvalue.^[59] For a skew-symmetric matrix S (S^T = -S), the transform is defined as Q = (I - S)(I + S)^{-1}, yielding Q^T Q = I and det(Q) = 1, provided I + S is invertible.^[59] This mapping is bijective onto the special orthogonal matrices without -1 eigenvalues and is useful for optimization and geodesic computations on the orthogonal manifold due to its rational form and avoidance of trigonometric functions.^[59]

Generalizations

Rectangular Orthogonal Matrices

A rectangular orthogonal matrix, also known as a semi-orthogonal matrix, extends the notion of orthogonality to non-square matrices. Consider an m \times n matrix Q where m \neq n. If m \geq n (a tall matrix), Q satisfies Q^\top Q = I_n, ensuring its columns form an orthonormal set in \mathbb{R}^m. Conversely, if m \leq n (a fat matrix), Q satisfies Q Q^\top = I_m, ensuring its rows form an orthonormal set in \mathbb{R}^n.^[60] For the tall case, the defining equation is

Q^\top Q = I_n,

which implies that the columns of Q are orthonormal, as the (i,j)-th entry of Q^\top Q equals the inner product of the i-th and j-th columns, yielding the Kronecker delta \delta_{ij}. This orthonormality extends to the fat case analogously for rows.^[60] These matrices exhibit partial isometry properties, preserving norms within their defined subspaces. For a tall Q, left-multiplication by Q acts as an isometry from \mathbb{R}^n to the column space of Q: for any \mathbf{x} \in \mathbb{R}^n,

\| Q \mathbf{x} \|_2 = \| \mathbf{x} \|_2,

since

\| Q \mathbf{x} \|_2^2 = \mathbf{x}^\top Q^\top Q \mathbf{x} = \mathbf{x}^\top I_n \mathbf{x} = \| \mathbf{x} \|_2^2.

A similar preservation holds for the row space in the fat case.^[61] Rectangular orthogonal matrices find key applications in numerical methods for overdetermined systems and data analysis. In solving least squares problems \min_{\mathbf{x}} \| A \mathbf{x} - \mathbf{b} \|_2 for an m \times n matrix A with m > n and full column rank, the thin QR decomposition A = Q R—where Q is m \times n rectangular orthogonal and R is n \times n upper triangular—yields the solution by forward substitution in R \mathbf{x} = Q^\top \mathbf{b}. This factorization leverages the orthogonality of Q for numerical stability and efficiency.^[62] They also play a central role in data reduction techniques, such as dimensionality reduction, by providing an orthonormal basis for projecting high-dimensional data onto a lower-dimensional subspace while maintaining inner products within that subspace. In principal component analysis (PCA), the leading principal components form the columns of a tall rectangular orthogonal matrix that captures the maximum variance when projecting the data.^[63] Unitary matrices represent the complex generalization of orthogonal matrices, extending the preservation of the Euclidean inner product from real vector spaces to complex ones. Specifically, a complex square matrix U is unitary if it satisfies U^\dagger U = I, where U^\dagger denotes the conjugate transpose, analogous to how an orthogonal matrix Q satisfies Q^T Q = I for real entries.^[64] This generalization allows unitary matrices to preserve the Hermitian inner product in \mathbb{C}^n, mirroring the length-preserving property of orthogonal matrices in \mathbb{R}^n.^[65] Symplectic matrices, in contrast, preserve a specific type of nondegenerate alternating bilinear form rather than the standard inner product. Defined over even-dimensional real or complex vector spaces, a matrix S is symplectic if S^T J S = J, where J is the standard symplectic matrix, ensuring the invariance of the symplectic form \omega(u, v) = u^T J v.^[66] This structure arises in Hamiltonian mechanics and phase space transformations, distinguishing symplectic matrices from orthogonal ones by their focus on area-preserving rather than length-preserving maps.^[67] The Stiefel manifold provides a geometric framework encompassing orthogonal matrices as the space of orthonormal frames in Euclidean space. Denoted V_{k,n}(\mathbb{R}), it consists of all ordered sets of k orthonormal vectors in \mathbb{R}^n (with k \leq n), which can be represented by n \times k matrices with orthonormal columns; when k = n, this reduces to the orthogonal group O(n).^[68] This manifold generalizes the concept of orthogonal matrices to rectangular cases, serving as a foundational structure in differential geometry and optimization.^[69] Orthogonal matrices form a special case of normal matrices, which are square matrices A satisfying A^\dagger A = A A^\dagger, ensuring simultaneous diagonalization by a unitary matrix. Since real orthogonal matrices are unitary (with conjugate transpose reducing to transpose) and thus commute with their adjoint, they inherit the spectral properties of normal matrices, such as having orthogonal eigenvectors.^[70]^[71]

References

[1]
[PDF] RES.18-011 (Fall 2021) Lecture 12: Orthogonal Matrices
Orthogonal matrices preserve lengths, as well as preserving angles up to sign. In general, a set of matrices satisfying some well-behaved properties of a set ...
[2]
[PDF] orthogonal matrices - UTK Math
Basic properties. (1) A matrix is orthogonal exactly when its column vectors have length one, and are pairwise orthogonal; likewise for the row vectors. In ...
[3]
[PDF] Math 224 Properties of Orthogonal Matrices - Kenyon College
The vectors formed by the first and last rows of an orthogonal matrix must be orthogonal. Mathematics Department. 1. Math 224: Linear Algebra. Page 2. Kenyon ...
[4]
Orthogonal Transformations and Orthogonal Matrices - UTSA
Jan 29, 2022 · Numerical analysis takes advantage of many of the properties of orthogonal matrices for numerical linear algebra, and they arise naturally.
[5]
[PDF] The orthogonal matrix and its applications
Orthogonal matrices are related to vectors and vector spaces, and are used in analytic geometry, including rotation of axes and reduction of quadratic forms.<|control11|><|separator|>
[6]
Orthogonal Matrix -- from Wolfram MathWorld
A n×n matrix A is an orthogonal matrix if AA^(T)=I, (1) where A^(T) is the transpose of A and I is the identity matrix. In particular, an orthogonal matrix ...
[7]
[PDF] Recitation 1: Mathematical Preliminaries 1.1 Linear Algebra
An orthogonal matrix is a real square matrix whose rows and columns are orthonormal vectors (orthogonal and unit norm). Formally, let Q ∈ Rn×n, then Q is ...
[8]
[PDF] Orthogonal matrices and Gram-Schmidt - MIT OpenCourseWare
A square orthonormal matrix Q is called an orthogonal matrix. If Q is square, then QTQ = I tells us that QT = Q-1. 0 0 1. 0 1 0. For example, if Q = 1 0 0.Missing: definition | Show results with:definition
[9]
[PDF] Linear Algebra 1: Matrices and Least Squares - MIT OpenCourseWare
Note that, unlike an orthogonal matrix, we do not require the matrix to be square. Just like orthogonal matrices, we have. QTQ = I , where I is an n × n ...
[10]
Unitary Matrix -- from Wolfram MathWorld
Unitary matrices leave the length of a complex vector unchanged. For real matrices, unitary is the same as orthogonal. In fact, there are some similarities ...
[11]
Lecture 17: Orthogonal matrices and Gram-Schmidt | Linear Algebra
These video lectures of Professor Gilbert Strang teaching 18.06 were ... 11:03So there's a -- there now I have got an orthogonal matrix,. 11:08in fact ...
[12]
[PDF] Lecture 4
All eigenvalues lie on the unit circle. Determinant has a magnitude that is one. |detA| = 1. For orthogonal matrix with real elements , this means that detA = ± ...
[13]
[PDF] 5. Orthogonal matrices
𝐴𝑥 is a permutation of the elements of 𝑥: 𝐴𝑥 = (𝑥𝜋1. ,𝑥𝜋2. ,...,𝑥𝜋𝑛). • 𝐴 has exactly one element equal to 1 in each row and each column.
[14]
[PDF] Rotations and reflections in the plane - Purdue Math
Recall that the transpose of an n ⇥ n matrix A is the n ⇥ n matrix with entries aji. A matrix is called orthogonal if AT A = I = AAT (the second equation is ...
[15]
[PDF] KEITH CONRAD - 1. Introduction An isometry of Rn is a function h ...
We return to orthogonal n × n matrices for any n ≥ 1. The geometric meaning of the condition A>A = In is that the columns of A are mutually perpendicular unit ...Missing: interpretation | Show results with:interpretation
[16]
[PDF] The Cartan–Dieudonné Theorem - CIS UPenn
In this chapter the structure of the orthogonal group is studied in more depth. In particular, we prove that every isometry in O(n) is the compo-.
[17]
Determinants and linear transformations - Math Insight
On the other hand, if det(A) is negative, the associated linear transformation reverses orientation by also reflecting the object (taking its mirror image).
[18]
Lecture 26 Orthogonal Matrices
The (counterclockwise) rotation matrices are the orthogonal matrices of determinant 1. The matrices corresponding to reflection followed by rotation are the ...<|separator|>
[19]
Householder matrix - StatLect
Learn how a Houselder matrix (or elementary reflector) is defined, constructed and used. With detailed explanations, proofs and solved exercises.
[20]
[PDF] Householder transformations - Cornell: Computer Science
Householder transformations are orthogonal reflections, expressed as H = I − 2vvT, used to introduce subdiagonal zeros in Gaussian elimination.
[21]
Givens rotation matrix - StatLect
Learn how a Givens rotation matrix is defined, constructed and used. With detailed explanations, proofs, examples and solved exercises.Definition · Orthogonality · Equivalent transformations · Annihilation of entries
[22]
Rotation Matrix -- from Wolfram MathWorld
Orthogonal matrices have special properties which allow them to be manipulated and identified with particular ease. Let A and B be two orthogonal matrices. By ...
[23]
Orthogonal Group -- from Wolfram MathWorld
The orthogonal group O(n) is the group of n×n orthogonal matrices. These matrices form a group because they are closed under multiplication and taking inverses.
[24]
Euler Angles -- from Wolfram MathWorld
153). If the coordinates of two sets of n points x_i and x_i^' are known, one rotated with respect to the other, then the Euler rotation matrix can be obtained ...
[25]
Quaternion -- from Wolfram MathWorld
### Summary: How Quaternions Represent 3D Rotations
[26]
Rodrigues' Rotation Formula -- from Wolfram MathWorld
Rodrigues' rotation formula gives an efficient method for computing the rotation matrix R in SO(3) corresponding to a rotation by an angle theta ...
[27]
Householder Matrix -- from Wolfram MathWorld
### Summary of Householder Matrix Formula in 3D (Using Normal Vector)
[28]
https://mathworld.wolfram.com/HouseholderMatrix.html
[29]
[PDF] Clifford Algebras and Spin Groups - Columbia Math Department
This group has two components, with the component of the identity SO(n, R), the orthogonal matrices of determinant 1. We'll mostly restrict attention to. SO(n, ...
[30]
The generation of all rational orthogonal matrices in R p,q
May 1, 2016 · We propose a method to find all rational orthogonal matrices in real indefinite inner product spaces, which is based on a proof of the Cartan–Dieudonné theorem.
[31]
[PDF] Similarity transforms Eigenvalue perturbations: a 2-by-2 illustration
Similarity transformations preserve a matrix's eigenvalue structure. Small perturbations to a matrix can split eigenvalues, and very small perturbations can ...
[32]
[PDF] §9.2 Orthogonal Matrices and Similarity Transformations
▷ PageRank vector x is eigenvector for G: G x = 1 · x, where 1 is always a simple eigenvalue of G. ▷ Power Method ...
[33]
[PDF] 3. Symmetric eigendecomposition
Symmetric eigendecomposition involves a symmetric matrix where all eigenvalues are real, and eigenvectors corresponding to different eigenvalues are orthogonal ...
[34]
[PDF] Eigenvalues - UC Berkeley math
the trance is invariant under similarity transformations. 235. Prove that tr A = − P λi, where λi are the roots of det(λI −A) = 0. Page 11. 1. The Spectral ...
[35]
[PDF] Lecture 2 Orthogonal Vectors and Matrices, Norms - DSpace@MIT
Sep 12, 2006 · Orthogonal Vectors and Matrices, Norms. MIT 18.335J / 6.337J ... • Therefore, lengths of vectors and angles between vectors are preserved.
[36]
[PDF] 9 Linear algebra
This is called a similarity transformation of the matrix B. To retain: • similarity transformations preserve traces and determinants: Tr M = Tr M0,. detM = ...
[37]
[PDF] Chapter 7 The Singular Value Decomposition (SVD)
I will first describe the SVD in terms of those basis vectors. Then I can also describe the SVD in terms of the orthogonal matrices U and V . (using vectors) ...
[38]
[PDF] Orthogonal Matrices and the Singular Value Decomposition
An m × n matrix A of rank r maps the r-dimensional unit hypersphere in rowspace(A) into an r-dimensional hyperellipse in range(A). Thus, a hypersphere is ...
[39]
[PDF] Lie Groups. Representation Theory and Symmetric Spaces
In this Chapter we discuss elementary properties of Lie groups, Lie algebras and their relationship. We will assume a good knowledge of manifolds, vector.
[40]
applications of the theory of matrices
The canonical forms are block diagonal, and each block is a 1 x 1 or a ... formation (Append I) then into real block diagonal form as shown. The product ...
[41]
[PDF] Canonical Forms (recall Jordan, why we wan - People @EECS
Theorem (Real Schur Canonical Form) Given any real square A, there is. a a real orthogonal Q such that Q*A*Q^T is block upper triangular with 1x1 and. 2x2 ...
[42]
[PDF] Further Properties of Involutory and Idempotent Matrices
Sep 17, 2025 · A is involutory if and only if A is similar to a diagonal matrix of the form diag(1,...,1,−1,...,−1). Proposition 2. Let A and B be real square ...
[43]
[PDF] The Exponential Map, Lie Groups, and Lie Algebras - UPenn CIS
The Lie algebra so(n,R) consisting of real skew symmet- ric n×n matrices is the corresponding set of infinitesimal rotations. The geometric link between a Lie ...
[44]
[PDF] Topics in Representation Theory: The Adjoint Representation 1 The ...
So, for any Lie group, we have a distinguished representation with dimension of the group, given by linear transformations on the Lie algebra. Later we will.
[45]
[PDF] Representations of the Rotation Groups SO(N)
Dec 14, 2019 · The adjoint representation is so useful because it is a characterization of how the generators of the Lie Algebra transform, and can be ...
[46]
[PDF] Clifford Algebras, Clifford Groups, and a Generalization ... - CIS UPenn
Jan 31, 2012 · Intuitively speaking,. SO(n) is more twisted than Spin(n). In fact, we will see that Spin(n) is a double cover of. SO(n). Since the spinor ...
[47]
[PDF] Topics in Representation Theory: The Spinor Representation
The fundamental irreducible representations that are missed by this con- struction are called the spinor representations. They are true representations of Spin( ...
[48]
https://www.math.columbia.edu/~woit/notes19.pdf
[49]
January 1928: The Dirac equation unifies quantum mechanics and ...
Nov 19, 2024 · A seminal paper by Paul Dirac, who relied on mathematical intuition, laid the foundation for quantum electrodynamics.
[50]
[PDF] QR decomposition: History and its Applications - Duke People
Dec 17, 2010 · There are several methods for computing the QR decomposition: • GS or modified GS,. • Givens rotations (real A), or. • Householder reflections.
[51]
Unitary Triangularization of a Nonsymmetric Matrix | Journal of the ...
A. S. HOUSEHOLDER (1958), A class of methods for inverting matrices. J ... Unitary Triangularization of a Nonsymmetric Matrix. Computing methodologies.
[52]
[PDF] Early History of the Singular Value Decomposition - UC Davis Math
Eugenio. Beltrami (1835-1899), Camille Jordan (1838-1921), ...
[53]
[PDF] Singular value decomposition and least squares solutions
To compute the singular value decomposition of a given matrix A, Forsythe and Henrici [2], Hestenes [8], and Kogbetliantz [9] proposed methods based on plane ...<|control11|><|separator|>
[54]
[PDF] Orthogonal Random Features - arXiv
Oct 28, 2016 · E(KORF(x, y)) = e−||x−y||2/2σ2 . 1We first generate the random Gaussian matrix G in (1). Q is the orthogonal matrix obtained from the QR.
[55]
Random matrices from QR - MathOverflow
Aug 11, 2019 · Let M=QR be the QR decomposition of M (Q is orthogonal, R is upper-triangular, with positive diagonal entries). Then Q is Haar-uniform.How to correctly generate uniformly distibuted random elements ...Intuition for Haar measure of random matrix - MathOverflowMore results from mathoverflow.net
[56]
[PDF] general solution of the Orthogonal Procrustes problem
The least-squares problem of transforming a given matrix A into a given matrix B by an orthogonal transformation matrix T so that the sums of squares of the ...Missing: seminal | Show results with:seminal
[57]
[PDF] Computing the Polar Decomposition—with Applications
been proposed for computing the orthogonal polar factor ofa nearly-orthogonal matrix. Bjorck and Bowie [7] derive a family of iterative methods with orders ...
[58]
[PDF] Random orthogonal matrices and the Cayley transform - arXiv
Oct 5, 2018 · The Cayley transform, as introduced in Cayley [5], is a map from skew-symmetric matrices to special orthogonal matrices. Given X in Skew(p) = {X ...
[59]
[PDF] On Orthogonalities in Matrices - arXiv
The paper discusses orthogonal, quasi-orthogonal, semi-orthogonal, and non-orthogonal matrices. An orthogonal matrix satisfies A'A = AA' = I.
[60]
OrthogonalMatrix: Construct an Orthogonal Matrix—Wolfram ...
A rectangular orthogonal matrix: Copy to clipboard. In[1]:=1. ✖. https ... Properties & Relations (1)Properties of the function, and connections to ...Missing: definition | Show results with:definition
[61]
[PDF] Golub and Van Loan - EE IIT Bombay
Theory of Matrices in Numerical Analysis, Blaisdell, New York. Reprinted in ... Over the complex field the unitary matrices correspond to the orthogonal matrices.
[62]
[PDF] A Tutorial on Principal Component Analysis
PCA assumes that all basis vectors {p1,..., pm} are orthonormal (i.e. pi · pj = δij). In the language of linear algebra, PCA assumes P is an orthonormal matrix.
[63]
[PDF] 8.7 Complex Matrices
Jul 8, 2020 · A square complex matrix U is called unitary if U−1 =UH. Thus a real matrix is unitary if and only if it is orthogonal. Example 8.7.5. The ...
[64]
Unitary matrices - Ximera - The Ohio State University
An complex matrix is unitary if , or equivalently if . Just as orthogonal matrices are exactly those that preserve the dot product, we have. A complex matrix is ...
[65]
[PDF] bilinear forms - keith conrad
Orthogonal bases for symmetric bilinear forms are the subject of Section 4. Symplectic bases for alternating bilinear forms are discussed in Section 5.
[66]
[PDF] Computing Generators of Groups preserving a Bilinear Form over ...
Aug 7, 2012 · Consider a matrix J ∈ Mn(Z) which is non-singular modulo p and for which JT = ±J. This matrix J then defines a bilinear form over the field Fp.
[67]
[PDF] Obstruction-theoretic definition of characteristic classes
Stiefel manifolds and their low-dimensional homotopy groups. The Stiefel manifold VkRn is the space of orthonormal k-frames at the origin in Rn. Such a frame ...
[68]
[PDF] Stiefel Manifolds and Polygons - Colorado State University
Pairs of vectors satisfying the above conditions are called orthonormal 2-frames for Rn, and the space of all such pairs is the Stiefel manifold St2(Rn).
[69]
[PDF] Lecture 3.26. Hermitian, unitary and normal matrices - Purdue Math
A Hermitian matrix satisfies A* = A, a unitary matrix satisfies U*U = I, and a normal matrix satisfies A*A = AA*.
[70]
[PDF] Unit 17: Spectral theorem
Orthogonal and unitary matrices are all normal. 17.2. Theorem: Symmetric matrices have only real eigenvalues. Proof. We extend the dot product to complex ...<|control11|><|separator|>