2016-07-12 52 views
2

Ich bin neu in Mathematica, und ich versuche, wieMathematica Problem: Lösung Matrix Gleichung AX = lambdaBX Symbolisch

AX = \lambda BX 

Hier eine Matrix-Gleichung in einer Form zu lösen, A und B sind 4*4 Matrizen im Folgenden ist \lambda ein Wert, X ist die Eigenvektor-4*1 Matrix.

A = {{a1 + b1, c, d, f}, 
    {c, a2 + b2 , f , e}, 
    {d , f , a3 + b1 , c}, 
    { f, e , c, a4 + b2}} 

B = {{1, 0, 0 , 0}, 
    {0, 1 , 0 , 0}, 
    {0 , 0 , -1 , 0}, 
    {0, 0 , 0, -1}} 

Ich mag würde diese Matrix-Gleichung zu lösen und die symbolische Lösung für \lambda mit a1,a2,a3,a4,b1,b2,c,d,e,f usw.

Es wäre viel dankbar, wenn mir jemand sagen kann.

Mit freundlichen Grüßen,

mike

+0

Gleichungen der Form (A-lambda * I) X = 0 Eigenwertprobleme genannt. Wenn Sie Ihr Formular in diesem Formular neu anordnen, kann es einfacher sein, eine allgemeine Lösung zu finden. Premultiply beide Seiten von B (invers) und Sie haben es. – duffymo

+0

Danke, aber Matrix B ist keine Identitätsmatrix. Und es wäre dankbar, wenn Sie einige Befehle auf Mathematica schreiben könnten. – Mike22LFC

+0

Ich weiß, B ist keine Identitätsmatrix; genau deshalb habe ich vorgeschlagen, beide Seiten mit B (invers) zu multiplizieren, um es in die Eigenwertform zu bringen. Ich bin kein Mathematica-Benutzer, daher kann ich keine Befehle posten. – duffymo

Antwort

1

Wolfram: Matrix Computations See - speziell im Abschnitt 'Generalized Eigenwert'.

Für n × n Matrizen A, B die verallgemeinerte Eigenwerte sind die n Wurzeln seiner charakteristischen Polynoms, p() = det (A - B). Für jeden verallgemeinerten Eigenwert, λ ε λ (A, B), die Vektoren, , dass

A χ = χ λ B

beschrieben als verallgemeinerte Eigenvektoren entsprechen.

Beispiel unter Verwendung von symbolischen Werten:

matA = {{a11, a12}, {a21, a22}}; 
matB = {{b11, b12}, {b21, b22}}; 

Eigenvalues[{matA, matB}] 

{(1/(2 (B21 -B12 + b11 b22))) (a22 b11-b12 a21-a12 a11 + b21 b22-Sqrt [(-a22 b11 + a21 b12 + a12 b21-a11 b22)^2-4 (-a12 a21 + a11 a22) (-b12 b21 + b11 b22)]), (1/(2 (-b12 b21 + b11 b22))) (a22 b11-a21 b12-a12 b21 + a11 b22 + Sqrt [(-a22 b11 + a21 b12 + a12 b21-a11 b22)^2-4 (-a12 a21 + a11 a22) (-b12 b21 + b11 b22)])}

Eigenvectors[{matA, matB}] 

...

+0

Danke @Chris. Das hilft mir sehr. Für A χ = λ B χ, in Mathematica, kann ich Eigenwerte [{A, B}] verwenden. Ähnlich fand ich in Matlab [V, D] = eig (A, B). Nochmals vielen Dank, Prost. – Mike22LFC