Linear Algebra
8.610 min read

Symmetric Matrices

A square matrix AA is symmetric when it equals its own transpose: A=AA^\top = A. Entry-by-entry, this means aij=ajia_{ij} = a_{ji} — the matrix is a mirror image across its main diagonal. Familiar examples include the identity matrix and any diagonal matrix.

The key insight: every quadratic function f(x)=xAxf(\mathbf{x}) = \mathbf{x}^\top A \mathbf{x} has a unique symmetric matrix representation. Given any matrix AA, you can replace it with 12(A+A)\frac{1}{2}(A + A^\top) without changing the quadratic form at all. This works because for a scalar s=xAxs = \mathbf{x}^\top A \mathbf{x}, we have s=s=xAxs = s^\top = \mathbf{x}^\top A^\top \mathbf{x}, so you can average the two representations.

The symmetric version has off-diagonal entry (i,j)(i,j) equal to the average of the original aija_{ij} and ajia_{ji}. For example, the non-symmetric matrix (320070401)\begin{pmatrix}3 & 2 & 0\\0 & 7 & 0\\4 & 0 & 1\end{pmatrix} defines the same quadratic as (312170201)\begin{pmatrix}3 & 1 & 2\\1 & 7 & 0\\2 & 0 & 1\end{pmatrix}. From now on, we always work with symmetric matrices for quadratic forms.

Formal View

Definition 8.5 — Symmetric Matrix
An n×nn \times n matrix AA is symmetric if A=AA^\top = A, or equivalently aij=ajia_{ij} = a_{ji} for all 1i,jn1 \leq i, j \leq n.
Lemma 8.1 — Unique Symmetric Representation
For any n×nn \times n matrix AA, we have xAx=x(A+A2)x\mathbf{x}^\top A \mathbf{x} = \mathbf{x}^\top \left(\frac{A + A^\top}{2}\right) \mathbf{x} for all xRn\mathbf{x} \in \mathbb{R}^n. The matrix 12(A+A)\frac{1}{2}(A + A^\top) is symmetric and is the unique symmetric matrix representing the same quadratic form.

Without loss of generality, the matrix of a quadratic form is always assumed to be symmetric.

Interactive Visualization

Transpose Visualizer

Why This Matters

Symmetric matrices arise naturally whenever a relationship between two things is mutual — covariance, stiffness, and Hessians are all symmetric.

  • Covariance matrices in statistics: Σij=Cov(Xi,Xj)=Cov(Xj,Xi)\Sigma_{ij} = \text{Cov}(X_i, X_j) = \text{Cov}(X_j, X_i).
  • Stiffness matrices in structural engineering: force-displacement relationships are mutual.
  • Hessian matrices of second derivatives are symmetric by Clairaut's theorem.

Quiz

Question 1

Every diagonal matrix is symmetric.

Question 2

What is the symmetric matrix associated with the quadratic form defined by A=(1402)A = \begin{pmatrix}1 & 4\\0 & 2\end{pmatrix}?

Common Mistakes

  • Thinking symmetry is just about the diagonal entries — it is about all off-diagonal pairs: aij=ajia_{ij} = a_{ji}.
  • Forgetting to symmetrize before computing eigenvalues — technically xAx\mathbf{x}^\top A \mathbf{x} works for any AA, but spectral theorem requires AA to be symmetric.