Matriz exponencial de una matriz hamiltoniana


10

Deje ser matrices reales, cuadradas y densas. G y Q son simétricos. DejarA,G,QGQ

H=[AGQAT]

ser una matriz hamiltoniana. Quiero calcular la exponencial de la matriz . Necesito la matriz exponencial completa, e t H , no solo el producto matriz-vector. ¿Hay algún algoritmo especializado o biblioteca disponible para calcular la exponencial de una matriz hamiltoniana?HetH


2
¿Desea la matriz exponencial en sí misma, o realmente solo quiere resolver la ODE ? z˙=Hz
Daniel Shapero

Necesito la matriz exponencial en sí misma. Pero lo que es equivalente i puede resolver el ODE . Z˙=HZ, Z(0)=I
Max Behr

2
La estructura de Benner que preserva los eigensolvers puede manejar la transformación de similitud para facilitar el cálculo exponencial de la matriz.
percusión

1
@RichardZhang La forma brutal es la descomposición QZ. Consulte, por ejemplo, a partir de link.springer.com/article/10.1007/s002110050315 para obtener más detalles.
percusión

2
El documento 19 Formas dudosas de calcular el exponencial de una matriz, 25 años después, cubre muchas formas malas (y algunas buenas) de calcular el exponencial de la matriz. No es específico de los problemas de Hamilton pero es muy valioso si está trabajando en este tipo de problemas.
Daniel Shapero

Respuestas:


6

Respuesta muy rápida ...

La exponencial de una matriz hamiltoniana es simpléctica, una propiedad que probablemente desee preservar, de lo contrario simplemente usaría un método que no preserva la estructura. De hecho, no existe una ventaja de velocidad real al usar un método estructurado, solo la preservación de la estructura.

Una posible forma de resolver su problema es la siguiente. En primer lugar encontrar una matriz simpléctica tal que H = M - 1 H M = [ A - G 0 - A T ] es hamiltoniano y bloquear triangular superior, y A tiene valores propios en el semiplano de la izquierda. Obtiene esta matriz, por ejemplo, tomando [ I 0 X I ] , donde X resuelve la ecuación de Riccati asociada a HH^=M1HM=[A^G^0A^T]A^[I0XI]XH, o (más estable ya que es ortogonal) reordenando la descomposición de Schur de y aplicando el truco de Laub (es decir, reemplazando el factor de Schur unitario [ U 11 U 12 U 21 U 22 ] con [ U 11 - U 12 U 12 U 11 ] ). Es posible que tenga problemas para hacerlo si el Hamiltoniano tiene valores propios en el eje imaginario, pero esa es una larga historia y por ahora supongo que no sucede en su problema.H[U11U12U21U22][U11U12U12U11]

Mexp(H)=Mexp(H^)M1

exp(H^)=[exp(A^)X0exp(A^T)],
X
A^X+XA^T=exp(A^)G^G^exp(A^T)
exp(H^)H^=H^exp(H^)

Entonces los tres factores son exactamente simplécticos. Simplemente utilícelos por separado: no calcule el producto o perderá esta propiedad numéricamente.


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

1

H

AGQHHHAGQAGQ provienen de una ecuación integral que también explicará su estructura densa y el potencial de compresión (dependiendo del núcleo).

(HλI)1HHAGQ

H

HH

Desventajas de este enfoque:

  • AGQ
  • no aprovecha la estructura hamiltoniana

Positivos:

  • la representación comprimida de la matriz exponencial, aunque sigue siendo una matriz, no solo una forma de hacer un MVP
  • Complejidad logarítmica lineal (siempre que exista el supuesto de bajo rango)
  • la biblioteca puede aprovechar la transposición y la simetría en los bloques
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.