Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| documentation:language_reference:objects:matrix:functions:eigensystem [2024/12/12 17:27] – Maurits W. Haverkort | documentation:language_reference:objects:matrix:functions:eigensystem [2025/05/06 14:11] (current) – Maurits W. Haverkort | ||
|---|---|---|---|
| Line 1: | Line 1: | ||
| {{indexmenu_n> | {{indexmenu_n> | ||
| - | ====== Eigensystem ====== | + | ====== |
| ### | ### | ||
| Line 90: | Line 90: | ||
| <code Quanty Example.Quanty> | <code Quanty Example.Quanty> | ||
| A = Matrix.New({{1, | A = Matrix.New({{1, | ||
| - | {1-2E-10,3,7}, | + | {1-1E-7,3,5}, |
| {3,5,1}}) | {3,5,1}}) | ||
| val, funL, funR = Eigensystem(A) | val, funL, funR = Eigensystem(A) | ||
| Line 98: | Line 98: | ||
| print(" | print(" | ||
| - | print(" | + | print(" |
| </ | </ | ||
| Line 107: | Line 107: | ||
| With potential complex eigenvalues | With potential complex eigenvalues | ||
| The eigenvalues are | The eigenvalues are | ||
| - | { -4.5383431004156 | + | { -3.7873476689872 |
| The left eigenfunctions are | The left eigenfunctions are | ||
| - | { { -0.3497 , -0.4901 , 0.8089 } , | + | { { -0.3763 , -0.5129 , 0.7715 } , |
| - | { -0.924 , | + | { -0.853 , |
| - | { -0.3429 , -0.6415 , -0.6945 } } | + | { -0.3616 , -0.6854 , -0.632 |
| The right eigenfunctions are | The right eigenfunctions are | ||
| - | { { -0.2791 , -0.6392 , 0.7283 } , | + | { { -0.3763 , -0.5129 , 0.7715 } , |
| - | { -0.8588 , 0.5199 , -0.0562 } , | + | { -0.853 |
| - | { -0.3174 , -0.7492 , -0.5912 } } | + | { -0.3616 , -0.6854 , -0.632 |
| The matrix transformed to a diagonal matrix by its eigenfunctions is | The matrix transformed to a diagonal matrix by its eigenfunctions is | ||
| - | { { -4.5383 , 0 , 0 } , | + | { { -3.7873 , 0 , 0 } , |
| - | { 0 , 0.591 | + | { 0 , 0.6489 , 0 } , |
| - | { 0 , 0 , 8.9473 } } | + | { 0 , 0 , 8.1385 } } |
| - | Be carefull, using only the left handed eigenvectors is not sufficient to diagonalize the matrix | + | The left and right hand vectors are now only marginally different |
| - | { { -4.6149 , -0.562 | + | { { -3.7873 , -5.4e-8 , |
| - | { | + | { |
| - | { -1.1403 , 0.598 | + | { -4.9e-9 , -8.4e-8 , |
| </ | </ | ||