View as the matrix of a linear map of matrix rank . Choose vectors whose images form a basis of , extend them by a basis of to a basis of , and extend to a basis of . In these two bases the matrix of is the rank normal form
The two changes of basis give invertible with the displayed matrix equal to .
For the block upper triangular matrix , every nonzero term in the Leibniz formula for determinants sends the rows belonging to into the columns belonging to ; bijectivity then sends the remaining rows into the columns. The permutation sum consequently factors into the determinant sums for and , proving
Finally suppose , so : the map is an intertwining operator. If lies in the generalized eigenspace of for , then for some ,
Because and have no common eigenvalue, is invertible, hence . The generalized eigenspaces of span , so . Thus the Sylvester equation operator is injective.
Since and ,
Thus the first column of is , and is a block upper triangular matrix:
where is its bottom-right submatrix. Therefore
A similarity transformation preserves the characteristic polynomial, so
with algebraic multiplicities included.
To construct , use a Householder transformation. With , choose the sign of to avoid cancellation and set
Then is an orthogonal matrix and maps to up to the chosen sign. If already lies on the first coordinate axis, take a suitable diagonal sign matrix. This is the one-vector case of orthogonal coordinate reduction of a subspace.