Sistema de Ecuaciones Métodos Directos

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

Sistema de Ecuaciones Métodos Directos

4.2 MÉTODOS DIRECTOS

4.2.5 FACTORIZACIÓN DIRECTA DE MATRICES

1. TEOREMA. Las siguientes afirmaciones son equivalentes

a) Ax = 0 tiene una única solución x = 0

b) Ax = b tiene una única solución x=0 , b n-dimensional

c) A-1

d) det (A) 0

e) Eliminación Gaussiana se puede aplicar a Ax = b

2. TEOREMA. Si el procedimiento de eliminación Gaussiana se puede aplicar


a Ax = b sin

intercambio de filas Y la matriz A puede factorizarse como.

A = LU

Donde

U: Es una matriz triangular superior.

L: Es una matriz triangular inferior

3. TEOREMA. Si A es una matriz positiva definida Y existe A-1

4. TEOREMA. Si A es una matriz estrictamente diagonal dominante Y existe A-1

5. ALGORITMO DE FACTORIZACION DIRECTA

Sea Ax = b

Si A = LU  LU x = b

Dra. Pilar Nuñez Blas Página 1


Sistema de Ecuaciones Métodos Directos
donde L = [lij] y U = [uij] son matrices triangulares inferior y superior
respectivamente.

Sea Ux = z  Lz = b

El método de eliminación Gaussiana con sustitución hacia adelante resuelve el


sistema y

con sustitución hacia atrás resuelve el sistema

OBSERVACIÓN: Si lii = 1, el método se llama MÉTODO DE DOOLITLE.

Si uii = 1, el método se llama MÉTODO DE CROUT.

6. ALGORITMO CHOLESKI

Teorema: Si la matriz simétrica A es positiva definida 

A=LL'

donde L es una matriz triangular inferior

7. ALGORITMO DE FACTORIZACION DE CROUT

Sea A una matriz tridiagonal.

A=LU

donde

Dra. Pilar Nuñez Blas Página 2


Sistema de Ecuaciones Métodos Directos

Efectuando A = L U, se obtiene

Ejemplos

1.- Factorización Directa de matrices A=LU

2.- Factorización de Choleski A= L.L'

3.- Facotrización de Métod de Crout

Ejemplos
1) FACTORIZACIÓN DIRECTA DE MATRICES A=LU

Dra. Pilar Nuñez Blas Página 3


Sistema de Ecuaciones Métodos Directos

l11 un = a11 u11 = a11 / l11

l11 u12 = a12 u12 = a12 / l11

l11 u13 = a13 u13 = a13 / l11

l11 u14 = a14 u14 = a14 / l11

l21 u11 = a21 l21 = a21 / u11

l31 u11 = a31 l31 = a31 / u11

l41 u11 = a41 l41 = a41 / u11

l21 u12 + l22 u22 = a22 u22 = (1 / l22) [a22 - l21 u12 ]

l21 u13 + l22 u23 = a23 u23 = (1 / l22) [a23 - l21 u13 ]

l21 u14 + l22 u24 = a24 u24 = (1 / l22) [a24 - l21 u14 ]

l31 u12 + l32 u22 = a32 l22 = (1 / u22) [a32 - l31 u12 ]

l41 u12 + l42 u22 = a42 l42 = (1 / u22) [a42 - l41 u12 ]

l31 u13 + l32 u23 + l33 u33 = a33 u33 = (1 / l33) [a33 - (l31 u13 + l32 u23)]

l31 u14 + l32 u24 + l33 u34 = a34 u34 = (1 / l33) [a43 - (l31 u14 + l32 u24)]

Dra. Pilar Nuñez Blas Página 4


Sistema de Ecuaciones Métodos Directos

Se aplica el método de Doolitle, haciendo lii = 1

Ejemplo

4) Aplicando factorización directa de matrices, resolver:

6.4375 x1 + 2.1849 x2 - 3.7474 x3 + 1.8822 x4 = 4.6351

2.1356 x1 + 5.2101 x2 + 1.5220 x3 - 1.1234x4 = 5.2131

-3.7362 x1 + 1.4998 x2 + 7.6461 x3 + 1.2324 x4 = 5.8665

1.8666 x1 - 1.1104 x2 + 1.2460 x3 + 8.3312 x4 = 4.1322

Solucion:

Lo solucionaremos aplicando el algoritmo:

Paso 01: Tomar l11 = 1

u11 = a11 / l11 = 6.4375

Paso 02: Para j = 2,3,4

u12 = a12 / l11 = 2.1849

u13 = a13 / l11 = -3.7474

u14 = a14 / l11 = 1.8822

l21 = a21 / u11 = 0.33174

l31 = a31 / u11 = -0.5804

l41 = a41 / u11 = 0.28996

Dra. Pilar Nuñez Blas Página 5


Sistema de Ecuaciones Métodos Directos
Paso 03: Para i = 2

Paso 04: Tomar l22 = 1

u22 = a22 - l21 u12 = 4.48528

Paso 05: Para j = 3,4

u23 = (1 / l22) [a23 - l21 u13 ] = 2.7652

u24 = (1 / l22) [a24 - l21 u14 ] = 1.7478

l32 = (1 / u22) [a32 - l31 u12 ] = 0.6171

l42 = (1 / u22) [a42 - l41 u12 ] = -0.388

Paso 03: Para i = 3

Paso 04: Tomar l33 = 1

u33 = (1 / l33) [a34 - (l31 u13 + l32 u23)] = 3.7647

Paso 05: Para j = 4

u34 = (1 / l33) [a43 - (l31 u14 + l32u24)] =30.40340

l34 = (1 / U 33) [a43 - (l41 u13 + l42 u23)] = 0.9052

Paso 06: Para l44 = 3

Paso 07: Tomar z1 = a15 / l11 = 4.6351

Paso 08:

z2 = (1 / l22) [a25 - l21 z1] = 3.6745

z3 = (1 / l 33) [a35 - (l31z1 + l32z2)] = 6.28859

z4 = (1 / l 44) [a45 - (l41z1 + l42z2+ l43z3)] = -1.4752

Paso 09:

x4 = z4 / u44 = -0.3665x3 = (1 / u33) [z3 - u34 x4] = 2.00174

x2 = (1 / u22) [z2 -( u23 x3+u24 x4)] = -0.5575

Dra. Pilar Nuñez Blas Página 6


Sistema de Ecuaciones Métodos Directos
x1 = (1 / u11) [z1 -( u12 x2+u13 x3+u14 x4)] =2.1816

Además L y U son

Dra. Pilar Nuñez Blas Página 7

También podría gustarte