Gaussjord Solucion en Maple

3.1.3 Eliminación Gaussiana básica
Ilustraremos el método de Gauss aplicando el procedimiento a un sistema de cuatro

ecuaciones con cuatro incógnitas:
6 2 4 x1 12
][ ][ ]
−2
[ 12 −8
3 −13
−6 4 1 −18
x
6 10 = 2 = 34
9 3 x3 27
x 4 −38
El sistema esta representado en forma matricial dando forma a una función tipo Ax =b
Donde;
A = matriz de coeficientes de las variables

x = vector de variables del sistema
b= vector de términos dependientes
El primer paso es saber si el sistema tiene o no solución, para ello existen varios
procedimientos, en nuestro caso nos enfocaremos en el determinante de la matriz. Si al
calcular el determinante este es distinto de cero entonces el sistema tiene solución. Es decir
| A|≠ 0 , es un sistema compatible deteminado
En el paso 1, es obtener el factor por el cual se va a multiplicar la primera fila para que el
−a21 −12
elemento a 21 sea cero, por lo que vamos a multiplicar la primera fila por = =2 y
a11 2
−3 −1
la sumamos a la segunda, después multiplicamos la primera fila por = y la
6 2
6
sumamos a la tercera y finalmente multiplicamos la primera fila por =1 y la sumamos a
6
la cuarta.
1
Los números −2 ,− y 1 son los factores multiplicadores del primer paso del proceso de
2
eliminación.
El número 6 es el elemento pivote de este primer paso y la primera fila, que no sufre
modificación alguna, se denomina fila pivote. El sistema resultante es equivalente al
primero con lo que el (SEL) tiene el siguiente aspecto:
6 −2 2 4 x1 12
( 0 −4
0 −12
0 2
8 1
3 −14
x3
x4
)()( )
2 2 = x 2 = 10
21
−26
Paso 2. En este paso, la segunda fila se emplea como fila pivote y -4 como elemento
pivote. Aplicamos del nuevo el proceso: multiplicamos la segunda fila por
−12
=−3y la sumamos a la tercera fila y después multiplicamos la segunda fila por
4
12 2 1 1
= = y la sumamos a la cuarta. Los multiplicadores son en esta ocasión 3 y y el
4 4 2 2
sistema de ecuaciones se reduce a:
6 −2 2 4 x1 12
( 0 −4
0 0 2 −5
0 0 4 −13
x3
x4
)( )( )
2 2 = x 2 = 10
−9
−21
Paso 3. Este último paso tomamos las fila 3 como pivote y el elemento 2 será el pivote por
−4
lo vamos ahora a multiplicar la tercera fila por ecuación por =2 y la sumamos a la
2
cuarta fila. El sistema resultante resulta ser:
6 −2 2 4 x1 12
( 0 0 0 −3
x3
x4
)()( )
0 −4 2 2 = x 2 = 10
0 0 2 −5 −9
−3
El sistema resultante es triangular superior y equivalente al sistema original (las soluciones

de ambos sistemas coinciden). Sin embargo, este sistema es fácilmente solucionado
aplicando el algoritmo de sustitución regresiva o hacia atrás explicado en el
apartado 6.1.1. La solución del sistema de ecuaciones resulta ser:
1
x=
()
−3
−2
1
Si colocamos los multiplicadores utilizados al transformar el sistema en una matriz

triangular inferior unitaria (L) ocupando cada uno de ellos la posición del cero que
contribuyó a producir, obtenemos la siguiente matriz:
Por otra parte, la matriz triangular superior (U) formada por los coeficientes resultantes tras
aplicar el algoritmo de Gauss (ecuación 61), es:
Estas dos matrices nos dan la factorización LU de la matriz inicial de coeficientes, A,

expresada por la ecuación (58):
3.1.4 Método de Gauss-Jordan
Como hemos visto, el método de Gauss transforma la matriz de coeficientes en una matriz
triangular superior. El método de Gauss-Jordan continúa el proceso de transformación hasta
obtener una matriz diagonal unitaria (aij=0 para cualquier ).
Veamos el método de Gauss-Jordan siguiendo con el ejemplo empleado en el apartado

anterior. Aplicando el método de Gauss habíamos llegado a la siguiente ecuación:
Ahora seguiremos un procedimiento similar al empleado en el método de Gauss.
Tomaremos como pivote el elemento a44=-3; multiplicamos la cuarta ecuación por y la

restamos a la primera:
Realizamos la misma operación con la segunda y tercera fila, obteniendo:
Ahora tomamos como pivote el elemento a33=2, multiplicamos la tercera ecuación por
y la restamos a la primera:
Repetimos la operación con la segunda fila:
Finalmente, tomamos como pivote a22=-4, multiplicamos la segunda ecuación por y la

sumamos a la primera:
El sistema de ecuaciones anterior es, como hemos visto, fácil de resolver. Empleando la
ecuación (46) obtenemos las soluciones:
Solución del sistema anterior en Maple

> restart:
> with(linalg):
> a:=matrix([[6,-2,2,4],[12,-8,6,10],[3,-13,9,3],[-6,4,1,-
18]]);
 6 -2 2 4
 
12 -8 6 10

a :=  
 3 -13 9 3

 -6 4 1 -18 

> b:=vector([12,34,27,-38]);
b := [ 12, 34, 27, -38 ]
> res:=evalf(linsolve(a,b));
res := [ 1., -3. , -2. , 1. ]
> aum:=augment(a,b);
 6 -2 2 4 12
 
12 -8 6 10 34
aum :=  
 3 -13 9 3 27

 -6 4 1 -18 -38 

> res1:=gaussjord(aum);
1 0 0 0 1
 
0 1 0 0 -3
res1 :=  
0 0 1 0 -2

0 0 0 1 1

> aum:=addrow(aum,1,2,(-aum[2,1]/aum[1,1]));
6 -2 2 4 12
 
0 -4 2 2 10
aum :=  
3 -13 9 3 27

-6 4 1 -18 -38 

6 -2 2 4 12
 
0 -4 2 2 10

aum :=  
 0 -12 8 1 21

-6 4 1 -18 -38 

6 -2 2 4 12
 
0 -4 2 2 10

aum :=  
0 -12 8 1 21

0 2 3 -14 -26 

6 0 1 3 7
 
0 -4 2 2 10

aum :=  
0 -12 8 1 21

0 2 3 -14 -26 

6 0 1 3 7
 
0 -4 2 2 10

aum :=  
0 0 2 -5 -9

0 2 3 -14 -26 
 
6 0 1 3 7
 
0 -4 2 2 10

aum :=  
0 0 2 -5 -9

0 0 4 -13 -21 
 
6 0 1 3 7
 
0 -4 2 2 10

aum :=  

0 0 2 -5 -9
 
0 0 0 -3 -3
 
6 0 1 3 7
 
0 -4 0 7 19
aum :=  
0 0 2 -5 -9

0 0 0 -3 -3

6 0 0 11 23 
 
 2 2 

aum := 0 -4 0 7 19 
 
0 0 2 -5 -9 
 
0 0 0 -3 -3 
6 0 0 0 6
 
0 -4 0 7 19

aum :=  

0 0 2 -5 -9
 
0 0 0 -3 -3
 
6 0 0 0 6
 
0 -4 0 0 12

aum :=  

0 0 2 -5 -9
 
0 0 0 -3 -3
 
6 0 0 0 6
 
0 -4 0 0 12

aum :=  

0 0 2 0 -4
 
0 0 0 -3 -3
 
> x4:=aum[4,5]/aum[4,4];
x4 := 1
> x3:=aum[3,5]/aum[3,3];
x3 := -2
> x2:=aum[2,5]/aum[2,2];
x2 := -3
> x1:=aum[1,5]/aum[1,1];
x1 := 1
>

Gaussjord Solucion en Maple

Cargado por

Copyright:

Formatos disponibles

Gaussjord Solucion en Maple

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Gaussjord Solucion en Maple

Cargado por

Copyright:

Formatos disponibles

3.1.3 Eliminación Gaussiana básica

Ilustraremos el método de Gauss aplicando el procedimiento a un sistema de cuatro

A = matriz de coeficientes de las variables

| A|≠ 0 , es un sistema compatible deteminado

El sistema resultante es triangular superior y equivalente al sistema original (las soluciones

Si colocamos los multiplicadores utilizados al transformar el sistema en una matriz

Estas dos matrices nos dan la factorización LU de la matriz inicial de coeficientes, A,

3.1.4 Método de Gauss-Jordan

obtener una matriz diagonal unitaria (aij=0 para cualquier ).

Veamos el método de Gauss-Jordan siguiendo con el ejemplo empleado en el apartado

Tomaremos como pivote el elemento a44=-3; multiplicamos la cuarta ecuación por y la

Realizamos la misma operación con la segunda y tercera fila, obteniendo:

Repetimos la operación con la segunda fila:

Finalmente, tomamos como pivote a22=-4, multiplicamos la segunda ecuación por y la

Solución del sistema anterior en Maple

También podría gustarte