Hier ich erklär dir mal, wie man das richtig macht. Die lassen dich doch dumm sterben.
Die wichtigste Erkenntnis; Diktat für Formelsammlung, Regelheft und Spickzettel:
" Jede Matrix löst ihre eigene Säkulardeterminante ( SD ) "
Fürdiagonalisierbare ===> halbeinfache Matrizen lässt sich das ja trivial einsehen; aber es gilt eben allgemein.
Und wie findet man die SD? Komisch; in den BNüchern ist das immer so Mega kompliziert erklärt. Du machst erst mal den quadratischen Ansatz
p_M ( x ) = x ² - p x + q ( 1 )
p und q folgen über Vieta das geschmähte Stiefkind.
p = E1 + E2 = Sp ( M ) ( 2a )
q = E1 E2 = det ( M ) ( 2b )
p_M ( x ) = x ² - x Sp ( M ) + det ( M ) ( 2c )
Jetzt tust du die Matrix M einsetzen in ihr Eigenpolynom ( 2c )
M ² - M Sp ( M ) + det ( M ) * 1| = 0 | * M ^ - 1 ( 3a )
M - 1| * Sp ( M ) + M ^ - 1 det ( M ) = 0 ( 3b )
M ^ - 1 = [ 1 / det ( M ) ] [ Sp ( M ) * 1| - M ] ( 3c ) ( 3c )