4 svar
135 visningar
Dani163 1009
Postad: 21 apr 2023 19:38 Redigerad: 21 apr 2023 19:44

LU-faktorisering av rektangulär (m×n) matris

Jag försöker just nu faktorisera en matris AA genom att skriva den som produkten av en nedre triangulär matris $L$ och en övre triangulär matris UU: A=LUA = LU. Jag har dock fastnat på följande problem:

A=123121423122-111A=\begin{bmatrix}1&2&3&1&2\ 1&4&2&3&1\ 2&2&-1&1&1\end{bmatrix}


Jag har försökt att subtrahera första raden från andra och tredje rader, för att erhålla:

1231202-12-10-2-7-1-3\left[\begin{array}{ccccc}1 & 2 & 3 & 1 & 2 \0 & 2 & -1 & 2 & -1 \0 & -2 & -7 & -1 & -3\end{array}\right]


Därefter bytte jag plats på den andra och tredje raden, vilket gav:

123120-2-7-1-302-12-1\left[\begin{array}{ccccc}1 & 2 & 3 & 1 & 2 \0 & -2 & -7 & -1 & -3 \0 & 2 & -1 & 2 & -1\end{array}\right]


Sedan subtraherade jag den andra raden från den tredje raden, vilket gav:

123120-2-7-1-300632\left[\begin{array}{ccccc}1 & 2 & 3 & 1 & 2 \0 & -2 & -7 & -1 & -3 \0 & 0 & 6 & 3 & 2\end{array}\right]


Här fastnade jag och kunde inte fortsätta. Jag skulle verkligen uppskatta all hjälp och vägledning jag kan få för att fortsätta lösa detta problem. Tack på förhand!

SaintVenant 3845
Postad: 22 apr 2023 00:55 Redigerad: 22 apr 2023 00:59

Du är färdig.  Edit: Ser nu att du gör fel i sista steget. Dubbelkolla vad du gjort så ser du det säkert.

Sedan fortsätt med nedan.

Gör en undertriangulär och verifiera att multiplikationen med den övertriangulära blir ursprungsmatrisen.

Dani163 1009
Postad: 22 apr 2023 01:18
SaintVenant skrev:

Du är färdig.  Edit: Ser nu att du gör fel i sista steget. Dubbelkolla vad du gjort så ser du det säkert.

Sedan fortsätt med nedan.

Gör en undertriangulär och verifiera att multiplikationen med den övertriangulära blir ursprungsmatrisen.

Tack för svar. Jag fick veta senare att man inte får byta rader när man ska göra LU-faktorisering, enligt min bok:

Men det jag kom fram till slutligen är dessa:

Och jag vet inte om jag har gjort rätt. här. Vet inte hur man skulle kunna bekräfta med Wolfram Alpha då jag får detta:

Via 

LUDecomposition[{{1,2,3,1,2},{1,4,2,3,1},{2,2,-1,1,1}}]

SaintVenant 3845
Postad: 22 apr 2023 15:36

Du har gjort rätt om multiplikationen blir AA

Det blir den inte i detta fall. Tänk på att LL ska agera på UU och göra dina radoperationer ogjorda för att få tillbaka AA

Du har alltså gjort något fel när du tog fram UU. Du ska inte hålla på och dividera dina rader med någon faktor utan att bokföra det i din LL-matris. Visa dina steg.

Kolla här för lösningsförslag:

https://www.emathhelp.net/calculators/linear-algebra/lu-decomposition-calculator/?i=%5B%5B1%2C2%2C3%2C1%2C2%5D%2C%5B1%2C4%2C2%2C3%2C1%5D%2C%5B2%2C2%2C-1%2C1%2C1%5D%5D

Dani163 1009
Postad: 22 apr 2023 19:35 Redigerad: 22 apr 2023 19:36
SaintVenant skrev:

Du har gjort rätt om multiplikationen blir AA

Det blir den inte i detta fall. Tänk på att LL ska agera på UU och göra dina radoperationer ogjorda för att få tillbaka AA

Du har alltså gjort något fel när du tog fram UU. Du ska inte hålla på och dividera dina rader med någon faktor utan att bokföra det i din LL-matris. Visa dina steg.

Kolla här för lösningsförslag:

https://www.emathhelp.net/calculators/linear-algebra/lu-decomposition-calculator/?i=%5B%5B1%2C2%2C3%2C1%2C2%5D%2C%5B1%2C4%2C2%2C3%2C1%5D%2C%5B2%2C2%2C-1%2C1%2C1%5D%5D

Tack för länken, jag drog inspiration ifrån lösningen men jag vet inte om det går att förbättra. Jag bifogar bild längre ner.

 

Från vad jag har sett ska lösningen skrivas på följande sätt, är min uppfattning korrekt att varje steg i lösningen som bokförs till L-matrisen skrivs som inversen av de elementära matriserna, E1, ..., Ep? Inte klart för mig hur man inverterar dessa eller hur man kommer fram till de. Enligt detta utdrag ur math stackexchange har de skrivit såhär, där diagonalen verkar ha bevarats men att man byter tecken på de siffror som läggs in i matriserna. 



Och detta är ifrån s. 158 i boken Linear Algebra and Its Applications, 6:e utgåvan:



Och lösningen jag skrev. Något du skulle förbättra?


Svara Avbryt
Close