Matriks eksponensial dari matriks Hamilton

10

Biarkan menjadi nyata, kotak, matriks padat. G dan Q simetris. MembiarkanA,G,QGQ

H=[AGQAT]

menjadi matriks Hamilton. Saya ingin menghitung eksponensial matriks . Saya membutuhkan eksponensial matriks penuh, e t H , tidak hanya produk matriks-vektor. Apakah ada algoritma atau perpustakaan khusus yang tersedia untuk menghitung eksponensial dari matriks Hamilton?HetH

Max Behr
sumber
2
Apakah Anda ingin matriks eksponensial itu sendiri, atau apakah Anda benar-benar hanya ingin memecahkan ODE ? z˙=Hz
Daniel Shapero
Saya Membutuhkan Matriks Eksponensial itu sendiri. Tapi ekuivalen saya bisa memecahkan ODE . Z˙=HZ, Z(0)=I
Max Behr
2
Struktur yang mempertahankan eigensolvers Benner dapat menangani transformasi kemiripan untuk memudahkan komputasi eksponensial matriks.
perkusi
1
@RichardZhang Cara brutal adalah dekomposisi QZ. Periksa misalnya mulai dari link.springer.com/article/10.1007/s002110050315 untuk detail lebih lanjut.
perkusi
2
Makalah 19 Cara Meragukan untuk Menghitung Eksponensial Matriks, 25 Tahun Kemudian mencakup banyak cara buruk (dan beberapa cara yang baik) untuk menghitung eksponensial matriks. Ini tidak spesifik untuk masalah Hamiltonian tetapi tetap sangat berharga jika Anda bekerja pada masalah semacam ini.
Daniel Shapero

Jawaban:

6

Jawaban yang sangat cepat ...

Eksponensial dari matriks Hamilton adalah symplectic, properti yang mungkin ingin Anda pertahankan, jika tidak, Anda hanya akan menggunakan metode pelestarian non-struktur. Memang, tidak ada keuntungan kecepatan nyata dalam menggunakan metode terstruktur, hanya pelestarian struktur.

Cara yang mungkin untuk menyelesaikan masalah Anda adalah sebagai berikut. Pertama menemukan matriks symplectic sehingga H = M - 1 H M = [ A - G 0 - Sebuah T ] adalah Hamiltonian dan memblokir segitiga atas, dan A memiliki nilai eigen di babak-pesawat kiri. Anda mendapatkan matriks ini misalnya dengan mengambil [ I 0 X I ] , di mana X memecahkan persamaan Riccati yang terkait dengan HH^=M1HM=[A^G^0A^T]A^[I0XI]XH, Atau (lebih stabil karena itu ortogonal) oleh penataan kembali Schur dekomposisi dan menerapkan trik Laub (yaitu, menggantikan faktor Schur kesatuan [ U 11 U 12 U 21 U 22 ] dengan [ U 11 - U 12 U 12 U 11 ] ). Anda mungkin mengalami kesulitan melakukannya jika Hamiltonian memiliki nilai eigen pada sumbu imajiner, tetapi itu adalah cerita yang panjang dan untuk saat ini saya kira itu tidak terjadi pada masalah Anda.H[U11U12U21U22][U11U12U12U11]

Setelah Anda memiliki , Anda memiliki exp ( H ) = M exp ( H ) M - 1 , dan Anda dapat menghitung exp ( H ) = [ exp ( A ) X 0 exp ( - A T ) ] , di mana X memecahkan persamaan Lyapunov tertentu, saya percaya sesuatu seperti A X + X A TMexp(H)=Mexp(H^)M1

exp(H^)=[exp(A^)X0exp(A^T)],
X (tanda-tanda mungkin salah; memaksakan exp ( H ) H = H exp ( H ) dan memperluas blok untuk mendapatkan persamaan yang benar. Cari "metode Schur-Parlett" untuk referensi trik ini).
A^X+XA^T=exp(A^)G^G^exp(A^T)
exp(H^)H^=H^exp(H^)

Kemudian ketiga faktor tersebut bersifat simplektik. Cukup gunakan secara terpisah: jangan menghitung produk atau Anda akan kehilangan properti ini secara numerik.

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

H

AGQHHHAGQAGQ berasal dari persamaan integral yang juga akan menjelaskan struktur padat dan potensi kompresi (tergantung pada kernel).

(HλI)1HHAGQ

H

HH

Kelemahan dari pendekatan ini:

  • AGQ
  • tidak mengambil keuntungan dari struktur Hamilton

Positif:

  • representasi terkompresi dari matriks eksponensial, meskipun masih berupa matriks, bukan hanya cara untuk melakukan MVP
  • kompleksitas linear-logaritmik (asalkan ada asumsi peringkat rendah)
  • perpustakaan mungkin mengambil keuntungan dari transposisi dan simetri dalam blok
Anton Menshov
sumber