Visual Tools
Calculators
Tables
Mathematical Keyboard
Converters
Other Tools


Complex Eigenvalues






When No Real Direction Survives

A real matrix can have complex eigenvalues — this happens when the characteristic polynomial has no real roots. Geometrically, no real direction is mapped to a scalar multiple of itself; the transformation involves rotation. Complex eigenvalues of real matrices always come in conjugate pairs, and each pair corresponds to a rotation-scaling action on a two-dimensional subspace.

For step-by-step results on your own numbers, use the eigenvalue and eigenvector calculators on CalculateMatrix.



When Complex Eigenvalues Appear

An eigenvector is a direction that a matrix only stretches. Look for one in a rotation.
Take R=(0−110)R = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}, the rotation by 90°90°, and test a few vectors.
Every vector comes out at right angles to itself, so no real direction survives.
vRva multiple of v?
(1, 0)(0, 1)no
(1, 1)(−1, 1)no
(1, 2)(−2, 1)no
(a, b)(−b, a), at right angles to vnever, unless v = 0
RvRv is always perpendicular to vv, so it is never a multiple of vv. Yet the equation det⁡(R−λI)=λ2+1=0\det(R - \lambda I) = \lambda^2 + 1 = 0 still has solutions: λ=±i\lambda = \pm i.
The characteristic polynomial has no real roots, so the eigenvalues leave the real numbers.
For a 2×22 \times 2 matrix, the discriminant of the characteristic polynomial determines the eigenvalue type:
Discriminant Classification 2x2
Δ=tr(A)2−4det⁡(A){Δ>0:two distinct real eigenvaluesΔ=0:one repeated real eigenvalueΔ<0:complex conjugate pair\Delta = \text{tr}(A)^2 - 4\det(A) \quad \begin{cases} \Delta > 0: & \text{two distinct real eigenvalues} \\ \Delta = 0: & \text{one repeated real eigenvalue} \\ \Delta < 0: & \text{complex conjugate pair} \end{cases}
Complex eigenvalues appear in the Δ<0\Delta < 0 case: the quadratic formula produces λ=tr(A)±negative2\lambda = \frac{\text{tr}(A) \pm \sqrt{\text{negative}}}{2}, which involves −1=i\sqrt{-1} = i.
The definition
A complex eigenvalue of a real matrix AA is a non-real root λ\lambda of its characteristic polynomial; it comes with a complex eigenvector vv satisfying Av=λvAv = \lambda v.
A real matrix with real entries can have complex eigenvalues. This occurs when the characteristic polynomial — a polynomial with real coefficients — has roots that are not real.
So the 90°90° rotation has eigenvalues ±i\pm i, and v=(1,−i)v = (1, -i) satisfies Rv=ivRv = iv.
Complex eigenvalues are the algebra's signal of rotation.
The geometric interpretation is clear: no real direction in R2\mathbb{R}^2 is mapped to a scalar multiple of itself. Every vector is rotated, not just stretched or compressed. The simplest example is rotation by 90°90°, with matrix (0−110)\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} and characteristic polynomial λ2+1=0\lambda^2 + 1 = 0. The eigenvalues are ±i\pm i — purely imaginary, reflecting pure rotation with no scaling.
A2×20−110→eigenvalues λi−i
A rotation has no real eigenvector

The rotation preset: a real 2 × 2 matrix whose characteristic polynomial has a negative discriminant. Its roots are a complex conjugate pair, and (A − λI)v = 0 has no nonzero real solution — every real vector is turned off its own line, so no real direction is merely stretched. Compare it with the real cases on the eigenvalues and eigenvectors visualizer.

Complex eigenvalues are the algebra's way of reporting that the transformation turns every real direction.
They come in pairs: a real matrix never has just one.
The characteristic polynomial has real coefficients, so if λ=a+bi\lambda = a + bi is a root, so is λˉ=a−bi\bar{\lambda} = a - bi. For the rotation these are ii and −i-i. The pairing is the subject of conjugate pairs.
The rest of the page covers notation, conjugate pairs, the 2×2 case in detail, complex eigenvectors, the real canonical form, larger matrices, dynamical systems, the link with the fundamental theorem of algebra, and a summary.

Complex Eigenvalue Notation

Notation

Complex Eigenvalue Notation

