Wie kann ich eine Basis für eine Matrix-Lie-Algebra bei einer endlichen Menge von Generatoren berechnen?

Gegeben eine beliebige Menge von (numerisch) quadratisch komplexer Matrices , ich bin an der Berechnung der reellen Matrix Liealgebra von erzeugten , rufen sie . Das heißt, würde ich eine Grundlage für gleiche wobei definiert ist als rekursiv $\mathcal{A}=\{A_1,A_2,\cdots,A_m\}$ $\mathcal{A}$ $\mathcal{L_\mathcal{A}}$

L_{A} = s p a n_{R} {B : B \in \cup_{k = 1}^{\infty} C_{k}}

$\mathcal{L_\mathcal{A}} = \mathbb{span_R}\{B:B\in\cup_{k=1}^{\infty}\mathcal{C}_k\}$

C_{k}

$\mathcal{C}_k$

, und

für

C_{1} = A

$\mathcal{C_1}=\mathcal{A}$

C_{k + 1} = {[X, Y] : X, Y \in \cup_{j = 1}^{k} C_{j}}

$\mathcal{C_{k+1}}=\{[X,Y]:X,Y\in\cup_{j=1}^k\mathcal{C_j}\}$

k \geq 1

$k\geq 1$

Diese Berechnung kommt in der (Quanten-) Steuerungstheorie vor.

Derzeit verwende ich eine hier gefundene Methode, die nur wiederholte Lie-Klammern durchsucht (dh solche der Form ) und wird garantiert beendet. Ich bin jedoch interessiert zu wissen, ob es andere (schnellere) Methoden gibt. Vielleicht mit P. Hall Basen? Vielleicht ein rekursiver Algorithmus? Meine derzeitige Standardsprache ist Matlab. $[A_{j_1},[A_{j_2},[A_{j_3},\cdots[A_{j_{n-1}},A_{j_n}]\cdots]]]$

linear-algebra matlab basis-set Ian Hincks
quelle

Ich vermute, dass Ihre ursprünglichen Generatoren hermitisch sind. Ist das wahr? Wenn ja, würde ich mir vorstellen, dass der erste Schritt darin besteht, die Eigenräume der Generatoren zu vergleichen, da Kommutatoren nur dann ungleich Null sind, wenn sich die Eigenräume unterscheiden.

Jack Poulson

@JackPoulson Ja, die A's stammen von Hamiltonianern und sind daher schief-hermitisch (nicht hermitisch, weil sie mit dem i in Schrödingers Gleichung multipliziert werden). Ich bin mir nicht sicher, warum das ein guter erster Schritt wäre. Wäre es nicht schneller, die Kommutatoren zu berechnen und zu überprüfen, ob sie nicht Null sind, als mit Eigenräumen herumzuspielen?

Ian Hincks

Für eine einzelne Ebene von Kommutatoren wahrscheinlich ja. Es gibt jedoch eine kombinatorische Explosion, wenn Sie mehrere Kommutatorebenen in Betracht ziehen. Ich kenne keinen Algorithmus, aber normalerweise ist es eine gute Idee, so viel Struktur wie möglich auszunutzen. Ich würde sorgfältig darüber nachdenken, ob Sie andere Eigenschaften kennen, die auch Ihre Generatoren betreffen.

Jack Poulson

Wie kann ich eine Basis für eine Matrix-Lie-Algebra bei einer endlichen Menge von Generatoren berechnen?

Antworten: