Matrix eingeschränkt auf einen Kern/eine Basis/einen Raum, der durch eine Basis aufgespannt wird |
31.07.2011, 19:12 | stevewilson | Auf diesen Beitrag antworten » | ||||
Matrix eingeschränkt auf einen Kern/eine Basis/einen Raum, der durch eine Basis aufgespannt wird in meinem Skript zu DGL steht folgende Definition: Ak und A sollen Matrizen sein. Dieses N ist im Endeffekt ein Kern einer Matrix. Also ist wohl gemeint: Eine Einschränkung der Matrix auf eine Basis bzw. auf einen Raum, der durch eine Basis eingespannt wird. Was kann das bedeuten? Wir können ja mal ein Beispiel nehmen: eingeschränkt auf: Kann mir jemand sagen, wie ich das ausrechne? |
||||||
01.08.2011, 15:51 | Mazze | Auf diesen Beitrag antworten » | ||||
Was ich mir vorstellen könnte ist, dass man sich das Bild der Matrix anschaut, wenn man das Urbild auf die entsprechende Menge einschränkt. Man kann ja jede Matrix mit einer linearen Abbildung beschreiben, und wenn man eine Abbildung auf eine Menge einschränkt wissen wir ja, wie man damit umgeht. |
||||||
01.08.2011, 18:38 | stevewilson | Auf diesen Beitrag antworten » | ||||
Mittels Spektralzerlegung soll das gehen: http://de.wikipedia.org/wiki/Schurzerlegung Aber wie geht das dann bei einer Matrix, bei der gilt? Z.b. bei der Matrix |
||||||
01.08.2011, 23:50 | Mazze | Auf diesen Beitrag antworten » | ||||
Möchtest Du jetzt die Schurform einmal durchrechnen? Ich sehe nicht wie das mit dem Ausgangsproblem zu tun hat. Kannst Du eventuel noch etwa smehr dazu sagen (in welchem Bereich Du dich befindest, was die eigentliche Aufgabe ist etc.) |
||||||
02.08.2011, 01:29 | stevewilson | Auf diesen Beitrag antworten » | ||||
So, jetzt muss ich erstmal gedankliches Backtracing machen... Also... Es gibt eine Lösungsformel für DGL-Systeme 1.Ordnung, die da lautet: Sofern man ein inhomogenes AWP hat: Mein DGL-System lautet: Anfangswerte: Nun muss ich ja exp(tA) berechnen. Doch A ist nicht nilpotent. Daher zerlege ich mittels Spektralzerlegung die Matrix und bekomme einen nilpotenten Teil und kann somit exp(tA) berechnen... So machen die das in diesem Skript: http://www.ifam.uni-hannover.de/~ehrnstr...ta/DE_short.pdf Seite 11 unten, bis Seite 13 Mitte... |
||||||
03.08.2011, 09:46 | Mazze | Auf diesen Beitrag antworten » | ||||
In dem Fall führt die Schurzerlegung zum Ziel. Du spaltestest dann das Ergebnis in eine Summe aus Diagonalmatrix und nilpotenter Matrix auf. So, der einzige doppelte Eigenwert ist die 1. Ein zugehörige Eigenvektor wäre (1,1). Als erstes geben wir einen Vektor an, so dass {(1,1),v} eine orthogonale Basis ist. Danach werden die Vektoren normiert und man erhält tatsächlich wobei D eine obere Dreiecksmarix ist. Rechne mal deinen Weg vor, weil ich komme auf eine obere Dreiecksmatrix. |
||||||
Anzeige | ||||||
|
||||||
03.08.2011, 10:33 | stevewilson | Auf diesen Beitrag antworten » | ||||
Du meinst wohl Wie geht es dann weiter? Richtig? |
||||||
03.08.2011, 10:56 | stevewilson | Auf diesen Beitrag antworten » | ||||
oder und ??? |
||||||
03.08.2011, 11:19 | stevewilson | Auf diesen Beitrag antworten » | ||||
ich meinte: und |
||||||
03.08.2011, 12:49 | Mazze | Auf diesen Beitrag antworten » | ||||
Der algorithmus für die Schurform sieht vor, dass der Eigenvektor zu einer orthonormalen Basis ergänzt wird. Die Vektoren (1,1) und (1,0) sind zwar eine Basis, aber nicht orthogonal. Nach dem Algorithmus hätte man zum Beispiel (1,1) und (1,-1) betrachtet (und normiert). Unabhängig davon funktioniert es für deine Matrix. Wichtig zu wissen ist : Das heißt, wenn D eine obere Dreiecksmatrix ist mit , kannst Du leicht berechnen, in dem Du berechnest. hast Du mit richtig bestimmt. |
||||||
03.08.2011, 13:06 | stevewilson | Auf diesen Beitrag antworten » | ||||
Achso, jetzt weiß ich was eine orthonormale Basis ist. Für mein habe ich raus: Richtig? Wie setze ich das jetzt in die u(t)-Gleichung mit dem Integral ein, die ich oben gepostet habe? |
||||||
03.08.2011, 13:29 | Mazze | Auf diesen Beitrag antworten » | ||||
In dem Integral hast Du ja ein Matrix-Vektor-produkt. Heraus kommt also ein Vektor als Integrand. Das Integral eines Vektors ist definiert als das Integral der Komponenten des Vektors (sofern Du nicht über einen Vektor integrierst, was hier nicht der Fall ist). Deine Matrix habe ich jetzt nicht überprüft, aber das Matrixprodukt kriegst Du ja gut hin |
||||||
04.08.2011, 04:20 | stevewilson | Auf diesen Beitrag antworten » | ||||
Okay, habe eine Lösung rausbekommen, die leider falsch ist. Habe ich beim Einsetzen in die DGL gemerkt. Vermutlich irgendwo nen Vorzeichenfehler, der sich fortgepflanzt hat. Habe in die Musterlösung geschaut und deren Lösungsweg ist kürzer, weil die das exp(tA) mithilfe des sogenannten Spektrums von A lösen. ES gibt wohl ein Polynom q, für das gilt: Ich frage mich was dieses D^j ist. Jedenfalls kann man ein Polynom vom Grad = der algebraischen Vielfachheit des Eigenvektors minus 1 bestimmen, mit dem man dann exp(tA) ausrechnen kann.
Wie wäre denn das Integral definiert, wenn ich über einen Vektor integrieren würde? Das hat mich bei Berechnen eines elektrischen Feldes immer interessiert... |
||||||
04.08.2011, 12:18 | Mazze | Auf diesen Beitrag antworten » | ||||
Im Normalfall ist das eine allgemeine Notation für eine partielle Ableitung (j ist dabei ein Multiindex). Allerdings weiß ich nicht ob das hier wirklich gemeint ist.
Das geht dann in die Integrationstheorie und würde an dieser Stelle zu weit führen. |
|
Verwandte Themen
Die Beliebtesten » |
|
Die Größten » |
|
Die Neuesten » |
|