Microsoft Excel en La SoluciÓn de Problemas de Álgebra Lineal

Revista Pedagogía Universitaria Vol. 8 No. 3 2003

MICROSOFT EXCEL EN LA SOLUCIÓN DE PROBLEMAS DE ÁLGEBRA LINEAL

Lic. Mirelis Rasúa López

Profesor Asistente Dpto Matemáticas y Computación ISMM

[email protected]

RESUMEN

El objetivo del presente trabajo es mostrar las posibilidades que brinda el Microsoft Excel para la solución de problemas típicos del Álgebra Lineal. Parece importante que el estudiante conozca que el universo de problemas que se pueden resolver con este software no se reduce al manejo de tablas y a la Estadística Descriptiva, sino que aún sin usar herramientas avanzadas de Excel (como las macros), es posible resolver Sistemas de Ecuaciones Lineales, calcular determinantes, determinar el rango de una matriz, resolver problemas de Programación Lineal y Cuadrática, e incluso hallar raíces de ecuaciones polinómicas o trascendentes. Todos los ejemplos que se muestran están confeccionados en Microsoft Excel 2000 aunque pueden funcionar perfectamente en el Microsoft Excel 97. El trabajo puede constituir un medio de superación y de consulta para profesores y estudiantes, así como otros interesados. Se ha escrito con fines didácticos, de modo que resulte asequible para aquellos a quienes va destinado.

INTRODUCCIÓN

Muchos profesores de Matemáticas y/o Computación cuando necesitan abordar los procedimientos automatizados del Álgebra Lineal utilizan paquetes especializados de Métodos Numéricos (Ej. MatNum) o de Cálculo Simbólico (Ej: Matemática, Derive, etc.). Es posible que algunos de ellos resuelvan de manera más cómoda o eficiente los problemas que serán objeto de análisis en este trabajo, pero lo cierto es que muchos de esos sistemas no tienen ni remotamente la amplia difusión que tiene el Excel. Es casi posible asegurar que en el 90% de los casos, las computadoras con que un recién graduado se encontrará en su vida laboral, tendrán (mientras Microsoft mantenga su posición dominante en el mercado) instalado el Microsoft Office y por ende el Excel. Por otro lado desde el punto de vista del profesor el uso del Excel como herramienta auxiliar en la solución de los problemas mencionados permitiría disminuir el tiempo de familiarización que necesitan los estudiantes cada vez que se enfrentan con un nuevo paquete especializado (Ej: Derive, MathCad, MatNum, Matemática) y sistematizar los conocimientos que los mismos reciben en el primer año de la carrera sobre el Microsoft Office en la asignatura Introducción a la Computación logrando una mayor integración horizontal y vertical en la carrera.

26

mailto:[email protected]


Debe tenerse en cuenta que al resolver problemas en el ordenador no siempre se obtiene la solución exacta del problema de partida, sino cierta solución aproximada, ¿A qué se debe el error que surge?. Pueden ser indicadas tres razones principales a consecuencia de las cuales surgen errores en la resolución numérica del problema matemático de partida. Ante todo, los datos de entrada del problema de partida (condiciones iniciales y de frontera, coeficientes y segundos miembros de las ecuaciones) se dan siempre con cierta inexactitud. Un error del método numérico condicionado por la prefijación inexacta de los datos de entrada suele denominarse error inevitable. Luego, al sustituir el problema de partida por otro problema discreto aparece un error que se llama error de discretización o, de otra forma, error del método. Finalmente, el orden finito de los números que se suministran al ordenador lleva a errores de redondeo que pueden acumularse en el transcurso de los cálculos. Todo esto debe siempre tenerlo presente el estudiante y el docente, puesto que en ocasiones un problema cuya solución teórica es, por ejemplo, x = 0 puede que a traves del Excel se obtenga un resultado pequeño pero diferente de cero (por ejemplo x = 0,9x10-6 ).

DESARROLLO

Solución de Sistemas de Ecuaciones Lineales.

El problema fundamental del Álgebra Lineal consiste en la resolución del sistema de ecuaciones

(1) fAu =

donde u = (ul, . . ., uN) es el vector buscado, f = (fl, f2,..., fN) es un vector conocido de dimensión N, A = (aij) (i, j = 1, 2, . . ., N) es una matriz cuadrada de dimensión NxN con elementos aij.

Se supondrá inicialmente que la matriz A es regular, es decir que su determinante es distinto de cero, de modo que la ecuación Au = 0 tiene sólo una solución trivial, y el sistema (1) tiene la única solución (2) fAu 1−=

