Let $\mathcal{F}=\{A_1,A_2,\ldots,A_r\}$ be a triangulable commuting family of $n\times n$ matrices (that is, each $A_i$ is triangulable and $A_iA_j=A_jA_i$ for every $i,j$). I know that $\mathcal{F}$ can be simultaneous triangularization, but what is the algorithm of finding the invertible matrix $P$ such that $P^{-1}A_iP$ is triangular? As a working example consider the matrices $$ A= \begin{pmatrix} -3 & 2 & -4 \\ -1 & 0 & -1\\ 2 & -2 & 3 \end{pmatrix}\qquad B= \begin{pmatrix} 3 & -2 & 2 \\ -1 & 2 & -1\\ -2 & 2 & -1 \end{pmatrix}\qquad C= \begin{pmatrix} 1 & 0 & 1 \\ 1 & 0 & 1\\ 0 & 0 & 0 \end{pmatrix} $$ Here $AB=BA$, $AC=CA$ and $BC=CB$. In addition, the characteristic polynomials of these matrices are $$ f_A(x)=x(x-1)(x+1)\\ f_B(x)=(x-2)(x-1)^2\\ f_C(x)=x^2(x-1) $$ so each one of them is triangulable. Thanks!
-
The characteristic polynomial for $A$ seems wrong. – Emilio Novati May 24 '17 at 07:32
-
A typo..I fixed that ($-4$ in the 13 position). Thanks @Emilio Novati. – boaz May 24 '17 at 07:34
-
Do you know how to find the matrix $P$ for a single triangularizable matrix $A$? – Greg Martin May 24 '17 at 08:22
-
yes, i know how to do the trianglazation process with a single matrix. – boaz May 24 '17 at 09:03
2 Answers
Here, it is easy. Since $A$ has $3$ distinct eigenvalues and $B,C$ commute with $A$, we can deduce that $B,C$ are polynomials in $A$ and it suffices to triangularize $A$.
In the general case.
Step 1. Find a common eigenvector of $A,B,C$.
Step 2. Proceed by recurrence.
Moreover, you can choose $P$ as an orthogonal matrix.
-
Thanks @loup blanc. Two questions: (1) Does the matrices $A_i$ must have a common eigenvectors? (2) What do you mean by "Proceed by recurrence"? – boaz May 24 '17 at 09:08
-
Yes, they do have a common eigenvector $u$. Consider a basis in the form $u,\cdots$ and use the same reasoning about the "right-lower" blocks of the new matrices. – May 24 '17 at 09:16
-
More explicitly: find a common eigenvector v1 among your given matrices. Now complete a basis arbitrarily *u2,...,un. Let Q = [v,u1,...,un ] a matrix. Now Q-1AQ, Q-1BQ, and Q-1CQ should all have an (n-1) x (n-1) submatrix in the bottom right with only zeros in the entries of first column directly to the left of it. You then want to find a common eigenvector w among all of these submatrices. Let U =[u1,...,un]. Then v2 =Uw. With this update *Q = [v1,v2,u3,...,un ] * where u3,...,un may have been rechosen in order to complete the space. You rinse and repeat with n-2, n-3 and so on.
If you need more information check out the algorithm in the paper by Dubi:
