unterbestimmte Lineare Gl < Numerik < Hochschule < Mathe < Vorhilfe
|
Status: |
(Frage) reagiert/warte auf Reaktion | Datum: | 17:09 Mi 24.03.2010 | Autor: | melchior2k |
Aufgabe | Lösung von eines Unter bestimmten Linearen Gleichungsystem Ax=b mittels Lu Zerlegung |
geg: Ax=b, A=(18,250), b=(1,18)
ich möchte das Gleichungssystem in Matlab lösen...ich dachte mir mittels Lu zerlegung A=LU
Lux=b --> Problem dabei dann immer Invertieren muss und das nicht funktioniert
Lc=b c=L(inv)*b
c=ux---> x=u(inv)*c
x=u(inv)*L(inv)*b
Problem ist das die Dimensionen nicht mehr übereinstimmen
wie kann ich das anders lösen?
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
|
|
Status: |
(Antwort) fertig | Datum: | 18:27 Mi 24.03.2010 | Autor: | leduart |
Hallo
ist A eine Zeilenmatrix?
dann hast du links ne Zahl, rechts einen Vektor stehen, also versteh ich deine Aufgabe nicht.
Gruss leduart
|
|
|
|
|
ich meine die Abmessungen also A=8*250 gross und b=1*18 sorry bei rand(,) macht man das so...
|
|
|
|
|
Status: |
(Antwort) fertig | Datum: | 21:16 Mi 24.03.2010 | Autor: | Blech |
Hi,
> Korrektur
> ich meine die Abmessungen also A=8*250 gross und b=1*18
8*250=2000. Wie habe ich mir eine 2000 große Matrix vorzustellen?
1. Gibt's unten einen Formeleditor. Dort ist das Zeichen, das Du suchst, direkt klickibunti zum Anklicken da.
2. Muß bei Ax=b die Zeilenanzahl von A und B übereinstimmen. Das tut sie bei Dir nicht,
3. Ist [mm] $A=18\times [/mm] 250$ immer noch genauso falsch. Die Matrix A ist nicht gleich [mm] "$18\times [/mm] 250$"; das ist eine völlig sinnlose Aussage.
4. Ist es wirklich zuviel verlangt, daß Du Deine Frage nicht in einer Form stellst, bei der sich selbst ein Erstsemester schämen müßte? Immerhin sollen wir Dir hier kostenlos helfen. Ein bißchen Einsatz Deinerseits sollte man da erwarten können.
Setzen, schämen. =)
ciao
Stefan
|
|
|
|
|
Status: |
(Antwort) fertig | Datum: | 21:06 Mi 24.03.2010 | Autor: | Blech |
Hi,
> Lösung von eines Unter bestimmten Linearen Gleichungsystem
> Ax=b mittels Lu Zerlegung
> geg: Ax=b, A=(18,250), b=(1,18)
Wenn das die Dimensionen sein sollen, dann muß es b=(18,1) heißen, sonst stimmen sie hier schon nicht.
>
> ich möchte das Gleichungssystem in Matlab lösen...ich
> dachte mir mittels Lu zerlegung A=LU
> Lux=b --> Problem dabei dann immer Invertieren muss und
> das nicht funktioniert
Wieso in Gottes Namen solltest Du invertieren wollen? Sinn, Zweck und Daseinsberechtigung der LU Zerlegung ist doch, daß man nicht mehr invertieren muß.
Lc=b, hier kannst Du doch einfach von oben nach unten eine Zeile nach der anderen lösen.
Dann das gleiche für c=Ux, nur daß x drastisch unterbestimmt ist, also kannst Du jede Menge frei wählen.
> Lc=b c=L(inv)*b
> c=ux---> x=u(inv)*c
> x=u(inv)*L(inv)*b
>
> Problem ist das die Dimensionen nicht mehr übereinstimmen
Überleg Dir mal die Dimensionen von L und U. Du kannst nur quadratische Matrizen invertieren.
> wie kann ich das anders lösen?
Zeilenweise. Schau Dir irgendeine Einführung zur LU Zerlegung an. Irgendwo mußt Du die Aufgabe ja her haben. =)
ciao
Stefan
|
|
|
|
|
Hab mal gegoogelt...aber hab es nicht so ganz verstanden...kannst du mir es an einem kleinen Beispiel erklären bitte...bzw..ich mach das auch mit iterativen verfahren nur bekomme ich bei jedem verfahren genau das selbe raus...also scheinbar gibts nur unterschiede in iterationschritten und genauigkeit?
Kennst du eine Seite wo zu den einzelnen verfahren etwas Theorie steht also zu lsqr und so gibts bei wiki nicht:-(
gibts eine Funktion die Rundungsfehler herausrechnet oder so?
bzw eine der man sagen kann das ein element kleiner 1 gleich 0 sein soll?
Danke schon mal
basti
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 18:20 Sa 27.03.2010 | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|