Veamos con un ejemplo concreto cómo resolver un sistema de ecuaciones lineales a través del Excel. En esencia utilizaremos la fórmula (2), es decir, primero calcularemos la matriz inversa y luego multiplicaremos la matriz inversa por el término libre. Este esquema se basa en que en Excel existen dos funciones apropiadas para esto y que comentaremos mas adelante que son MINVERSA, que devuelve la matriz inversa de una matriz y la función MMULT, que devuelve el producto matricial de dos matrices. Estas funciones son ejemplos de las llamadas fórmulas matriciales del ExUna fórmula matricial puede ejecutar varias operaciones y devolver un único resultado o varios resultados. Las fórmulas matriciales actúan en uno, dos o más conjuntos de valores denominados argumentos matriciales. Las fórmulas matriciales se crean del mismo modo que las demás fórmulas, excepto que al concluir su introducción se debe presionar la combinación de teclas CTRL+MAYÚS+ENTRAR.

Cuando se introduce una fórmula matricial, Microsoft Excel inserta de forma automática la fórmula entre llaves ({}).

Supongamos se desea resolver el sistema:

27


=

−−

−−

−

1311975

1111111111111111111111111

5

4

3

2

1

xxxxx

(3)

cuya solución es obviamente:

5,4,3,2,1 54321 ===== xyxxxx (4)

Para ello introduzcamos la matriz del sistema (3) y el vector libre como un bloque rectangular de 5 filas y 6 columnas en una hoja de Excel tal y como se muestra a continuación:

Ahora procederemos a calcular la matriz inversa de la matriz del sistema. Para ello: seleccionaremos un bloque vacío de 5x5 donde se ubicará la matriz inversa y con el rango de 5x5 celdas seleccionado se oprime el signo “=” (significa que se introducirá una fórmula);

28


Luego mediante el asistente de funciones seleccionamos la función MINVERSA;

Esta fórmula tiene un argumento que es la matriz que se desea invertir. Por ello el Excel nos pide que introduzcamos dicho argumento en el diálogo que se muestra a continuación. Observe que aquí podemos sencillamente introducir el rango de celdas correspondientes a la matriz original u oprimir el botón que permite seleccionar este rango de celdas mediante el ratón.

Puede apreciarse que escribimos el rango A1:E5 que corresponde a la matriz original del Sistema (3). Una vez establecido el rango se oprime el botón aceptar y se obtiene el siguiente resultado:

29


Aparentemente, no se ha obtenido la respuesta esperada pues sólo se ha calculado el valor de la primera celda de la matriz inversa. Realmente el problema es que queda un paso importante que es convertir la fórmula a una fórmula matricial. Para ello basta con hacer clic con el ratón en la barra de fórmulas y oprimir CTRL+MAYÚS+ENTRAR. Una vez que se hace esto se obtiene:

Nótese que la fórmula aparece ahora encerrada entre llaves, lo que indica que es una fórmula matricial. En el rango G1:K5 aparece la matriz inversa a la matriz del sistema y basta con multiplicar esta matriz por el vector libre para obtener la solución del sistema (3). Para ello seleccionemos el rango de celdas donde depositaremos el producto (que debe ser un vector columna), por ejemplo podríamos seleccionar el rango L1:L5 y repitiendo los pasos anteriores mediante el asistente de funciones seleccionamos la función MMULT luego de lo cual aparece un diálogo para que determinemos los argumentos de esta fórmula:

30


Aquí evidentemente debemos introducir en Matriz1 el rango G1:K5 que es donde se aloja la matriz inversa ya calculada y en Matriz 2 el rango F1:F5 que es donde está almacenado el vector libre del sistema (3). Una vez aceptados los parámetros y convertida la fórmula a matricial obtendremos:

Observe que ya se obtuvo en el rango L1:L5 la solución correcta (4).

Si no existe seguridad de que la matriz del sistema sea no degenerada es conveniente chequear esta condición primero. Para ello el Excel soporta la función MDETERM que calcula el determinante de una Matriz. Si el determinante es distinto de cero, podemos proseguir con el proceso descrito anteriormente, si es igual cero se necesitarían otros métodos para calcular las soluciones (si existen) del sistema. La función MDETERM no es una función matricial. Ella tiene un solo argumento que es la matriz a la que se le desea calcular el determinante. A modo de ejemplo supongamos deseamos verificar que la matriz:

