Matrix exponentiell einer Hamiltonschen Matrix

10

Sei reelle, quadratische, dichte Matrizen. G und Q sind symmetrisch. LassenA,G,QGQ

H=[AGQAT]

sei eine Hamiltonsche Matrix. Ich möchte die Matrix Exponential von berechnen . Ich brauche das Exponential der Vollmatrix, e t H , nicht nur das Matrixvektorprodukt. Gibt es spezielle Algorithmen oder Bibliotheken, um das Exponential einer Hamilton-Matrix zu berechnen?HetH

Max Behr
quelle
2
Möchten Sie die Exponentialmatrix selbst oder möchten Sie wirklich nur die ODE lösen ? z˙=Hz
Daniel Shapero
Ich brauche das Matrix Exponential selbst. Aber äquivalent kann ich die ODE lösen . Z˙=HZ, Z(0)=I
Max Behr
2
Benners strukturerhaltende Eigensolver können sich mit der Ähnlichkeitstransformation befassen, um die exponentielle Matrixberechnung zu vereinfachen.
Percusse
1
@RichardZhang Der brutale Weg ist die QZ-Zerlegung. Weitere Informationen finden Sie beispielsweise unter link.springer.com/article/10.1007/s002110050315 .
Percusse
2
Das Papier 19 Zweifelhafte Methoden zur Berechnung des Exponentials einer Matrix, 25 Jahre später, behandelt viele schlechte (und einige gute) Methoden zur Berechnung des Exponentials der Matrix. Es ist nicht spezifisch für Hamilton-Probleme, aber dennoch sehr wertvoll, wenn Sie an solchen Problemen arbeiten.
Daniel Shapero

Antworten:

6

Sehr schnelle Antwort ...

Das Exponential einer Hamilton-Matrix ist symplektisch, eine Eigenschaft, die Sie wahrscheinlich beibehalten möchten, andernfalls würden Sie einfach eine nicht strukturerhaltende Methode verwenden. In der Tat gibt es keinen wirklichen Geschwindigkeitsvorteil bei der Verwendung einer strukturierten Methode, sondern nur die Erhaltung der Struktur.

Eine Möglichkeit, Ihr Problem zu lösen, ist die folgende. Zuerst findet eine symplektischer Matrix , so daß H = M - 1 H M = [ A - G 0 - A T ] Hamilton - Operator ist und oberes Dreieck Block und A hat Eigenwert in der linken Halbebene. Sie erhalten diese Matrix zum Beispiel, indem Sie [ I 0 X I ] nehmen , wobei X die mit H verbundene Riccati-Gleichung löstH^=M1HM=[A^G^0A^T]A^[I0XI]XHoder (stabiler, da es orthogonal ist) durch Neuordnung der Schur-Zerlegung von und Anwenden des Laub-Tricks (dh Ersetzen des einheitlichen Schur-Faktors [ U 11 U 12 U 21 U 22 ] durch [ U 11 - U 12 U 12 U 11 ] ). Möglicherweise haben Sie Probleme damit, wenn der Hamilton-Operator Eigenwerte auf der imaginären Achse hat, aber das ist eine lange Geschichte, und ich gehe davon aus, dass dies bei Ihrem Problem nicht der Fall ist.H[U11U12U21U22][U11U12U12U11]

Sobald Sie haben , haben Sie exp ( H ) = M exp ( H ) M - 1 , und kann man berechnen exp ( H ) = [ exp ( A ) X 0 exp ( - A T ) ] , wobei X löst eine bestimmte Lyapunov - Gleichung, ich so etwas wie glauben A X + X A TMexp(H)=Mexp(H^)M1

exp(H^)=[exp(A^)X0exp(A^T)],
X (Zeichen können falsch sein; verhängen exp ( H ) H = H exp ( H ) und erweiternblockieren die richtige Gleichung zu erhalten. Suchen Sie nach "Schur-Parlett-Methode", um einen Hinweis auf diesen Trick zu erhalten.
A^X+XA^T=exp(A^)G^G^exp(A^T)
exp(H^)H^=H^exp(H^)

Dann sind die drei Faktoren genau symplektisch. Verwenden Sie sie einfach separat: Berechnen Sie das Produkt nicht, sonst verlieren Sie diese Eigenschaft numerisch.

Federico Poloni
quelle
1
HH~=[A^G0A^T]XLA^XL+XLA^T=GM2=[IXL0I]H^H^^A^AT^
1

H

AGQHHHAGQAGQ stammen aus einer Integralgleichung, die auch ihre dichte Struktur und ihr Kompressionspotential (abhängig vom Kernel) erklärt.

(HλI)1HHAGQ

H

HH

Nachteile dieses Ansatzes:

  • AGQ
  • nutzt die Hamiltonsche Struktur nicht aus

Positiv:

  • Die komprimierte Darstellung der Matrix ist exponentiell, obwohl es sich immer noch um eine Matrix handelt, und nicht nur um einen MVP
  • linear-logarithmische Komplexität (vorausgesetzt, die Annahme mit niedrigem Rang ist vorhanden)
  • Die Bibliothek kann die Transposition und Symmetrie in den Blöcken nutzen
Anton Menschow
quelle