LR-Zerlegung; Matlab < Numerik < Hochschule < Mathe < Vorhilfe
|
| Status: |
(Frage) überfällig  | | Datum: | 19:44 Do 18.11.2010 | | Autor: | dennis2 |
| Aufgabe | Schreiben Sie mit MATLAB ein Programm zur Lösung eines linearen Gleichungssystems Ax=b. Zerlegen Sie dazu die Matrix A mittels einer LR-Zerlegung mit Spaltenpivotsuche und lösen Sie das System mittels Vorwärts- und Rückwärtseinsetzen. Achten Sie dabei auf eine effiziente Speicherung der Matrizen L und R.
Testen Sie Ihr Programm für n=8 an der Matrix
[mm] A_{i,j}=1/(i+j-1), [/mm] i,j=1,...,n und dem Vektor [mm] b_i=i+3, [/mm] i=1,...n. |
Zunächst muss ich mitteilen, dass ich ein absoluter Nichtskönner auf dem Gebiet der Programmierung bin. Ich habe nur sehr kleine MATLAB-Kenntnisse.
Ich bin daher für diese Aufgabe zu folgendem uneinheitlichen Ansatz gekommen:
Ich schreibe drei einzelne function-files:
einen zur Herstellung der LR-Zerlegung, einen für das Vorwärtseinsetzen und schließlich einen fürs Rückwärtseinsetzen.
Meine Fragen sind nun:
1.) Wie kann man das zu EINEM Programm machen?
2.) Wie baut man die Spaltenpivotsuche hierin ein?
Ich habe bisher immer nur eine LR-Zerlegung ohne Zeilentausch und/oder Pivoting zustande bekommen.
Bei der Spaltenpivotsuche sucht man sich ja aber den jeweils betragsgrößten Eintrag in einer Spalte. Und dann tauscht man die Zeilen, falls der betragsgrößte Eintrag eben in einer anderen Zeile steht.
Ich habe keine Ahnung, wie ich das implementieren kann.
Ich stelle hier meinen Code lieber nicht ein; denn ich habe ein bisschen Sorge, dass ich dann irgendwelche Urheberrechte verletze, da ich mich an fremden Lösungen orientiert habe.
Hat jemand die Geduld und die Motivation mir bei der Lösung dieser Aufgabe zu helfen?
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 19:55 Do 18.11.2010 | | Autor: | dennis2 |
Ich bin mir bewusst, dass es sehr schwer ist, mir zu helfen, wenn ich keinen Quelltext angebe.
Vielleicht ist eine Hilfe dennoch möglich.
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 20:20 Sa 20.11.2010 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