31


=

1234523456765436543254321

A

es no degenerada. Para ello calculamos el determinante en el Excel. Introduzcamos la matriz A como un bloque rectangular de 5 filas y 5 columnas en una hoja de Excel tal y como se muestra a continuación:

luego seleccionamos una celda en la cual colocaremos el resultado del determinante (por ejemplo F1) y mediante el asistente de funciones introducimos la fórmula MDETERM luego de lo cual aparece un diálogo para que precisemos el argumento de esta fórmula:

En nuestro ejemplo la matriz se encuentra en el rango A1:E5 y ese es el rango que debemos introducir en el argumento. Una vez que se oprime aceptar se obtiene el resultado siguiente:

32


con lo que obtenemos que el determinante de la matriz A es igual a cero y ella es por tanto degenerada.

Usando el mismo esquema para la matriz:

=

1543221543321544321554321

B

se obtendría

de lo cual se desprende que la matriz B es no degenerada ya que su determinante es igual a 1875.

Solución de Sistemas de Ecuaciones Lineales por el método de Gauss.

Hay varias variantes de cálculo del método de Gauss basado en la idea de eliminación sucesiva. El proceso de resolución del sistema de ecuaciones algebraicas lineales (1) o

33


∑=

==N

jijij Nifua

1,,3,2,1, Κ

por el método de Gauss, consta de dos etapas. En la primera etapa (procedimiento directo) el sistema (1) se reduce a un sistema con matriz triangular y en la segunda etapa (procedimiento inverso) se resuelve el sistema triangular.

No es difícil resolver el sistema (1) en el caso particular en que la matriz del sistema A es una matriz triangular superior, es decir si aij = 0 para j<i (i,j = 1, 2, . . ., N),

,

00

0 222

11211

=

NN

N

N

a

aaaaa

A

ΚΚΚΚΚ

ΚΚ

en efecto en este caso el sistema (1) tiene por expresión:

(2)

NNNN

NNNNNNN

NN

NN

fuafuaua

fuauafuauaua

==+

=++=++

−−−−− 1,111,1

22222

11212111

ΚΚΚΚ

Κ

y los componentes del vector u = (u1, u2, ..., uN) se pueden determinar sucesivamente, según las fórmulas:

