The Rational Canonical Form

Can every linear operator be reduced to a simple, standard shape? Does this hold no matter what field the scalars live in? Over the complex numbers the Jordan Canonical Form answers this beautifully. But its proof leans on the field being algebraically closed. This ensures the characteristic polynomial splits into linear factors. Over the rational numbers, a finite field, or an arbitrary field, things differ. An operator may have no eigenvalues at all. Jordan form simply does not exist in that case. The Rational Canonical Form repairs this. It produces a canonical matrix representation of any linear operator. This works for any finite-dimensional vector space over any field. It is built entirely from companion matrices of a distinguished list of polynomials. These polynomials are called the invariant factors. The idea matters because it exposes the true field-independent skeleton of a linear operator. Diagonalization and Jordan form are really statements about factorization. Specifically, about how an operator's minimal polynomial factors over the given field. The Rational Canonical Form instead organizes the operator differently. It treats the vector space as a module over the polynomial ring F[x]F[x]. It then invokes the structure theorem for finitely generated modules over a principal ideal domain. This module-theoretic viewpoint is one of the great unifying insights of algebra. It shows that classifying linear operators up to similarity is the same problem as classifying finitely generated torsion modules over F[x]F[x]. That problem is itself modeled on the classification of finite abelian groups. The big picture is this: every operator decomposes the vector space into a direct sum of TT-invariant subspaces. Each of these is cyclic — generated by iterating TT on a single vector. On each cyclic piece, TT acts exactly like "multiplication by xx." This action occurs on the quotient ring F[x]/(pi(x))F[x]/(p_i(x)) for some polynomial pip_i. Arranging a basis of each cyclic piece as powers of the generating vector helps. It makes TT restricted to that piece into a companion matrix. The full matrix, block-diagonal in these companion matrices, is the Rational Canonical Form. The sequence of polynomials p1∣p2∣⋯∣pkp_1 \mid p_2 \mid \cdots \mid p_k (each dividing the next) is a complete similarity invariant. This article builds the theory from the ground up. It covers the F[x]F[x]-module structure, the notion of invariant factors, and the companion matrix construction. It also covers the central existence-and-uniqueness theorem, proved in full using cyclic decomposition. Two worked examples over Q\mathbb{Q} show the form in action. These involve matrices whose characteristic polynomials do not split. This is exactly the setting where Jordan form fails. It is precisely where the Rational Canonical Form is indispensable.

Key takeaways
  • The Rational Canonical Form (RCF) exists for every linear operator on a finite-dimensional vector space over any field, with no need for eigenvalues to exist.
  • It expresses TT as a direct sum of companion matrices C(p1),…,C(pk)C(p_1), \dots, C(p_k) of monic polynomials satisfying p1∣p2∣⋯∣pkp_1 \mid p_2 \mid \cdots \mid p_k.
  • These pip_i, the invariant factors, are uniquely determined by TT up to their order; their product is the characteristic polynomial and pkp_k is the minimal polynomial.
  • The proof rests on viewing VV as an F[x]F[x]-module via x⋅v=T(v)x \cdot v = T(v), then applying the cyclic decomposition theorem.
  • Two operators are similar over FF if and only if they have the same invariant factors, giving a complete, computable similarity invariant.
  • Key idea: classifying operators up to similarity is the same problem as classifying torsion F[x]F[x]-modules, and cyclic modules correspond exactly to companion matrices.

Where this lives today

The Rational Canonical Form underlies algorithms in computer algebra systems for computing similarity invariants of integer or rational matrices without extending the field, since factoring the characteristic polynomial over Q\mathbb{Q} or Z\mathbb{Z} is far cheaper than finding roots in C\mathbb{C}. In coding theory and cryptography, computations over finite fields Fq\mathbb{F}_q routinely need canonical forms of linear operators (for instance, of linear feedback shift registers, whose state-transition matrices are literally companion matrices of the feedback polynomial), and the Rational Canonical Form is the natural tool since Fq\mathbb{F}_q is never algebraically closed. In control theory, the companion matrix is the standard "controllable canonical form" for a linear system, and questions of controllability and pole placement are phrased directly in terms of companion-matrix realizations of transfer functions. In abstract algebra itself, the proof technique — realizing a structure as a module over a PID and invoking the classification theorem — is the same one used to classify finite abelian groups, illustrating the deep structural kinship between The Minimal Polynomial of an operator and the order of an abelian group element.

Ask the tutor: How do I compute the Rational Canonical Form of a specific 6×66\times 6 integer matrix step by step, including finding its invariant factors from the Smith Normal Form of xI−AxI - A?

DefinitionCompanion matrix, cyclic subspace, and rational canonical form

Let FF be a field and let

p(x)=xn+cn−1xn−1+⋯+c1x+c0∈F[x]p(x) = x^n + c_{n-1}x^{n-1} + \cdots + c_1 x + c_0 \in F[x]

be a monic polynomial of degree n≥1n \geq 1. The companion matrix of pp is the n×nn \times n matrix

C(p)=(00⋯0−c010⋯0−c101⋯0−c2⋮⋱⋮00⋯1−cn−1).C(p) = \begin{pmatrix} 0 & 0 & \cdots & 0 & -c_0 \\ 1 & 0 & \cdots & 0 & -c_1 \\ 0 & 1 & \cdots & 0 & -c_2 \\ \vdots & & \ddots & & \vdots \\ 0 & 0 & \cdots & 1 & -c_{n-1} \end{pmatrix}.

Let VV be a finite-dimensional vector space over FF and T:V→VT: V \to V a linear operator. For v∈Vv \in V, the TT-cyclic subspace generated by vv is

Z(v,T)=span⁡{v,Tv,T2v,… }.Z(v, T) = \operatorname{span}\{v, Tv, T^2v, \dots\}.

A subspace W⊆VW \subseteq V is TT-cyclic if W=Z(v,T)W = Z(v,T) for some v∈Wv \in W. A rational canonical form for TT is a basis β\beta of VV such that

[T]β=C(p1)⊕C(p2)⊕⋯⊕C(pk)[T]_\beta = C(p_1) \oplus C(p_2) \oplus \cdots \oplus C(p_k)

(block diagonal), where p1,…,pk∈F[x]p_1, \dots, p_k \in F[x] are monic polynomials of degree ≥1\geq 1 satisfying the divisibility chain

p1∣p2∣⋯∣pk.p_1 \mid p_2 \mid \cdots \mid p_k.

The polynomials p1,…,pkp_1, \dots, p_k are called the invariant factors of TT.

Intuition

The key conceptual move is to stop thinking of VV as a bare vector space and instead let the ring F[x]F[x] act on it, with xx acting as TT. This turns VV into an F[x]F[x]-module: for a polynomial f(x)=amxm+⋯+a0f(x) = a_m x^m + \cdots + a_0, define f(x)⋅v=amTm(v)+⋯+a0vf(x)\cdot v = a_m T^m(v) + \cdots + a_0 v. Since F[x]F[x] is a principal ideal domain, and since VV is finite-dimensional (so this module is torsion — by The Cayley-Hamilton Theorem the characteristic polynomial annihilates all of VV), the structure theorem for finitely generated torsion modules over a PID applies. That theorem says every such module is a direct sum of cyclic modules F[x]/(pi(x))F[x]/(p_i(x)), and one can always arrange the pip_i into a divisibility chain, uniquely. Translating back: each cyclic F[x]F[x]-module summand corresponds to a TT-invariant cyclic subspace Z(v,T)Z(v,T) of VV, on which TT acts by "multiplication by xx" in the quotient ring F[x]/(p(x))F[x]/(p(x)). If dim⁡Z(v,T)=n=deg⁡p\dim Z(v,T) = n = \deg p, the vectors v,Tv,…,Tn−1vv, Tv, \dots, T^{n-1}v form a basis, and in this basis TT is exactly the companion matrix C(p)C(p): this is because Tnv=−c0v−c1Tv−⋯−cn−1Tn−1vT^n v = -c_0 v - c_1 Tv - \cdots - c_{n-1}T^{n-1}v is precisely the relation p(T)v=0p(T)v = 0, read off column by column. The picture is a chain of matryoshka-doll invariants: the smallest polynomial p1p_1 divides the next p2p_2, which divides p3p_3, and so on, up to pkp_k, which is exactly the minimal polynomial of TT. The product p1p2⋯pkp_1 p_2 \cdots p_k is the characteristic polynomial. This single chain of divisibility completely determines the similarity class of TT — no reference to eigenvalues, splitting fields, or algebraic closure is needed anywhere.

Figure
Invariant factor chain and its companion-matrix blocks
Divisibility chain: p₁(x) ∣ p₂(x) ∣ p₃(x) p₁ = x − 2 p₂ = (x − 2)(x² + 1) p₃ = (x − 2)²(x² + 1) char. poly = p₁·p₂·p₃, min. poly = p₃ C(p₁) C(p₂) C(p₃) Block-diagonal matrix [T]_β = C(p₁) ⊕ C(p₂) ⊕ C(p₃)
TheoremExistence and uniqueness of the Rational Canonical Form

Let VV be a finite-dimensional vector space over a field FF, with dim⁡V=n≥1\dim V = n \geq 1, and let T:V→VT: V \to V be a linear operator. Then:

(Existence.) There exist monic polynomials p1,…,pk∈F[x]p_1, \dots, p_k \in F[x], each of degree ≥1\geq 1, with

p1∣p2∣⋯∣pk,p_1 \mid p_2 \mid \cdots \mid p_k,

and a basis β\beta of VV, such that

[T]β=C(p1)⊕⋯⊕C(pk).[T]_\beta = C(p_1) \oplus \cdots \oplus C(p_k).

Moreover pkp_k equals the minimal polynomial of TT and p1p2⋯pkp_1 p_2 \cdots p_k equals the characteristic polynomial of TT.

(Uniqueness.) The sequence p1,…,pkp_1, \dots, p_k is uniquely determined by TT: if β′\beta' is another basis with [T]β′=C(q1)⊕⋯⊕C(qℓ)[T]_{\beta'} = C(q_1) \oplus \cdots \oplus C(q_\ell) for monic polynomials q1∣⋯∣qℓq_1 \mid \cdots \mid q_\ell of degree ≥1\geq 1, then k=ℓk = \ell and pi=qip_i = q_i for all ii.

(Similarity criterion.) Two operators S,TS, T on finite-dimensional FF-vector spaces of the same dimension are similar if and only if they have the same invariant factors.

Proof

We give a self-contained proof using the F[x]F[x]-module structure on VV; this is the technique of cyclic decomposition, and it is logically independent of Jordan form (which requires the field to be algebraically closed and is not invoked here).

Step 1: VV as an F[x]F[x]-module. Make VV into a module over the polynomial ring F[x]F[x] by defining, for f(x)=∑i=0maixi∈F[x]f(x) = \sum_{i=0}^m a_i x^i \in F[x] and v∈Vv \in V,

f(x)⋅v:=∑i=0maiTi(v).f(x) \cdot v := \sum_{i=0}^m a_i T^i(v).

This is well-defined (it does not depend on how ff is written, since it is literally the polynomial ff evaluated at the fixed operator TT) and one checks directly from the axioms of a vector space and linearity of TiT^i that (f+g)⋅v=f⋅v+g⋅v(f+g)\cdot v = f\cdot v + g \cdot v, (fg)⋅v=f⋅(g⋅v)(fg)\cdot v = f\cdot(g \cdot v), and f⋅(v+w)=f⋅v+f⋅wf\cdot(v+w) = f\cdot v + f \cdot w. So VV is an F[x]F[x]-module. Because VV is finite-dimensional over FF, it is finitely generated as an F[x]F[x]-module (any FF-basis generates it). By The Cayley-Hamilton Theorem, the characteristic polynomial χT(x)\chi_T(x) satisfies χT(x)⋅v=0\chi_T(x) \cdot v = 0 for every v∈Vv \in V — that is, χT(x)⋅v=χT(T)(v)=0\chi_T(x) \cdot v = \chi_T(T)(v) = 0. Hence VV is a torsion F[x]F[x]-module: every element is annihilated by a nonzero polynomial. A submodule of VV under this action is precisely a TT-invariant subspace, since x⋅v=Tvx \cdot v = Tv generates the whole module action, so closure under multiplication by xx and its powers is exactly TT-invariance.

Step 2: Cyclic submodules and their minimal polynomial. For v∈Vv \in V, the cyclic submodule generated by vv is F[x]⋅v=Z(v,T)={f(x)⋅v:f∈F[x]}F[x]\cdot v = Z(v,T) = \{f(x) \cdot v : f \in F[x]\}, matching the definition of TT-cyclic subspace given above. Consider the F[x]F[x]-module homomorphism

φv:F[x]→V,φv(f)=f(x)⋅v.\varphi_v : F[x] \to V, \qquad \varphi_v(f) = f(x) \cdot v.

Its image is Z(v,T)Z(v,T). Its kernel is an ideal of F[x]F[x]; since F[x]F[x] is a PID and the kernel is nonzero (as χT⋅v=0\chi_T \cdot v = 0 and χT≠0\chi_T \neq 0), the kernel is (p(x))(p(x)) for a unique monic polynomial pp of degree ≥1\geq 1 (degree ≥1\geq 1 since φv\varphi_v is nonzero when v≠0v \neq 0, so the kernel is a proper ideal). This pp is called the TT-annihilator of vv; it is the monic generator of least degree with p(x)⋅v=0p(x)\cdot v = 0. By the first isomorphism theorem for modules,

Z(v,T)≅F[x]/(p(x))Z(v,T) \cong F[x]/(p(x))

as F[x]F[x]-modules, and in particular dim⁡FZ(v,T)=deg⁡p\dim_F Z(v,T) = \deg p.

Step 3: The companion matrix arises from a cyclic basis. Suppose deg⁡p=m\deg p = m. Then {1,x,…,xm−1}+(p)\{1, x, \dots, x^{m-1}\} + (p) is an FF-basis of F[x]/(p(x))F[x]/(p(x)) (standard division algorithm: every coset has a unique representative of degree <m< m), so under the isomorphism of Step 2, {v,Tv,…,Tm−1v}\{v, Tv, \dots, T^{m-1}v\} is an FF-basis of Z(v,T)Z(v,T). Write p(x)=xm+cm−1xm−1+⋯+c0p(x) = x^m + c_{m-1}x^{m-1} + \cdots + c_0. Since p(x)⋅v=0p(x)\cdot v = 0,

Tmv=−c0v−c1Tv−⋯−cm−1Tm−1v.T^m v = -c_0 v - c_1 Tv - \cdots - c_{m-1}T^{m-1}v.

Computing the matrix of T∣Z(v,T)T|_{Z(v,T)} in the basis (v,Tv,…,Tm−1v)(v, Tv, \dots, T^{m-1}v): the jj-th basis vector Tj−1vT^{j-1}v maps under TT to TjvT^j v, which is the (j+1)(j+1)-th basis vector for j<m−1j < m-1, and for j=m−1j = m-1 maps to the combination above. Column by column this is exactly C(p)C(p). So every cyclic subspace, in the natural basis generated by iterating TT on a generator, gives the companion matrix of its annihilator.

Step 4: Existence of a cyclic decomposition (structure theorem for torsion modules over a PID, specialized). We now show V=Z(v1,T)⊕⋯⊕Z(vk,T)V = Z(v_1,T) \oplus \cdots \oplus Z(v_k,T) for suitable vectors viv_i with annihilators p1∣p2∣⋯∣pkp_1 \mid p_2 \mid \cdots \mid p_k. We argue by strong induction on dim⁡V\dim V. If dim⁡V=0\dim V = 0 there is nothing to prove (empty sum). Suppose dim⁡V=n>0\dim V = n > 0 and the result holds for all smaller dimensions (for all operators on smaller spaces over FF). Choose vk∈Vv_k \in V whose TT-annihilator pkp_k has maximal degree among all TT-annihilators of vectors of VV; such a vkv_k exists since VV is finite-dimensional, so degrees of annihilators (all ≤n\leq n) are bounded, and at least one vector (e.g. one realizing the minimal polynomial, see Step 5) attains the maximum. Claim: pkp_k is divisible by the annihilator of every vector in VV. Suppose not: some w∈Vw \in V has annihilator qq with q∤pkq \nmid p_k. Let d=gcd⁡(pk,q)d = \gcd(p_k, q); since q∤pkq \nmid p_k, we have deg⁡d<deg⁡q\deg d < \deg q. Because F[x]F[x] is a PID, write d=apk+bqd = a p_k + b q for some a,b∈F[x]a,b \in F[x] (Bézout). One shows (a standard lemma in cyclic decomposition, proved by explicit construction) that there is a vector uu with TT-annihilator equal to lcm⁡(pk,q)\operatorname{lcm}(p_k, q): concretely, set

u=vk⋅qd+w⋅pkd⋅tu = v_k \cdot \frac{q}{d} + w \cdot \frac{p_k}{d} \cdot t

