Ejemplo LU Con Pivoteo Parcial 3x3

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 2

Estado inicial:

A es la matriz cuya descomposición LU con pivoteo


parcial deseamos obtener. Inicializamos a L y a P con la matriz identidad.
- 3 10 - 2 1 0 0 1 0 0
A= 6 5 - 10 , L= 0 1 0 , P= 0 1 0
12 - 7 4 0 0 1 0 0 1

Paso k=1:
Busco en A el número de mayor valor absoluto de la columna
k=1 que esté en o debajo de la diagonal principal. Resulta estar en la fila l=3.
- 3 10 - 2
A= 6 5 - 10
12 - 7 4

Intercambio las filas k=1 y l=3 de A y P . Hago lo mismo


con L pero solo para las entradas entre las columnas 1 y k-1=0 (0 columnas).
12 - 7 4 1 0 0 0 0 1
A= 6 5 - 10 , L= 0 1 0 , P= 0 1 0
- 3 10 - 2 0 0 1 1 0 0

En las entradas de L en la columna k=1 y en las filas entre la k+1=2 y la


n=3 ponemos el multiplicador: La razón entre el correspondiente elemento de la
matriz A y el elemento en la posición (k,k) de A (enmarcado en líneas punteadas).
1 0 0
1
L= 2
1 0
- 14 0 1

Efectuamos el paso de eliminación (creación de ceros) en la


matriz A. A las filas de A de índice j con j entre k+1=2 y n=3 se les resta la fila k=1
de la misma matriz A multiplicada por la actual entrada en la posición (j,k) de L.
12 - 7 4
17
A= 0 2
- 12
33
0 4
-1

Paso k=2:
Busco en A el número de mayor valor absoluto de la columna
k=2 que esté en o debajo de la diagonal principal. Resulta estar en la fila l=2.
12 - 7 4
17
A= 0 2
- 12
33
0 4
-1

Intercambio las filas k=2 y l=2 de A y P (en otras palabras, dejo a la fila k donde está). Hago
lo mismo con L pero solo para las entradas entre las columnas 1 y k-1=1 (1 columna).
12 - 7 4 1 0 0 0 0 1
A= 0 172
- 12 , L= 1
2
1 0 , P= 0 1 0
0 33
-1 - 14 0 1 1 0 0
4
2

En las entradas de L en la columna k=2 y en las filas entre la k+1=3 y la


n=3 ponemos el multiplicador: La razón entre el correspondiente elemento de la
matriz A y el elemento en la posición (k,k) de A (enmarcado en líneas punteadas).
1 0 0
1
L= 2
1 0
- 14 33
34
1

Efectuamos el paso de eliminación (creación de ceros) en la


matriz A. A las filas de A de índice j con j entre k+1=3 y n=3 se les resta la fila k=2
de la misma matriz A multiplicada por la actual entrada en la posición (j,k) de L.
12 - 7 4
17
A= 0 2
- 12
0 0 181 17

Estado final:
Las matrices de salida L y P ya están listas. La matriz de salida U es el estado final de la matriz A.
1 0 0 12 - 7 4 0 0 1
1
L= 2
1 0 , U= 0 172
- 12 , P= 0 1 0
- 14 33
34
1 0 0 181 17
1 0 0

Por último, P multiplicada por la matriz de entrada original debe ser igual a L*U.
- 3 10 - 2 12 - 7 4 12 - 7 4
P* 6 5 - 10 = 6 5 - 10 , L*U= 6 5 - 10
12 - 7 4 - 3 10 - 2 - 3 10 - 2

También podría gustarte