Compute the Smith reduction of a 2×2 polynomial matrix:

Compute the Smith reduction of a 2×3 polynomial matrix:

The matrix can have rational coefficients:

The matrix can have complex coefficients:

The matrix can have approximate coefficients:

The matrix can have symbolic coefficients:

When a matrix has rational function entries, PolynomialSmithReduce will compute the Smith–McMillan form, where elements of s have a denominator:

The second element of the Smith decomposition is a matrix with diagonal elements that are scalar multiples of the Smith normal form diagonal elements given by PolynomialSmithReduce.

Compute the Smith reduction of a 2×2 polynomial matrix:

Compute the decomposition and extract the second part:

The last element in the diagonal of the Smith normal form of a characteristic matrix gives the minimal polynomial:

Compare to MatrixMinimalPolynomial:

The polynomial Smith reduction can be used to compute the rational canonical form of a matrix; this is the second element in a Frobenius decomposition.

Compute the Smith form of a characteristic matrix and extract the diagonal:

Use this to create companion matrix blocks:

Use SparseArray and Band to reconstruct the companion matrix for m as a diagonal block matrix:

Check the result: