Apndice

lgebralinealconwxMaxima

Objetivos
1. Definir matrices con wxMaxima.
2. Aplicar con wxMaxima operaciones con matrices.
3. Aplicar transformaciones elementales de matrices.
4. Calcular el determinante de una matriz cuadrada.
5. Interpretar las propiedades del determinante de una matriz cuadrada.
6. Plantear y resolver sistemas de ecuaciones lineales con wxMaxima.
7. Discutir y resolver sistemas de ecuaciones lineales.
Contenidos
A-1. Matrices. Operaciones con matrices.
A-2 Transformaciones elementales de matrices.
A-3 Rango de una matriz. Inversa de una matriz cuadrada regular. Factorizacin
LU.
A-4. Determinante de una matriz cuadrada. Propiedades. Clculo
A-3. Sistemas de ecuaciones lineales.
Referencias
AEM11

ALANINOS PRATS, J; EXTREMERA LIZANA, J; MUOZ RIVAS, P.

(2011)
Clculo con wxMaxima.

G07

GLASNER, MOSES (2007)

A Maxima Guide for Calculus Students

RR05

REDONDO NEBLE, M. VICTORIA; RODRGUEZ GALVN, J.

RAFAEL (2005)
Introduccin a Maxima

RG07

RODRGUEZ GALVN, J. RAFAEL (2007)

Maxima con wxMaxima: software libre en el aula de matemticas

RR08

RODRGUEZ RIOTORTO, MARIO (2008)

Primeros pasos en Maxima

RV09

RODRGUEZ, MARIO; VILLATE, JAIME (2009)

Manual de Maxima.

VR09

VALLEJO RODRGUEZ, JOS ANTONIO (2009)

Clculo diferencial con Maxima

Apndice:lgebralinealconwxMaxima

A-1.- Matrices. Operaciones con matrices

Los contenidos de este apartado se desarrollan en el archivo Apendice-1.wxm.
En este apartado se ver la definicin de matrices con wxMaxima y las diferentes
operaciones con matrices.

A.1.1. Definicin de matrices submatrices.

Definicin de matrices:

Explicitar una fila o una columna de una matriz:

Obtencin de submatrices por eliminacin de una fila:

Apndice:lgebralinealconwxMaxima

Obtencin de submatrices por eliminacin de una columna:

Obtencin de submatrices por eliminacin de una fila y una columna:

Obtencin de submatrices por eliminacin de ms de una fila y una columna:

Definicin de un vector columna (obsrvese que previamente se ha de cargar el paquete

eigen):

Definicin de una matriz por una funcin de los ndices de fila y de columna:

Apndice:lgebralinealconwxMaxima

Definicin de una matriz diagonal:

Definicin de las matrices unidad (o identidad):

Definicin de una matriz nula de un cierto orden:

Apndice:lgebralinealconwxMaxima

Clculo de la traza de una matriz cuadrada:

A.1.2 Operaciones algebraicas con matrices.

Suma de matrices de la misma estructura:

Producto de una matriz por un escalar:

Combinacin lineal de matrices:

Para el producto de matrices hay que observar que en algunos programas el smbolo es
el asterisco (*). Veamos que pasa con wxMaxima:

Apndice:lgebralinealconwxMaxima

Obviamente esta operacin no es el producto de matrices, para las que hace falta que el
nmero de columnas de la primera sea igual al nmero de filas de la segunda. Veamos
la sintaxis con wxMaxima:

Sirva esto para recordar que el producto de matrices no es conmutativo.

Producto de una matriz por un vector columna:

Transposicin de matrices:

Apndice:lgebralinealconwxMaxima

A.1.3 Concatenacin de matrices.

A partir de unas matrices iniciales se pueden construir nuevas matrices por
concatenacin de filas, es decir, escribir a continuacin de cada fila de una matriz las
filas de la segunda:

Tambin se pueden construir nuevas matrices por concatenacin de columnas, es decir,

escribir a continuacin de cada columna de una matriz las columnas de la segunda:

A.1.4 Matrices por bloques.

A partir de matrices cuadradas se pueden construir nuevas matrices con la estructura
llamada diagonal por bloques o bloque-diagonal:

Apndice:lgebralinealconwxMaxima

Apndice:lgebralinealconwxMaxima

A-2.- Transformaciones elementales de matrices

Los contenidos de este apartado se desarrollan en el archivo Apendice-2.wxm.

A.2.1. Transformaciones elementales de filas y columnas

Se llaman transformaciones elementales de matrices las siguientes:
1. Intercambio de dos filas o dos columnas.
2. Multiplicacin de una fila o columna por un escalar no nulo.
3. Sumar a una fila o columna un mltiple escalar (no nulo) de otra fila o columna.
Veamos la sintaxis con wxMaxima.
1. Intercambio de dos filas o dos columnas:

Cabe decir que esta transformacin se obtiene premultiplicando la matriz por la matriz
unidad del mismo nmero de filas en la que se haya aplicado la transformacin que se
quiere aplicar en la matriz. En efecto:

10

Apndice:lgebralinealconwxMaxima

Intercambio de columnas:

Esta transformacin se obtiene postmultiplicando la matriz por la matriz unidad del

mismo nmero de columnas en la que se haya aplicado la transformacin que se quiere
aplicar en la matriz. En efecto:

2. Multiplicacin de una fila o columna por un escalar no nulo.

No hay una instruccin de wxMaxima que realice esta transformacin. Para esto se ha
de recurrir a la metodologa de premultiplicacin (filas) o postmultiplicacin (columnas)
comentada en el caso anterior:

Apndice:lgebralinealconwxMaxima

11

3. Sumar a una fila o columna un mltiple escalar (no nulo) de otra fila o columna.

Veamos tambin su implementacin con la premultiplicacin por la transformacin

hecha en la matriz unidad:

Y ahora con una columna:

12

Apndice:lgebralinealconwxMaxima

A.2.2. Obtencin de matrices escalonadas y matrices escalonadas reducidas.

Estas matrices se obtienen con las instrucciones triangularize, resultando una matriz
triangular superior, y echelon que da una matriz triangular superior reducida. Veamos
la sintaxis mediante unos ejemplos.

Apndice:lgebralinealconwxMaxima

13

14

Apndice:lgebralinealconwxMaxima

A-3.- Rango de una matriz. Inversa de una matriz. Factorizacin LU.

Los contenidos de este apartado se desarrollan en el archivo Apendice-3.wxm.

A.3.1. Rango de una matriz

Para calcular el rango de una matriz hay que aplicar la instruccin rank:

A.3.2. Inversa de una matriz

Para calcular la inversa de una matriz cuadrada regular hay que verificar primero que su
rango es el adecuado y aplicar despus la instruccin invert:

Como se puede ver a continuacin el clculo mediante los mtodos manuales

tradicionales puede resultar complicado:

Apndice:lgebralinealconwxMaxima

15

A.3.3. Factorizacin LU
Para calcular la descomposicin LU de una matriz cuadrada regular, es decir, dos
matrices la primera de les cuales (L) es triangular inferior y la segunda de las cuales (U)
es triangular superior, que multiplicadas en este orden dan la matriz A, hay
implementada en wxMaxima una instruccin que da el resultado correspondiente.
Veamos un ejemplo.

Veamos que la matriz es regular:

Y ahora apliquemos la instruccin que calcula estos factores:

16

Apndice:lgebralinealconwxMaxima

La segunda matriz es la triangular inferior (L) y la tercera es la triangular superior (U).

En efecto:

Apndice:lgebralinealconwxMaxima

17

A-4.- Determinante de una matriz cuadrada.

Los contenidos de este apartado se desarrollan en el archivo Apendice-4.wxm.
Veamos el desarrollo del clculo del determinante de una matriz cuadrada de orden n,
en los casos n = 2 y n = 3:

Como es sabido, el determinante de una matriz triangular es igual al producto de los

elementos de la diagonal. En efecto:

18

Apndice:lgebralinealconwxMaxima

Veamos ahora un ejemplo numrico:

Como es conocido, la inversa de una matriz cuadrada se puede calcular aplicando el

concepto de determinante, mediante la matriz de adjuntos o cofactores transpuesta
multiplicada por el inverso del determinante de la matriz. As, por ejemplo:

Apndice:lgebralinealconwxMaxima

19

Esta matriz es la inversa de la matriz M1; en efecto:

20

Apndice:lgebralinealconwxMaxima

A-5.- Sistemas de ecuaciones lineales.

Los contenidos de este apartado se desarrollan en el archivo Apendice -5.wxm.

A.5.1. Discusin y resolucin de sistemas lineales.

La mecnica de resolucin de sistemas de ecuaciones lineales sencillas con wxMaxima
ya es conocida. Recordmosla con un ejemplo.
Ejemplo A.5.1. Resolver el sistema de ecuaciones lineales:
2 x - 4 y = -2

-3 x + 5 y = 1
Introducimos las ecuaciones y asignemos a cada ecuacin una referencia:

Para resolver el sistema aplicamos la instruccin linsolve indicando las ecuaciones y

las variables:

Interesa ahora ilustrar el procedimiento de discusin y resolucin de sistemas de

ecuaciones lineales aplicando el teorema de Rouch-Frobenius. Lo hacemos con dos
ejemplos.
Ejemplo A.5.2. Discutir y resolver el sistema de ecuaciones lineales:
x + 2 y + 3z = 1

2 x - y + z = -1

-x + 3 y - 2 z = 2

Como se ha dicho antes, hay que introducir las ecuaciones y asignar a cada ecuacin
una referencia:

Apndice:lgebralinealconwxMaxima

21

Ahora calcularemos con wxMaxima las dos matrices asociadas al sistema lineal: la
matriz del sistema formada por los coeficientes de las incgnitas o variables y a
continuacin la matriz ampliada, resultante de concatenar la matriz del sistema con la
columna de los trminos independientes. Matriz del sistema:

Matriz ampliada: hay una instruccin de wxMaxima que da esta matriz, pero en una
forma ligeramente diferente ya que lo hace pasando al primer miembro los trminos
independientes y, por lo tanto, cambiando su signo. En efecto:

Para calcular la matriz ampliada en la forma estndar, hay que definir el vector de
trminos independientes y hacer la concatenacin con la matriz del sistema. Por lo
tanto:

Ahora se ha de calcular el rango de cada una de estas matrices:

Por lo tanto el sistema es compatible determinado. Para resolverlo por la metodologa

clsica se calcula la matriz escalonada reducida de la matriz ampliada:

22

Apndice:lgebralinealconwxMaxima


De aqu se puede obtener la solucin, as como con la instruccin de wxMaxima a estos
efectos:

Ejemplo A.5.3. Discutir y resolver el sistema de ecuaciones lineales:

x + 2 y + 3 z = 1

2 x - y + z = -1
Introducimos las ecuaciones y asignamos a cada ecuacin una referencia:

Ahora calcularemos con wxMaxima las dos matrices asociadas al sistema lineal. Matriz
del sistema:

Matriz ampliada:

Ahora calculemos el rango de cada una de estas matrices:

Apndice:lgebralinealconwxMaxima

23

Por lo tanto el sistema es compatible y simplemente indeterminado. Para resolverlo

aplicamos la instruccin de wxMaxima:

Es decir la solucin es:

5 z +1
x =
5

y = 3 - 5 z

z = z (parametro)

Ejemplo A.5.4. Discutir y resolver el sistema de ecuaciones lineales:

x + 2 y + 3 z = 1

2 x + 4 y + 6 z = -1
Introducimos las ecuaciones y asignamos a cada ecuacin una referencia:

Ahora calcularemos con wxMaxima las dos matrices asociadas al sistema lineal. Matriz
del sistema:

Matriz ampliada:

24

Apndice:lgebralinealconwxMaxima

Ahora calculemos el rango de cada una de estas matrices:

Por lo tanto el sistema es incompatible. Veamos qu sucede si se quiere resolver con

wxMaxima:

A.5.2. Sistemas lineales con parmetros.

La discusin de un sistema de ecuaciones lineales a menudo involucra la existencia de
uno o ms parmetros en las ecuaciones; la tipologa y la solucin del sistema depende
del valor de estos parmetros y hay que hacer la discusin correspondiente. Ilustraremos
el procedimiento con un par de ejemplos.
Ejemplo A.5.5. Discutir y resolver en funcin del parmetro a el sistema de
ecuaciones lineales:
3 x - 2 y + z = 1

4 x + y - 2 z = 2

2 x - 5 y - az = 3
Introducimos las ecuaciones y asignamos a cada ecuacin una referencia:

Matriz del sistema y matriz ampliada:

Apndice:lgebralinealconwxMaxima

25

Triangularizamos la matriz ampliada:

Ahora tenemos dos situaciones:

1) Si -11a - 44 = 0 , es decir, si a = - 4 , entonces rang ( A) = 2; rang ( A0 ) = 3 y
por lo tanto el sistema es incompatible.
2) Si - 11a - 44 0 , es decir, si a - 4 , entonces rang ( A) = rang ( A0 ) = 3 y
por lo tanto el sistema es compatible determinado.
En efecto:

Solucin del sistema en el segundo caso ( a - 4 ):

26

Apndice:lgebralinealconwxMaxima

Ejemplo A.5.6. Discutir y resolver en funcin de los parmetros a, b el sistema de

ecuaciones lineales:
ax + y + z = b

x + ay + z = b

x + y + az = b
Introducimos las ecuaciones y asignamos a cada ecuacin una referencia:

Matriz del sistema y matriz ampliada:

Discutiremos el rango de la matriz del sistema ya que slo depende de un parmetro y lo

haremos con su determinante, ya que es cuadrada:

Ahora tenemos tres situaciones:

1) a = 1 ;
2) a = -2 ;
3) a 1 y a - 2 .
Primer caso. La matriz del sistema y la matriz ampliada son:

Apndice:lgebralinealconwxMaxima

27

Rango de estas matrices:

Por lo tanto el sistema es compatible doblemente indeterminado. Solucin del sistema:

Segundo caso. La matriz del sistema y su rango son:

Matriz ampliada:

Triangularizamos la matriz:

28

Apndice:lgebralinealconwxMaxima


Por lo tanto si b0, el sistema es incompatible y si b=0, el sistema es compatible
simplemente indeterminado. Solucin del sistema:

Tercer caso. Se sabe que el sistema es compatible determinado y entonces se puede

calcular ya la solucin:

Apndice:lgebralinealconwxMaxima

29