I am trying to understand how do we form a Jordan Normal Basis and more specifically why the chains are linearly independent. I saw some relevant question, but still, I can't figure that out.
We start with a highest order eigenvector and form its Jordan chain. Then, if we haven't had yet a basis, we take another eigenvector of highest possible order and form its Jordan chain and we continue like this. Why will two jordan chains we obtain be linearly independent?
Thanks in advance
edit:In this post, he says:
The other method is to first determine the lengths of the chains (there are formulas you can use to find them). If the longest chains have length kk, then you find linearly independent vectors in $(T−\lambda I)^{k−1}$ that are not in $(T−\lambda I)^{k−2}$. Then you apply $(T−\lambda I)$ to them, and if necessary find linearly independent vectors from the ones you have so far that are in $(T−\lambda I)^{k−2}$ but not in $(T−\lambda I)^{k−3}$. You continue this way, always selecting linearly independent vectors, until you are down to finding enough linearly independent vectors in $(T−\lambda I)$ besides the ones you already have; you will be guaranteed they are linearly independent because at each stage we are always selecting vectors that are linearly independent.
My question is actually, how do we know that after applying $(T−\lambda I)$ to the vectors we already have selected, we get a linearly independent set of vectors?