Eigenraumbasis kontinuierlich abhängig von Parametern

8

Ich habe eine hermitische Matrix die von zwei Parametern abhängt, z. B. x und y . Wenn ich es an zwei engen Punkten ( x 1 , y 1 ) und ( x 2 , y 2 ) diagonalisiere, erhalte ich zwei enge Eigenwerte ( ε 1 und ε 2 ) und zwei entsprechende Eigenräume ( S 1 und S 2 ) derselben Dimension .Hxy(x1,y1)(x2,y2)ε1ε2S1S2

Beachten Sie, dass es sich nicht um Eigenwerte derselben Matrix handelt. Es gibt zwei verschiedene Matrizen: und H 2 = H ( x 2 , y 2 ) .H1=H(x1,y1)H2=H(x2,y2)

Ich habe ein Netz von Punkten und möchte den Eigenwert und den Eigenraum an jedem Punkt durch Interpolation finden. Das Problem ist, dass die Basen von S 1 und S 2 völlig unabhängig sind , da die Matrizen numerisch diagonalisiert sind. Selbst wenn ( x 1 , y 1 ) und ( x 2 , y 2 ) sehr nahe beieinander liegen, können die Basisvektoren sehr unterschiedliche Komponenten haben.(xi,yi)S1S2(x1,y1)(x2,y2)

Für die Interpolation benötige ich eine Basis, die kontinuierlich von und y abhängt , dh je näher die Eigenräume S 1 und S 2 sind, desto näher sollten die Basisvektoren sein.xyS1S2

Wenn und S 2 Ebenen im dreidimensionalen euklidischen Raum sind, besteht eine gute Möglichkeit, eine Basis in S2 auszuwählen, darin, die Basis von S1 um die Linie zu drehen, die den Schnittpunkt der Ebenen darstellt. Gibt es etwas Analoges dazu im komplexen mehrdimensionalen Raum?S1S2

Maksim Zholudev
quelle

Antworten:

7

t

Damit Sie kontinuierliche Eigenräume haben können, müssen Sie davon ausgehen, dass sich die zugehörigen Eigenwerte nicht annähernd kreuzen. (Bei nahezu sich kreuzenden Eigenwerten kann es sehr gut vorkommen, dass die Eigenräume im Wesentlichen ausgetauscht werden, obwohl sich die Eigenwertkurven nicht berühren. Dies tritt regelmäßig bei symmetrischen parametrischen Eigenwertproblemen auf und wird als Phänomen vermiedener Kreuzungen bezeichnet.)

H(t)λ(t)Qkt=tkQ(t)=Q0+jΦj(t)ZjΦj(t)Q0QkZjYkQkYkjΦj(tk)ZjQ0

Ersetzen Sie mit 2 Parametern die B-Splines durch FEM-Formfunktionen. Ihr Problem der kleinsten Quadrate könnte jetzt groß werden, und es sind geeignete zusätzliche Tricks erforderlich, um das Problem lösbar zu machen, wenn die direkte Lösung nicht durchführbar ist.

Arnold Neumaier
quelle
ΦjYkQk
Ykk