A bar that has to be applied twice — once to a number and once to a whole vector — plus a ±\pm that is a statement about pairing rather than a shorthand for two options, and a normal form whose block is chosen so that no complex number has to be written down at all. All catalogued among the linear algebra symbols.
λ\lambda, the spectrum and the characteristic polynomial come from eigenvalue notation; the overbar and a+bia + bi from conjugate and complex number notation; Δ\Delta as discriminant from quadratic notation; P−1APP^{-1}AP from matrix representation notation.
Av=λv  ⇒  Avˉ=λˉvˉA\mathbf{v} = \lambda\mathbf{v} \;\Rightarrow\; A\bar{\mathbf{v}} = \bar{\lambda}\bar{\mathbf{v}}
conjugating both the eigenvalue and the eigenvector gives another eigenpair
The overbar works at two levels here: over λ\lambda it conjugates a single number, over v\mathbf{v} it conjugates every entry of a vector. The same mark, applied to objects of different type — and the implication only holds because AA is real, as Conjugate Pairs below proves.
CasesPhysics and much of numerical linear algebra prefer A∗A^{*} or AHA^{H} for the conjugate transpose of a matrix, reserving the plain bar for entrywise conjugation without the transpose; the two differ, and on a real matrix the distinction disappears, which is why it goes unnoticed until complex entries arrive.
Do not confuseConjugating only the eigenvalue. λˉ\bar{\lambda} pairs with vˉ\bar{\mathbf{v}}, not with v\mathbf{v} — the original eigenvector does not belong to the conjugate eigenvalue, and dropping the bar on the vector produces a statement that is simply false.
λ=a±bi\lambda = a \pm bi
a pair of eigenvalues, not a choice between two
Elsewhere ±\pm offers alternatives; here it asserts that both values occur, together, as a package — a real matrix cannot have one without the other. The compression is honest because the pair always appears in full, which is what makes the single line a complete description of two eigenvalues.
CasesThe sign convention runs through the geometry: bb is the rotation rate and its sign fixes the direction of turning, so writing a±bia \pm bi records a rotation together with its mirror image; b=0b = 0 collapses the pair to a single real eigenvalue and the notation degenerates gracefully.
Do not confuseThe ±\pm of the quadratic formula, where the two signs are two candidate roots you may pick between. Here selecting one is not an option — the eigenvalues are both present in the matrix.
(a−bba)\begin{pmatrix} a & -b \\ b & a \end{pmatrix}
the real canonical block standing in for the pair a plus or minus b i
A 2×22 \times 2 block that represents a conjugate pair without writing any complex number — the whole point of Real Canonical Form below. The off-diagonal sign asymmetry is the notation: −b-b above the diagonal and bb below is what encodes rotation.
CasesScaled out, the block reads r(cos⁡θ−sin⁡θsin⁡θcos⁡θ)r\begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta\end{pmatrix} with r=∣λ∣r = |\lambda| and θ=arg⁡λ\theta = \arg\lambda — the modulus and argument of the eigenvalue reappearing as a scaling factor and a rotation angle.
Do not confuseA diagonal matrix. The whole reason this form exists is that a real matrix with complex eigenvalues cannot be diagonalised over R\mathbb{R}; the block is the closest real substitute, and treating it as diagonal loses the rotation it was built to carry.
DownloadPNG2×PDF
Notation · Cheat sheet
Eigenvalue Notation at a Glance
The marks from the notation cards across the eigenvalue pages: how each is read and what it means. Every group links to the full card on its page.
Characteristic equation
Av=λvA\mathbf{v} = \lambda\mathbf{v}
A v equals lambda v
λ\lambda is the eigenvalue, v≠0\mathbf{v} \neq \mathbf{0} the eigenvector.
det⁡(A−λI)=0\det(A - \lambda I) = 0
the characteristic equation
The II is needed: a number cannot be taken from a matrix.
EλE_{\lambda}
the eigenspace for lambda
Null⁡(A−λI)\operatorname{Null}(A - \lambda I): the eigenvectors plus 0\mathbf{0}.
Diagonalization
A=PDP−1A = PDP^{-1}
P D P inverse
Eigenvectors as the columns of PP, eigenvalues on the diagonal of DD.
P:,j↔DjjP_{:,j} \leftrightarrow D_{jj}
matching order
Column jj of PP goes with entry jj of DD.
defective
too few eigenvectors
Cannot be diagonalised.
Properties
f(A)v=f(λ)vf(A)\mathbf{v} = f(\lambda)\mathbf{v}
same eigenvector, new eigenvalue
Powers, inverses and polynomials of AA keep v\mathbf{v}.
q(A)=c0I+⋯q(A) = c_0 I + \cdots
a polynomial evaluated at A
The constant term becomes c0Ic_0 I.
wTA=λwT\mathbf{w}^{T} A = \lambda \mathbf{w}^{T}
a left eigenvector
Multiplies from the left; plain eigenvector means the right one.
Complex eigenvaluesThis page
λ=a±bi\lambda = a \pm bi
a conjugate pair
Both values occur together for a real matrix.
Avˉ=λˉvˉA\bar{\mathbf{v}} = \bar{\lambda}\bar{\mathbf{v}}
conjugate eigenpair
Conjugating both gives another eigenpair.
(a−bba)\begin{pmatrix} a & -b \\ b & a \end{pmatrix}
the real canonical block
Stands for the pair a±bia \pm bi with real entries.
Same marks, different job
±\pmBoth values together for eigenvalues; elsewhere, a choice between two.
λ\lambda vs λI\lambda IA number against the matrix it becomes inside A−λIA - \lambda I.
vˉ\bar{\mathbf{v}}Every entry conjugated, not a mean.
Eigenvalue Notation at a Glance·/linear-algebra/eigen/characteristic-equationLearn Math Class

Conjugate Pairs

For a real matrix, complex eigenvalues always come in conjugate pairs. If λ=a+bi\lambda = a + bi is an eigenvalue (with b≠0b \neq 0), then λˉ=a−bi\bar{\lambda} = a - bi is also an eigenvalue. The corresponding eigenvectors are also conjugates:

Complex Conjugate Pairs
A∈Rn×n,    Av=λv⇒Avˉ=λˉvˉA \in \mathbb{R}^{n \times n}, \;\; A\mathbf{v} = \lambda\mathbf{v} \Rightarrow A\bar{\mathbf{v}} = \bar{\lambda}\bar{\mathbf{v}}

The proof uses the fact that the characteristic polynomial has real coefficients. If p(λ)=0p(\lambda) = 0 and every coefficient of pp is real, then p(λˉ)=p(λ)‾=0‾=0p(\bar{\lambda}) = \overline{p(\lambda)} = \overline{0} = 0. Complex conjugation passes through the polynomial because the coefficients are their own conjugates.

For the eigenvector identity: applying complex conjugation to Av=λvA\mathbf{v} = \lambda\mathbf{v} and using that AA has real entries gives Avˉ=λˉvˉA\bar{\mathbf{v}} = \bar{\lambda}\bar{\mathbf{v}}.

One consequence for odd-dimensional real matrices: at least one eigenvalue must be real. Complex roots pair up, consuming an even number of the nn roots. When nn is odd, at least one root is left unpaired, and an unpaired root of a real polynomial must be real.
Real entries force complex eigenvalues, and their eigenvectors, to arrive in conjugate pairs.

The 2×2 Case in Detail

Fact: No crystal has fivefold symmetry
A rotation that maps a plane lattice onto itself has integer entries in a basis of lattice vectors, so its trace 2cos⁡θ2\cos\theta is an integer (in space, 1+2cos⁡θ1 + 2\cos\theta). For a complex pair that leaves only θ=60°\theta = 60°, 90°90° or 120°120°: crystals allow 6-, 4- and 3-fold turns, never the 72°72° turn of 5-fold symmetry.

A 2×22 \times 2 real matrix with complex eigenvalues λ=a±bi\lambda = a \pm bi (where b≠0b \neq 0) acts, in a suitable real basis, as a rotation composed with a scaling. The eigenvalues have modulus r=a2+b2r = \sqrt{a^2 + b^2} and argument θ\theta with cos⁡θ=a/r\cos\theta = a/r, sin⁡θ=b/r\sin\theta = b/r (this is arctan⁡(b/a)\arctan(b/a) only when a>0a > 0).

When r=1r = 1 (equivalently, det⁡(A)=a2+b2=1\det(A) = a^2 + b^2 = 1), the transformation is, in that basis, a pure rotation by angle θ\theta. When r>1r > 1, it is a rotation with outward spiraling. When r<1r < 1, it is a rotation with inward spiraling.

Worked Example


For A=(3−211)A = \begin{pmatrix} 3 & -2 \\ 1 & 1 \end{pmatrix}: tr(A)=4\text{tr}(A) = 4, det⁡(A)=5\det(A) = 5, discriminant =16−20=−4= 16 - 20 = -4. The eigenvalues are

The 2×2 Case in Detail
λ=4±−42=2±i\lambda = \frac{4 \pm \sqrt{-4}}{2} = 2 \pm i

Pitfall: Taking the determinant as the modulus
Reading det⁡(A)=5\det(A) = 5 as the modulus says distances are scaled by 55. But det⁡(A)=a2+b2=r2\det(A) = a^2 + b^2 = r^2, so r=5≈2.24r = \sqrt{5} \approx 2.24; the factor 55 belongs to areas.

The modulus is r=4+1=5r = \sqrt{4 + 1} = \sqrt{5} and the argument is θ=arctan⁡(1/2)\theta = \arctan(1/2). In the basis of the real canonical form, the transformation rotates by arctan⁡(1/2)≈26.6°\arctan(1/2) \approx 26.6° while scaling distances by 5≈2.24\sqrt{5} \approx 2.24. Since r>1r > 1, repeated application spirals outward.
DownloadPNG2×PDF
Eigen · complex pairs

Three regimes, set by the modulus

The argument θ\theta fixes how fast the plane turns and never affects whether it grows. That is decided entirely by the modulus, and the three cases below are the only ones there are.

3regimes
Neither growing nor decaying1
geometric actionpure rotation by θ\theta in the basis of the real canonical form, no scaling
det of the 2×2 blockdet⁡A=1\det A = 1
orbit of a pointpoints on a closed ellipse (a circle in the basis of the real canonical form)
In the basis of the real canonical form, lengths are preserved exactly and the block is orthogonal; AA itself need not be. Every orbit closes if θ\theta is a rational multiple of 2π2\pi and fills the circle densely otherwise — the boundary case in every stability argument.
Growing or decaying2
geometric actionrotation with outward scaling by rr
det of the 2×2 blockdet⁡A>1\det A > 1
orbit of a pointa spiral outward
Each application turns by θ\theta and stretches by rr, so iterates diverge while circling. In a discrete system this is instability — and the divergence is geometric, so it is fast.
geometric actionrotation with inward scaling by rr
det of the 2×2 blockdet⁡A<1\det A < 1
orbit of a pointa spiral inward to the origin
The stable case for iteration: every orbit converges to the origin, still turning as it goes. This is what ∣λ∣<1|\lambda| < 1 means geometrically, and why the condition is stated on the modulus rather than on the real part.
The determinant of the 2×22 \times 2 block is a2+b2=r2a^2 + b^2 = r^2, so the area scaling and the modulus are the same fact stated twice — which is why the determinant can be read as a stability marker without computing an eigenvalue at all.
Three regimes, set by the modulus·/linear-algebra/eigen/complexLearn Math Class
The three cases exhaust the possibilities because rr is a positive real and can only be less than, equal to, or greater than one. The argument θ\theta plays no part in the classification — it sets the speed of rotation and nothing else, which is why two systems can spiral at completely different rates and be equally stable.
Re vIm v12345678A⁸x₀x₀
Every direction turned, none left fixed

No arrow survives this matrix pointing where it started: every direction is rotated. That is what a complex eigenvalue looks like from the real plane — not an exotic number but the absence of any invariant line. The rotation angle is the argument of the eigenvalue. Turn the matrix and watch the fixed directions vanish on the complex eigenvalue explorer.

The conjugate partner of this eigenvalue describes the same rotation turning the other way, which is why the pair always arrives together.

Complex Eigenvectors

The eigenvectors corresponding to complex eigenvalues have complex entries. To find them, solve (A−λI)v=0(A - \lambda I)\mathbf{v} = \mathbf{0} using complex arithmetic.

Continuing the example with A=(3−211)A = \begin{pmatrix} 3 & -2 \\ 1 & 1 \end{pmatrix} and λ=2+i\lambda = 2 + i:

Complex Eigenvectors
A−(2+i)I=(1−i−21−1−i)A - (2 + i)I = \begin{pmatrix} 1 - i & -2 \\ 1 & -1 - i \end{pmatrix}

The second row gives v1+(−1−i)v2=0v_1 + (-1 - i)v_2 = 0, so v1=(1+i)v2v_1 = (1 + i)v_2. Setting v2=1v_2 = 1: v=(1+i,1)T\mathbf{v} = (1 + i, 1)^T.

This eigenvector splits into real and imaginary parts: v=(1,1)T⏟u+i(1,0)T⏟w\mathbf{v} = \underbrace{(1, 1)^T}_{\mathbf{u}} + i\underbrace{(1, 0)^T}_{\mathbf{w}}. The real vectors u\mathbf{u} and w\mathbf{w} encode the rotation — they span the two-dimensional subspace on which the rotation-scaling acts. The conjugate eigenvector for λˉ=2−i\bar{\lambda} = 2 - i is vˉ=(1−i,1)T\bar{\mathbf{v}} = (1 - i, 1)^T.
A complex eigenvector is two real vectors in disguise, spanning the plane on which the matrix rotates and scales.

Real Canonical Form

A 2×22 \times 2 real matrix with eigenvalues a±bia \pm bi is similar (over R\mathbb{R}) to a rotation-scaling block:

Real Canonical Form 2x2
λ=a±bi⇒P−1AP=(a−bba)=r(cos⁡θ−sin⁡θsin⁡θcos⁡θ)\lambda = a \pm bi \Rightarrow P^{-1}AP = \begin{pmatrix} a & -b \\ b & a \end{pmatrix} = r\begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}

where r=a2+b2r = \sqrt{a^2 + b^2} and θ\theta satisfies cos⁡θ=a/r\cos\theta = a/r, sin⁡θ=b/r\sin\theta = b/r. This is a rotation by θ\theta scaled by rr.

The similarity is achieved by the real matrix P=[w  u]P = [\mathbf{w} \; \mathbf{u}] — imaginary part first, then real part — where v=u+iw\mathbf{v} = \mathbf{u} + i\mathbf{w} is the complex eigenvector for λ=a+bi\lambda = a + bi. Then P−1AP=(a−bba)P^{-1}AP = \begin{pmatrix} a & -b \\ b & a \end{pmatrix}.

Pitfall: Putting the real part first
Taking P=[u  w]=(1110)P = [\mathbf{u} \; \mathbf{w}] = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} gives P−1AP=(21−12)P^{-1}AP = \begin{pmatrix} 2 & 1 \\ -1 & 2 \end{pmatrix}, a rotation the other way. With P=[w  u]=(1101)P = [\mathbf{w} \; \mathbf{u}] = \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix} the result is (2−112)\begin{pmatrix} 2 & -1 \\ 1 & 2 \end{pmatrix}.

For the running example, λ=2+i\lambda = 2 + i and v=(1+i,1)T\mathbf{v} = (1 + i, 1)^T from the previous section, so u=(1,1)T\mathbf{u} = (1, 1)^T and w=(1,0)T\mathbf{w} = (1, 0)^T. These two real vectors become the columns of PP, and P−1APP^{-1}AP should come out as the block with a=2a = 2 and b=1b = 1. Computing it is a direct check on the eigenvector found above.

This is the real alternative to diagonalization: instead of a diagonal matrix with complex entries (which is valid over C\mathbb{C} but not over R\mathbb{R}), we get a real 2×22 \times 2 rotation-scaling block. The transformation is expressed in its simplest real form.
The check the section sets up: the columns w and u of P turn A into its rotation-scaling block.
DownloadPNG2×PDF
Real canonical form

A rotation-scaling block

Instead of a diagonal matrix with complex entries, a real 2 × 2 block: a rotation by θ scaled by r.
Before
Standard basis
A=(3−211)A = \begin{pmatrix} 3 & -2 \\ 1 & 1 \end{pmatrix}
Eigenvalues 2 ± i
Change basis
P = [w u]: w = (1, 0), then u = (1, 1)
After
Basis w, u
(2−112)\begin{pmatrix} 2 & -1 \\ 1 & 2 \end{pmatrix}
a = 2, b = 1, so r = √5, θ ≈ 26.6°
Taking P=[u  w]=(1110)P = [\mathbf{u} \; \mathbf{w}] = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} gives P−1AP=(21−12)P^{-1}AP = \begin{pmatrix} 2 & 1 \\ -1 & 2 \end{pmatrix}, a rotation the other way. With P=[w  u]=(1101)P = [\mathbf{w} \; \mathbf{u}] = \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix} the result is (2−112)\begin{pmatrix} 2 & -1 \\ 1 & 2 \end{pmatrix}.
What changed
The basis: in the basis w, u the map reads as a rotation by θ ≈ 26.6° scaled by r = √5.
What stayed
The map and its eigenvalues 2 ± i; every entry stays real.
The running example in real canonical form·/linear-algebra/eigen/complexLearn Math Class
P⁻¹AP comes out as the block with a = 2 and b = 1, which confirms the eigenvector (1 + i, 1) found above.

Complex Eigenvalues in Larger Matrices

For an n×nn \times n real matrix, complex eigenvalues still come in conjugate pairs. When AA is diagonalizable over C\mathbb{C}, each pair contributes a 2×22 \times 2 rotation-scaling block to the real canonical form, and each real eigenvalue contributes a 1×11 \times 1 block (just the eigenvalue itself).

A 4×44 \times 4 matrix with eigenvalues 2±3i2 \pm 3i and −1±i-1 \pm i has real canonical form

Complex Eigenvalues in Larger Matrices
(2−300320000−1−1001−1)\begin{pmatrix} 2 & -3 & 0 & 0 \\ 3 & 2 & 0 & 0 \\ 0 & 0 & -1 & -1 \\ 0 & 0 & 1 & -1 \end{pmatrix}

A 3×33 \times 3 matrix with eigenvalues 55 and 1±2i1 \pm 2i has real canonical form

Complex Eigenvalues in Larger Matrices
(50001−2021)\begin{pmatrix} 5 & 0 & 0 \\ 0 & 1 & -2 \\ 0 & 2 & 1 \end{pmatrix}

Fact: Only 8 of 100 eigenvalues are real
A 100×100100 \times 100 matrix with independent, normally distributed entries has on average only about 8.48.4 real eigenvalues; the other 9090 or so pair off into 2×22 \times 2 blocks. The expected number of real eigenvalues grows only like 2n/π\sqrt{2n/\pi} (Edelman, Kostlan and Shub, 1994).

The real canonical form is the real analogue of diagonalization: for a matrix diagonalizable over C\mathbb{C} it achieves the simplest possible real matrix that is similar to AA, with the block structure directly reflecting the eigenvalue spectrum; otherwise larger real Jordan blocks appear.

How the nn eigenvalues split between the two kinds of block depends on the matrix. A symmetric matrix has only real eigenvalues, so its real canonical form is made of 1×11 \times 1 blocks alone, while a rotation of R4\mathbb{R}^4 in two separate planes is made of two 2×22 \times 2 blocks and nothing else. A matrix whose entries are chosen at random falls somewhere between these extremes.

Dynamical Systems Interpretation

Complex eigenvalues produce oscillatory behavior in dynamical systems.

In the discrete system xn+1=Axn\mathbf{x}_{n+1} = A\mathbf{x}_n, an eigenvalue λ=a+bi\lambda = a + bi with modulus r=∣λ∣=a2+b2r = |\lambda| = \sqrt{a^2 + b^2} generates a spiraling trajectory. When r>1r > 1, the spiral grows outward — the system is unstable. When r<1r < 1, the spiral decays inward — the system converges to the origin. When r=1r = 1, the points stay on a closed curve — perpetual oscillation without growth or decay.

Pitfall: Using the real-part test on a discrete system
Judging xn+1=Axn\mathbf{x}_{n+1} = A\mathbf{x}_n by the sign of the real part calls λ=12±12i\lambda = \frac{1}{2} \pm \frac{1}{2}i unstable, since Re⁡(λ)>0\operatorname{Re}(\lambda) > 0. Here the modulus decides: ∣λ∣=22≈0.71<1|\lambda| = \frac{\sqrt{2}}{2} \approx 0.71 < 1, so the iterates spiral in to 0\mathbf{0}.

In the continuous system x′=Ax\mathbf{x}' = A\mathbf{x}, an eigenvalue λ=a+bi\lambda = a + bi contributes a term eat(cos⁡bt,sin⁡bt)e^{at}(\cos bt, \sin bt) to the solution. The real part aa determines exponential growth (a>0a > 0) or decay (a<0a < 0). The imaginary part bb determines the oscillation frequency. When a=0a = 0, the oscillation is sustained. When a<0a < 0, it is damped. When a>0a > 0, it grows without bound.

THEOREM
Stability of a linear system reduces to eigenvalue analysis: the system is stable if and only if every eigenvalue has negative real part (continuous) or modulus less than 11 (discrete).
Regime Discrete: xn+1 = A xn Continuous: x′ = A x
Stable (converging) |λ| < 1 → trajectory spirals inward to the origin Re(λ) < 0 → damped oscillation decaying to zero
Marginally stable |λ| = 1 → trajectory stays on a closed curve Re(λ) = 0 → sustained oscillation at constant amplitude
Unstable (diverging) |λ| > 1 → trajectory spirals outward without bound Re(λ) > 0 → oscillation with exponentially growing amplitude
Re vIm v123456789A¹²x₀x₀Re vIm v12345678A⁸x₀x₀
Spiralling in, then spiralling out

The rotation is the same in both; what differs is the modulus of the eigenvalue. Below one, repeated application pulls everything towards the origin; above one, it throws everything outward. The argument sets how fast it turns and the modulus sets whether it survives, which is the whole stability story for a discrete system. Tune both on the complex eigenvalue explorer.

Whether a system settles, oscillates forever or blows up is therefore read off a single number: the modulus.

The Fundamental Theorem of Algebra Connection

THEOREMFundamental Theorem of Algebra
Over C\mathbb{C}, every polynomial of degree nn factors completely into nn linear factors.

This means every n×nn \times n matrix — real or complex — has exactly nn eigenvalues counted with algebraic multiplicity when the characteristic polynomial is factored over C\mathbb{C}.

Over R\mathbb{R}, some factors may be irreducible quadratics rather than linear, corresponding to conjugate pairs of complex eigenvalues. The characteristic polynomial of a real matrix factors into linear terms (real eigenvalues) and irreducible quadratic terms (complex conjugate pairs).

Working over C\mathbb{C} simplifies the theory considerably. Every matrix is triangularizable over C\mathbb{C} (Schur decomposition): A=UTU∗A = UTU^* where UU is unitary and TT is upper triangular with eigenvalues on the diagonal. Diagonalizability depends only on whether the geometric multiplicities match the algebraic multiplicities, with no additional complications from irreducible quadratics.

The choice between working over R\mathbb{R} and C\mathbb{C} is a recurring theme. Real matrices are the natural objects of computation, but complex eigenvalues are the natural objects of spectral theory. Both perspectives are needed.

Summary: Complex Eigenvalues at a Glance

The behavior of complex eigenvalues for a real matrix is governed by a handful of formulas and statements that connect the discriminant, the modulus, the argument, and the real part to specific geometric and dynamical consequences. The table below collects each one alongside what it produces and what it means, providing a single reference for the page's key facts.
A complex eigenvalue is often met as a nuisance — the characteristic polynomial refuses to factor over the reals and the calculation appears to have failed. It has not. The pair is describing a rotation, and once that is granted the two numbers it carries are exactly the two things a rotation has: an angle and a scale. The grouping below is recognition first, then those two numbers, then what they predict.
DownloadPNG2×PDF
Eigen · complex eigenvalues

Reading a conjugate pair

A complex pair is not a failure of the real case — it is a rotation, and its two parts say by how much and through what angle. The modulus and argument are the whole content, and everything below is one of them applied.

8readings
Recognising the pair3
1
2×22 \times 2 real matrices
tr⁡(A)2−4det⁡(A)<0\operatorname{tr}(A)^2 - 4\det(A) < 0
The discriminant settles it before any solving. A negative value means no real direction is mapped to a multiple of itself — the transformation turns every vector, so no real eigenvector can exist.
2
AA real — not for complex entries
λ=a+bi⇒λˉ=a−bi\lambda = a + bi \Rightarrow \bar{\lambda} = a - bi
Complex roots of a real polynomial arrive in pairs, and the eigenvectors pair with them: v\mathbf{v} and vˉ\bar{\mathbf{v}}. Nothing has to be computed twice — finding one half of the pair gives the other by conjugation.
3
real matrices only
nn odd ⇒\Rightarrow at least one real λ\lambda
Complex roots come in pairs, so an odd number of them cannot pair off entirely and at least one must be real. This is why a rotation in R3\mathbb{R}^3 has an axis while a rotation in R2\mathbb{R}^2 has none — the eigenvector of the leftover real eigenvalue spans the axis.
Reading the pair3
4
the scaling factor
r=a2+b2r = \sqrt{a^2 + b^2}
How much the transformation stretches the plane per application. Since r2=det⁡r^2 = \det of the 2×22 \times 2 block, the determinant is the area scaling and the modulus is its square root — consistent with the geometric reading.
5
the rotation angle
θ=arg⁡(a+bi)\theta = \arg(a + bi)
How far the plane turns per application. A pair with a=0a = 0 is a quarter turn; the sign of bb picks the direction, and the two conjugates give the two orientations of the same rotation.
6
similar over R\mathbb{R}, no complex arithmetic
[a−bba]=r⋅R(θ)\begin{bmatrix} a & -b \\ b & a \end{bmatrix} = r \cdot R(\theta)
The pair does not have to be handled in C\mathbb{C} at all. Over R\mathbb{R} the matrix is similar to a scaled rotation block, so the whole analysis stays real — which is what makes complex eigenvalues usable in applied work.
What the pair predicts2
7
for xk+1=Axk\mathbf{x}_{k+1} = A\mathbf{x}_k
∣λ∣<1|\lambda| < 1 for every eigenvalue
The modulus alone decides whether iteration spirals inward or outward; the argument only sets how fast it turns. Note ∣λ∣=1|\lambda| = 1 is the boundary case — in the basis of the real canonical form, a pure rotation that neither grows nor decays.
8
for x′=Ax\mathbf{x}' = A\mathbf{x}
Re⁡(λ)<0\operatorname{Re}(\lambda) < 0 for every eigenvalue
A different criterion for a different system, and confusing the two is the usual error. Discrete systems compare the modulus against one; continuous systems compare the real part against zero — because the solution is eλte^{\lambda t} and its growth is governed by Re⁡(λ)\operatorname{Re}(\lambda).
The single idea underneath: a real matrix with eigenvalues re±iθre^{\pm i\theta} acts on the corresponding real plane, in the basis of the real canonical form, as a rotation by θ\theta scaled by rr. Recognition, canonical form and stability are all that statement read at different moments — which is why the modulus and argument are worth computing before anything else.
Reading a conjugate pair·/linear-algebra/eigen/complexLearn Math Class
The real canonical form is what keeps this practical. A conjugate pair can be handled entirely in R\mathbb{R}: the matrix is similar to a block [a−bba]\begin{bmatrix} a & -b \\ b & a \end{bmatrix}, which is a rotation by θ\theta scaled by rr, and no complex arithmetic is required at any point. The complex numbers are a route to the answer rather than a property of it.
The two stability criteria are worth separating carefully, because they look alike and are not. A discrete system xk+1=Axk\mathbf{x}_{k+1} = A\mathbf{x}_k is stable when every ∣λ∣<1|\lambda| < 1 — the modulus decides, since each step multiplies by it. A continuous system x′=Ax\mathbf{x}' = A\mathbf{x} is stable when every Re⁡(λ)<0\operatorname{Re}(\lambda) < 0 — the real part decides, since the solution behaves as eλte^{\lambda t}. Applying the wrong test gives a confidently wrong answer.

Common Pitfalls

Four mistakes come up again and again when working with complex eigenvalues. Three are marked with a box where they happen in the sections above, and the first, the formula for the angle, is flagged in the text of the 2×2 section. Here each is worked through in full, with the wrong and the right working side by side.

