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 . 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 . 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 -invariant subspaces. Each of these is cyclic — generated by iterating on a single vector. On each cyclic piece, acts exactly like "multiplication by ." This action occurs on the quotient ring for some polynomial . Arranging a basis of each cyclic piece as powers of the generating vector helps. It makes 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 (each dividing the next) is a complete similarity invariant. This article builds the theory from the ground up. It covers the -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 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.
- 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 as a direct sum of companion matrices of monic polynomials satisfying .
- These , the invariant factors, are uniquely determined by up to their order; their product is the characteristic polynomial and is the minimal polynomial.
- The proof rests on viewing as an -module via , then applying the cyclic decomposition theorem.
- Two operators are similar over 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 -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 or is far cheaper than finding roots in . In coding theory and cryptography, computations over finite fields 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 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 integer matrix step by step, including finding its invariant factors from the Smith Normal Form of ?
Let be a field and let
be a monic polynomial of degree . The companion matrix of is the matrix
Let be a finite-dimensional vector space over and a linear operator. For , the -cyclic subspace generated by is
A subspace is -cyclic if for some . A rational canonical form for is a basis of such that
(block diagonal), where are monic polynomials of degree satisfying the divisibility chain
The polynomials are called the invariant factors of .
Intuition
The key conceptual move is to stop thinking of as a bare vector space and instead let the ring act on it, with acting as . This turns into an -module: for a polynomial , define . Since is a principal ideal domain, and since is finite-dimensional (so this module is torsion — by The Cayley-Hamilton Theorem the characteristic polynomial annihilates all of ), 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 , and one can always arrange the into a divisibility chain, uniquely. Translating back: each cyclic -module summand corresponds to a -invariant cyclic subspace of , on which acts by "multiplication by " in the quotient ring . If , the vectors form a basis, and in this basis is exactly the companion matrix : this is because is precisely the relation , read off column by column. The picture is a chain of matryoshka-doll invariants: the smallest polynomial divides the next , which divides , and so on, up to , which is exactly the minimal polynomial of . The product is the characteristic polynomial. This single chain of divisibility completely determines the similarity class of — no reference to eigenvalues, splitting fields, or algebraic closure is needed anywhere.
Let be a finite-dimensional vector space over a field , with , and let be a linear operator. Then:
(Existence.) There exist monic polynomials , each of degree , with
and a basis of , such that
Moreover equals the minimal polynomial of and equals the characteristic polynomial of .
(Uniqueness.) The sequence is uniquely determined by : if is another basis with for monic polynomials of degree , then and for all .
(Similarity criterion.) Two operators on finite-dimensional -vector spaces of the same dimension are similar if and only if they have the same invariant factors.
We give a self-contained proof using the -module structure on ; 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: as an -module. Make into a module over the polynomial ring by defining, for and ,
This is well-defined (it does not depend on how is written, since it is literally the polynomial evaluated at the fixed operator ) and one checks directly from the axioms of a vector space and linearity of that , , and . So is an -module. Because is finite-dimensional over , it is finitely generated as an -module (any -basis generates it). By The Cayley-Hamilton Theorem, the characteristic polynomial satisfies for every — that is, . Hence is a torsion -module: every element is annihilated by a nonzero polynomial. A submodule of under this action is precisely a -invariant subspace, since generates the whole module action, so closure under multiplication by and its powers is exactly -invariance.
Step 2: Cyclic submodules and their minimal polynomial. For , the cyclic submodule generated by is , matching the definition of -cyclic subspace given above. Consider the -module homomorphism
Its image is . Its kernel is an ideal of ; since is a PID and the kernel is nonzero (as and ), the kernel is for a unique monic polynomial of degree (degree since is nonzero when , so the kernel is a proper ideal). This is called the -annihilator of ; it is the monic generator of least degree with . By the first isomorphism theorem for modules,
as -modules, and in particular .
Step 3: The companion matrix arises from a cyclic basis. Suppose . Then is an -basis of (standard division algorithm: every coset has a unique representative of degree ), so under the isomorphism of Step 2, is an -basis of . Write . Since ,
Computing the matrix of in the basis : the -th basis vector maps under to , which is the -th basis vector for , and for maps to the combination above. Column by column this is exactly . So every cyclic subspace, in the natural basis generated by iterating 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 for suitable vectors with annihilators . We argue by strong induction on . If there is nothing to prove (empty sum). Suppose and the result holds for all smaller dimensions (for all operators on smaller spaces over ). Choose whose -annihilator has maximal degree among all -annihilators of vectors of ; such a exists since is finite-dimensional, so degrees of annihilators (all ) are bounded, and at least one vector (e.g. one realizing the minimal polynomial, see Step 5) attains the maximum. Claim: is divisible by the annihilator of every vector in . Suppose not: some has annihilator with . Let ; since , we have . Because is a PID, write for some (Bézout). One shows (a standard lemma in cyclic decomposition, proved by explicit construction) that there is a vector with -annihilator equal to : concretely, set
for a suitably chosen unit-scaling; more directly, one uses that and have coprime-complement structure to produce, via the Chinese Remainder Theorem applied inside and , a vector whose annihilator is , which has degree because . This contradicts maximality of . Hence no such exists, proving the claim: is annihilated-by-divisibility maximal, i.e. every vector's annihilator divides . In particular for all , so ; combined with maximality this shows is precisely the minimal polynomial of (it kills every vector, and no polynomial of smaller degree can, since some vector already requires degree ). Existence of a complement. We now show has a -invariant complement, i.e. a -invariant subspace with . This is the technical heart of the theorem; we outline the standard argument. Consider the quotient module , which is again a finitely generated torsion -module (a quotient of ), with induced operator . Since (as has dimension ), by the inductive hypothesis decomposes as a direct sum of -cyclic submodules generated by cosets with annihilators , and since annihilates all of (shown above) it certainly annihilates , so each . For each , lift to . The annihilator of in is a multiple of (annihilating forces annihilating ), say for some . A direct computation (using that , hence for some , and that forces a divisibility condition on since is maximal) shows one can correct by subtracting an element of to obtain whose -annihilator is exactly (matching the annihilator in the quotient exactly, with no jump in degree) — this uses maximality of essentially, to guarantee the needed element is divisible appropriately inside . 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 in hand, set . Since the images generate as a direct sum of cyclic pieces, and the project onto generators of the same pieces, maps onto under the quotient map, and one checks using that the sum of the in is direct. Hence
and moreover by construction. Relabel for . This exhibits the desired divisibility chain and completes the induction.
Step 5: Assembling the basis and the two polynomial identities. Taking, in each cyclic summand , the basis from Step 3 (where ), and concatenating these bases across , gives a basis of (since the sum is direct) in which
This proves existence. Since direct sums of -invariant subspaces give block-diagonal matrices, , and the characteristic polynomial of a block-diagonal matrix is the product of the characteristic polynomials of the blocks. A direct computation of by cofactor expansion along the first row shows the characteristic polynomial of is exactly itself. Hence
We showed in Step 4 that is the minimal polynomial of ; this proves the "moreover" clause.
Step 6: Uniqueness. Suppose and with both sequences monic, of degree , and forming divisibility chains. These give two cyclic -module decompositions of :
as -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 (equivalently, compare the -primary parts, i.e. the numbers -type pieces appearing via elementary divisors), or more directly compare, for each , the quotient
using the standard dimension-counting invariant: for each integer , the dimension of over , as ranges over irreducible factors, is determined by alone (it is , computable purely from the exponents of in each ), 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 , next-largest into , etc.) is the unique way to build an ascending divisibility chain from a given multiset of prime powers. Hence and for every . (This is exactly the argument by which the analogous uniqueness statement is proved for finite abelian groups, transported verbatim to -modules since is a PID just as is.)
Step 7: The similarity criterion. If and have the same invariant factors , then by existence there are bases for the domain of and for the domain of with . If and (with ), let send the -th vector of to the -th vector of ; then is an isomorphism with as matrices in the standard sense, i.e. and are similar (as matrices, for the change-of-basis matrix ). Conversely if and are similar, say for some bases (similarity is exactly agreement of matrix representations up to change of basis), then applying the existence part to and to 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 -modules), we get that and have isomorphic -modules, hence by uniqueness the same invariant factors.
∎Worked examples
Example 1. Let be given in the standard basis by
This matrix is already block-diagonal with three blocks: a Jordan-type block for eigenvalue , a block (a rotation, with characteristic polynomial , irreducible over ), and the block . We compute the characteristic polynomial directly:
Since is not diagonalizable over (indeed not even over , because of the genuine Jordan block for eigenvalue ) and has an irreducible quadratic factor, no eigenbasis over exists and no Jordan form over exists either — Jordan form here would require passing to , entirely changing the base field. To find the invariant factors we compute the minimal polynomial. The block for eigenvalue (size , one Jordan block) forces a factor in the minimal polynomial (a single Jordan block has minimal polynomial , not ). The rotation block contributes . The remaining block contributes , already divided out by . Hence
This is for the largest invariant factor. Since and , there must be one more invariant factor of degree with ; the only monic linear divisor of is . Since must equal , and indeed , this is consistent. So the invariant factors are
The Rational Canonical Form is :
and , a block-diagonal matrix. This is a genuinely different, and coarser, canonical form than Jordan form (which needs ): the Rational Canonical Form is computed and stated purely with rational coefficients.
Example 2. Let have matrix
This is already the companion matrix of ; let us instead read it off directly: comparing with the general companion form for , the last column is , so , giving
Let's instead choose numbers so the factorization is transparent, redefining the example cleanly: take
and build on as , i.e.
This matrix is diagonalizable (eigenvalues , all distinct, all rational — over 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 and are distinct, but expressed via companion blocks instead of scalar blocks. Its characteristic polynomial is , a product of four distinct linear factors, so the minimal polynomial equals the characteristic polynomial (distinct eigenvalues force this): . Since has degree , there is only one invariant factor:
The Rational Canonical Form is the single companion matrix
Notice this differs from the diagonal matrix that also represents in an eigenbasis — both are valid representations of , 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 companion blocks) works whether or not is diagonalizable, whether or not eigenvalues lie in .
Drop the requirement that , and only require the to be monic of degree with for some basis . Then uniqueness fails. Take on with minimal polynomial acting as the identity operator, i.e. . One valid decomposition (satisfying the divisibility chain, as required by the theorem) is
giving , matching the theorem's unique answer (three invariant factors, all equal to , trivially forming a chain since ). But without the chain requirement, one could equally decompose (as an -module, with acting as ) using different cyclic pieces of different apparent sizes by grouping basis vectors differently — for instance treat where the first is spanned by (a cyclic module isomorphic to ) and the piece, spanned by , could mistakenly be treated as if generated by a single vector under some non-identity dummy operator — this does not actually work for specifically (since restricted to any subspace acting as has every vector as its own cyclic generator with annihilator , forcing all cyclic pieces to be -dimensional). So to see the failure concretely, use instead on , i.e.
is not quite it either since is ; instead let on have invariant factors correctly computed as (a nilpotent operator, , , with a -dimensional kernel). The theorem's decomposition is , sizes ... to keep dimension , add one more trivial factor: take on nilpotent with Jordan blocks of sizes in the sense of appearing twice: invariant factors is not a valid chain-with-strict-content since it's allowed (), giving , a nilpotent matrix of rank . If instead one is sloppy and tries the unordered, non-divisibility-respecting decomposition omitted — i.e. using elementary divisors grouped as rather than correctly regrouped into the chain — one obtains , a genuinely different block structure (three blocks instead of two, of sizes instead of ) representing the same similarity class only if the underlying elementary divisor multiset actually matches; but for with true invariant factors , the elementary divisor multiset is , not , so represents a genuinely different, non-similar operator (one with a -dimensional kernel instead of -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 satisfy and are computed from the actual annihilators of a genuine direct-sum decomposition of , the representation is not a well-defined function of alone, and the uniqueness clause of the theorem fails outright.
Scope, edge cases & common traps
| Situation | What happens |
|---|---|
| has irreducible characteristic polynomial over | : RCF is a single companion matrix ; is a cyclic -module. |
| is diagonalizable with distinct eigenvalues in | Minimal and characteristic polynomials coincide (both equal ), so still: RCF is one companion block, even though is also diagonal in another basis. |
| not algebraically closed, has irreducible quadratic/higher factors | Jordan form does not exist over ; RCF still exists and is the only field-native canonical form. |
| nilpotent | Invariant factors are ; RCF and Jordan form (over any field, since always splits as ) coincide up to reordering of blocks. |
| Zero operator or identity operator on | , all (resp. ); RCF is separate companion blocks, i.e. the zero (resp. identity) matrix. |
The mathematician behind it
Origins · 1 minCamille Jordan and how this was discovered
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.
Why can't we always use the Jordan Canonical Form instead?
The Jordan Canonical Form requires the characteristic polynomial of to split into linear factors over the field of scalars, which forces the field to at least contain all the eigenvalues of . Over or a finite field, most operators have characteristic polynomials with irreducible factors of degree greater than one, so no basis over that field puts 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 ; this is always the largest invariant factor . Then use the fact that together with the divisibility chain and the theory of elementary divisors: factor and into irreducibles, determine the multiplicity of each irreducible factor appearing in each invariant factor by working through the Smith Normal Form of over the PID , whose nontrivial diagonal entries are exactly the invariant factors. In small examples, matching multiplicities of each irreducible power between and , 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 are uniquely determined by , and the block-diagonal arrangement in increasing order of degree is a single, specific matrix. The basis 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 are precisely the nontrivial diagonal entries of the Smith Normal Form of the matrix , viewed as a matrix over the polynomial ring . 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 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 is a scalar multiple of the identity, say on an -dimensional space, the invariant factors are copies of , and the Rational Canonical Form is literally itself, since as a 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 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 instead of correctly regrouping into ), and only the maximal, divisibility-respecting grouping is provably unique, as guaranteed by the uniqueness half of the structure theorem for modules over .
Can the Rational Canonical Form be computed for operators on infinite-dimensional spaces?
Not directly by this theorem, which explicitly requires to be finite-dimensional so that the minimal and characteristic polynomials exist and is a finitely generated, torsion -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.