,),(1; ,111,1

1 ΚNNNNNN

NNN

NN uaf

au

afu −−

−−− −==

.1,2,,2,1),(11

Κ−−=−= ∑+=

NNnuafa

uN

nkknkn

nnn (3)

lo que requiere aproximadamente N2 operaciones. Por ello el paso crucial en el método de Gauss es la transformación del sistema en un sistema con matriz triangular superior. Una vez realizado este paso el sistema resultante puede resolverse por las fórmulas (3) si resulta ser no degenerado.

Pasemos a la exposición detallada del método. El primer paso del método de Gauss consiste en la eliminación de la incógnita u1 de todas las ecuaciones con excepción de la primera. Supongamos que a11 ≠ 0, y escribamos el sistema (1) en la forma

,11212111 fuauaua NN =+++ Κ (4)

Multipliquemos la ecuación (4) por ail / a1l , donde i = 2, 3, . . ., N, y sustrayendo la ecuación obtenida de la i-ésima ecuación (5):

34


,,,3,2,)()( 111

1

11

112

11

1212 Nif

aafu

aaaau

aaaa i

iNNi

iNi

i ΚΚ =−=−++−

Introduciendo las designaciones

,,,3,2,,111

11

11

11)1( Njifaaff

aaa

aa iii

jiijij Κ=−=−= (6)

rescribamos el sistema de ecuaciones obtenido (que es equivalente al sistema (1)) en la forma

,11212111 fuauaua NN =+++ Κ

.,,3,2,)1()1(2

)1(2 Nifuaua iNiNi ΚΚ ==++

La primera columna de la matriz de este sistema se compone de ceros, a excepción del primer elemento para i = 1, j = 1.

El paso segundo consiste en la eliminación de u2 del sistema

1)1(2

)1(2

12

)1(22

)1(22

NNNNN

NN

fuaua

fuaua

=++

=++

ΚΚ

Κ (7)

Para ello, multipliquemos la primera ecuación por y sumemos con la ecuación )/( 22)1(

2 aai−

,,,4,3,)1()1(3

)1(32

)1(2 Nifuauaua iNiNii ΚΚ ==+++

Continuando los razonamientos, obtendremos tras el (N-1)-ésimo paso (es decir, al haber excluido u1, u2, . . . . ., uN-1)

)1()1( −− = N

NNN

NN fua (8)

Llegamos en fin al sistema (2) con la matriz triangular superior

,11313212111 ϕ=++++ NN ubububub Κ

,22323222 ϕ=+++ NN ububub Κ

Κ (9)

,1,111,1 −−−−− =+ NNNNNNN ubub ϕ

.NNNN ub ϕ=

35


que puede ser resuelto como se explicó anteriormente.

Usualmente para el cálculo manual es cómodo realizar las transformaciones elementales por el método de Gauss sobre la matriz ampliada del sistema que es la siguiente matriz:

NNNNN

N

N

faaa

faaafaaa

ΚΚΚΚΚΚ

ΚΚ

11

222221

111211

(10)

El concepto esencial para utilizar el Excel en la resolución de sistemas por el método de Gauss parte de que las transformaciones elementales que se realizan para llevar la matriz a una matriz escalonada equivalen a la multiplicación de ésta por la izquierda por matrices de forma especial, en cierta medida “elementales”. Por ejemplo para hacer cero los elementos de la matriz (10) que están por debajo de la diagonal principal basta con multiplicarla por la izquierda por la matriz “elemental”:

−

−

10

01001

11

1

11

21

ΚΚΚΚΚ

ΚΚ

aa

aa

N

Note que los elementos de la matriz resultante se determinan también por las fórmulas (6).

Veamos un ejemplo concreto. Supongamos se desea escalonar la matriz:

=

523132131312

A

que es la matriz ampliada del sistema:

=++=++=++

523323132

321

321

321

xxxxxxxxx

en el primer paso multiplicaremos la matriz A por la izquierda por la matriz elemental

36


−

−

1021

0123

001

es fácil ver que:

−

−

1021

0123

001

523132131312

=

−−

29

21

250

23

25

210

1312

Para continuar el proceso de escalonamiento basta con multiplicar la matriz resultante con la matriz elemental:

150010001

Esta matriz elemental coincide con la matriz unitaria en todas las celdas a excepción de aquellas en que deseamos hacer cero los elementos de la matriz a escalonar. Como ahora necesitamos hacer cero el elemento a32 de la matriz, colocamos en esa posición de la matriz

elemental el valor 22

32

aa

− que como se puede apreciar es igual a 5.

Efectuando el producto

150010001

−−

29

21

250

23

25

210

1312

=

−

−−

12120023

25

210

1312

de donde fácilmente se obtiene . 1;2;1 123 ==−= xxx

La realización de este esquema en Excel es muy apropiado debido a las facilidades que posee para la multiplicación de matrices.

El primer paso es colocar la matriz ampliada del sistema en un bloque de celdas:

37


en este caso se tomó el rango A2:D4. Luego, a la derecha de la matriz ampliada (dejando una columna de separación) escribimos la matriz elemental necesaria para eliminar los elementos de la primera columna de la matriz ampliada que están por debajo de la diagonal principal:

Debe recordarse que la matriz elemental coincide con la unitaria a excepción de las celdas en las cuales deben anularse los elementos de la Matriz ampliada. Los valores en estas celdas son

correspondientemente 11

12

aa

− y 11

13

aa

− . Ahora efectuamos la multiplicación de la matriz

elemental por la Matriz Ampliada y colocaremos el resultado en el rango A6:D8, es decir, dejando una fila de separación directamente debajo de la matriz ampliada. Para ello debemos utilizar la función matricial MMULT pasando como primer argumento el rango F2:H4 (matriz elemental) y como segundo argumento la matriz a transformar (rango A2:D4). Después de efectuar la operación matricial obtendremos:

38


donde se aprecia claramente el primer paso del escalonamiento. Debemos ahora escribir la matriz elemental del segundo paso y efectuar de nuevo la multiplicación por ella de donde se obtiene:

Al terminar el segundo paso ya se obtiene la matriz escalonada y puede resolverse fácilmente el sistema. Claro está que el método de la matriz inversa es menos trabajoso en Excel que el Método de Gauss, pero esta última variante ofrece ciertas ventajas ya que permite al estudiante ver y comprender los pasos que se siguen, además de que posibilita enfrentar otros problemas como la determinación del rango de una matriz y la solución de sistemas degenerados.

En la siguiente tabla se muestra como determinar el rango de la matriz:

39


=

1232123432345432343212321

A

de donde se aprecia que el número de filas no nulas en la matriz escalón resultante es igual a 2 y por ello rang . 2)( =A

40


Solución de Sistemas de Ecuaciones Lineales por el método de Gauss-Jordán.

El método de Gauss-Jordán es similar al método de Gauss sólo que en vez de aplicar la eliminación sucesiva para reducir el sistema a uno triangular superior se trata de reducir el sistema a otro cuya matriz sea la matriz unitaria.

Pasemos a la exposición detallada del método. El primer paso del método de Gauss-Jordán consiste en la eliminación de la incógnita u1 de todas las ecuaciones con excepción de la primera. Supongamos que a11 ≠0, y escribamos el sistema (1) en la forma

,11212111 fuauaua NN =+++ Κ

,2,2211 Nifuauaua iNiNii ≤≤=+++ Κ

Multipliquemos la primera ecuación por ail / a1l , donde i es cualquiera de los números i = 2, 3, . . ., N, y sustrayendo la ecuación obtenida de la i-ésima ecuación se obtiene :

,,,3,2,)()( 111

1

11

112

11

1212 Nif

aafu

aaaau

aaaa i

iNNi

iNi

i ΚΚ =−=−++−

por lo que introduciendo las designaciones

,,,3,2,,111

11

11

11)1( Njifaaff

aaa

aa iii

jiijij Κ=−=−=

podemos rescribir el sistema de ecuaciones obtenido (que es equivalente al sistema (1)) en la forma

,11212111 fuauaua NN =+++ Κ

.,,3,2,)1()1(2

)1(2 Nifuaua iNiNi ΚΚ ==++

La primera columna de la matriz de este sistema se compone de ceros, a excepción del primer elemento para i = 1, j = 1.

El segundo paso consiste en la eliminación de u2 de todas las ecuaciones del sistema (incluyendo la primera, a diferencia del método de Gauss en que sólo se elimina en las filas j > 2). Si escribimos el sistema después de la primera eliminación en la forma:

iiNN aafuauaua 1)1(

11)1(

12)1(

121)1(

11 , ==+++ Κ

41


1)1(

2)1(2

12

)1(22

)1(22

NNNNN

NN

fuaua

fuaua

=++

=++

ΚΚ

Κ

entonces el próximo paso sería multiplicar la segunda ecuación por y sumarla con cada una de las ecuaciones

)/( )1(22

)1(2 aai−

,,,4,3,1,)1()1(3

)1(32

)1(21

)1(1 Nifuauauaua iNiNiii ΚΚ ==++++

Continuando los razonamientos, obtendremos tras el (N-1)-ésimo paso (es decir, al haber excluido u1, u2, . . . . ., uN-1)

)1()1( −− = NNN

NNN fua

y llegaremos en fin a un sistema con matriz diagonal

,1111 ϕ=ub

,2222 ϕ=ub

Κ

.NNNN ub ϕ=

multiplicando ahora la ecuación i por 1 se reduce el sistema a un sistema con matriz unitaria:

iib/

,/ 1111 bu ϕ=

,/ 2222 bu ϕ=

Κ

NNNN bu /.ϕ=

en el cual la solución es sencillamente el término libre.

Usualmente para el cálculo manual es cómodo realizar las transformaciones elementales por el método de Gauss-Jordán sobre la matriz ampliada del sistema que es la siguiente matriz:

NNNNN

N

N

faaa

faaafaaa

ΚΚΚΚΚΚ

ΚΚ

11

222221

111211

El concepto esencial para utilizar el Excel para resolver sistemas por el método de Gauss-Jordán parte de que las transformaciones elementales que se realizan para llevar la matriz a una matriz diagonal equivalen a la multiplicación de ésta por la izquierda por matrices

42


“elementales”. Por ejemplo, para hacer cero los elementos de la matriz (11) que están por debajo de la diagonal principal basta con multiplicarla por la izquierda por la matriz “elemental”:

−

−

10

01

001

11

1

11

21

Κ

ΚΚΚΚ

Κ

Κ

aa

aa

N

Veamos un ejemplo concreto. Supongamos se desea resolver el sistema, resuelto por el método de Gauss anteriormente, cuya matriz es:

=

523132131312

A

Obviamente esta es la matriz ampliada del sistema:

=++=++=++

523323132

321

321

321

xxxxxxxxx

en el primer paso multiplicaremos la matriz A por la izquierda por la matriz elemental

−

−

1021

0123

001

es fácil ver que:

−

−

1021

0123

001

523132131312

=

−−

29

21

250

23

25

210

1312

Para continuar el proceso de conversión del sistema a un sistema diagonal basta con multiplicar la matriz resultante con la matriz elemental:

43


150010021

Esta matriz elemental coincide con la matriz unitaria en todas las celdas a excepción de aquellas en que deseamos hacer cero los elementos de la matriz a diagonalizar. Como ahora necesitamos hacer cero el elemento a12 y el elemento a32 de la matriz colocamos en esa

posición de la matriz elemental los valores 22

12

aa

− y 22

32

aa

− que como se puede apreciar son

iguales a 2 y a 5.

Efectuando el producto

150010021

−−

29

21

250

23

25

210

1312

=

−

−−

−

12120023

25

210

4202

Ahora debemos hacer cero los elementos que no están en la diagonal principal de la tercera columna. Para ello sencillamente multiplicamos por la matriz elemental:

−

−

100245106101

Efectuemos el producto

−

−

100

24510

6101

−

−−

−

121200

23

25

210

4202

=

−

−−

121200

10210

2002

por último multiplicaremos la matriz ampliada por la homotecia:

−−

12/100020002/1

que corresponde a la división de cada ecuación por 1 iib/

44


−

−

12100

020

0021

−

−−

121200

10210

2002

=

−110020101001

Obsérvese que la matriz del sistema es ahora la matriz unitaria y de ahí fácilmente se obtiene.

La realización de este esquema en Excel es muy apropiado debido a las facilidades que posee para la multiplicación de matrices.

El primer paso es colocar la matriz ampliada del sistema en un bloque de celdas:

En este caso se tomó el rango A2:D4. Luego a la derecha de la matriz ampliada (dejando una columna de separación) escribimos la matriz elemental necesaria para eliminar los elementos de la primera columna de la matriz ampliada que estén fuera de la diagonal principal:

45


debe recordarse que la matriz elemental coincide con la unitaria a excepción de las celdas en las cuales deben anularse los elementos de la Matriz ampliada. Los valores en estas celdas son

correspondientemente 11

12

aa

− y 11

13

aa

− . Ahora efectuamos la multiplicación de la matriz

elemental por la Matriz Ampliada y colocaremos el resultado en el rango A6:D8 es decir dejando una fila de separación directamente debajo de la matriz ampliada. Para ello debemos utilizar la función matricial MMULT pasando como primer argumento el rango F2:H4 (matriz elemental) y como segundo argumento la matriz a transformar (rango A2:D4). Después de efectuar la operación matricial obtendremos:

donde se aprecia claramente el primer paso del proceso de conversión a una matriz diagonal. Debemos ahora escribir la matriz elemental del segundo paso y efectuar de nuevo la multiplicación por ella de donde se obtiene:

46


Al terminar el segundo paso ya se han hecho cero los elementos que están fuera de la diagonal principal en las dos primeras columnas. Escribiendo la matriz elemental para la tercera columna y efectuando el producto se obtiene:

La matriz resultante ya es diagonal y sólo se necesita realizar la homotecia para transformarla en unitaria. La matriz elemental en ese caso es una matriz diagonal cuyos elementos en la diagonal principal son los recíprocos de los correspondientes elementos de la matriz obtenida. Efectuando la operación obtenemos:

47


Por tanto la solución del sistema es la última columna de la matriz obtenida. Aunque el método de Gauss-Jordán necesita de más operaciones para resolver un sistema lineal ofrece la ventaja que a través de él se puede determinar la matriz inversa.

Solución de Sistemas de ecuaciones lineales por métodos iterativos.

Pasemos a la descripción general del método de iteraciones para un sistema de ecuaciones algebraicas lineales

(1) fAu =

Con el fin de resolverlo se elige cierta aproximación inicial y0∈Rn y se hallan sucesivamente soluciones aproximadas (iteraciones) de la ecuación. El valor de una iteración yk+l se expresa en términos de las iteraciones precedentes conocidas yk, yk-1, .... Si, al calcular yk+l, se utiliza sólo una iteración precedente yk, entonces el método iterativo se denomina de un paso (o de dos capas); si, en cambio, yk+1 se expresa en términos de dos iteraciones, yk y yk-l, el método se llama de dos pasos (o de tres capas). Aquí se analizarán solo los métodos de un paso.

En los métodos explícitos la iteración yk+1 se obtiene por la fórmula:

),(11 fAyyy kkkk −−= ++ τ

En los métodos iterativos implícitos para determinar yk+1 hace falta resolver la ecuación:

48


Κ,2,1,0)( ,11 ==−−= ++ kFfAyByBy kkkkk τ

Es natural exigir que el volumen de los cálculos para resolver el sistema Byk+1 = Fk sea sustancialmente inferior al volumen de los cálculos para la resolución directa del sistema

. fAu =

La exactitud del método iterativo se caracteriza por la magnitud del error zk = yk - u, es decir, por la diferencia entre la solución aproximada yk y la solución exacta u del sistema inicial de ecuaciones algebraicas lineales. La sustitución yk = zk + u lleva a una ecuación homogénea para el error.

Suele decirse que un método iterativo converge si

,),(,0 zzzdondezlim kk==

∞→

En el caso general se prefija cierto error (relativo) ε > 0 con el que se debe hallar la solución aproximada yk, los cálculos se dan por terminados cuando queda cumplida la condición:

,0 uyuyn −≤− ε

Si n = n(ε) es el mínimo de los números, para los cuales se verifica la condición anterior, entonces el número total de operaciones aritméticas que han de realizarse para hallar la solución aproximada de la ecuación (1) es igual a Qn(ε) = n(ε)q0, donde q0 es el número de operaciones que se realizan para hallar una iteración. E1 problema consiste en minimizar Qn(ε) eligiendo de modo adecuado B y los parámetros {τk}. Comencemos por analizar los métodos iterativos más simples.

Método de la iteración simple.

Para la resolución del sistema de ecuaciones (1) puede emplearse el método de la iteración simple

)()( )(

1

)(1

iN

j

jkijkkkk fyayfAyyy −−=−−= ∑

=+ ττ

donde τ es el parámetro de la iteración. Al comparar con (2) vemos que el método de la iteración simple se da mediante un esquema explícito de dos capas con el parámetro constante τk ≡ τ.

Método de Seidel Es de amplio uso (en particular cuando es insuficiente la información sobre la matriz A) el método iterativo de Seidel en una de las siguientes formas:

,,,2,1,0,1

)()(

1

)(1 Niafyaya ii

N

ij

ijkij

i

j

jkij Κ=≠=+∑∑

+==+

49


,,,2,1,0,1

)()(1

1

)( Niafyaya ii

N

ij

ijkij

i

j

jkij Κ=≠=+∑∑

+=+

=

Los componentes del vector yk+1 se hallan sucesivamente de ambas fórmulas. Así por ejemplo de ellas determinamos sucesivamente :,,, )(

1)2(1

)1(1

Nkkk yyy +++ Κ

),(12

)()1(

11

)1(1 ∑

=+ −=

N

j

jkijk yaf

ay

),(1 1

1

)(1

1

)()()(1 ∑∑

+

=+

+=+ −−=

i

j

jkij

N

ij

jkij

i

ii

ik yayaf

ay

Con respecto a la convergencia de estos dos métodos tienen lugar los siguientes teoremas:

Teorema: El método de Seidel converge con la velocidad de una progresión geométrica de razón q<1, si A=(aij) es una matriz positiva y

1,,,2,1, <=≤∑≠

qNiaqaij

iiij Κ

Esto significa que para el error zk = yk - u tenemos

01

1 zqzqz kkk

++ ≤≤

Teorema: El método de iteración simple converge con la velocidad de una progresión geométrica de razón q<1, si A=(aij) es una matriz no negativa, con valores propios mínimos y máximos γ1 , γ2 respectivamente, tomando como parámetro τ al valor:

21

2γγ

τ+

=

y en este caso:

2

1,11

γγ

ξξξ

=+−

=q

En caso de que no se conozcan los valores propios de la matriz A (como sucede en muchas ocasiones) se puede tomar τ=1, pero en este caso no se garantiza la velocidad de convergencia que se expresó en el teorema anterior.

Ilustremos el método de iteración simple para resolver el sistema:

50


=+−−=++=−+

74224

13

321

321

321

xxxxxxxxx

Tomemos como aproximación inicial el vector . Dividiendo cada ecuación por el elemento a

)1,1,1(0 =yii obtenemos:

=+−

−=++

=−+

47

21

41

21

41

41

31

31

31

321

321

321

xxx

xxx

xxx

que en notación matricial puede escribirse como:

−=

−

−

472131

121

41

411

41

31

311

3

2

1

x

x

x

Para la primera iteración calculamos primero el vector residuo z0

−

=

−−

−

−

=

−−

−

−

=

1

2

32

472131

1

1

1

121

41

411

41

31

311

472131

121

41

411

41

31

311

00 yz

y luego la aproximación 1y

−=

−

−

=−=

2

1

31

1

2

32

1

1

1

001 zyy

repitiendo el proceso obtenemos:

)

)

