Unidad III: Resolución Numérica de Sistemas de Ecuaciones

lineales

Tema 01 : Métodos Directos:

 Métodos de Gauss-Jordán

 Métodos de Eliminación Gaussiana: Sin pivoteo

 Métodos de Eliminación Gaussiana: Con pivoteo

parcial

 Método de Factorización de Doolitle

 Método de Factorización de Crout

 Factorización con Pivoteo parcial

 Método de factorización QR

I) Dados los siguientes Sistemas de Ecuaciones Lineales: (En forma

manual)
 x  y  2z - w  -13
 x  y  z  -4

Ejercicio N 0 72: 
 x - y  4z  3 w  -13
2x - 2y  3z  3w  -25

Resolver mediante el:

a) Método de Gauss-Jordan;
b) Método de Eliminación Gaussiana; sin pivoteo
c) Método de Eliminación Gaussiana, con pivoteo parcial;
d) Método de Factorización A = LU de Doolitle;
e) Método de Factorización A = LU de Crout;
f) Método de factorización QR

Solución
En forma matricial: AX=B
1  1 2  1  x   - 13 
1 1 1 0   y   - 4 
 
1  1 4 3   z    13 
    
2  2 3  3  w   25
Métodos de Eliminación:
a) Método de Gauss-Jordan:
Matriz ampliada
1  1 2  1  13 
 
1 0 4 
A B  11 11 4 3  13 
 
2  2 3  3  25

Transformación en la matriz ampliada tal que A es matriz diagonal

1  1 2  1  13  1  1 2  1  13
  - F1  F2  
1 0 4  0 2 1 1 9 
A B  11 11 4 3  13 
- F1  F3 
0 0 2 4 0 
  - 2F1  F4  
2  2 3  3  25 0 0  1  1 1 

1  1 2  1  13
 
1/2F3  F 40 2  1 1 9 
0 0 2 4 0 
 
0 0 0 1 1 

1  1 2 0  12 1  1 0 0 8
- 4F4  F3    
0 2 1 0 8  1 / 2 F 3  F 2 0 2 0 0 6 
- F4  F2  1 / 2 F 2  F1
0 0 2 0  4   F 3  F 1 0 0 2 0  4
F4  F1    
0 0 0 1 1  0 0 0 1 1 

Solución de DX=C:

1 0 0 0  x    5 
0 2 0 0  y   6 
 
0 0 2 0  z    4
    
0 0 0 1  w  1 
Retornando al sistema de ecuaciones:
 x  5
 6
 2y

 2z  4
 w 1

Solución del sistema de ecuaciones:

x  -5  y  3  z  -2  w  1

b) Método de Eliminación Gaussiana; sin pivoteo

Matriz ampliada
1  1 2  1  13 
 
1 0 4 
A B  11 11 4 3  13 
 
2  2 3  3  25

Transformación de la matriz ampliada en su triangular superior: U

1  1 2  1  13  1  1 2  1  13
  - F1  F2  
1 0 4  0 2 1 1 9 
A B  11 11 4 3  13 
- F1  F3 
0 0 2 4 0 
  - 2F1  F4  
2  2 3  3  25 0 0  1  1 1 

1  1 2  1  13
 
1/2F3  F 40 2  1 1 9 
0 0 2 4 0 
 
0 0 0 1 1 

Solución de UX=C:

1  1 2  1  x   13
0 2  1 1   y   9 
     
0 0 2 4  z   0 
    
 0 0 0 1   w  1 
Retornando al sistema de ecuaciones:
 x  y  2z - w  -13
 2y - z  w  9


 2z  4w  0
 w 1

Solución del sistema de ecuaciones por el proceso regresivo:

 x  5

y  3

z  -2
w  1
c) Método de Eliminación Gaussiana, con pivoteo parcial:
Matriz ampliada
1  1 2  1  13 
 
1 0 4 
A B  11 11 4 3  13 
 
2  2 3  3  25

Transformación de la matriz A ampliada en su forma triangular superior: U

1  1 2  1  13  2  2 3  3  25
   
1 0 4 1 1 1 0 4 
A B  11 11 4 3  13 
F1 xF4 
1  1 4 3  13 
   
2  2 3  3  25 1  1 2  1  13 

2  2 3  3  25  2  2 3  3  25 
- 1/2F1  F2    
 0 2  1 / 2 3 / 2 17 / 2  1  0 2  1 / 2 3 / 2 17 / 2 
- 1/2F1  F3 - F3  F4
0 0 5 / 2 9 / 2  1 / 2 5 0 0 5/ 2 9 / 2  1/ 2 
- 1/2F1  F4    
0 0 1 / 2 1 / 2  1 / 2 0 0 0  2 / 5  2 / 5

Solución de DX=C:

Retornando al sistema de ecuaciones:

2  2 3  3   x    25 
0 2  1 / 2 3 / 2   y   17 / 2 
     
0 0 5/ 2 9 / 2  z    1/ 2 
    
0 0 0  2 / 5  w  2 / 5
Retornando al sistema de ecuaciones:
 2x  2 y  3z - 3w  -25
 2y - 1/2z  3/2w  17 / 2


 5/2z  9/2w  1 / 2
  2 / 5w  2 / 5

Solución del sistema de ecuaciones por el proceso regresivo:

 x  5

y  3

z  -2
w  1

Métodos de Factorización de

Matrices:
d) Método de Factorización A = LU de Doolitle
Factorización de la matriz A=LU:
Cálculo de U:
Transformación de la matriz A en su forma triangular superior: U
1  1 2 1  1  1 2  1 
  - F1  F2 
0 2  1 1 
A B  11 11 1 0 
4 3 
- F1  F3 
0 0 2 4 
  - 2F1  F4  
2  2 3 3  0 0  1  1 
1  1 2  1 
1/2F3  F 40 2  1 1 
 U
0 0 2 4 
 
0 0 0 1 

Cálculo de L: (Matriz triangular inferior, cuyos elementos de su diagonal

son todos iguales a 1 y los restantes son los multiplicadores con el signo
contrario)
1 0 0 0
1 1 0 0
L
1 0 1 0
 
2 0  1/ 2 1
Comprobación:
1 0 0 0 1  1 2  1 1  1 2  1
1 1 0 0 0 2  1 1  1 1 1 0 
LU     A  A  LU
1 0 1 0  0 0 2 4  1  1 4 3
    
2 0  1/ 2 1  0 0 0 1   2  2 3  3
Solución del sistema de ecuaciones:
L(Y )  B
Si AX  B  A  LU  ( LU ) X  B  L(UX )  B  
UX  Y
Cálculo de L(Y )  B

1 0 00  y1    13 
1
 1 00  y 2    4 

1 0 10  y3    13 
    
2 0  1 / 2 1  y 4   25

Retorno al sistema de ecuaciones:

 y1  13
 y1  y 2  4


 y1  y 3  13
 2y1  1 / 2 y 3  y 4  25

Resolviendo por el proceso de sustitución progresiva:

 y1  13
y  9
 2

y 3  0
 y 4  1

Cálculo de UX  Y
1  1 2  1  x   13
0 2  1 1   y   9 
     
0 0 2 4 z    0 
    
 0 0 0 1   w  1 
Retornando al sistema de ecuaciones lineales:
 x - y  2z - w  4
 2y - z  w  9


 2z  4w  3 / 2
 w 1

Solución del sistema de ecuaciones por el proceso regresivo:

 x  5

y  3

z  -2
w  1

d) Método de Factorización A = LU de Crout

Factorización de la matriz A=LU:
1  1 2  1 l11 0 0 0  1 u12 u13 u14 
1 1
 1 0  l 21 l 22 0 0  0 1 u 23 u 24 

1  1 4 3  l31 l32 l33 0  0 0 1 u 34 
    
2  2 3  3 l 41 l 42 l 43 l 44  0 0 0 1

l11 l11u12 l11u13 l11u14 

l l u l l 21u13  l 22 u 23 l 21u14  l 22 u 24 
  21 21 12 22 
l 31 l 31u12  l 32 l 31u13  l 32 u 23  l 33 l 31u14  l 32 u 24  l 33u 34 
 
l 41 l 41u12  l 42 l 41u13  l 42 u 23  l 43 l 41u14  l 42 u 24  l 43u 34  l 44 

Por igualdad de matrices:

l11  1  l 21  1  l31  1  l 41  2
l u  1  u  1
 11 12 12

l11u13  2  u13  2

l11u14  1  u14  1
l u  l  1  l  2
 21 12 22 22

l31u12  l32  1  l32  0


l 41u12  l 42  2  l 42  0
l u  l u  1  u  1 / 2
 21 13 22 23 23

l31u13  l32u 23  l33  4  l33  2


l 41u13  l 42u 23  l 43  3  l 43  1
l 21u14  l 22u 24  0  u 24  1 / 2

l31u14  l32u 24  l33u 34  3  u 34  2
l u  l u  l u  l  3  l  1
 41 14 42 24 43 34 44 44

Luego:
  l11 0 0 0  1 0 0 0
1 l 0 0  1 2 0 0
L   21 22 
 l31 l32 l33 0  1 0 2 0
   
 l 41 l 42 l 43 l 44  2 0  1 1

1 u12 u13 u14  1  1 2 1

0 1 u 23  
u 24  0 1  1 / 2 1 / 2
U  
0 0 1 u 34  0 0 1 2 
   
