LR Zerlegung ohne Pivotisierung in Matlab |
19.05.2012, 12:50 | steviehawk | Auf diesen Beitrag antworten » | |||||||||||||||
LR Zerlegung ohne Pivotisierung in Matlab Hallo Leute, ich habe ein Programm geschrieben, dass mir zu einer Matrix A die LR Zerlegung vornimmt. Ich möchte es dahingehend verbessern, dass es mir einen Fehler ala: Zerlegung nicht möglich ausgibt wenn es eben nicht möglich ist. Hier mal das Programm:
So ich dachte erst, an den Ansatz über die Determinante der Hauptminore (also der Unterhauptmatrizen) oder wie man die nennt. Die muss ja immer verschiedenen von Null sein, dass es ohne Pivotisierung geht. Nur habe ich das Problem, dass Matlab anscheinend keine Determinante bestimmen kann, zumindest nicht besonders gut, vielleicht schnell, aber nicht genau, statt 0 kommt in Matlab z.B.: 1,5*10-8 oder so raus.. In unserem Skript steht, dass die Pivotelemente nicht null werden dürfen, das habe ich dann so implementiert:
am Ende des Programms. Aber für die Matrix
wird eine LR Zerlegung bestimmt obwohl die Determinante des 3x3 Hauptminors = 0 ist. Also kann es über die Pivotelemente in Programm wohl nicht gesteuert werden oder was mache ich da falsch?? Gibt es einen anderen Weg?? Meine Ideen: Danke für die Hilfe!! Sorry, das ist falsch hier in Analysis, kann das jemand in Software-Abteilung bzw, in Numerik verschieben?? Danke |
|
Verwandte Themen
Die Beliebtesten » |
Die Größten » |
|
Die Neuesten » |
|