Reading the angle from the arctangent alone

The argument θ\theta of a+bia + bi is computed as arctan⁡(b/a)\arctan(b/a) whatever the sign of aa. The arctangent only returns angles between −90°-90° and 90°90°, so when a<0a < 0 the angle lands in the wrong half of the plane.

Wrong working: for the pair −1±i-1 \pm i of the 4×44 \times 4 example, r=2r = \sqrt{2} and
θ=arctan⁡(1−1)=−45°\theta = \arctan\left(\frac{1}{-1}\right) = -45°

2(cos⁡(−45°)−sin⁡(−45°)sin⁡(−45°)cos⁡(−45°))=(11−11)\sqrt{2}\begin{pmatrix} \cos(-45°) & -\sin(-45°) \\ \sin(-45°) & \cos(-45°) \end{pmatrix} = \begin{pmatrix} 1 & 1 \\ -1 & 1 \end{pmatrix}

which is not the block (−1−11−1)\begin{pmatrix} -1 & -1 \\ 1 & -1 \end{pmatrix} of that pair.

Right working: use cos⁡θ=a/r\cos\theta = a/r and sin⁡θ=b/r\sin\theta = b/r together.
cos⁡θ=−12,sin⁡θ=12⇒θ=135°\cos\theta = -\frac{1}{\sqrt{2}}, \quad \sin\theta = \frac{1}{\sqrt{2}} \quad\Rightarrow\quad \theta = 135°

2(cos⁡135°−sin⁡135°sin⁡135°cos⁡135°)=(−1−11−1)\sqrt{2}\begin{pmatrix} \cos 135° & -\sin 135° \\ \sin 135° & \cos 135° \end{pmatrix} = \begin{pmatrix} -1 & -1 \\ 1 & -1 \end{pmatrix}


Why it happens: most worked examples, including 2±i2 \pm i on this page, have a>0a > 0, where arctan⁡(b/a)\arctan(b/a) is correct, so the shortcut is remembered as the rule.

How to check: the signs of cos⁡θ\cos\theta and sin⁡θ\sin\theta must match the signs of aa and bb. Here a<0a < 0 and b>0b > 0, so θ\theta lies between 90°90° and 180°180°.

Taking the determinant as the modulus

The determinant det⁡(A)=a2+b2\det(A) = a^2 + b^2 is read as the modulus rr of the eigenvalues, so the scaling factor comes out squared.

Wrong working: for A=(3−211)A = \begin{pmatrix} 3 & -2 \\ 1 & 1 \end{pmatrix},
det⁡(A)=3⋅1−(−2)⋅1=5⇒r=5\det(A) = 3 \cdot 1 - (-2) \cdot 1 = 5 \quad\Rightarrow\quad r = 5

and the rotation-scaling is said to stretch distances by 55.

Right working: the eigenvalues are 2±i2 \pm i, and
r=∣2+i∣=22+12=5≈2.24,det⁡(A)=(2+i)(2−i)=5=r2r = |2 + i| = \sqrt{2^2 + 1^2} = \sqrt{5} \approx 2.24, \qquad \det(A) = (2 + i)(2 - i) = 5 = r^2

In the basis of the real canonical form, distances are scaled by 5\sqrt{5}; the factor 55 is how much areas grow.

Why it happens: the determinant is often described as the scaling factor of a transformation. It is the factor for areas, and in the plane that is the square of the factor for lengths.

How to check: for a real 2×22 \times 2 matrix with a complex pair, r=det⁡(A)r = \sqrt{\det(A)}. The verdict on growth is the same either way, since r>1r > 1 exactly when det⁡(A)>1\det(A) > 1, so only the rate gives the slip away.

Putting the real part first

The matrix PP is built with the real part u\mathbf{u} of the eigenvector as its first column and the imaginary part w\mathbf{w} second. The result is still a rotation-scaling block, but the sign of bb is flipped, so it turns the other way.

Wrong working: with v=(1+i,1)T\mathbf{v} = (1 + i, 1)^T, u=(1,1)T\mathbf{u} = (1, 1)^T and w=(1,0)T\mathbf{w} = (1, 0)^T,
P=[u  w]=(1110),P−1=(011−1)P = [\mathbf{u} \; \mathbf{w}] = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}, \qquad P^{-1} = \begin{pmatrix} 0 & 1 \\ 1 & -1 \end{pmatrix}

P−1AP=(011−1)(1321)=(21−12)P^{-1}AP = \begin{pmatrix} 0 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 1 & 3 \\ 2 & 1 \end{pmatrix} = \begin{pmatrix} 2 & 1 \\ -1 & 2 \end{pmatrix}

the block for a=2a = 2 and b=−1b = -1, a rotation by −θ-\theta.

Right working: imaginary part first.
P=[w  u]=(1101),P−1=(1−101)P = [\mathbf{w} \; \mathbf{u}] = \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix}, \qquad P^{-1} = \begin{pmatrix} 1 & -1 \\ 0 & 1 \end{pmatrix}

P−1AP=(1−101)(3112)=(2−112)P^{-1}AP = \begin{pmatrix} 1 & -1 \\ 0 & 1 \end{pmatrix}\begin{pmatrix} 3 & 1 \\ 1 & 2 \end{pmatrix} = \begin{pmatrix} 2 & -1 \\ 1 & 2 \end{pmatrix}


