Jordan Normal Form physical applications

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 4K views
Maybe_Memorie
Messages
346
Reaction score
0
Are there any applications of the Jordan Normal Form of a matrix in physics?
If so, please explain?
 
Physics news on Phys.org
The Jordan Normal Form is very often used to solve linear differential equations. Such an equations can often be represented in the form

[tex]y^\prime=Ay,~y(0)=y_0[/tex]

with A a matrix. The solution to such an equation is [tex]y(t)=e^{At}y_0[/tex]. The problem is to calculate the matrix exponential [tex]e^{At}[/tex], this can be very difficult. To ease the problem, we calculate the Jordan Normal Form of A. In this from, calculating the matrix exponential becomes very easy!
The fun part is of course to calculate the Jordan Normal Form, as this can be quite difficult...
 
A slight variation on what micromass said:
If you have the d.e. y'= Ay, and B is such that [itex]B^{-1}AB= J[/itex], the "Jordan form" for A, then we can multiply both sides of the equation by [itex]B^{-1}[/itex]:
[tex](B^{-1}y)'= B^{-1}Ay= B^{-1}A(BB^{-1})y= (B^{-1}AB)(B^{-1}y)[/tex]
and, letting [itex]x= B^{-1}y[/itex] write the equation as
[tex]x'= Jx[/itex].<br /> <br /> Now, since J is a Jordan form matrix, that gives a system of "almost uncoupled" equations. If we write x as a column matrix<br /> [tex]\begin{bmatrix}x_1 \\ x_2 \\ ---\\ x_n\end{bmatrix}[/tex]<br /> <br /> Then every equation is of the form [itex]x_k'= \lambda_kx_k[/itex] or [itex]x_k'= \lambda_nx_k+ x_{k+1}[/itex]. The very last equation is, of course, [itex]x_n'= \lambda_n x_n[itex], because there is no "[itex]x_{n+1}[/itex]" and is easy to solve. Then, if [itex]x_{n-1}'= \lambda_{n-1}x_{n-1}+ x_n[/itex], because we already have [itex]x_n[/itex], that is easy to solve, etc.<br /> <br /> Finally, since [itex]x= B^{-1}y[/itex], [itex]y= Bx[/itex].<br /> <br /> If a matrix happens to be "diagonalizable" (all those "1" above the diagonal are unecessary), we can completely "uncouple" all those equations, whether differential equations or not, and have, basically, n separate, simpler, problems. If the matrix is NOT diagonalizable, we can still simplify as much as possible using the "Jordan Normal Form".[/itex][/itex][/tex]