Ejemplo de factorización PLU - Egor...

Post on 22-Sep-2020

13 views 0 download

Transcript of Ejemplo de factorización PLU - Egor...

Ejemplo de factorizacion PLU

Egor Maximenkohttp://www.egormaximenko.com

Instituto Politecnico Nacional,Escuela Superior de Fısica y Matematicas

Mexico, D.F.

5 de enero de 2015

PrerrequisitosAlgoritmo de factorizacion LU

(guardando L y U juntas)Matrices de permutacion

Factorizacion PLU

Problema de factorizacion PLU

Hallar una factorizacion PLU de la matriz

A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

.

Se buscan matrices P, L, U tales quePA = LU,P es una matriz de permutacion,L es unitriangular inferior,U es triangular superior.

Una matriz puede tener varias factorizaciones PLU.En este ejemplo una respuesta posible es

P =

1 0 0 00 0 1 00 0 0 10 1 0 0

, L =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

, U =

5 2 4 30 3 −1 30 0 1 00 0 0 −2

.

Problema de factorizacion PLU

Hallar una factorizacion PLU de la matriz

A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

.

Se buscan matrices P, L, U tales quePA = LU,P es una matriz de permutacion,L es unitriangular inferior,U es triangular superior.

Una matriz puede tener varias factorizaciones PLU.En este ejemplo una respuesta posible es

P =

1 0 0 00 0 1 00 0 0 10 1 0 0

, L =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

, U =

5 2 4 30 3 −1 30 0 1 00 0 0 −2

.

Problema de factorizacion PLU

Hallar una factorizacion PLU de la matriz

A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

.

Se buscan matrices P, L, U tales quePA = LU,P es una matriz de permutacion,L es unitriangular inferior,U es triangular superior.

Una matriz puede tener varias factorizaciones PLU.

En este ejemplo una respuesta posible es

P =

1 0 0 00 0 1 00 0 0 10 1 0 0

, L =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

, U =

5 2 4 30 3 −1 30 0 1 00 0 0 −2

.

Problema de factorizacion PLU

Hallar una factorizacion PLU de la matriz

A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

.

Se buscan matrices P, L, U tales quePA = LU,P es una matriz de permutacion,L es unitriangular inferior,U es triangular superior.

Una matriz puede tener varias factorizaciones PLU.En este ejemplo una respuesta posible es

P =

1 0 0 00 0 1 00 0 0 10 1 0 0

, L =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

, U =

5 2 4 30 3 −1 30 0 1 00 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1

=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1

=========

5 2 4 3

3 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1

R3 +=−6R1R4 += 5R1

=========

5 2 4 33

0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1

R4 += 5R1

=========

5 2 4 33

0 0 −2

6

3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33

0 0 −2

6

3 −1 −3

−5

−12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0

0 −2

6

3 −1 −3

−5

−12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0

−2

6

3 −1 −3

−5

−12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26

3 −1 −3

−5

−12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3

−1 −3

−5

−12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1

−3

−5

−12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5

−12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12

5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5

12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2

=========

5 2 4 36 3 −1 −33

0 0 −2

−5

−4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2

R4 += 4R2

=========

5 2 4 36 3 −1 −33 0

0 −2

−5

−4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0

0 −2

−5 −4

1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0 0

−2

−5 −4

1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4

1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1

0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Solucion (en cada paso PA = LU)

P1,2,3,4A =

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

R2 +=−3R1R3 +=−6R1R4 += 5R1=========

5 2 4 33 0 0 −26 3 −1 −3−5 −12 5 12

R2↔R3−−−−→

P1,3,2,4 A =

5 2 4 36 3 −1 −33 0 0 −2−5 −12 5 12

R3 += 0R2R4 += 4R2=========

5 2 4 36 3 −1 −33 0 0 −2−5 −4 1 0

R3↔R4−−−−→

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

R4 += 0R3=========

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

El ultimo estado de la solucion es

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Respuesta:1 0 0 00 0 1 00 0 0 10 1 0 0

︸ ︷︷ ︸

P

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

︸ ︷︷ ︸

A

=

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

︸ ︷︷ ︸

L

5 2 4 30 3 −1 30 0 1 00 0 0 −2

︸ ︷︷ ︸

U

.

El ultimo estado de la solucion es

P1,3,4,2 A =

5 2 4 36 3 −1 −3−5 −4 1 0

3 0 0 −2

.

Respuesta:1 0 0 00 0 1 00 0 0 10 1 0 0

︸ ︷︷ ︸

P

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

︸ ︷︷ ︸

A

=

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

︸ ︷︷ ︸

L

5 2 4 30 3 −1 30 0 1 00 0 0 −2

︸ ︷︷ ︸

U

.

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 330 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 330 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5

2 4 330 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2

4 330 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4

330 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30

12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3

24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1

18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3

−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25

− 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12

− 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1

− 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 0

15 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015

6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0

12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0

9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 330 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15

−25 −22 −15 −315 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15

−25 −22 −15 −315 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

. X

Comprobacion

LU =

1 0 0 06 1 0 0−5 −4 1 0

3 0 0 1

5 2 4 30 3 −1 30 0 1 00 0 0 −2

=

5 2 4 3

30 12 + 3 24− 1 18− 3−25 − 10− 12 − 20 + 4 + 1 − 15 + 20 + 015 6 + 0 12 + 0 + 0 9 + 0 + 0− 2

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

.

P1,3,4,2 A = P1,3,4,2

5 2 4 3

15 6 12 730 15 23 15−25 −22 −15 −3

=

5 2 4 3

30 15 23 15−25 −22 −15 −3

15 6 12 7

. X

EjercicioConstruir una factorizacion PLU de la siguiente matriz y hacer la comprobacion.

A =

0 −2 1 −23 1 5 66 0 11 66 8 1 20

.

¡Ejercer para aprender!

Una respuesta posible:

P =

0 1 0 01 0 0 00 0 0 10 0 1 0

, L =

1 0 0 00 1 0 02 −3 1 02 1 0 1

, U =

3 1 5 60 −2 1 −20 0 −6 20 0 0 −4

.

EjercicioConstruir una factorizacion PLU de la siguiente matriz y hacer la comprobacion.

A =

0 −2 1 −23 1 5 66 0 11 66 8 1 20

.

¡Ejercer para aprender!

Una respuesta posible:

P =

0 1 0 01 0 0 00 0 0 10 0 1 0

, L =

1 0 0 00 1 0 02 −3 1 02 1 0 1

, U =

3 1 5 60 −2 1 −20 0 −6 20 0 0 −4

.

EjercicioConstruir una factorizacion PLU de la siguiente matriz y hacer la comprobacion.

A =

0 −2 1 −23 1 5 66 0 11 66 8 1 20

.

¡Ejercer para aprender!

Una respuesta posible:

P =

0 1 0 01 0 0 00 0 0 10 0 1 0

, L =

1 0 0 00 1 0 02 −3 1 02 1 0 1

, U =

3 1 5 60 −2 1 −20 0 −6 20 0 0 −4

.