Why it happens: v=u+iw\mathbf{v} = \mathbf{u} + i\mathbf{w} lists the real part first, and that order is carried into PP. Some textbooks do use [u  w][\mathbf{u} \; \mathbf{w}], but pair it with the eigenvalue a−bia - bi, and the two conventions get mixed.

How to check: AwA\mathbf{w} must equal aw+bua\mathbf{w} + b\mathbf{u}. Here Aw=(3,1)T=2(1,0)T+1⋅(1,1)TA\mathbf{w} = (3, 1)^T = 2(1, 0)^T + 1 \cdot (1, 1)^T, so w\mathbf{w} belongs in the first column, and the entry below the diagonal of P−1APP^{-1}AP must be b=1b = 1.

Using the real-part test on a discrete system

The stability of xn+1=Axn\mathbf{x}_{n+1} = A\mathbf{x}_n is judged by the sign of Re⁡(λ)\operatorname{Re}(\lambda), the test for continuous systems, or the other way round: x′=Ax\mathbf{x}' = A\mathbf{x} is judged by ∣λ∣|\lambda|.

Wrong working: A=12(1−111)A = \frac{1}{2}\begin{pmatrix} 1 & -1 \\ 1 & 1 \end{pmatrix} has tr(A)=1\text{tr}(A) = 1 and det⁡(A)=12\det(A) = \frac{1}{2}, so
λ=1±1−22=12±12i,Re⁡(λ)=12>0\lambda = \frac{1 \pm \sqrt{1 - 2}}{2} = \frac{1}{2} \pm \frac{1}{2}i, \qquad \operatorname{Re}(\lambda) = \frac{1}{2} > 0

and the discrete system is declared unstable.

Right working: for a discrete system the modulus decides.
∣λ∣=14+14=22≈0.71<1|\lambda| = \sqrt{\frac{1}{4} + \frac{1}{4}} = \frac{\sqrt{2}}{2} \approx 0.71 < 1

Starting from x0=(1,0)T\mathbf{x}_0 = (1, 0)^T, the iterates (12,12)T(\frac{1}{2}, \frac{1}{2})^T and (0,12)T(0, \frac{1}{2})^T have lengths 0.710.71 and 0.50.5: the trajectory spirals in to 0\mathbf{0}. For x′=Ax\mathbf{x}' = A\mathbf{x} with the same matrix the verdict reverses: ∣λ∣<1|\lambda| < 1, but Re⁡(λ)>0\operatorname{Re}(\lambda) > 0, and the solutions grow like et/2e^{t/2}.

Why it happens: both tests compare an eigenvalue with a boundary, the unit circle for discrete systems and the imaginary axis for continuous ones, and the two boundaries are easily swapped.

How to check: ask what the system does to the part of x\mathbf{x} along an eigenvector. A discrete step multiplies it by λ\lambda, so after nn steps its size is ∣λ∣n|\lambda|^n; a continuous solution behaves like eλte^{\lambda t}, whose size is eate^{at}.

Complex Eigenvalues FAQ

Does λˉ\bar{\lambda} pair with v\mathbf{v} or with vˉ\bar{\mathbf{v}}?

+
With vˉ\bar{\mathbf{v}}. Conjugating the whole eigenvalue equation conjugates the vector as well, so the conjugate eigenvalue belongs to the conjugate eigenvector. The original v\mathbf{v} does not satisfy the equation for λˉ\bar{\lambda}, and writing it that way produces a statement that is simply false rather than merely imprecise.Read more →

Is the ±\pm in a±bia \pm bi like the one in the quadratic formula?

+
It looks identical but does a different job. In the quadratic formula the two signs offer two candidate roots and you pick whichever the problem needs. Here selecting one is not an option: both eigenvalues genuinely belong to the matrix at the same time, and the mark is recording a pair rather than a choice.Read more →

Is the real canonical block just a diagonal matrix?

+
No, and that is exactly why it exists. A real matrix with complex eigenvalues cannot be diagonalised over the reals at all, so the block is the closest real substitute. Written out, it is a scaled rotation, with the eigenvalue's modulus setting the scale and its argument the angle. Treating it as diagonal discards the rotation.Read more →

Does every real matrix have at least one real eigenvalue?

+
No, but every real matrix of odd size does. Complex eigenvalues of a real matrix come in conjugate pairs, so they use up an even number of the nn roots; when nn is odd one root is left over, and it must be real. That is why a rotation of R3\mathbb{R}^3 always has an axis, while the 90°90° rotation of the plane has eigenvalues ±i\pm i and no real eigenvector at all.Read more →

How do you find the rotation angle of a complex eigenvalue?

+
Use cos⁡θ=a/r\cos\theta = a/r and sin⁡θ=b/r\sin\theta = b/r together, with r=a2+b2r = \sqrt{a^2 + b^2}. The shortcut θ=arctan⁡(b/a)\theta = \arctan(b/a) only returns angles between −90°-90° and 90°90°, so it is right only when a>0a > 0. For λ=−1+i\lambda = -1 + i it gives −45°-45°, but the cosine is negative and the sine positive, so the angle is 135°135°, and only 135°135° rebuilds the matrix's rotation-scaling block.Read more →

Why can crystals not have fivefold symmetry?

+
Because a rotation that maps a lattice onto itself has integer entries in a basis of lattice vectors, so its trace is an integer. In the plane its eigenvalues are e±iθe^{\pm i\theta} and the trace is 2cos⁡θ2\cos\theta, which must be one of −2,−1,0,1,2-2, -1, 0, 1, 2. That leaves only 2-, 3-, 4- and 6-fold turns. A 72°72° turn would need 2cos⁡72°≈0.6182\cos 72° \approx 0.618, which is not an integer.Read more →