0 0 0 1  0 0 0 1 
Comprobación
1 0 0 0 1  1 2  1  1  1 2  1
1 2 0 0 0 1  1 / 2 1 / 2 1 1 1 0 
LU     A  A  LU
1 0 2 0  0 0 1 2  1  1 4 3
    
2 0 1 1  0 0 0 1  2  2 3  3

Solución del sistema de ecuaciones:

L(Y )  B
Si AX  B  A  LU  ( LU ) X  B  L(UX )  B  
UX  Y
Cálculo de L(Y )  B

1 0 00  y1    13 
1
 2 00  y 2    4 

1 0 20  y3    13 
    
2 0  1 1  y 4   25

Retorno al sistema de ecuaciones:

 y1  13
 y1  2 y 2  4


 y1  2 y 3  13
 2y1 - y 3  y 4  25

Resolviendo por el proceso de sustitución progresiva:

 y1  13
y  9 / 2
 2

y 3  0
 y 4  1

Cálculo de UX  Y
1  1 2 1   x   13
0 1  1 / 2 1 / 2   y   9 / 2 
     
0 0 1 2  z   0 
    
0 0 0 1  w   1 
Retornando al sistema de ecuaciones lineales:

 x  1y  2z  w  -13
 y - 1/2z  1/2w  9 / 2


 z  2w  0
 w 1

Resolviendo por el proceso de sustitución regresiva:

 x  5
y  3


 z  2
 w  1

Método de Factorización QR

Identificación de los vectores columna de A:

1  1 2  1
1 1 1 0 
Si A    A  [a 1 a 2 a 3 a 4 ]
1  1 4 3
 
2  2 3  3

Donde:
1 
1 
a1  1 1 1 2   
t

1 
 
 2
  1
1
a 2   1 1  1  2   
t

  1
 
  2
  2
1 
a3  2 1 4 3   
t

 4
 
 3
  1
0
a 4   1 0 3  3   
t

3
 
 3

Cálculo de los vectores unitarios ortogonales: q 1 , q 2 , q 3 y q 4

u1
u 1  a1 ; q 1 
u1
k 1


uj
u k  ak  (a k .q j )q j ; k  2,3,4,..., n Donde: q j  ; j  2,3,4,.., k  1
uj
j 1

Entonces para el ejemplo se tendría:

u1
u 1  a1 ; q 1 
u1
k 1

u k  ak 
 (a .q )q ; k  2,3,4
j 1
k j j

uj
qj  ; j  1,2,3
uj

Procediendo al cálculo:
1  1 
1   
  u1 1 1 
Si u 1  a1   q1  
1  u1 7 1 
   
 2  2

  1 1  1  0.2857
1 1     
1   1 1  1.7143 
Si u 2  a 2  (a 2 q1 )q1     ( 1 1  1  2 
t
)
  1 7 1 7 1  2.2857
       
  2  2 2   0.5714
   0.1543
 0.9258 
u2
Entonces: q 2   
u2   0.1543
 
 0.3086

Si u 3  a3  (a3 q2 )q2  (a3 q1 )q1

t t

 2   0.1543   0.1543  0.378   0.378 

1   0.9258   0.9258   0.378   0.378 
    (2 1 4 3  )    (2 1 4 3 ) 
 4  0.1542  0.1542  0.378   0.378 
         
 3  0.3086  0.3086 0.7559 0.7559
0
0
 
2
 
 1

 0 
 0 
u3  
Entonces: q3  
u 3  0.8944 
 
 0.4472
Si u 4  a4  (a4 q3 )q3  (a4 q2 )q2  (a4 q1 )q1
t t t

  1  0   0    0.1543   0.1543
0  0   0     
    ( 1 0 3  3 )   ( 1 0 3  3 0.9258  )  0.9258 
3  0.8944   0.8944    0.1543   0.1543
         
 3  0.4472  0.4472  0.3086  0.3086
 0.378   0.378   0.3333
 0.378   0.378   0 
 (1 0 3  3 )  
 0.378   0.378   0.0669 
     
0.7559 0.7559  0.1332 
 0.9128
 0 
u4  
Entonces: q 4  
u4  0.1832 
 
 0.3648 
Cálculo de la matriz ortogonal Q:

Q  [q 1 q 2 q 3 q4 ]

 0.378  0.1543 0  0.9128

 0.378 0.9258 0 0 
Q 
 0.378  0.1543 0.8944 0.1832 
 
0.7559  0.3086  0.4472 0.3648 

Cálculo de la matriz ortogonal R:

2.6458  1.8898 4.9137  1.5117

 0 1.8516  0.9258 0.6172 
R  QT * A  
 0 0 2.2360 4.0248 
 
 0 0 0.0016 0.368 
  13 
4
b 
  13 
 
 25
CALCULO DE X
X=R\Q T *b
  5
3
X  
  2
 
1
Por tanto: x  5^ y  3^ z  2^ w  1