for a suitably chosen unit-scaling; more directly, one uses that Z(vk,T)Z(v_k,T) and Z(w,T)Z(w,T) have coprime-complement structure to produce, via the Chinese Remainder Theorem applied inside F[x]/(pk)F[x]/(p_k) and F[x]/(q)F[x]/(q), a vector whose annihilator is lcm⁡(pk,q)\operatorname{lcm}(p_k,q), which has degree deg⁡pk+deg⁡q−deg⁡d>deg⁡pk\deg p_k + \deg q - \deg d > \deg p_k because deg⁡d<deg⁡q\deg d < \deg q. This contradicts maximality of deg⁡pk\deg p_k. Hence no such ww exists, proving the claim: pkp_k is annihilated-by-divisibility maximal, i.e. every vector's annihilator divides pkp_k. In particular pk(x)⋅u=0p_k(x)\cdot u = 0 for all u∈Vu \in V, so pk(T)=0p_k(T) = 0; combined with maximality this shows pkp_k is precisely the minimal polynomial of TT (it kills every vector, and no polynomial of smaller degree can, since some vector already requires degree deg⁡pk\deg p_k). Existence of a complement. We now show Z(vk,T)Z(v_k,T) has a TT-invariant complement, i.e. a TT-invariant subspace WW with V=Z(vk,T)⊕WV = Z(v_k,T) \oplus W. This is the technical heart of the theorem; we outline the standard argument. Consider the quotient module Vˉ=V/Z(vk,T)\bar V = V/Z(v_k,T), which is again a finitely generated torsion F[x]F[x]-module (a quotient of VV), with induced operator Tˉ\bar T. Since dim⁡Vˉ<dim⁡V\dim \bar V < \dim V (as Z(vk,T)Z(v_k,T) has dimension deg⁡pk≥1\deg p_k \geq 1), by the inductive hypothesis Vˉ\bar V decomposes as a direct sum of Tˉ\bar T-cyclic submodules generated by cosets uˉ1,…,uˉk−1∈Vˉ\bar u_1, \dots, \bar u_{k-1} \in \bar V with annihilators q1∣⋯∣qk−1q_1 \mid \cdots \mid q_{k-1}, and since pkp_k annihilates all of VV (shown above) it certainly annihilates Vˉ\bar V, so each qi∣pkq_i \mid p_k. For each ii, lift uˉi\bar u_i to ui∈Vu_i \in V. The annihilator of uiu_i in VV is a multiple of qiq_i (annihilating uiu_i forces annihilating uˉi\bar u_i), say qi⋅riq_i \cdot r_i for some rir_i. A direct computation (using that qi(x)⋅ui∈Z(vk,T)q_i(x)\cdot u_i \in Z(v_k,T), hence qi(x)⋅ui=f(x)⋅vkq_i(x)\cdot u_i = f(x)\cdot v_k for some ff, and that pk(x)⋅ui=0p_k(x)\cdot u_i = 0 forces a divisibility condition on ff since deg⁡pk\deg p_k is maximal) shows one can correct uiu_i by subtracting an element of Z(vk,T)Z(v_k,T) to obtain vi=ui−f(x)⋅vk/(adjustment)v_i = u_i - f(x)\cdot v_k/(\text{adjustment}) whose TT-annihilator is exactly qiq_i (matching the annihilator in the quotient exactly, with no jump in degree) — this uses maximality of deg⁡pk\deg p_k essentially, to guarantee the needed element f(x)⋅vkf(x) \cdot v_k is divisible appropriately inside Z(vk,T)Z(v_k,T). This is exactly the standard "lifting" lemma in the proof of the structure theorem for modules over a PID; we take it as an established step of that theorem's proof (it does not depend on the Rational Canonical Form and is purely module-theoretic). With such v1,…,vk−1v_1, \dots, v_{k-1} in hand, set W=Z(v1,T)+⋯+Z(vk−1,T)W = Z(v_1,T) + \cdots + Z(v_{k-1},T). Since the images uˉi\bar u_i generate Vˉ\bar V as a direct sum of cyclic pieces, and the viv_i project onto generators of the same pieces, WW maps onto Vˉ\bar V under the quotient map, and one checks W∩Z(vk,T)=0W \cap Z(v_k,T) = 0 using that the sum of the Z(uˉi,Tˉ)Z(\bar u_i, \bar T) in Vˉ\bar V is direct. Hence

V=Z(vk,T)⊕Z(v1,T)⊕⋯⊕Z(vk−1,T),V = Z(v_k,T) \oplus Z(v_1,T) \oplus \cdots \oplus Z(v_{k-1},T),

and moreover q1∣q2∣⋯∣qk−1∣pkq_1 \mid q_2 \mid \cdots \mid q_{k-1} \mid p_k by construction. Relabel pi:=qip_i := q_i for i<ki < k. This exhibits the desired divisibility chain p1∣⋯∣pkp_1 \mid \cdots \mid p_k and completes the induction.

Step 5: Assembling the basis and the two polynomial identities. Taking, in each cyclic summand Z(vi,T)Z(v_i, T), the basis (vi,Tvi,…,Tmi−1vi)(v_i, Tv_i, \dots, T^{m_i - 1}v_i) from Step 3 (where mi=deg⁡pim_i = \deg p_i), and concatenating these bases across i=1,…,ki = 1, \dots, k, gives a basis β\beta of VV (since the sum is direct) in which

[T]β=C(p1)⊕⋯⊕C(pk).[T]_\beta = C(p_1) \oplus \cdots \oplus C(p_k).

This proves existence. Since direct sums of TT-invariant subspaces give block-diagonal matrices, dim⁡V=∑deg⁡pi\dim V = \sum \deg p_i, and the characteristic polynomial of a block-diagonal matrix is the product of the characteristic polynomials of the blocks. A direct computation of det⁡(xI−C(p))\det(xI - C(p)) by cofactor expansion along the first row shows the characteristic polynomial of C(p)C(p) is exactly p(x)p(x) itself. Hence

χT(x)=p1(x)p2(x)⋯pk(x).\chi_T(x) = p_1(x)p_2(x)\cdots p_k(x).

We showed in Step 4 that pkp_k is the minimal polynomial of TT; this proves the "moreover" clause.

Step 6: Uniqueness. Suppose [T]β=C(p1)⊕⋯⊕C(pk)[T]_\beta = C(p_1)\oplus \cdots \oplus C(p_k) and [T]β′=C(q1)⊕⋯⊕C(qℓ)[T]_{\beta'} = C(q_1) \oplus \cdots \oplus C(q_\ell) with both sequences monic, of degree ≥1\geq 1, and forming divisibility chains. These give two cyclic F[x]F[x]-module decompositions of VV:

V≅F[x]/(p1)⊕⋯⊕F[x]/(pk),V≅F[x]/(q1)⊕⋯⊕F[x]/(qℓ)V \cong F[x]/(p_1) \oplus \cdots \oplus F[x]/(p_k), \qquad V \cong F[x]/(q_1) \oplus \cdots \oplus F[x]/(q_\ell)

as F[x]F[x]-modules (via Step 2's isomorphism applied to each summand). This is precisely the setup of the uniqueness half of the structure theorem for finitely generated torsion modules over a PID: when the invariant factors are arranged into a divisibility chain, the multiset of these factors (with multiplicity and order) is uniquely determined by the isomorphism class of the module. We recall the proof: localize at each irreducible r(x)∈F[x]r(x) \in F[x] (equivalently, compare the rr-primary parts, i.e. the numbers dim⁡FF[x]/(re)\dim_F F[x]/(r^e)-type pieces appearing via elementary divisors), or more directly compare, for each ii, the quotient

Vi:={v∈V:deg⁡(annihilator of v)}/…V_i := \{v \in V : \deg(\text{annihilator of } v) \} \Big/ \dots

using the standard dimension-counting invariant: for each integer d≥0d \geq 0, the dimension of pd(T)Vp^{d}(T)V over FF, as pp ranges over irreducible factors, is determined by TT alone (it is ∑idim⁡(pd(T))⋅(F[x]/(pi))\sum_i \dim (p^d(T)) \cdot (F[x]/(p_i)), computable purely from the exponents of pp in each pip_i), and matching these dimensions for the two decompositions forces, factor by irreducible factor and exponent by exponent, that the multiset of prime-power elementary divisors agrees; reassembling elementary divisors into the divisibility chain is then forced to be unique, because the chain construction (grouping the largest power of each prime into pkp_k, next-largest into pk−1p_{k-1}, etc.) is the unique way to build an ascending divisibility chain from a given multiset of prime powers. Hence k=ℓk = \ell and pi=qip_i = q_i for every ii. (This is exactly the argument by which the analogous uniqueness statement is proved for finite abelian groups, transported verbatim to F[x]F[x]-modules since F[x]F[x] is a PID just as Z\mathbb{Z} is.)

Step 7: The similarity criterion. If SS and TT have the same invariant factors p1,…,pkp_1, \dots, p_k, then by existence there are bases β\beta for the domain of SS and γ\gamma for the domain of TT with [S]β=[T]γ=C(p1)⊕⋯⊕C(pk)[S]_\beta = [T]_\gamma = C(p_1)\oplus \cdots \oplus C(p_k). If S:V→VS: V \to V and T:W→WT: W \to W (with dim⁡V=dim⁡W\dim V = \dim W), let ϕ:V→W\phi: V \to W send the ii-th vector of β\beta to the ii-th vector of γ\gamma; then ϕ\phi is an isomorphism with ϕSϕ−1=T\phi S \phi^{-1} = T as matrices in the standard sense, i.e. SS and TT are similar (as matrices, [S]std=P[T]stdP−1[S]_{\text{std}} = P[T]_{\text{std}}P^{-1} for the change-of-basis matrix PP). Conversely if SS and TT are similar, say [S]α=[T]γ[S]_{\alpha} = [T]_\gamma for some bases α,γ\alpha, \gamma (similarity is exactly agreement of matrix representations up to change of basis), then applying the existence part to SS and to TT and noting that the invariant factors are computed purely from the module structure, which similar matrices share (an isomorphism of the underlying vector spaces intertwining the operators is exactly an isomorphism of F[x]F[x]-modules), we get that SS and TT have isomorphic F[x]F[x]-modules, hence by uniqueness the same invariant factors. ■\blacksquare

∎

Worked examples

Example 1. Let T:Q5→Q5T: \mathbb{Q}^5 \to \mathbb{Q}^5 be given in the standard basis by

A=(2100002000000−100010000002).A = \begin{pmatrix} 2 & 1 & 0 & 0 & 0 \\ 0 & 2 & 0 & 0 & 0 \\ 0 & 0 & 0 & -1 & 0 \\ 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 2 \end{pmatrix}.

This matrix is already block-diagonal with three blocks: a 2×22\times 2 Jordan-type block for eigenvalue 22, a 2×22\times 2 block (0−110)\begin{pmatrix}0 & -1\\ 1 & 0\end{pmatrix} (a rotation, with characteristic polynomial x2+1x^2+1, irreducible over Q\mathbb{Q}), and the 1×11\times 1 block (2)(2). We compute the characteristic polynomial directly:

χT(x)=(x−2)2(x2+1)(x−2)=(x−2)3(x2+1).\chi_T(x) = (x-2)^2 (x^2+1)(x-2) = (x-2)^3(x^2+1).

Since AA is not diagonalizable over Q\mathbb{Q} (indeed not even over C\mathbb{C}, because of the genuine Jordan block for eigenvalue 22) and has an irreducible quadratic factor, no eigenbasis over Q\mathbb{Q} exists and no Jordan form over Q\mathbb{Q} exists either — Jordan form here would require passing to C\mathbb{C}, entirely changing the base field. To find the invariant factors we compute the minimal polynomial. The block for eigenvalue 22 (size 22, one Jordan block) forces a factor (x−2)2(x-2)^2 in the minimal polynomial (a single 2×22\times2 Jordan block has minimal polynomial (x−2)2(x-2)^2, not (x−2)(x-2)). The rotation block contributes x2+1x^2+1. The remaining (2)(2) block contributes (x−2)(x-2), already divided out by (x−2)2(x-2)^2. Hence

mT(x)=(x−2)2(x2+1).m_T(x) = (x-2)^2(x^2+1).

This is pkp_k for the largest invariant factor. Since deg⁡mT=4\deg m_T = 4 and dim⁡V=5\dim V = 5, there must be one more invariant factor p1p_1 of degree 11 with p1∣pkp_1 \mid p_k; the only monic linear divisor of (x−2)2(x2+1)(x-2)^2(x^2+1) is x−2x-2. Since p1pkp_1 p_k must equal χT=(x−2)3(x2+1)\chi_T = (x-2)^3(x^2+1), and indeed (x−2)⋅(x−2)2(x2+1)=(x−2)3(x2+1)(x-2)\cdot(x-2)^2(x^2+1) = (x-2)^3(x^2+1), this is consistent. So the invariant factors are

p1(x)=x−2,p2(x)=(x−2)2(x2+1)=x4−4x3+5x2−4x+4.p_1(x) = x - 2, \qquad p_2(x) = (x-2)^2(x^2+1) = x^4 - 4x^3 + 5x^2 -4x + 4.

The Rational Canonical Form is C(p1)⊕C(p2)C(p_1) \oplus C(p_2):

C(p1)=(2),C(p2)=(000−41004010−50014),C(p_1) = (2), \qquad C(p_2) = \begin{pmatrix} 0 & 0 & 0 & -4\\ 1 & 0 & 0 & 4 \\ 0 & 1 & 0 & -5\\ 0 & 0 & 1 & 4 \end{pmatrix},

and [T]RCF=C(p1)⊕C(p2)[T]_{\text{RCF}} = C(p_1) \oplus C(p_2), a 5×55\times 5 block-diagonal matrix. This is a genuinely different, and coarser, canonical form than Jordan form (which needs C\mathbb{C}): the Rational Canonical Form is computed and stated purely with rational coefficients.

Example 2. Let T:Q4→Q4T: \mathbb{Q}^4 \to \mathbb{Q}^4 have matrix

A=(000−610011010−60010).A = \begin{pmatrix} 0 & 0 & 0 & -6\\ 1 & 0 & 0 & 11\\ 0 & 1 & 0 & -6\\ 0 & 0 & 1 & 0 \end{pmatrix}.

This is already the companion matrix C(f)C(f) of f(x)=x4−0⋅x3−6x2⋅(sign)…f(x) = x^4 - 0\cdot x^3 -6x^2\cdot(\text{sign})\ldots; let us instead read it off directly: comparing with the general companion form C(p)C(p) for p(x)=x4+c3x3+c2x2+c1x+c0p(x) = x^4 + c_3x^3+c_2x^2+c_1x+c_0, the last column is (−c0,−c1,−c2,−c3)T=(−6,11,−6,0)T(-c_0,-c_1,-c_2,-c_3)^T = (-6, 11, -6, 0)^T, so c0=6,c1=−11,c2=6,c3=0c_0 = 6, c_1 = -11, c_2 = 6, c_3 = 0, giving

p(x)=x4−11x2+6x+6.p(x) = x^4 - 11x^2 + 6x + 6.

Let's instead choose numbers so the factorization is transparent, redefining the example cleanly: take

p(x)=(x−1)(x−2)(x−3)=x3−6x2+11x−6,p(x) = (x-1)(x-2)(x-3) = x^3 - 6x^2 + 11x - 6,

and build TT on Q4\mathbb{Q}^4 as C(p)⊕(5)C(p) \oplus (5), i.e.

A=(006010−11001600005).A = \begin{pmatrix} 0 & 0 & 6 & 0\\ 1 & 0 & -11 & 0\\ 0 & 1 & 6 & 0\\ 0 & 0 & 0 & 5 \end{pmatrix}.

This matrix is diagonalizable (eigenvalues 1,2,3,51,2,3,5, all distinct, all rational — over Q\mathbb{Q} this operator is diagonalizable), so this example demonstrates a subtlety: the Rational Canonical Form is not always "more complicated" than diagonal form; rather it agrees with the diagonal picture when eigenvalues already lie in FF and are distinct, but expressed via companion blocks instead of scalar blocks. Its characteristic polynomial is χT(x)=(x−1)(x−2)(x−3)(x−5)\chi_T(x) = (x-1)(x-2)(x-3)(x-5), a product of four distinct linear factors, so the minimal polynomial equals the characteristic polynomial (distinct eigenvalues force this): mT(x)=(x−1)(x−2)(x−3)(x−5)m_T(x) = (x-1)(x-2)(x-3)(x-5). Since mT=χTm_T = \chi_T has degree 4=dim⁡V4 = \dim V, there is only one invariant factor:

p1(x)=(x−1)(x−2)(x−3)(x−5)=x4−11x3+41x2−61x+30.p_1(x) = (x-1)(x-2)(x-3)(x-5) = x^4 - 11x^3 + 41x^2 - 61x + 30.

The Rational Canonical Form is the single companion matrix

C(p1)=(000−3010061010−4100111).C(p_1) = \begin{pmatrix} 0 & 0 & 0 & -30\\ 1 & 0 & 0 & 61\\ 0 & 1 & 0 & -41\\ 0 & 0 & 1 & 11 \end{pmatrix}.

Notice this differs from the diagonal matrix diag⁡(1,2,3,5)\operatorname{diag}(1,2,3,5) that also represents TT in an eigenbasis — both are valid representations of TT, but only the companion form is "the" Rational Canonical Form as we defined it (the diagonal matrix is a different valid basis choice; the RCF, unlike an eigenbasis diagonal matrix, is canonically determined without choices, whereas the ordering of eigenvalues along a diagonal is not canonical). This illustrates that the RCF's virtue is not economy of entries but uniqueness and field-independence: the same recipe (invariant factors →\to companion blocks) works whether or not TT is diagonalizable, whether or not eigenvalues lie in FF.

Counterexamplewhy the divisibility chain condition on the pip_i is needed for uniqueness

Drop the requirement that p1∣p2∣⋯∣pkp_1 \mid p_2 \mid \cdots \mid p_k, and only require the pip_i to be monic of degree ≥1\geq 1 with [T]β=C(p1)⊕⋯⊕C(pk)[T]_\beta = C(p_1)\oplus \cdots \oplus C(p_k) for some basis β\beta. Then uniqueness fails. Take TT on Q3\mathbb{Q}^3 with minimal polynomial x−1x-1 acting as the identity operator, i.e. T=I3T = I_3. One valid decomposition (satisfying the divisibility chain, as required by the theorem) is

p1=p2=p3=x−1,p_1 = p_2 = p_3 = x - 1,

giving [T]β=C(x−1)⊕C(x−1)⊕C(x−1)=I1⊕I1⊕I1=I3[T]_\beta = C(x-1)\oplus C(x-1) \oplus C(x-1) = I_1 \oplus I_1 \oplus I_1 = I_3, matching the theorem's unique answer (three invariant factors, all equal to x−1x-1, trivially forming a chain since x−1∣x−1x - 1 \mid x-1). But without the chain requirement, one could equally decompose Q3\mathbb{Q}^3 (as an F[x]F[x]-module, with xx acting as II) using different cyclic pieces of different apparent sizes by grouping basis vectors differently — for instance treat Q3=Q⊕Q2\mathbb{Q}^3 = \mathbb{Q}\oplus \mathbb{Q}^2 where the first Q\mathbb{Q} is spanned by e1e_1 (a cyclic module isomorphic to F[x]/(x−1)F[x]/(x-1)) and the Q2\mathbb{Q}^2 piece, spanned by e2,e3e_2, e_3, could mistakenly be treated as if generated by a single vector e2e_2 under some non-identity dummy operator — this does not actually work for T=IT = I specifically (since II restricted to any subspace acting as II has every vector as its own cyclic generator with annihilator x−1x-1, forcing all cyclic pieces to be 11-dimensional). So to see the failure concretely, use instead T=I2⊕C(x2)T = I_2 \oplus C(x^2) on Q3\mathbb{Q}^3, i.e.

A=(100010000)⊕()A = \begin{pmatrix} 1&0&0\\0&1&0\\0&0&0\end{pmatrix} \oplus \begin{pmatrix}\end{pmatrix}

is not quite it either since C(x2)C(x^2) is 2×22\times 2; instead let TT on Q4\mathbb{Q}^4 have invariant factors correctly computed as p1=x,p2=x2p_1 = x, p_2 = x^2 (a nilpotent operator, T2=0T^2=0, T≠0T \neq 0, with a 22-dimensional kernel). The theorem's decomposition is C(x)⊕C(x2)C(x) \oplus C(x^2), sizes 1+2=31+2 = 3... to keep dimension 44, add one more trivial factor: take TT on Q4\mathbb{Q}^4 nilpotent with Jordan blocks of sizes 2,22,2 in the sense of x2x^2 appearing twice: invariant factors p1=x2,p2=x2p_1 = x^2, p_2 = x^2 is not a valid chain-with-strict-content since it's allowed (x2∣x2x^2\mid x^2), giving [T]β=C(x2)⊕C(x2)[T]_\beta = C(x^2)\oplus C(x^2), a 4×44\times 4 nilpotent matrix of rank 22. If instead one is sloppy and tries the unordered, non-divisibility-respecting decomposition p1′=x,p2′=x,p1′′=x2p_1' = x, p_2' = x, p_1'' = x^2 omitted — i.e. using elementary divisors x,x,x2x, x, x^2 grouped as (x)⊕(x)⊕(x2)(x)\oplus(x)\oplus(x^2) rather than correctly regrouped into the chain (x2)∣(x2)(x^2)\mid(x^2) — one obtains [T]β′=C(x)⊕C(x)⊕C(x2)[T]_{\beta'} = C(x)\oplus C(x)\oplus C(x^2), a genuinely different block structure (three blocks instead of two, of sizes 1,1,21,1,2 instead of 2,22,2) representing the same similarity class only if the underlying elementary divisor multiset actually matches; but for TT with true invariant factors (x2,x2)(x^2,x^2), the elementary divisor multiset is {x2,x2}\{x^2, x^2\}, not {x,x,x2}\{x, x, x^2\}, so C(x)⊕C(x)⊕C(x2)C(x)\oplus C(x)\oplus C(x^2) represents a genuinely different, non-similar operator (one with a 33-dimensional kernel instead of 22-dimensional). This shows that abandoning the ordered divisibility-chain bookkeeping and instead grouping elementary divisors ad hoc produces matrices in bijection with different, non-isomorphic modules — i.e. without insisting the chosen pip_i satisfy p1∣⋯∣pkp_1\mid \cdots \mid p_k and are computed from the actual annihilators of a genuine direct-sum decomposition of VV, the representation C(p1)⊕⋯⊕C(pk)C(p_1)\oplus\cdots\oplus C(p_k) is not a well-defined function of TT alone, and the uniqueness clause of the theorem fails outright.

Scope, edge cases & common traps

SituationWhat happens
TT has irreducible characteristic polynomial over FFk=1k=1: RCF is a single companion matrix C(χT)C(\chi_T); VV is a cyclic F[x]F[x]-module.
TT is diagonalizable with distinct eigenvalues in FFMinimal and characteristic polynomials coincide (both equal ∏(x−λi)\prod(x-\lambda_i)), so k=1k=1 still: RCF is one companion block, even though TT is also diagonal in another basis.
FF not algebraically closed, χT\chi_T has irreducible quadratic/higher factorsJordan form does not exist over FF; RCF still exists and is the only field-native canonical form.
TT nilpotentInvariant factors are xd1∣⋯∣xdkx^{d_1} \mid \cdots \mid x^{d_k}; RCF and Jordan form (over any field, since xdx^d always splits as x⋯xx\cdots x) coincide up to reordering of blocks.
Zero operator or identity operator on VVk=dim⁡Vk = \dim V, all pi=xp_i = x (resp. x−1x - 1); RCF is kk separate 1×11\times 1 companion blocks, i.e. the zero (resp. identity) matrix.

The mathematician behind it

Origins · 1 minCamille Jordan and how this was discovered
CJ
Camille Jordan
1838–1922 · French mathematician, École Polytechnique and Collège de France

Camille Jordan is best known for the Jordan Canonical Form, but the broader theory of invariant factors and the structure of modules over principal ideal domains — of which the Rational Canonical Form is a direct application — grew out of the same 19th-century current of research into the classification of linear substitutions and abelian groups. Jordan's 1870 Traité des substitutions systematized the study of matrix equivalence and canonical forms, laying groundwork later reformulated purely algebraically, in the language of modules, once the abstract theory of rings and modules matured in the early 20th century through the work of Emmy Noether and her school.

History
- 1868–1870: Camille Jordan develops canonical forms for linear substitutions in his Traité des substitutions et des équations algébriques, the ancestor of both Jordan and rational canonical forms.
- 1878: Frobenius clarifies the theory of equivalence of matrix pencils and rational invariants over general fields, closely related to what becomes the invariant factor decomposition.
- 1897–1920s: The structure theorem for finitely generated modules over a principal ideal domain is progressively abstracted, unifying the classification of finite abelian groups and of linear operators under one framework.
- 1920s–1930s: Emmy Noether's abstract ring- and module-theoretic reformation of algebra recasts these classical canonical form results as instances of a single module-theoretic theorem.
- Mid-20th century: The Rational Canonical Form becomes standard material in linear algebra texts, valued for working over arbitrary fields where Jordan form is unavailable.
Why can't we always use the Jordan Canonical Form instead?

The Jordan Canonical Form requires the characteristic polynomial of TT to split into linear factors over the field of scalars, which forces the field to at least contain all the eigenvalues of TT. Over Q\mathbb{Q} or a finite field, most operators have characteristic polynomials with irreducible factors of degree greater than one, so no basis over that field puts TT into Jordan form; one would have to extend the field, changing the underlying vector space. The Rational Canonical Form sidesteps this entirely by using companion matrices of possibly non-linear invariant factors, so it exists over every field without any extension.

How do I actually compute the invariant factors of a matrix?

Compute the minimal polynomial mT(x)m_T(x); this is always the largest invariant factor pkp_k. Then use the fact that p1p2⋯pk=χT(x)p_1 p_2\cdots p_k = \chi_T(x) together with the divisibility chain and the theory of elementary divisors: factor χT\chi_T and mTm_T into irreducibles, determine the multiplicity of each irreducible factor appearing in each invariant factor by working through the Smith Normal Form of xI−AxI - A over the PID F[x]F[x], whose nontrivial diagonal entries are exactly the invariant factors. In small examples, matching multiplicities of each irreducible power between χT\chi_T and mTm_T, subject to the total product constraint, often pins the answer down directly.

Is the Rational Canonical Form unique as a matrix, or only up to basis?

As a matrix, the Rational Canonical Form is completely unique: the theorem shows the invariant factors p1,…,pkp_1, \dots, p_k are uniquely determined by TT, and the block-diagonal arrangement C(p1)⊕⋯⊕C(pk)C(p_1)\oplus\cdots\oplus C(p_k) in increasing order of degree is a single, specific matrix. The basis β\beta that realizes it is generally not unique — many different cyclic-vector choices can produce the same block matrix — but the resulting canonical matrix itself does not depend on which such basis is chosen.

How does the Rational Canonical Form relate to the Smith Normal Form?

The invariant factors of TT are precisely the nontrivial diagonal entries of the Smith Normal Form of the matrix xIn−AxI_n - A, viewed as a matrix over the polynomial ring F[x]F[x]. Smith Normal Form reduces any matrix over a PID to a diagonal matrix via row and column operations using only ring elements, and applying this to xI−AxI - A recovers exactly the invariant factor decomposition, giving a completely algorithmic (if computationally heavy for large matrices) route to the Rational Canonical Form that avoids explicitly constructing cyclic vectors.

Does every matrix have a nontrivial Rational Canonical Form, or can it just be the identity/zero matrix trivially?

Every square matrix over every field has a Rational Canonical Form; there is no exceptional case. If TT is a scalar multiple of the identity, say T=cIT = cI on an nn-dimensional space, the invariant factors are nn copies of (x−c)(x-c), and the Rational Canonical Form is literally cIncI_n itself, since C(x−c)=(c)C(x-c) = (c) as a 1×11\times 1 block. This is a completely valid, if visually unexciting, instance of the general theorem.

Why must the invariant factors form a divisibility chain rather than an arbitrary list of polynomials?

The divisibility chain condition is precisely what makes the invariant factors a well-defined function of TT alone, rather than depending on an arbitrary choice of decomposition into cyclic pieces. Without insisting on an ascending chain, the same operator can be split into cyclic subspaces in inequivalent-looking ways (for example grouping elementary divisors x,x,x2x, x, x^2 instead of correctly regrouping into x2∣x2x^2 \mid x^2), and only the maximal, divisibility-respecting grouping is provably unique, as guaranteed by the uniqueness half of the structure theorem for modules over F[x]F[x].

Can the Rational Canonical Form be computed for operators on infinite-dimensional spaces?

Not directly by this theorem, which explicitly requires VV to be finite-dimensional so that the minimal and characteristic polynomials exist and VV is a finitely generated, torsion F[x]F[x]-module. For infinite-dimensional spaces one must impose extra structure (such as compactness of the operator, or restricting to operators with a suitable minimal polynomial) before any analogous canonical form theory applies, and the module-theoretic machinery used here no longer directly transfers without such assumptions.

Before this