0,08333333,0,08333333,1(1 −=z

1,16666667,1,08333333-,1,33333333(2 =y

)0,29166667,0,04166667,0,25(2 =z

51


0,875)1,125,-,1,08333333(3 =y

)0,04166667-,0,13541667-,0,08333333(3 =z

)

)

)

)

0,91666667,0,98958333-1,(4 =y

)0,08854167-,0,01041667-,0,03125(4 =z

1,00520833,0,97916667-0,96875,(5 =y

)0,01302083-,0,01432292,7-0,0260416(5 =z

1,01822917,0,99348958-,0,99479167(6 =y

Como puede apreciarse, el proceso converge a la solución exacta del sistema que es evidentemente . La organización de los cálculos del proceso de iteración simple es en Excel extremadamente sencilla. En primer lugar reservemos un bloque para la matriz del sistema y para el vector libre tal como se aprecia en la siguiente tabla:

( 1,1,1 −=u

Luego, en la columna contigua al término libre colocamos la aproximación inicial y0 y el primer residuo z0 . Para calcular el primer residuo utilizamos una formula matricial tal y como se muestra a continuación:

Observe que la fórmula matricial para el residuo es la diferencia entre el producto de la matriz del sistema por el vector y y el término libre. Esta fórmula es conveniente modificarla para que pueda usarse en otras iteraciones. La modificación consiste en poner referencias absolutas a la matriz del sistema (que siempre estará en el bloque A2:C4) y al término libre D2:D4. Esto significa que debemos sustituir en la fórmula A2:C4 por $A$2:$C$4 y D2:D4 por $D$2:$D$4. Una vez efectuado este cambio obtenemos:

52


Al hacer la modificación no podemos olvidar que debemos concluir la introducción de la fórmula con la combinación de teclas CTRL+MAYÚS+ENTRAR.

Ahora debemos introducir la fórmula para el cálculo de la nueva iteración. Para el cálculo de la nueva aproximación usaremos otra fórmula matricial:

El cálculo del residuo para ésta aproximación y la nueva iteración se obtienen fácilmente si copiamos las 2 últimas columnas en la próxima columna libre (en este caso la columna H). Al hacer esto se obtiene:

Cada nueva iteración se reduce a copiar las últimas dos columnas a la próxima columna libre, por lo que se obtienen extremadamente rápido y de manera sencilla.

53


CONCLUSIONES

El trabajo aborda problemas del Álgebra Lineal muy relacionados entre sí.

Desde el punto de vista metodológico, los métodos analizados permiten apoyar y profundizar en la solución de los problemas.

El uso del Excel permitiría la comprobación rápida de los resultados obtenidos por otros métodos y la sistematización de habilidades y conceptos tales como: la conformidad de matrices para el producto, la obtención del rango de matrices, la relación entre la existencia y unicidad de la solución de un Sistema de Ecuaciones Lineales y el determinante de la misma, etc.

Si logramos que este material contribuya en algo al desarrollo del proceso docente y constituya una motivación que promueva el interés hacia la Computación, por la verdadera importancia que en la actualidad le corresponde a esta disciplina, nos sentiremos satisfechos.

BIBLIOGRAFÍA

1. A. A. Samarski. ¨Introducción a los métodos numéricos¨. Editorial ¨MIR¨. Moscú

2. George Forsythe, Cleve B. Moler. ¨Computer Solution of Linear Algebraic Systems¨.

3. I. S. Berezin and N.P.Zhidkov. ¨Computing Methods¨. Volumen II

4. Manuel Alvarez Blanco y otros. ¨Matemática Numérica¨. Editorial ¨Félix Varela¨. 1998.

5. María Virginia Varela Marcelo y otros.. “Álgebra Lineal”. Editorial Pueblo y educación. 1986.

6. D. Faddieev, I. Sominski. ¨Problemas del Álgebra Superior¨. Editorial ¨MIR¨. Moscú 1971.

54

Microsoft Excel en La SoluciÓn de Problemas de Álgebra Lineal

Documents

Transcript of Microsoft Excel en La SoluciÓn de Problemas de Álgebra Lineal