82RJGUJH3PTBUYT7SNH2
82RJGUJH3PTBUYT7SNH2
82RJGUJH3PTBUYT7SNH2
y sistemas
de ecuaciones
lineales
Antonio Montes Lozano
1,5 créditos
P00/75004/00191
⎛ a 11 a 21 … a 1n ⎞
⎜ ⎟
A = ⎜ a 21 a 22 … a 2n ⎟
⎜ ⎟
⎝ am 1 am 2 … am n ⎠
© FUOC • P00/75004/00191 Matrices, vectores y sistemas de ecuaciones lineales
Índice
Introducción............................................................................................... 5
Objetivos ...................................................................................................... 6
1. Preliminares .......................................................................................... 7
2. Matrices .................................................................................................. 8
2.1. Producto de matrices ......................................................................... 9
2.2. Suma de matrices............................................................................... 10
2.3. Matrices cuadradas ............................................................................ 12
2.4. La transpuesta de una matriz ............................................................ 14
6. Sistemas homogéneos.......................................................................... 38
Resumen....................................................................................................... 59
Solucionario................................................................................................ 64
Glosario ........................................................................................................ 76
Bibliografía................................................................................................. 77
© FUOC • P00/75004/00191 5 Matrices, vectores y sistemas de ecuaciones lineales
Introducción
Objetivos
Se pretende que, estudiando los conceptos de este módulo y con los ejemplos
y actividades que se incluyen, se alcancen los objetivos siguientes:
1. Preliminares
2x + 3y = 7 ⎫ 2x +3y = 4 ⎫ 2x + 3y = 6 ⎫
⎬ ⎬ ⎬
2x – 3y = 1 ⎭ 2x + 3y = 6 ⎭ 4x + 6y = 12 ⎭
El primer sistema representa dos rectas que se cortan en el punto (2, 1) y tiene una
solución única. El segundo caso representa dos rectas paralelas y no tiene solu-
ción. El tercero corresponde a dos rectas coincidentes y tiene infinitas soluciones.
En estos términos, los tres casos del ejemplo quedan clasificados de la forma
siguiente: el primero es compatible y determinado; el segundo es incompati-
ble y el tercero, compatible e indeterminado.
2. Matrices
Un sistema de m ecuaciones lineales y n incógnitas x1 , x 2 , ..., x n o variables con- … fue introducida en 1857
por los matemáticos británicos
siste en un conjunto de m ecuaciones de la forma siguiente: W. Hamilton y A. Cayley. El pri-
mero que utilizó la palabra
matriz fue J. J Sylvester para
a 11 x 1 + a 12 x 2 + … + a 1 n x n = b 1 ⎫ indicar una ordenación deter-
⎪ minada de números.
a 21 x 1 + a 22 x 2 + … + a 2 n x n = b 2 ⎪
⎬, (1)
⎪
⎪
a m1 x 1 + a m2 x 2 + … + a mn x n = b m ⎭
donde las aij y las b i son escalares conocidos, y el rango de variación de los índices
es 1 ≤ i ≤ m, 1 ≤ j ≤ n. El coeficiente aij es el factor que multiplica x j en la i-ésima
ecuación, y el coeficiente b i es el término independiente de la i-ésima ecuación.
Arthur Cayley
Una matriz m × n es una tabla de números aij donde se tiene 1 ≤ i ≤ m, (1821 - 1895)
De este modo, los coeficientesa ij del sistema de ecuaciones (1) forman una ma-
Advertencia
triz m × n, que denotamos como A y que llamamos matriz del sistema .
Para no confundirnos denota-
mos las matrices con negrita,
Pongamos A = ( a ij)1 ≤ i ≤ m, 1 ≤ j ≤ n para indicar que A es la matriz m × n que tiene con el objetivo de distinguirlas
de los números o elementos.
por elementos los a ij . Cuando sean claros los rangos de variación de i y de j
escribimos simplemente A = (a ij).
⎛ x 1⎞ ⎛ b 1⎞
⎜ ⎟ ⎜ ⎟
x b
x = ⎜ 2⎟ b = ⎜ 2⎟ .
⎜ ⎟ ⎜ ⎟
...
...
⎜ ⎟ ⎜ ⎟
⎝ x n⎠ ⎝ b m⎠
© FUOC • P00/75004/00191 9 Matrices, vectores y sistemas de ecuaciones lineales
⎛ e 11 … e 1n ⎞ ⎛ c 11 … c 1k ⎞ ⎛ d 11 … d 1n ⎞
⎜ ⎟ = ⎜ ⎟⋅ ⎜ ⎟,
...
...
...
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝e m 1 … e mn⎠ ⎝c m 1 … c m k⎠ ⎝ dk 1 … d k n⎠
donde el elemento eij de la i-ésima fila y j-ésima columna del producto E se ob-
tiene multiplicando los elementos de la i-ésima fila de C por los elementos de
la j-ésima columna de D y sumando:
En el producto
de matrices...
k
e ij = c i1 d 1 j + c i2 d 2j + … + c i kd k j = ∑ c il d lj . ... observamos que el número de
l=1 columnas de la matriz de la iz-
quierda y el de las filas de la ma-
triz de la derecha es el mismo ( k).
De forma abreviada escribimos:
E=C ·D
·: L (m ,k ) × L (k ,n ) → L ( m ,n ).
y, por lo tanto, para poder multiplicar dos matrices es necesario que la matriz
de la izquierda tenga el mismo número de columnas que el número de filas de
la matriz de la derecha.
A · x = b. (3)
Ejemplo 1
⎛2 –1⎞
⎛3 –1 0 4⎞ ⎜ ⎟
⎜–1 3 2⎟
⎜ 0 3 1⎟⎟ ⋅ ⎜⎜ .
0⎟⎟
⎜4
–
⎝2 5 –1 0⎠
⎝1 3⎠
Para poder multiplicar A y B es necesario que el número de columnas de A sea igual al de las
filas de B. En el ejemplo dado el número de columnas de A y de filas de B es el mismo, 4, y
el producto es posible. Cada elemento del producto tendrá 4 sumandos. El resultado es una
matriz de 3 filas (número de filas de A) y 2 columnas (número de columnas de B) que pre-
sentamos a continuación:
⎛ 7 7 ⎞
⎜–13 4 ⎟⎟ .
⎜
⎝ 23 8⎠
© FUOC • P00/75004/00191 10 Matrices, vectores y sistemas de ecuaciones lineales
⎛ ⎞
⎜ a 11 a 12 … a 1 n b 1 ⎟
⎜ a ⎟
⎜ 21 a 22 … a 2n b 2 ⎟ . (4)
⎜ ⎟
...
...
...
⎜ ⎟
⎜ a m 1 a m 2 … a mn b m⎟
⎝ ⎠
Ejemplo 2
2x + 3y – 6z = 7 ⎫
⎪
– 3 y + z = 1 ⎬,
⎪
– x + 4y – z = 4 ⎭
escribimos la matriz del sistema y la matriz ampliada. Escribimos el sistema de la forma que
expresa la ecuación (3).
⎛ 2 3 –6 ⎞
⎜ 0 –3 1 ⎟.
⎜ ⎟
⎝ –1 4 –1 ⎠
Matriz ampliada:
⎛2 3 –6 7⎞
⎜ ⎟
⎜0 –3 1 1⎟ .
⎝–1 4 –1 4⎠
Sistema:
⎛ 2 3 –6 ⎞ ⎛ x⎞ ⎛ 7⎞
⎜ 0 –3 1 ⎟⎟ ⋅ ⎜ y⎟ = ⎜ 1⎟ .
⎜ ⎜ ⎟ ⎜ ⎟
⎝–1 4 –1 ⎠ ⎝ z⎠ ⎝ 4⎠
⎛ c 11 c 12 … c 1n ⎞ ⎛ d 11 d 12 … d 1n ⎞
⎜ ⎟ ⎜ ⎟
⎜ c 21 c 22 … c 2n ⎟ ⎜ d 21 d 22 … d 2n ⎟
⎜ ⎟+ ⎜ ⎟ =
...
...
...
...
...
...
⎜ ⎟ ⎜ ⎟
⎝c m 1 c m2 … cmn ⎠ ⎝ d m1 d m2 … d m n⎠
⎛ c 11 + d 11 c12 + d 12 … c 1n + d 1n ⎞
⎜ ⎟
c + d 21 c22 + d 22 … c 2n + d 2n ⎟
= ⎜ 21 = E,
⎜ ⎟
...
...
...
⎜ ⎟
⎝ c m1 + d m1 c m 2 + d m2 … c mn + d mn⎠
© FUOC • P00/75004/00191 11 Matrices, vectores y sistemas de ecuaciones lineales
e ij = cij + d ij .
Ejemplo 3
⎛ ⎞ ⎛ ⎞
⎜ 3 – 1 0 4 ⎟ ⎜ 8 1 0 4 –4 ⎟
⎜ –1 0 3 1 ⎟ + ⎜ 3 2 1 –3 ⎟ .
⎜ ⎟ ⎜ ⎟
⎝ 2 5 – 1 0 ⎠ ⎝ 5 –5 8 2 ⎠
Solución:
⎛ 1194 0⎞
⎜ ⎟
⎜ 2 2 4 –2 ⎟ .
⎜ ⎟
⎝ 7 0 7 2⎠
También podemos asociar a cada matriz A su opuesta, −A, que tiene por ele-
mentos los opuestos de los de A. La suma de A y de −A da la matriz 0. Final-
mente, la suma es conmutativa.
1) (A + B) + C = A + ( B + C) asociativa
2) A + 0 = 0 + A = A elemento neutro (5)
3) A + (−A) = 0 elemento opuesto
4) A + B = B + A conmutativa
Todo conjunto en el que hay definida una operación que tiene las pro-
William Rowman
piedades (1), (2) y (3) recibe el nombre de grupo. Si, además, tiene la Hamilton (1805-1865)
propiedad conmutativa (4), decimos que es un grupo conmutativo.
Matemático nacido en Dublín,
fue un niño precoz que a los
cinco años, además de inglés,
leía latín, griego y hebreo; a los
Así, (L(m,n),+) es un grupo conmutativo. diez, francés, italiano, árabe y
sánscrito, y a los catorce persa.
Abandonó el estudio de las len-
guas para consagrarse a las ma-
Propiedades inmediatas temáticas. Son importantes sus
trabajos sobre números irracio-
nales (Algebra as the Science of
1) Dado un conjunto dotado con una operación interna (que denota- Pure Time , 1933-1835).
remos por “+”), si existe un elemento neutro 0 tal que para todo a se ve-
rifica que a + 0 = 0 + a = a, entonces este elemento es único.
© FUOC • P00/75004/00191 12 Matrices, vectores y sistemas de ecuaciones lineales
4) En un grupo, −(−a) = a.
a+ c=b+c⇒ a = b.
Demostración
1) Supongamos que existe otro neutro 0’ con las mismas propiedades. Por el he-
cho de que el 0 es el elemento neutro, se tendrá 0 + 0’ = 0. Por el hecho de que 0´
también es neutro, se tendrá 0 + 0’ = 0’. Por lo tanto, 0 = 0’, y el neutro es único.
a ′ = 0 + a ′ = ( a″ + a ) + a ′ = a ″ + ( a + a ′ ) = a ″ + 0 = a ″.
5) Para probarlo sólo tenemos que añadir −c por la derecha a cada lado de la
ecuación y utilizar la propiedad asociativa.
⎛0 0 … 0⎞
⎜ ⎟
0 0 … 0⎟
O = ⎜
⎜ ⎟
...
...
...
⎜ ⎟
⎝0 0 … 0⎠
© FUOC • P00/75004/00191 13 Matrices, vectores y sistemas de ecuaciones lineales
⎛1 0 … 0⎞
⎜ ⎟
0 1 … 0⎟
I = ⎜⎜ ⎟.
...
...
...
⎜ ⎟
⎝0 0 … 1⎠
1) (A · B)· C = A · (B · C) asociativa
3) A · (B + C) = A · B + A · C
(A + B) · C = A · C + B · C distributiva
En general: A · B ≠ B · A.
Ejemplo 4
1 3⎞ 1 1⎞
A = ⎛⎝ B = ⎛⎝ .
–2 1⎠ 3 2⎠
10 7⎞ –1 4⎞
AB = ⎛⎝ BA = ⎛⎝ .
1 0⎠ –1 11 ⎠
Un conjunto A dotado con dos operaciones: una suma, con las propie-
dades de grupo conmutativo, y un producto con las propiedades (5) y
(6), es un anillo con unidad. Si, además, el producto es conmutativo,
diremos que A es un anillo conmutativo con unidad .
Propiedades inmediatas
1) a ⋅ 0 = 0 .
2) a ⋅ ( –b ) = –( a ⋅ b ) = ( –a ) ⋅ b .
3) ( –a ) ⋅ ( –b ) = a ⋅ b.
© FUOC • P00/75004/00191 14 Matrices, vectores y sistemas de ecuaciones lineales
Demostración
a ⋅ 0= a ⋅ (0 + 0)= a ⋅ 0 + a ⋅ 0 .
Y, simplificando, resulta a · 0 = 0.
–( a ⋅ b ) = – ( a ⋅ b ) + 0 = – ( a ⋅ b ) + a ⋅ 0 = – ( a ⋅ b ) + a ⋅ ( b + ( –b ) )=
= ( – ( a ⋅ b ) + a ⋅ b ) + a ⋅ ( –b ) = 0 + a ⋅ ( –b )= a ⋅ ( –b ).
Dada una matriz, en ocasiones interesa intercambiar sus filas y sus columnas. La
matriz resultante se llama matriz transpuesta, y la denotaremos como At.
Ejemplo 5
⎛ 3 –1 0 4⎞
A = ⎜⎜ –1 0 3 1⎟⎟
⎝ 2 5 –1 0⎠
⎛ 3 –1 2⎞
⎜ ⎟
–1 0 5⎟
A = ⎜⎜
t
.
0 3 – 1⎟
⎜ ⎟
⎝ 4 1 0⎠
(A + B)t = At + B t
(A ⋅ B)t = Bt ⋅ At
Diremos que una matriz cuadrada que sea igual a su transpuesta es simétrica .
Las matrices que son simétricas lo son respecto a la diagonal principal, y no
cambian al intercambiar filas por columnas.
Ejemplo 6
La matriz siguiente:
⎛ ⎞
⎜ 1 2 3 ⎟
A = ⎜ 2 –5 6 ⎟
⎜ ⎟
⎝ 3 6 7 ⎠
Actividades
⎛ ⎞ ⎛ ⎞
A = ⎜ 2 5⎟ B = ⎜ 3 –1 ⎟
⎝ –1 3 ⎠ ⎝ 4 8 ⎠
Comprobad que A · B ≠ B · A.
2 . El producto de una matriz por sí misma se denota en forma de potencia. De este modo,
pondremos A · A · A = A 3 , etc.
⎛ ⎞
Dada la matriz A = ⎜ 7 –2 ⎟ calculad:
⎝ –1 2 ⎠
a) 2 A3 – 5 A 2 + 4 A + 3 I .
t
b) A ⋅ A .
4. Si A es una matriz 4 × 6:
c ) Ídem con A · A t .
5 . Dada la matriz:
⎛ ⎞
A= ⎜ a b ⎟
⎝ c d⎠
1 ⎛ ⎞
At = --- ⋅ ⎜ d –b ⎟
δ ⎝ –c a ⎠
es inversa de A, es decir A · At = A t · A = I.
6 . Sean:
⎛ x⎞ ⎛ a a a ⎞
⎜ ⎟ ⎜ 11 12 1 ⎟
x = ⎜ y ⎟ A = ⎜ a 12 a 22 a 2 ⎟
⎜ ⎟ ⎜ ⎟
⎝ 1 ⎠ ⎝ a1 a2 a ⎠
Probad que:
x t ⋅ A ⋅ x = a 11 x 2 + 2 a 12 x y + a 22 y2 + 2 a 1 x + 2 a 2 y + a .
3. Espacios vectoriales
Para representar
La multiplicación por escalares también tiene propiedades características
los escalares...
y fáciles de demostrar. Para todo A, B, λ, µ se verifican las propiedades que
... utilizamos las letras griegas λ
presentamos a continuación: (lambda) y µ (mu).
3) (λ + µ) · A = ( λA) + (µA)
Un conjunto dado de una operación interna (que llamaremos suma) con es-
tructura de grupo conmutativo (5) y una multiplicación por escalares (elemen-
tos de un cuerpo) con las propiedades (7) decimos que tiene estructura de
espacio vectorial (sobre el cuerpo de los escalares).
Nos interesa destacar esta estructura de espacio vectorial que tienen las matri-
ces y que también tienen otros muchos conjuntos matemáticos, como por
ejemplo los polinomios.
Propiedades inmediatas
1) 0 ⋅ a = 0 .
2) λ ⋅ 0 = 0
3) ( –1 ) ⋅ a = –a
© FUOC • P00/75004/00191 17 Matrices, vectores y sistemas de ecuaciones lineales
Demostración
Para fijar ideas y por tratarse del espacio vectorial más natural, consideraremos el
conjunto de vectores-columna (o matrices de n filas y una sola columna).
⎛ v1 ⎞
⎜ ⎟
v
v = ⎜⎜ 2 ⎟⎟ .
...
⎜ ⎟
⎝ vn ⎠
⎛ v 1⎞ ⎛ w 1 ⎞ ⎛ v 1 + w 1⎞
⎜ ⎟ ⎜ ⎟
v + w = ⎜ ⎟ + ⎜ ⎟ = ⎜⎜ ⎟
...
...
...
⎟
⎝ v n⎠ ⎝ w n ⎠ ⎝ v n + w n⎠
⎛ v 1⎞ ⎛ λ v1⎞
λv = λ⎜ ⎟ = ⎜ ⎟
...
...
⎜ ⎟ ⎜ ⎟
⎝ v n⎠ ⎝ λ vn⎠
y cumplen las propiedades (5) y (7), que lo configuran como espacio vectorial.
Dado un conjunto de vectores {v1 ,v2 ..., v k} de En, decimos que son li-
nealmente independientes si y sólo si la igualdad:
λ 1 v 1 + λ 2 v 2 + … + λ kv k = 0 (9)
© FUOC • P00/75004/00191 18 Matrices, vectores y sistemas de ecuaciones lineales
λ1 = λ2 = … = λk = 0.
⎛ 1⎞ ⎛ 0⎞ ⎛ 0⎞
⎜ 0⎟ ⎜ 1⎟ ⎜ 0⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
u1 = ⎜ 0⎟ , u 2 = ⎜ 0⎟ , …, u n = ⎜ 0⎟ (10)
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ 0⎠ ⎝ 0⎠ ⎝ 1⎠
⎛ λ 1⎞
⎜ ⎟
λ
λ 1 u 1 + λ 2 u 2 + … + λn u n = ⎜⎜ 2⎟⎟ = 0
...
⎜ ⎟
⎝ λ n⎠
se deriva: λ 1 = λ 2 = ... = λ n = 0 .
⎛v 1⎞
v = ⎜⎜ ⎟⎟ = v 1 u 1 + … + vn un (11)
...
⎝v n⎠
que se expresa diciendo que v es combinación lineal de los vectores u1 , ..., un.
Ejemplo 7
⎛ 1⎞ ⎛ –2⎞ ⎛ –5⎞
v1 = ⎜ –2⎟ v2 = ⎜ 3⎟ v3 = ⎜ 7 ⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ 5⎠ ⎝ –1⎠ ⎝ 2⎠
v 1 + 3 v 2 – v3 = 0 ,
es decir, hay una combinación lineal de ellos igual a 0 donde no todos los escalares son 0.
λ1 λi – 1 λk
v i = – -----v – … – ----
------v – … – -----v .
λi 1 λi i –1 λi k
3.2. Bases
Los vectores (10) son una base de E n, ya que verifican las dos condiciones an-
teriores, como ya hemos visto. Este conjunto de vectores recibe el nombre de
base canónica de En.
Proposición 1
Demostración
En efecto, sea {e1 , ..., ek } una base del espacio. Por la segunda condición de la
definición de base, cualquier vector v de En puede ser expresado en esta base,
es decir, existen v1 , ..., vk escalares tales que:
v = v 1e 1 + … + v k e k.
© FUOC • P00/75004/00191 20 Matrices, vectores y sistemas de ecuaciones lineales
v = v′ 1 e 1 + … + v′ k e k
0 = ( v1 – v′1 )e 1 + … + ( v k – v ′ k) e k
Proposición 2
v = v 1u 1 + … + v nu n.
Demostración
Debemos ver que {v, u2 , ..., un } cumplen las dos condiciones necesarias para
ser una base:
De la ecuación:
λ 1 v + λ 2 u2 + … + λn u n = 0 (12)
se deriva:
λ 1 ( v 1 u 1 + … + v n u n) + λ 2 u 2 + … + λ n un = 0 ,
© FUOC • P00/75004/00191 21 Matrices, vectores y sistemas de ecuaciones lineales
que implica:
λ 1 v 1 u 1 + ( λ 1 v2 + λ 2 )u 2 + … + ( λ 1 v n + λ n ) un = 0 ,
λ 1 v 1 = 0, λ 1 v 2 + λ 2 = 0, …, λ1v n + λn = 0 ,
Por lo tanto, queda probado que (12) implica que λ i = 0 para todo i = 1, ..., n
y los vectores son linealmente independientes.
1 v2 v
u1 = -----v 1 – -----u – … – ----n-u n . (13)
v1 v1 2 v1
w = w 1 u 1 + … + wn u n .
1 v vn ⎞
w = w 1 ⎛ -----v – ----2-u2 – … – -----u + w2 u 2 + … + w n u n =
⎝v1 v1 v 1 n⎠
w w 1 v2 ⎞ w 1 v n⎞
= -----1-v + ⎛ w 2 – ------
------ u + … + ⎛ wn – ------
------ u
v1 ⎝ v1 ⎠ 2 ⎝ v1 ⎠ n
que prueba que cualquier vector de E n se puede escribir como combinación li-
neal de los vectores {v, u 2, ..., un }
Demostración:
e 1 = a 11 u1 + … + a 1n un .
e 2 = a 21 e 1 + a 22 u 2 + … + a 2n un .
Necesariamente, alguno de los a 22 , a 23 , ..., a2n debe ser diferente de cero, ya que
si no, la expresión anterior implicaría que {e 1 , e 2 } son linealmente depen-
dientes, contra la hipótesis. Supongamos también, para fijar ideas, que sea
a 2 2 ≠ 0. Aplicando nuevamente la proposición 2 pasaremos a la nueva base
{e 1,e 2,u 3 ..., u n }.
Teorema 2
Teorema 3
Así, pues, el teorema 2 nos dice que no podemos tener una cantidad de vec-
tores mayor que la dimensión del espacio y que sean linealmente indepen-
dientes.
Ejemplo 8
⎛ 0⎞ ⎛ 0⎞
e 1 = ⎜ 3⎟ e 2 = ⎜ –1⎟ .
⎜ ⎟ ⎜ ⎟
⎝ –2⎠ ⎝ 1⎠
© FUOC • P00/75004/00191 23 Matrices, vectores y sistemas de ecuaciones lineales
b) Formemos una base de E 3 completando {e1 , e2} con vectores de la base canónica.
Solución
Tenemos:
⎛ 0⎞
e 1 = ⎜ 3⎟ = 3 u 2 – 2 u 3 ,
⎜ ⎟
⎝ –2⎠
de donde resulta:
1 2
u 2 = ---e 1 + --u 3 .
3 3
⎛ 0⎞
e 2 = ⎜ –1⎟ = – u 2 + u 3 = – ⎛ 1 --u ⎞ + u 3 = – 1
---e + 2 ---e 1 + 1
---u 3 .
⎜ ⎟ ⎝ 3 1 3 3⎠ 3 3
⎝ 1⎠
u 3 = e1 + 3 e2
1 2
u 2 = ---e 1 + -- ( e 1 + 3 e 2 ) = e 1 + 2 e 2 .
3 3
Ejemplo 9
Escribamos la base canónica para el espacio vectorial de las matrices n × m. ¿Qué dimensión
tiene?
La base canónica estará formada por cada una de las matrices que tienen un 1 en una posi-
ción diferente y todo el resto 0.
⎛1 0 … 0⎞ ⎛0 1 … 0⎞ ⎛0 0 … 1⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜0 0 … 0⎟ ⎜0 0 … 0⎟ ⎜0 0 … 0⎟
⎟, …,
⎜ ⎜ ⎟, ⎜ ⎟
...
...
...
...
...
...
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝0 0 … 0⎠ ⎝0 0 … 0⎠ ⎝0 0 … 0⎠
⎛0 0 … 0⎞ ⎛0 0 … 0⎞ ⎛0 0 … 0⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜1 0 … 0⎟ ⎜0 1 … 0⎟ ⎜0 0 … 1⎟
…,
⎜ ⎟, ⎜ ⎟, ⎜ ⎟
...
...
...
...
...
...
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝0 0 … 0⎠ ⎝0 0 … 0⎠ ⎝0 0 … 0⎠
...
...
...
⎛0 0 … 0⎞ ⎛0 0 … 0⎞ ⎛0 0 … 0⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜0 0 … 0⎟ ⎜0 0 … 0⎟ ⎜0 0 … 0⎟
, , …,
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
...
...
...
...
...
...
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝1 0 … 0⎠ ⎝0 1 … 0⎠ ⎝0 0 … 1⎠
© FUOC • P00/75004/00191 24 Matrices, vectores y sistemas de ecuaciones lineales
Es inmediato comprobar que las matrices anteriores son linealmente independientes y que
toda matriz puede expresarse como combinación lineal de éstas. En consecuencia, forman
una base y la dimensión del espacio vectorial de las matrices n × m es n · m.
Ejemplo 10
F 1 = { λv λ ∈ R } .
En particular, siv es el vector 0 , el subespacio se reduce al vector 0, que por sí solo constituye
un subespacio vectorial trivial, que decimos que es de dimensión 0.
F = { λ 1 v 1 + … + λk v k λ 1 ∈ R , …, λ k ∈ R }
Actividades
b) { ( x, y, z) ∈ R 3 : 5 x + 2 y – 7 z = 3 } .
c) { ( x , y, z) ∈ R 3 : x ≥ 0 } .
11. Sea F(R,R) el conjunto de funciones reales de variable real. Decid si son subespacios
vectoriales los subconjuntos de F(R,R) siguientes:
a) { f ∈ F ( R , R ) : f ( 1 ) = k } .
b) { f ∈ F ( R, R ) : f ( 0 ) = 2 f ( 1) } .
c) { f ∈ F ( R , R ) : f ( –x ) = f 2 ( x ) }.
© FUOC • P00/75004/00191 26 Matrices, vectores y sistemas de ecuaciones lineales
4. El método de Gauss
Karl F. Gauss
3x + 2y + z = 1 ⎫ (1777-1855)...
⎪
y – 2z = 2 ⎬ ... llamado Princeps Mathemati-
⎪
2z = –1 ⎭ corum (príncipe de los matemá-
ticos), ha sido, tal vez, el mayor
matemático de la historia. Fue
también, como E. Galois, un
niño precoz: a los 3 años descu-
es inmediato de resolver por sustitución hacia atrás: aislamos la z en la tercera brió un error de cálculo en la
paga de los obreros que tenía su
ecuación, después la sustituimos en la segunda para aislar la y, y, finalmente, padre (no sabemos si a favor de
sustituimos la y y la z en la primera y aislamos la x. De este modo, tenemos: los obreros o de su padre).
z = –1
--- ⎫
2 ⎪
⎪
1 ⎪
y = 2 + 2z = 2 + 2 ⎛⎝ –---⎞⎠ = 1 ⎬.
2 ⎪
⎪
1 1⎛ 1⎞ 1 ⎪
x = ---( 1 – 2y – z ) = ---⎝ 1 – 2 + ---⎠ = –---
3 3 2 6 ⎭
Decimos que un sistema como el del ejemplo está en forma triangular, ya que
la matriz del sistema tiene ceros en la parte inferior izquierda de la diagonal
principal.
En el objetivo de reducir el sistema a otro con las mismas soluciones pero más
sencillo, el método de Gauss utiliza tres tipos de transformaciones que presen-
tamos a continuación.
© FUOC • P00/75004/00191 27 Matrices, vectores y sistemas de ecuaciones lineales
a 11 x 1 + … + a 1n x n = b 1 ⎫
⎪
...
...
...
⎪
a i1 x 1 + … + a in x n = b i ⎪
⎪
...
...
...
⎪
⎪
a j – 1, 1 x 1 + … + a j – 1 , n x n = b j – 1 ⎬ (14)
⎪
( a j1 + λ a i1 )x 1 + … + ( a jn + λ a in )x n = b j + λ b i ⎪
⎪
a j + 1, 1 x 1 + … + a j + 1, n x n = b j + 1 ⎪
⎪
...
...
...
⎪
a m 1 x 1 + … + a mnx n = b m ⎭
y recíprocamente.
Pero esto es obvio: toda solución de (1) verifica (14), ya que la única ecuación
diferente, la j-ésima, es una combinación lineal de la i-ésima y j-ésima ecua-
ciones de (1), y recíprocamente, toda solución de (14) lo es de (1), ya que la
ecuación diferente, la j-ésima ecuación de (1), es igual a la j-ésima ecuación de
(14) menos λ por la ecuación i-ésima.
a ′11 x 1 + a′ 12 x 2 + … + a ′1n x n = b ′1 . ⎫
⎬
...
...
...
...
pivote. De este modo sustituimos la fila i por la fila i menos la fila 1 multipli-
a
cada por ----i1
---. La nueva fila i será:
a 11
a i1 ⎞ a a
0 + ⎛⎝ a i2 – -----
--a 12 x 2 + ⎛⎝ a i3 – ----i1
---a 13⎞ x 3 + … + ⎛⎝ a in – ----i1
---a 1 n⎞⎠ x n =
a 11 ⎠ a 11 ⎠ a 11
a
= b i – ----i1
---b .
a 11 1
Ejemplo 11
2x + 3 y + z = 6 ⎫
⎪
x – 2y – 5 z = 5 ⎬ .
⎪
3x + 5y + 3z = 8 ⎭
2x + 3y + z = 6 ⎫
⎪
–7 ⁄ 2 y − 1 1 ⁄ 2 z = 2 ⎬.
⎪
1 ⁄ 2 y + 3 ⁄ 2 z = –1 ⎭
2x + 3y + z = 6 ⎫⎪
–7 ⁄ 2 y − 1 1 ⁄ 2z = 2 ⎬.
⎪
5 ⁄ 7 z = –5 ⁄ 7 ⎭
Cuando hayamos llegado a este punto, el sistema habrá quedado triangular. Podemos aislar
la z de la última ecuación y tendremos:
z = –1.
© FUOC • P00/75004/00191 29 Matrices, vectores y sistemas de ecuaciones lineales
--⎛ 2 + 11
y = –2 ------z⎞ = –2
---⎛2 – 11
------⎞ = 1 .
7⎝ 2 ⎠ 7⎝ 2 ⎠
1 1
x = ---( 6 – 3 y – z ) = ---( 6 – 3 + 1 ) = 2.
2 2
Una forma más clara de seguir el procedimiento es utilizar la notación matricial. La matriz
ampliada del sistema dado es:
⎛2 3 1 6⎞
⎜ ⎟
⎜1 –2 –5 5⎟ .
⎜ ⎟
⎝3 5 3 8⎠
En el primer paso hemos elegido como pivote el coeficiente a 11 = 2, y lo que hemos hecho
ha sido multiplicar la primera fila por –1/2 y sumar el resultado a la segunda, y multiplicar la
primera fila por –3/2 y sumar el resultado a la tercera. El resultado es ahora:
⎛ 2 3 1 6⎞
⎜ ⎟
⎜0 –7 ⁄ 2 –11 ⁄ 2 2 ⎟ .
⎝0 1 ⁄2 3 ⁄ 2 –1⎠
⎛2 3 1 6 ⎞
⎜ ⎟
⎜0 –7 ⁄ 2 –11 ⁄ 2 2 ⎟. (15)
⎝0 0 5 ⁄ 7 –5 ⁄ 7⎠
No todos los sistemas tendrán siempre una única solución, como sucedía en
el ejemplo anterior. Veamos un caso de ello en el ejemplo siguiente:
Ejemplo 12
x + y –z = 0 ⎫
⎪
– 2x + 3z = 5 ⎬.
⎪
x + 5 y + z = 1 0⎭
⎛ 1 1 –1 0 ⎞
⎜ ⎟.
⎜ –2 0 3 5 ⎟
⎝ 1 5 1 10⎠
Aplicándole el método de Gauss (segunda transformación tres veces) obtenemos la matriz si-
guiente:
⎛1 1 –1 0 ⎞ ⎛1 1 – 1 0⎞
⎜ ⎟ ⎜ ⎟
⎜0 2 1 5 ⎟→⎜ 0 2 1 5⎟ .
⎝0 4 2 10 ⎠ ⎝ 0 0 0 0⎠
La matriz final tiene una fila menos que incógnitas, y ahora resulta que tenemos un grado
de libertad para elegir la z, que puede tomar cualquier valor; esto se debe a que la tercera ecua-
ción, que es la que queda con esta variable, es 0 · z = 0, que se cumple para cualquier valor
real de z. Ponemos, entonces, z = t .
© FUOC • P00/75004/00191 30 Matrices, vectores y sistemas de ecuaciones lineales
Sustituyendo esta expresión en y, y las dos en x, obtenemos que la solución general del siste-
ma depende de un parámetro arbitrario t, y es la siguiente:
z = t ⎫
⎪
y = 1 ⁄ 2 ⋅ (5 – z ) = 1 ⁄ 2 ⋅ ( 5 – t ) ⎬.
⎪
x = z – y = t – 1 ⁄ 2 ⋅ ( 5 – t ) = 1 ⁄ 2 ⋅ ( 3t – 5 ) ⎭
Ejemplo 13
⎛ 1 –1 2 0 2 3⎞
⎜ ⎟
⎜ 1 –1 2 0 3 5⎟
⎜ 1 0 2 3 6 1 ⎟.
⎜ ⎟
⎜ –1 2 –2 4 4 –6⎟
⎜ ⎟
⎜ 2 –2 4 0 5 8⎟
⎝ 1 –1 2 0 1 1⎠
Tomamos como pivote el elemento de la primera fila y primera columna a11 , y reducimos a
cero los restantes de la primera columna. El resultado es:
⎛ 1 –1 2 0 2 3⎞
⎜ ⎟
⎜ 0 0 0 0 1 2⎟
⎜ 0 1 0 3 4 –2⎟ .
⎜ ⎟
⎜ 0 1 0 4 6 –3⎟
⎜ ⎟
⎜ 0 0 0 0 1 2⎟
⎝ ⎠
0 0 0 0 –1 –2
Ahora tomamos como pivote el elemento de la segunda columna y tercera fila, y permuta-
mos así la segunda fila y la tercera:
⎛ 1 –1 2 0 2 3⎞
⎜
⎜ 0 1 0 3 4 –2⎟⎟
⎜ 0 0 0 0 1 2⎟ .
⎜ ⎟
⎜ 0 0 0 1 2 –3⎟
⎜ ⎟
⎜ 0 0 0 0 1 2⎟
⎝ 0 0 0 0 1 –2⎠
⎛ 1 –1 2 0 2 3⎞
⎜
⎜ 0 1 0 3 4 –2⎟⎟
⎜ 0 0 0 0 1 2⎟
⎜ ⎟.
⎜ 0 0 0 1 2 –1⎟
⎜ ⎟
⎜ 0 0 0 0 1 2⎟
⎝ 0 0 0 0 –1 –2⎠
En este ejemplo ha resultado, además, que se han anulado también los restantes elementos
de la tercera columna.
© FUOC • P00/75004/00191 31 Matrices, vectores y sistemas de ecuaciones lineales
Tomamos como pivote el elemento de la cuarta fila y cuarta columna. Así, permutamos las
filas tercera y cuarta:
⎛ 1 –1 2 0 2 3⎞
⎜ ⎟
⎜ 0 1 0 3 4 –2⎟
⎜ 0 0 0 1 2 –1⎟ .
⎜ ⎟
⎜ 0 0 0 0 1 2⎟
⎜ ⎟
⎜ 0 0 0 0 1 2⎟
⎝ ⎠
0 0 0 0 – 1 –2
Finalmente, tomamos como pivote el elemento de la cuarta fila y quinta columna y el siste-
ma queda reducido completamente:
⎛ 1 –1 2 0 2 3⎞
⎜ ⎟
⎜ 0 1 0 3 4 –2⎟
⎜ 0 0 0 1 2 –1⎟ .
⎜ ⎟ (16)
⎜ 0 0 0 0 1 2⎟
⎜ ⎟
⎜ 0 0 0 0 0 0⎟
⎝ ⎠
0 0 0 0 0 0
x5 = 2
x 4 = – 1 – 2 x 5 = –5
x3 = t
x2 = – 2 – 3 x4 – 4 x5 = 5
x 1 = 3 + x 2 – 2 x3 – 2 x 5 = 4 – 2 t
que, como vemos, depende de un parámetro arbitrario t (si bien en este ejemplo algunas in-
cógnitas tienen valor fijo, independientemente de t).
En este ejemplo hemos visto que, al reducir a cero los elementos de una co-
lumna (la segunda), se nos han reducido “por azar” los elementos de la si-
guiente. Cuando esto sucede podemos hacer también una permutación de
columnas, con el objetivo de continuar colocando el nuevo pivote en la dia-
gonal principal de la matriz, siempre que permutemos también el orden de las
variables correspondientes. Cuando reducimos con este criterio hablamos de
pivotamiento total. Esto no aporta nada esencial a la solución del sistema,
aparte del motivo estético que supone que los pivotes estén en la diagonal
principal. Por esta razón generalmente no lo haremos.
⎛ x1 x2 x4 x3 x5 ⎞
⎜ ⎟
⎜ 1 –1 0 2 2 3⎟
⎜ –2⎟⎟
⎜ 0 1 3 0 4
⎜ 0 0 1 0 2 –1⎟ .
⎜ ⎟
⎜ 0 0 0 0 1 2⎟
⎜ ⎟
⎜ 0 0 0 0 0 0⎟
⎝ 0 0 0 0 0 0⎠
⎛ x1 x2 x4 x5 x3 ⎞
⎜ ⎟
⎜ 1 –1 0 2 2 3⎟
⎜ –2⎟⎟
⎜ 0 1 3 4 0
⎜ 0 0 1 2 0 –1⎟ .
⎜ ⎟
⎜ 0 0 0 1 0 2⎟
⎜ ⎟
⎜ 0 0 0 0 0 0⎟
⎝ 0 0 0 0 0 0⎠
x3 = t
x5 = 2
x 4 = – 1 – 2x 5 = –5
x 2 = – 2 – 3x 4 – 4x 5 = 5
x 1 = 3 + x 2 – 2x 3 = 4 – 2t
que, como vemos, coincide exactamente con la solución dada con pivota-
miento parcial.
Las soluciones del sistema pueden ponerse también en forma de vectores co-
lumna:
⎛ x 1⎞ ⎛4⎞
⎜ ⎟ ⎛ –2 ⎞
x
⎜ ⎟ 2 ⎜ 5 ⎟ ⎜0⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
x = ⎜ x 3⎟ = ⎜ 0 ⎟ + t ⎜ 1 ⎟ .
⎜ ⎟ ⎜ ⎟
⎜ x 4⎟ ⎜ –5 ⎟ ⎜0⎟
⎜ ⎟ ⎝ ⎠ ⎝0⎠
⎝ x 5⎠ 2
Ejemplo 14
⎛7 –3 –1 0⎞
⎜ ⎟
⎜4 2 8 5⎟ .
⎝3 –5 –9 7⎠
⎛ 7 –3 –1 0⎞
⎜ ⎟
⎜ 0 26 ⁄ 7 60 ⁄ 7 5⎟ .
⎝0 –26 ⁄ 7 –6 0 ⁄ 7 7⎠
© FUOC • P00/75004/00191 33 Matrices, vectores y sistemas de ecuaciones lineales
Tomamos como pivote el elemento de la segunda fila y segunda columna a22 = 26/7 y redu-
cimos la segunda columna:
⎛ 7 –3 –1 0⎞
⎜ ⎟
⎜ 0 26 ⁄ 7 60 ⁄ 7 5 ⎟ .
⎝0 0 0 12⎠
Ahora el resultado consiste en que la tercera ecuación es imposible y, por lo tanto, el sistema
no tiene ninguna solución. Es un sistema incompatible.
Actividades
12.
x + y + 5z = 0 ⎫
⎪
2 x – y + 3 z = 1 ⎬.
⎪
– 3x + 4y – z = 0 ⎭
13.
5x + 2y – z = 2 ⎫
⎪
x – y + 3 z = 8 ⎬.
– 2 x + 3 y + 2 z = –3 ⎪⎭
14.
x + z– t = 2
⎫
2x – y + z= 1 ⎪
⎪
x + 3y + 4z – 7t = 11⎬.
⎪
– 4x + 3y – z – 2t = 1 ⎪
– 2x + 2y – t = 2 ⎭
15.
4x – 2y – z = 2 ⎫
⎪
3x + 5y + 3z = 2 ⎬.
⎪
x – 7y – 4z = 2 ⎭
16.
2 x1 + 2 x2 + 2 x3 – 4 x4 – 4 x5 = 8 ⎫
⎪
2 x 4 + 4 x 5 = –2 ⎪
⎪
x 1 + 4 x 2 + 4 x 3 + 4 x 4 + 4 x 5 = 1 ⎬.
⎪
x 1 + 3 x 2 + 4x 3 x2 4 x3 = 3 ⎪
⎪
3 x1 + 3 x2 + 3 x3 + 4 x4 + 4 x5 = 1 ⎭
© FUOC • P00/75004/00191 34 Matrices, vectores y sistemas de ecuaciones lineales
⎛ v11⎞ ⎛ v21 ⎞ ⎛ vk 1⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
v12 v22 vk 2
v 1 = ⎜⎜ ⎟⎟ , v 2 = ⎜⎜ ⎟⎟ , … , v k = ⎜⎜ ⎟⎟ .
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ v 1n⎠ ⎝ v2n ⎠ ⎝ vk n⎠
⎛ v 11 v 12 … v 1n⎞
⎜ ⎟
⎜ v 21 v 22 … v 2n⎟
. (17)
⎜ ⎟
...
...
...
⎜ ⎟
⎝ v k1 vk 2 … v kn⎠
Teniendo en cuenta que los restantes vectores no cambian con una segunda
transformación, los subespacios generados por {v1 , ..., vi , ... v j, ..., v k } y por
{v 1, ..., v′i , ... v ′j , ..., v k } son los mismos, ya que toda combinación lineal de
vectores del primer conjunto puede ponerse como combinación lineal de los
vectores del segundo conjunto y recíprocamente.
⎛ z0 ... z
0z0 ... z
0z0 ... z
0z0 ... z
0⎞
⎜ ⎟
⎜ 0 0 ... z0 ... zz
0 0 ... zz
0 0 ... z
0⎟
⎜ 0 0 ... 0 0 ... z0 ... z
0z 0⎟
0 ... z
⎜ ⎟
⎜ 0 0 ... 0 0 ... 0 0 ... z ... z⎟
(18)
⎜ ⎟
⎜ 0 0 ... 0 0 ... 0 0 ... 0 0 ... 0 ⎟
⎜ ⎟
...
...
...
...
...
...
...
...
⎜ ... ⎟
⎝ 0 0 ... 0 0 ... 0 0 ... 0 0 ... 0⎠
donde los son los pivotes que son estrictamente diferentes de cero y z son
números que pueden ser diferentes de cero o no. Los restantes elementos ha-
brán quedado reducidos forzosamente a cero.
Los vectores de la fila que han quedado son, por construcción, un conjunto
de generadores del subespacio que genera {v1 , ..., vk}.
Además, los vectores-fila que quedan forman una base del subespacio genera-
do.
Ejemplo 15
⎛ 1⎞ ⎛1 ⎞ ⎛ 1⎞ ⎛ –1⎞ ⎛ 2⎞ ⎛1⎞
⎜ –1 ⎟ ⎜–1 ⎟ ⎜ 0⎟ ⎜ 2⎟ ⎜ –2 ⎟ ⎜–1⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜ 2⎟ ⎜2 ⎟ ⎜ 2⎟ ⎜ –2⎟ ⎜ 4⎟ ⎜2⎟
v1 = ⎜ ⎟, v2 = ⎜ ⎟, v3 = ⎜ ⎟ , v4 = ⎜ ⎟ , v5 = ⎜ ⎟, v6 = ⎜ ⎟
⎜ 0⎟ ⎜0 ⎟ ⎜ 3⎟ ⎜ 4⎟ ⎜ 0⎟ ⎜0⎟
⎜ 2⎟ ⎜3 ⎟ ⎜ 6⎟ ⎜ 4⎟ ⎜ 5⎟ ⎜1⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ 3⎠ ⎝5 ⎠ ⎝ 1⎠ ⎝ –6⎠ ⎝ 8⎠ ⎝1⎠
© FUOC • P00/75004/00191 36 Matrices, vectores y sistemas de ecuaciones lineales
⎛ 1 –1 2 0 2 3 ⎞
⎜ 1 –1 2 0 3 5 ⎟
⎜ ⎟
⎜ 1 0 2 3 6 1 ⎟
⎜ ⎟
⎜ –1 2 –2 4 4 –6 ⎟
⎜ ⎟
⎜ 2 –2 4 0 5 8 ⎟
⎝ 1 –1 2 0 1 1⎠
y aplicamos el método de Gauss. Esta matriz es la misma del ejemplo 13, donde hemos visto
que el resultado era:
⎛1 –1 2 0 2 3 ⎞
⎜0 1 0 3 4 –2 ⎟⎟
⎜
⎜0 0 0 1 2 –1 ⎟
⎜ ⎟
⎜0 0 0 0 1 2 ⎟
⎜ ⎟
⎜0 0 0 0 0 0 ⎟
⎝0 0 0 0 0 0 ⎠
Por lo tanto, el rango de {v1 , ...,v6 } es 4 y, además, la reducción gaussiana nos permite deter-
minar una base del subespacio que generan:
⎛ 1⎞ ⎛ 0⎞ ⎛ 0⎞ ⎛ 0⎞
⎜ –1⎟ ⎜ 1⎟ ⎜ 0⎟ ⎜ 0⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜ 2⎟ ⎜ 0⎟ ⎜ 0⎟ ⎜ 0⎟
v 1′ = ⎜ ⎟ , v 2 ′ = ⎜ ⎟ , v 3 ′ = ⎜ ⎟ , v 4 ′ = ⎜ ⎟ .
⎜ 0⎟ ⎜ 3⎟ ⎜ 1⎟ ⎜ 0⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜ ⎟2 ⎜ ⎟4 ⎜ ⎟2 ⎜ 1⎟
⎝ 3⎠ ⎝ –2⎠ ⎝ –1⎠ ⎝ 2⎠
...
...
...
...
...
...
...
...
⎜ ⎟
⎝ 0 0 ... 0 0... 0 0 ... 0 0 ... 0 0⎠
Recordemos que...
El rango de una matriz es el máximo número de filas linealmente in-
... para determinar el rango de
dependientes que contiene, y es igual, a la dimensión del subespacio ge- una matriz podemos emplear
nerado por las filas consideradas como vectores-fila. el método de eliminación de
Gauss.
© FUOC • P00/75004/00191 37 Matrices, vectores y sistemas de ecuaciones lineales
Actividad
17. Comentad las actividades 10, 11, 12 , 13 y 14 desde el punto de vista del teorema de
Rouché-Fröbenius.
© FUOC • P00/75004/00191 38 Matrices, vectores y sistemas de ecuaciones lineales
6. Sistemas homogéneos
a 11 x 1 + a 12 x 2 + … + a 1n x n = 0 ⎫
⎪
a 21 x 1 + a 22 x 2 + … + a 2n x n = 0 ⎪
⎬ (19)
⎪
.. .
.. .
.. . ⎪
a m 1 x 1 + a m2 x 2 + … + a mn x n = 0 ⎭
A · x = 0.
⎛x 1⎞ ⎛0 ⎞
⎜ ⎟ ⎜ ⎟
x
⎜ ⎟ = ⎜0 ⎟
2
⎜ ⎟ ⎜ ⎟
...
...
⎜ ⎟ ⎜ ⎟
⎝x n⎠ ⎝0 ⎠
o, brevemente, x = 0.
Para los sistemas homogéneos, la matriz ampliada contiene sólo una columna
de ceros más que la matriz del sistema, de forma que el rango de la matriz am-
pliada será igual al rango de la matriz del sistema. La primera parte del teorema
de Rouché, aplicado ahora, nos dice que siempre tendrá solución, algo eviden-
te ya que x = 0 es una solución trivial. Pero en el caso de los sistemas homo-
géneos, las soluciones que interesan son precisamente las que no son triviales.
⎛ x 1⎞ ⎛ x 11 ⎞
⎜ ⎟ ⎜ ⎟
x x
x = ⎜ 2⎟ = ⎜ 12 ⎟
⎜ ⎟ ⎜ ⎟
...
...
⎜ ⎟ ⎜ ⎟
⎝ x n⎠ ⎝ x 1n⎠
© FUOC • P00/75004/00191 39 Matrices, vectores y sistemas de ecuaciones lineales
⎛ x 1⎞ ⎛ x 11 ⎞
⎜ ⎟ ⎜ ⎟
⎜ x 2⎟ = t ⎜ x 12 ⎟ .
⎜ ⎟ ⎜ ⎟
...
...
⎜ ⎟ ⎜ ⎟
⎝ x n⎠ ⎝ x 1n⎠
x = t1 x 1 + t 2 x 2 + … + t s x s
⎛ x 1⎞ ⎛ x 11⎞ ⎛ x 21 ⎞ ⎛ x s1 ⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜ 2⎟ = t ⎜ 12⎟ + t ⎜ 22 ⎟ + … + t ⎜ x s2 ⎟ .
x x x
(20)
⎜ ⎟ 1 ⎜ ⎟ 2 ⎜ ⎟ s ⎜ ⎟
...
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ x n⎠ ⎝ x 1n⎠ ⎝ x 2n⎠ ⎝ x sn ⎠
Ejemplo 16
Resolvamos el sistema
5x + 2y – 2w = 0 ⎫
z+
⎪
x – y + 3 z – 7w = 0 ⎬
⎪
2 x + 5 y – 1 0z + 23 w = 0 ⎭
⎛ 5 2 –1 2 0⎞
⎜ ⎟
⎜1 –1 3 –7 0 ⎟ .
⎝ 2 5 –10 23 0 ⎠
Resulta más cómodo tomar el primer pivote de la segunda ecuación y permutar las filas una
y dos:
⎛1 –1 3 –7 0⎞
⎜ ⎟
⎜5 2 –1 2 0⎟ .
⎝2 5 –10 23 0⎠
⎛1 –1 3 –7 0⎞
⎜ ⎟
⎜0 7 –16 37 0⎟ .
⎝0 7 –16 37 0⎠
© FUOC • P00/75004/00191 40 Matrices, vectores y sistemas de ecuaciones lineales
⎛ 1 –1 3 –7 0⎞
⎜ ⎟
⎜ 0 7 –16 37 0⎟ .
⎝ 0 0 0 0 0⎠
w = t1 ⎫
⎪
z = t2 ⎪
⎬.
y = 16 ⁄ 7 ⋅ z – 37 ⁄ 7 ⋅ w = 16 ⁄ 7 ⋅ t 2 – 37 ⁄ 7 ⋅ t 1 ⎪
⎪
x = y – 3 z + 7 w = –5 ⁄ 7 ⋅ t 2 + 12 ⁄ 7 ⋅ t 1 ⎭
⎛x⎞ ⎛ 12 ⁄ 7 ⎞ ⎛ –5 ⁄ 7⎞ ⎛ 12 ⎞ ⎛ –5 ⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
y
⎜ ⎟ = t ⎜ – 37 ⁄ 7 ⎟+t ⎜ 16 ⁄ 7 ⎟ = t ′⎜ –37 ⎟ + t ′ ⎜ 16⎟ .
⎜z⎟ 1⎜
0 ⎟ 2⎜
1 ⎟ 1 ⎜
0 ⎟ 2 ⎜
7⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ w⎠ ⎝ 1 ⎠ ⎝ 0 ⎠ ⎝ 7 ⎠ ⎝ 0⎠
Actividades
2x – y + 3z = 0 ⎫
⎪
– x + 3y + 7z = 0 ⎬.
⎪
– x + 8 y + 24 z = 0 ⎭
x + y + 5z = 0 ⎫
⎪
2 x + 2 y + 10 z = 0 ⎬ .
⎪
– 3 x – 3 y – 15 z = 0 ⎭
x + y + 5z = 0 ⎫
⎪
4x + 7y – 3z = 0 ⎬.
⎪
– x – 3y + 2z = 0 ⎭
x + y + 5z – 3t = 0 ⎫
⎪
2x – 3y + z + 2t = 0 ⎪
⎬.
4 x – y + 11 z – 4 t = 0 ⎪
x – 4 y – 4 z – t = 0 ⎪⎭
t1 x 1 + t2 x 2 + … + ts x s
también es solución del sistema para todo conjunto de valores t 1, t 2, ..., t s reales.
© FUOC • P00/75004/00191 41 Matrices, vectores y sistemas de ecuaciones lineales
b1 a 12⎞
x 1 det A = det ⎛⎝ .
b2 a ⎠
22
a 11 b 1⎞
x 2 det A = det ⎛⎝ .
a 21 b⎠ 2
© FUOC • P00/75004/00191 42 Matrices, vectores y sistemas de ecuaciones lineales
Ahora, si:
a 11 a 12
det A = ≠ 0, (21)
a 21 a 22
b1 a 12 a 11 b1
b2 a 22 a 21 b2
x 1 = -------
---------------------
----- x 2 = -------
-----------------------
---.
a 11 a 12 a 11 a 12
a 21 a 22 a 21 a 22
Ejemplo 17
2x + 5y = 7 ⎫
⎬.
3x – 4y = 4 ⎭
∆ = 2 5 = –23 ≠ 0.
3 –4
Por lo tanto el sistema es compatible y determinado, y las soluciones son: El uso matemático de DD
Ejemplo 18
Dado el sistema:
3 kx + 8 y = 2 ⎫
⎬
4 x + 7 y = –5 ⎭
discutid para qué valores del parámetro k el sistema es compatible y determinado, y en este
caso dad la solución.
© FUOC • P00/75004/00191 43 Matrices, vectores y sistemas de ecuaciones lineales
3k 8
∆ = = 21 k – 32.
4 7
Por lo tanto, si:
32
k ≠ ------
21
2 8 3k 2
– 1 5k – 8
---5
--------------------7 54 4 5
x = --– –
--- = ----------
------------- y = ---------------------------- = ------------------------.
∆ 2 1k – 32 ∆ 2 1k – 3 2
⎛ c 11 c 12 c 13⎞
C = ⎜⎜ c 21 c 22 c 23⎟⎟ ,
⎝ c 31 c 32 c 33⎠
c11 c 12 c 13
det C = c21 c 22 c 23 = c11 γ 11 – c 21 γ21 + c 31 γ31 =
c31 c 32 c 33
c 22 c 23 c12 c 13 c12 c 13
= c11 – c 21 + c 31 . (22)
c 32 c 33 c32 c 33 c22 c 23
La definición (22) o (23), aunque parezca arbitraria, es de hecho consecuencia Ved los apartados 7.1 y 7.3 de este
mismo módulo didáctico.
de generalizar, al caso 3 × 3, las fórmulas para aislar las incógnitas que hemos
visto para sistemas 2 × 2.
Ejemplo 19
1 4 –3
∆ = 3 –2 1.
4 –1 5
Solución
∆ = 1 ( –2 ) 5 + 3 ( –1 ) ( –3 ) + 1 ⋅ 4 ⋅ 4– 4( – 2) ( –3 ) – 1 ( –1 ) ⋅ 1 –3 ⋅ 4 ⋅ 5 = –68.
⎛ a 11 a 12 … a 1n ⎞
⎜ ⎟
⎜ a 21 a 22 … a 2n ⎟
⎜ ⎟.
⎜ ⎟
⎝ an 1 a n2 … a nm⎠
σ
det A = ∑ ( –1 ) a 1i 1 a 2i 2 … a ni n (24)
σ
( 1, 2 , 3 )+ , ( 1, 3 , 2 )– , ( 2, 1 , 3 ) – , ( 2, 3 , 1) + , (3 , 1 , 2 ) + , ( 3, 2, 1 )–
det A = a 11a 22 a 23 – a 11 a 23 a 32 – a 12 a 21 a 33 +
+ a 12 a 23 a 31 – a 13 a 21 a 32 – a 13 a 22 a 31 ,
Las propiedades que ahora daremos de los determinantes son válidas para los
determinantes de cualquier orden. Las demostraciones de estas propiedades,
hechas a partir de la definición (24), son un poco técnicas y no las ponemos
en el texto para no alargarlo innecesariamente. Podéis consultarlas en la bi-
bliografía.
+ – + …
– + – ….
+ – + …
...
...
...
det C = c i1 C i1 + c i2 C i2 + … + c in C in = c1 i C 1i + c2 i C 2i + … + c ni C ni .
c11 c 12 … c 1n c 11 c21 … c n1
c21 c 22 … c 2n c 12 c22 … c n2
= .
...
...
...
...
...
...
c n1 c n2 … c nn c 1n c2 n … cn n
3) Permutar filas o columnas. Al intercambiar entre sí dos filas o bien dos co-
lumnas el valor de un determinante cambia de signo.
c11 … c1 i … c 1j … c 1n
c21 … c2 i … c 2j … c 2n
=
...
...
...
...
c n1 … c ni … c nj … c nn
c11 … c1 j … c 1i … c 1n
c21 … c2 j … c 2i … c 2n
= – .
...
...
...
...
c n1 … cn j … c ni … c nn
c 11 … c 1i … c 1i … c 1n
c 21 … c 2i … c 2i … c 2n
= 0.
...
...
...
...
cn 1 … c ni … c ni … c nn
a 11 a 12 … a1 n ′
a 11 ′
a 12 … a 1′ n
c 21 c 22 … c 2n c 21 c 22 … c 2n
+ =
...
...
...
...
...
...
c n1 c n2 … cn n c n1 c n2 … cn n
...
...
cn 1 c n2 … c nn
λc 11 λ c12 … λ c 1n c 11 c12 … c1 n
c 21 c22 … c2 n c 21 c22 … c2 n
= λ .
...
...
...
...
...
...
c n1 cn 2 … c nn c n1 cn 2 … c nn
Teniendo en cuenta esto (utilizando esta propiedad para cada una de las n filas
o bien n columnas de A) resultará:
det ( λ A ) = λn ⋅ d e tA .
7) Añadir una fila (o una columna) multiplicada por un escalar a una pa-
ralela. Un determinante no varía si añadimos a una fila (o una columna) los
elementos de una paralela multiplicados por una misma constante:
c 11 c12 … c1 n + λc 1j c 11 c 12 … c 1n
c 21 c22 … c2 n + λc 2j c 21 c 22 … c 2n
= si j ≠ n.
...
...
...
...
...
...
c n1 cn 2 … c nn + λ c nj c n1 c n2 … c nn
a i1 A j1 + a i2 A j2 + … + a in A jn = 0 si i ≠ j.
© FUOC • P00/75004/00191 48 Matrices, vectores y sistemas de ecuaciones lineales
Ejemplo 20
1 –2 4
∆ = 3 –5 1.
2 7 2
El resultado es:
1 –2 4
∆ = 0 1 –11 .
0 11 –6
∆ = 1 –1 1 = – 6 + 121 = 115.
11 –6
Este método es mucho más eficiente para calcular determinantes que (23). La
utilidad de las fórmulas (24) o (23) es más teórica que práctica. Sirve para de-
mostrar resultados u obtener expresiones compactas, pero se utiliza poco para
calcular efectivamente determinantes.
Las propiedades mencionadas también nos permiten asegurar que cierto de-
terminante es cero sin necesidad de desarrollarlo.
Ejemplo 21
Probemos que:
1 –2 –1
∆ = 3 –5 –2 = 0 .
2 7 9
Solución: en efecto, teniendo en cuenta que la tercera columna es la suma de las dos prime-
ras, y utilizando primero la propiedad (7) y después la (5) obtenemos que el resultado es 0,
porque se trata de la suma de dos determinantes nulos (es decir, que valen 0).
© FUOC • P00/75004/00191 49 Matrices, vectores y sistemas de ecuaciones lineales
a 11 x 1 + a 12 x 2 + … + a 1n x n = b 1 ⎫
⎪
a 21 x 1 + a 22 x 2 + … + a 2n x n = b 2 ⎪
⎬.
...
...
⎪
...
...
a n 1 x 1 + a n 2 x 2 + … + a nn x n = b n ⎪⎭
( a 11 A 11 + a 21 A 21 + … + a n1 A n 1 )x 1 +
+ ( a 12 A 11 + a 22 A 21 + … + a n2 A n1 )x 2 + … +
+ ( a 1n A 11 + a 2 n A 21 + … + a nn A n1 )x n =
= b 1 A 11 + b 2 A 21 + … + b n A n1 .
b1 a 12 … a1 n
b2 a 22 … a2 n
x 1 det A = .
...
...
...
bn a n2 … a nn
Podemos repetir el cálculo anterior con los adjuntos de las columnas segunda,
tercera, ..., i-ésima, respectivamente, y si:
det A ≠ 0 (25)
∆
x i = ----i (26)
∆
donde:
a 11 … a 1, i – 1 b1 a 1, i + 1 … a 1n
a 21 … a 2, i – 1 b2 a 2, i + 1 … a 2n
∆i =
...
...
...
...
...
a n1 … a n, i – 1 bn a n, i + 1 … a nn
© FUOC • P00/75004/00191 50 Matrices, vectores y sistemas de ecuaciones lineales
a 11 a 12 … a1 n
a 21 a 22 … a2 n
∆ = (27)
...
...
...
a n1 a n2 … a nn
Observad que:
Ejemplo 22
x + 2y – z = 1 ⎫
⎪
– 5x + 4y + 3z = 5 ⎬.
⎪
2x + y + 2z = 6 ⎭
1 2 –1
∆ = –5 4 3 = 1 ⋅ 4 ⋅ 2 + ( –5 ) ⋅ 1 ⋅ ( –1 ) + 2 ⋅ 3 ⋅ 2 – 2 ⋅ 4 ⋅ ( –1 ) –
2 1 2
–2 ⋅ ( –5 ) ⋅ 2 –1 ⋅ 3 ⋅ 1 = 50 ≠ 0
1 2 –1
x = --1 40
---- 5 4 ----- = 4
3 = - ---
50 50 5
6 1 2
1 1 –1
1 48 24
y = ------ –5 5 3 = -
----- = ------
50 50 25
2 6 2
1 2 1
1 86 43
z = ------ –5 4 ----- = ------.
5 = -
50 50 25
2 1 6
Veamos un ejemplo 4 × 4:
Ejemplo 23
x + y + 3z – 2t = 3 ⎫
⎪
2x – 4 y + 7 z + 2 t = –1 ⎪
⎬.
3x – 2y + 9z – t = 6 ⎪
x + 3 y – z – t = 0 ⎪⎭
© FUOC • P00/75004/00191 51 Matrices, vectores y sistemas de ecuaciones lineales
1 1 3 –2
2 –4 7 2
∆ = .
3 –2 9 –1
1 3 –1 –1
1 1 3 –2
–6 1 6
∆ = 0 –6 1 6 =
–5 0 5.
0 –5 0 5
2 –4 1
0 2 –4 1
Obtenemos:
0 1 6
∆ = 0 0 5 = 3 ⋅ 5 = 15 ≠ 0.
3 –4 1
3 1 3 –2 0 –1 1 24 4
∆ 1 = –1 –4 7 2 = –1 –4 7 2 =
6 –2 9 –1 0 –26 51 11
0 3 –1 –1 0 3 –1 –1
–11 24 4 61 24 –20
= –2 6 51 11 = 127 51 –40 =
3 –1 –1 0 –1 0
61 –20 = 100,
=
127 –4 0
100 20
x = ---------- = ------.
15 3
Finalmente,
16 17 89
y = –------, z = –------, t = –------.
3 5 15
Como puede comprobarse, habría resultado mucho más corto hacer los cálculos con el mé-
todo de Gauss.
Hemos visto que la teoría de los determinantes permite una formulación com-
pacta para resolver los sistemas n × n. Pero el precio que se paga por este mé-
todo es bastante elevado, ya que para resolver efectivamente un sistema es
mucho menos eficiente que el método de Gauss.
© FUOC • P00/75004/00191 52 Matrices, vectores y sistemas de ecuaciones lineales
Teorema 5
Para utilizar este teorema debemos encontrar un menor no nulo de orden tan
grande como seamos capaces, algo que nos indicará que los vectores de la fila
implicados son linealmente independientes. A continuación, vamos amplián-
dolo con una fila y una columna, de modo que el nuevo determinante menor
continúe siendo diferente de cero. Cuando todos los menores orlados del dado
sean nulos, habremos determinado el rango.
Actividades
23. Utilizando determinantes, ved qué podemos decir a propósito de un sistema homo-
géneo de n ecuaciones y n incógnitas:
a 11 x 1 + a 12 x 2 + … + a 1n x n = 0 ⎫
⎪
a 21 x 1 + a 22 x 2 + … + a 2n x n = 0 ⎪
⎬. (28)
⎪
...
...
...
a n1 x 1 + a n2 x 2 + … + a nn x n = 0 ⎪⎭
b) En el caso de que tenga, y que ésta tenga a su vez 1 grado de libertad, ¿qué expresión
podemos dar para la solución?
24. ¿Para qué valores de λ el sistema homogéneo siguiente tiene soluciones no triviales?
En estos casos, encontradlas.
( – 2 0 – λ) x + 9y – 5z = 0 ⎫
⎪
– 51 x + ( 2 3 – λ )y – 12 z = 0 ⎬.
⎪
3x – y + ( 2 – λ )z = 0 ⎭
2 x – 3y – z = 2 ⎫
⎪
x + 2y – 3 z = 4 ⎪
⎬.
x + 2y – z = 3 ⎪
λ x – 2 y + 2 z = 1 ⎪⎭
26. Discutid las soluciones del sistema siguiente, según los valores del parámetro λ:
λx + y + z = 1 ⎫
⎪
x + λy + z = λ ⎬.
⎪
x + y + λz = λ2⎭
27. Determinad λ para que el sistema tenga soluciones no triviales, y dad las soluciones
correspondientes en forma de vectores-columna.
λx + y + z = 0 ⎫
⎪
3 x + λ y – 2 z = 0 ⎬.
⎪
2x + 4 y + z = 0 ⎭
28.
a) Demostrad que:
1 x x2
1 y y
2 = ( x – y ) ( y – z ) ( z – x ).
2
1 z z Determinante
de Vandermonde
b) Generalizando el método anterior, calculad el determinante siguiente, conocido
como determinante de Vandermonde: Determinante que tiene unos
en la primera fila, mientras que
la segunda fila es arbitraria y las
2 n –1
1 x1 x1 … x1 siguientes están formadas por
2 n –1 las potencias sucesivas de los
1 x2 x2 … x2 . elementos de la segunda fila.
Lo mismo se puede aplicar en
...
...
...
...
a 11 a 12 a 13 … a1 n
0 a 22 a 23 … a2 n
0 0 a 33 … a 3n .
...
...
...
...
0 0 0 … an n
© FUOC • P00/75004/00191 55 Matrices, vectores y sistemas de ecuaciones lineales
8. Matriz inversa
Diremos que la matriz cuadrada A tiene inversa si hay una matriz, que
denotaremos como A –1, la cual verifica:
A ⋅ A–1 = A –1 ⋅ A = I.
Teniendo en cuenta que el determinante de un producto de matrices es igual Recordad la novena propiedad de
los determinantes
al producto de los determinantes y que el determinante de I es 1, para que A
tenga inversa es necesario que:
det A ≠ 0. (29)
⎛ A 11 A 21 … A n1 ⎞
⎜ ⎟
1 A 12 A 22 … A n2 ⎟
= ---------------⎜⎜
–1
A ⎟ (30)
det A
...
...
...
⎜ ⎟
⎝ A 1n A 2n … A nn ⎠
formada por los adjuntos de los elementos (transpuestos) divididos por el de-
terminante es inversa de A .
La fórmula (39) de la matriz inversa tiene, como siempre, más utilidad teórica
que práctica, ya que resulta muy costoso hacer el cálculo de la inversa a partir
de ésta. Podemos emplear también el método de Gauss.
⎛ 1⎞ ⎛ 0⎞ ⎛0 ⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎜ 0⎟ , ⎜ ⎟ , …, ⎜ 0 ⎟
1
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ 0⎠ ⎝ 0⎠ ⎝1 ⎠
⎛ a 11 … a 1n ⎞ ⎛ x 11 … x 1n ⎞ ⎛1 … 0⎞
⎜ ⎟ ⋅⎜ ⎟ = ⎜ ⎟
...
...
...
...
...
...
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ a n1 … a nn ⎠ ⎝ x n1 … x nn⎠ ⎝0 … 1⎠
o, en forma de matrices:
A ⋅X =I (31)
⎛ a 11 … a 1n 1 … 0⎞
⎜ ⎟.
...
...
...
⎜ ⎟
...
⎝a … a nn 0 … 1⎠
n1
⎛1 … 0 x 11 … x 1 n⎞
⎜ ⎟,
...
...
...
...
⎜ ⎟
⎝0 … 1 x n1 … x nn⎠
Ejemplo 24
⎛1 1 1⎞
A = ⎜5 4 5 ⎟⎟ .
⎜
⎝2 –1 1⎠
© FUOC • P00/75004/00191 57 Matrices, vectores y sistemas de ecuaciones lineales
⎛ 1 1 1 1 0 0⎞
⎜ 0 ⎟⎟ .
⎜ 5 4 5 0 1
⎝ 2 –1 1 0 0 1⎠
⎛ 1 1 1 1 0 0⎞
⎜ –5 1 0⎟⎟ .
⎜ 0 –1 0
⎝ 0 0 –1 13 –3 1⎠
⎛1 0 0 9 –2 1⎞
⎜ 5 –1 0⎟⎟ .
⎜0 1 0
⎝0 0 1 –13 3 –1⎠
⎛ 9 –2 1⎞
⎜ 5 –1 0⎟⎟ .
⎜
⎝ –13 3 –1⎠
Vemos así que el coste de encontrar la inversa de una matriz n × n por el mé-
todo de Gauss es resolver simultáneamente n sistemas por Gauss con la propia
matriz del sistema, que es muy inferior a n veces el coste de resolver un siste-
ma; esto se debe a que, para cada sistema, únicamente debemos añadir los cál-
culos relativos a la columna correspondiente de la matriz ampliada.
A·x=b
ya que multiplicando esta ecuación por A−1 por la izquierda obtenemos lo si-
guiente:
x = A−1 · b
Actividades
⎛ 1 –2 4⎞
⎜ 2 5 –3 ⎟⎟ .
⎜
⎝ –3 1 4⎠
31.
⎛ cos α –s e n α⎞
.
⎝ senα cos α ⎠
© FUOC • P00/75004/00191 58 Matrices, vectores y sistemas de ecuaciones lineales
b) Comprobad que A · A t = A t · A = I.
c ) Las matrices que verifican la propiedad anterior reciben el nombre de ortogonales. ¿Es
posible conseguir que la matriz siguiente sea ortogonal eligiendo a, b y c de forma adecuada?
⎛3 ⁄ 5 –4 ⁄ 13 a⎞
⎜4 ⁄ 5 3 ⁄ 13 b ⎟⎟ .
⎜
⎝ 0 12 ⁄ 13 c⎠
© FUOC • P00/75004/00191 59 Matrices, vectores y sistemas de ecuaciones lineales
Resumen
Por otro lado, la teoría de los determinantes permite una formulación teórica más
directa, pero los cálculos son mucho más complicados. Para un sistema n × n de-
terminado, las fórmulas de Cramer nos dan una expresión directa de las solucio-
nes. A pesar de ello, en el momento de calcular determinantes también utilizamos
el método de Gauss.
Los determinantes resultan también útiles para saber si una matriz tiene inver-
sa, y hemos llegado a la conclusión de que, para ello, es necesario únicamente
que su determinante sea diferente de cero.
© FUOC • P00/75004/00191 61 Matrices, vectores y sistemas de ecuaciones lineales
Ejercicios de autoevaluación
3x + y + 2z = 19 ⎫
⎪
x + 2 y + 3 z = 27 ⎬ .
⎪
4 x + 2 y + 5z = 42 ⎭
3x + y + 2z = 19 ⎫
⎪
x + 2 y + 3 z = 27 ⎬ .
⎪
4 x + 3 y + 5z = 42 ⎭
3 x + y + 2z = 19 ⎫
⎪
x + 2 y + 3 z = 27 ⎬ .
⎪
4 x + 3 y + 5z = 46 ⎭
3x + y + 2z = 0 ⎫
⎪
x + 2y + 3z = 0 ⎬.
⎪
4 x + 2 y + 5z = 0 ⎭
3x + y + 2z = 0 ⎫
⎪
x + 2y + 3z = 0 ⎬.
⎪
4 x + 3 y + 5z = 0 ⎭
2 . Discutid el sistema siguiente según el valor del parámetro m. Resolvedlo por los valores de
m que hacen que sea compatible:
–x + mz = m ⎫
⎪
x + y + 3 z = 5 ⎬.
⎪
2x + my + = 0 ⎭
a x + by + z = 1 ⎫
⎪
x + aby + z = b ⎬.
x + b y + a z = 1 ⎪⎭
4. Consideremos los vectores v1 = (1, 2, 0, –1), v2 = (2, –1, 1, 3), v3 = (–1, 3, –2, 2) y v4 = (4, 3, 0, 7).
⎛ 9 –2 –4 ⎞
A = 1
---⎜ 1 2 –1 ⎟
5⎜ ⎟
⎝ –12 1 7 ⎠
6 . Dada la matriz
2 –1 ⎞
A = ⎛⎝
4 –2 ⎠
© FUOC • P00/75004/00191 62 Matrices, vectores y sistemas de ecuaciones lineales
7 . Lo mismo que en el ejercicio anterior, suponiendo que ahora A es tal que det( A)≠0 .
8 . Una empresa tiene dos centros de producción, uno en Lérida y otro en Mataró. En cada
uno de los centros se han hecho las compras siguientes:
• En Lérida se han comprado 1 tonelada de tinta, 2 toneladas de papel, 3 toneladas de
disolvente y 5 toneladas de madera.
• En Mataró se han comprado 2 toneladas de tinta, 6 toneladas de papel, 4 toneladas de di-
solvente y 3 toneladas de madera. Sabiendo que los precios para cada artículo son: tinta 3
u.m. (unidades monetarias), papel 2 u.m., disolvente 5 u.m. y madera 2 u.m.
A B C
Mezclando las tres aleaciones podemos obtener diferentes productos con varias proporcio-
nes de plata, cobre y oro.
a) Encontrad una matriz tal que, multiplicándola por los kilogramos de A, de B y de C que
mezclemos, nos proporcione los kilogramos de plata, cobre y oro que hay en la mezcla fi-
nal. ¿Cuál es la suma de cada una de sus columnas?
b) Mezclando 50 kilogramos de A, 20 de B y 20 de C, cuántos kilogramos de plata, cobre y oro
habrá en la mezcla resultante? ¿Cuál será la proporción de plata, cobre y oro en la mezcla final?
c) Mezclando 0,5 kilogramos de A, 0,2 de B y 0,3 de C, ¿cuántos kilogramos de plata, cobre y
oro habrá en la mezcla resultante? ¿Cuál será la proporción de plata, cobre y oro en la mezcla
final?
d) ¿Podemos obtener una mezcla con el 20% de plata y el 80% de oro? ¿Qué interpretación
de la solución podemos dar cuando alguna incógnita toma valor negativo?
e) ¿Podemos obtener una mezcla con el 11,5% de plata, el 30,5% de cobre y el 58% de oro?
10. Supongamos que las composiciones de las aleaciones del ejercicio anterior son:
A B C
A B C
para obtener una mezcla que tenga el 12% de plata, n% de cobre y el resto de oro. Plantead
el sistema correspondiente y discutidlo según el valor de los parámetros m y n. Resolvedlo
cuando sea posible.
© FUOC • P00/75004/00191 64 Matrices, vectores y sistemas de ecuaciones lineales
Solucionario
Actividades
1 . Tenemos que:
2 5 ⎞ ⎛ 3 –1 ⎞ 26 38 ⎞
A ⋅ B = ⎛⎝ = ⎛⎝ .
–1 3 ⎠ ⎝ 4 8 ⎠ 9 25 ⎠
3 –1 ⎞ ⎛ 2 5 ⎞ 7 12 ⎞
B ⋅ A = ⎛⎝ = ⎛⎝ .
4 8 ⎠ ⎝ –1 3 ⎠ 0 44 ⎠
2.
526 –194 ⎞
2 A – 5 A + 4 A + 3 I = ⎛⎝
3 2
–97 41 ⎠
7 –1 ⎞ ⎛ 7 –2 ⎞ 50 –1 6 ⎞
A A = ⎛⎝ = ⎛⎝
t
.
–2 2 ⎠ ⎝ –1 2 ⎠ –16 8 ⎠
∑a
n
Si consideramos ahora la matriz (A · B) t, tenemos que ((A · B )t) ij = (A · B) ji = jk b k i.
k= 1
Sin embargo:
∑ (B ) ∑b ∑a
t t n t t n n
( B A ) ij = ik ( A )k j = ki a jk = jk b k i.
k = 1 k = 1 k= 1
con lo cual ((A · B )t) ij = (Bt · A t)ij y las matrices (A·B ) t y Bt·A t son iguales.
4 . Tenemos:
t t t
A · A = B1 y A · A = B2
6×4 4× 6 6× 6 4× 6 6×4 4×4
Por lo tanto, los dos productos se pueden calcular. Observemos que, por las propiedades
de la matriz transpuesta:
t t t t t t t
B 1 = ( A A) = A (( A) ) = A A = B 1
5.
⎛ a b ⎞ ⎛ -----------
1------------⎛ d –b ⎞⎞ -----------1
------------⎛
a b ⎞ ⎛ d –b ⎞
⎝ c d ⎠ ⎝ ( a d – b c ) ⎝–c a ⎠⎠ = ( ad – bc ) ⎝ c d ⎠ ⎝ –c a ⎠ =
1------------⎛ a d – b c 0 ⎞ 1 0⎞
= ----------- = ⎛ .
( a d – bc ) ⎝ 0 a d – b c⎠ ⎝0 1 ⎠
6.
⎛ a 11 a 12 a 1⎞ ⎛ x⎞
( x y 1 ) ⎜⎜ a 12 a 22 a 2⎟⎟ ⎜⎜ y⎟⎟ =
⎝ a 1 a 2 a ⎠ ⎝ 1⎠
⎛ a 11 x + a 12 y + a 1⎞
= ( x y 1 )⎜ a 12 x + a 22 y + a 2⎟ =
⎜ ⎟
⎝ a1 x + a2y + a ⎠
2 2
= a 11 x + 2 a 12 xy + a 22 y + 2 a 1 x + 2a 2 y + a .
© FUOC • P00/75004/00191 65 Matrices, vectores y sistemas de ecuaciones lineales
7 . Recordemos que ((AB )C)ij se obtiene haciendo el producto de la fila i de la matriz (AB ) por
la columna j de C y, por lo tanto, podemos escribir:
∑ ( A B)
n
( ( A B )C) ij = ( A B ) i1 C 1j + ... + ( A B ) in C nj = im C m j.
m = 1
∑a
n
( A B ) im = a i1 b 1m + a i2 b 2m + ... + a im b n m = ik bk m ;
k= 1
⎛ ⎞
( ( A B )C ) ij = ∑ (A B )
n
m = 1
im Cmj = ∑ ⎜⎝ ∑ a
n
m= 1 k = 1
n
ik b km ⎟ c m j .
⎠
Haciendo los productos y reorganizando los sumatorios (lo cual puede hacerse por la pro-
piedad distributiva del producto respecto a la suma), obtenemos:
⎛ ⎞ ⎛ ⎞
( ( A B )C ) ij = ∑ ⎜⎝ ∑ a
n
m = 1 k= 1
n
ik b km c m j⎟ =
⎠
∑ a ⎜⎝ ∑ b
n
k= 1
ik
n
m = 1
km c m ⎟j =
⎠
∑a
n
= ik ( B C) k j = ( A (B C ) ) ij.
k = 1
8.
a) La ecuación λ 1 e1 + λ 1 e1 = 0 se escribe:
λ 1 ( 2 , 0 , –3 ) + λ 2 ( –1 , 3 , 5 ) = 0
2 λ1 + λ2 = 0
3 λ2 = 0
– 3 λ1 + 5 λ2 = 0
e 1 = 2 u1 – 3 u 3
y así:
u1 = 1
---e 1 + 3
---u 3 .
2 2
1 3
v = v1 u 1 + v 2 u 2 + v 3 u 3 = v 1 ⎛⎝ ---e 1 + --u 3 ⎞⎠ + v 2 u 2 + v 3 u 3 =
2 2
1 3
= ---v1 e 1 + v 2 u 2 + ⎛⎝ ---v1 + v3 ⎞⎠ u 3 .
2 2
1 3 1 7
e 2 = ---( –1 )e 1 + 3u 2 + ⎛ ---( –1 ) + 5⎞ u 3 = –---e1 + 3 u 2 + ---u 3 .
2 ⎝2 ⎠ 2 2
1 1 7
u 2 = ---e 1 + ---e 2 – ---u 3 .
6 3 6
1--v e + v u + ⎛ 3
v= - --v + v 3⎞⎠ u 3 =
2
1 1 2 2 ⎝-2 1
= 1
---v 1 e 1 + v 2 ⎛ 1
---e + 1
--e – 7
---u ⎞ + ⎛ 3
---v + v 3⎞ u 3 =
2 ⎝ 6 1 3 2 6 3⎠ ⎝ 2 1 ⎠
1 1 1 ⎛3 7
= ⎛⎝ ---v 1 + ---v 2 ⎞⎠ e 1 + ---v2 e 2 + --v – --v + v 3 ⎞⎠ u 3 .
2 6 3 ⎝-
2 1 6 2
u1 = 1
---e 1 + 3
---u 3
2 2
1 1 1
u 2 = ---e 1 + ---e 2 – ---u 3 .
6 3 7
9.
a) El segundo vector debe pertenecer al subespacio generado por { e1,e 2}. Por lo tanto, debe
ser combinación lineal:
e ′2 = λ ( 2, 0, –3 ) + 2 λ ( –1 , 3, 5 ) = λ ( 0, 6, 7 )
b)
1 1 1
e 2 = ---( e ′2 –λe 1 ) = – ---e 1 + ---e ′2 .
µ 2 2
10.
a) Se toman dos elementos del conjunto A = {(x, y, z) ∈ R 3 : –3x + 2 y + z = 0}:(u 1,u2 ,u3) y
(v1,v2 ,v3).
• Debemos comprobar si se verifica que (u 1,u2 ,u3) + (v1 ,v2,v3) = (u 1 + v1, u2 + v′2 u 3+v3) ∈ A.
Para ello, es necesario que veamos que: –3 ( u 1 + v1 ) + 2( u 2 + v 2) + u 3 + v3 = 0 .
Efectivamente:
– 3 u 1 + 2 u 2 + u 3 – 3 v 1 + 2 v 2 + v3 = 0 porque –3 u 1 + 2 u 2 + u 3 = 0 y –3 v 1 + 2 v2 + v 3 = 0
• Debemos comprobar si k (u 1, u2, u 3) ∈ A. Para ello es necesario ver que –3 ku1 + 2ku2 + ku3 = 0.
Efectivamente:
– 3 k u1 + 2 k u 2 + k u3 = k ( – 3 u 1 + 2u 2 + u 3 ) = 0 .
5 ( u 1 + v1 ) + 2 ( u 2 + v 2 ) – 7 ( u 3 + v 3 ) = 3 .
5 u 1 + 2 u 2 – 7 u 3 + 5 v 1 + 2 v2 – 7 v3 = 3 + 3 = 6;
c) Se toman dos elementos del conjunto C = {( x, y, z) ∈ R 3 : x ≥ 0}: (u 1 , u2, u3) y (v1 , v2, v3).
• Debemos comprobar si se verifica que ( u1 + v 1, u2 + v 2u 3 + v3) ∈ C. Para ello, es necesario
ver que:
( u1 + v 1 ) ≥ 0 .
11.
a) Se toman dos elementos del conjunto A = {f ∈ F(R,R) : f(1) = k}:f1 y f2.
• Queremos ver si f1 + f 2 ∈ A. Para ello es necesario que se verifique que (f1 + f 2)(1) = k.
Tenemos
( f 1 + f 2 ) ( 1 ) = f 1 ( 1 ) + f 2 ( 1) = k + k = 2k .
Así pues, sólo se verifica en el caso k = 0. Por lo tanto, en general, A no es subespacio vec-
torial de F(R,R).
( f1 + f 2 ) ( 0 ) = f 1 ( 0 ) + f 2 ( 0 ) = 2 f 1 ( 1 ) + 2 f 2 ( 1) = 2 ( f 1 + f 2 ) ( 1 ).
• Queremos comprobar si kf 1 ∈ B. Para ello es necesario ver si (kf1)(0) = 2(kf 1)(1). Efecti-
vamente:
( k f 1 ) ( 0 ) = kf 1 ( 0 ) = k2 f 1 ( 1) = 2 ( kf 1 ) ( 1 ).
Tenemos
2 2 2
( f 1 + f 2 ) ( –x ) = f1 ( –x ) + f 2 ( –x ) = f 1 ( x ) + f 2 ( x ) ≠ ( f 1 + f 2 ) ( x ).
12.
⎛ 1 1 5 0⎞
⎛1 1 5 0⎞ ⎛ 1 1 5 0⎞
⎜ 2 –1 3 1 ⎟ → ⎜ 0 –3 –7 1 ⎟ → ⎜⎜ 0 –3 –7 1⎟⎟
⎜ ⎟ ⎜ ⎟
⎝ –3 4 –1 0 ⎠ ⎝ 0 7 1 4 0⎠ ⎜ 7 7⎟
⎝ 0 0 –-3-- -3-⎠
y, por lo tanto, el sistema tiene como solución única x = 3; y = 2; z = –1. Podemos compro-
bar que el resultado es correcto sustituyendo los valores encontrados por las incógnitas en
cada una de las ecuaciones, con lo que obtenemos:
x + y + 5 z = 3 + 2 + 5 ⋅ ( –1) = 0 (correcto),
2 x – y + 3 z = 2 ⋅ 3 – 2 + 3 ⋅ ( –1 ) = 1 (correcto),
– 3 x + 4 y – z = ( –3 ) ⋅ 3 + 4 ⋅ 2 –( –1) = 0 (correcto).
13. Observad las permutaciones entre las filas y cómo facilitan los cálculos:
© FUOC • P00/75004/00191 68 Matrices, vectores y sistemas de ecuaciones lineales
⎛ 5 2 –1 2⎞ ⎛ 1 –1 3 8⎞ ⎛ 1 –1 3 8⎞
⎜ 1 –1 3 8⎟ → ⎜ 5 2 –1 2⎟ → ⎜ 0 7 –16 –38⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ –2 3 2 –3⎠ ⎝ –2 3 2 – 3⎠ ⎝ 0 1 8 13⎠
⎛ 1 –1 8 ⎞⎟ ⎛ 1 –1 3 8 ⎞
⎜ 3 ⎜ ⎟
→⎜ 0 1 8 13 ⎟ → ⎜ 0 1 8 13 ⎟
⎜ ⎟ ⎜ ⎟
⎝ 0 7 – 1 6 – 38 ⎠ ⎝ 0 0 – 7 2 –129 ⎠
31 4 43
x = ------, y = –---, z = ------.
24 3 24
14.
⎛ 1 0 1 1 2⎞ ⎛1 0 1 –1 2⎞
⎜ 2 –1 1 0 1⎟ ⎜0 – 1–1 2 –3 ⎟
⎜ ⎟ ⎜ ⎟
⎜ 1 3 4 –7 11⎟ → ⎜ 0 3 3 –6 9 ⎟ →
⎜ ⎟ ⎜ ⎟
⎜ –4 3 –1 –2 1 ⎟ ⎜0 3 3 –6 9 ⎟
⎝ –2 2 0 –1 2 ⎠ ⎝0 2 2 –3 6 ⎠
⎛ 1 0 1 – 1 2⎞ ⎛ 1 0 1 0 2⎞
⎜ 0 –1 –1 2 –3⎟ ⎜ 0 1 1 0 3⎟
⎜ ⎟ ⎜ ⎟
→ ⎜ 0 0 0 0 0 ⎟ → ⎜0 0 0 1 0 ⎟
⎜ ⎟ ⎜ ⎟
⎜0 0 0 0 0 ⎟ ⎜0 0 0 0 0 ⎟
⎝0 0 0 1 0 ⎠ ⎝0 0 0 0 0 ⎠
x = 2 – z ; y = 3 – z ; t = 0;
15.
⎛ 4 –2 –1 2⎞ ⎛ 1 –7 –4 2 ⎞
⎜ 3 5 3 2⎟ → ⎜ 3 5 3 2⎟ →
⎜ ⎟ ⎜ ⎟
⎝ 1 –7 –4 2⎠ ⎝ 4 –2 –1 2 ⎠
⎛1 –7 –4 2⎞ ⎛1 –7 –4 2⎞
→ ⎜⎜ 0 26 15 –4⎟⎟ → ⎜⎜ 0 26 15 –4⎟⎟
⎝0 26 15 –6 ⎠ ⎝0 0 0 –2⎠
16.
⎛2 2 2 –1 –1 8⎞ ⎛1 1 1 1 1 1⎞
⎜0 0 0 –2 ⎟1 ⎜0 0 1 0 1 1 –2 ⎟
⎜ ⎟ ⎜ ⎟
⎜1 1 1 1 1 1⎟ → ⎜2 2 2 –1 –1 8⎟ →
⎜ ⎟ ⎜ ⎟
⎜1 1 1 0 0 3⎟ ⎜1 1 1 0 0 3⎟
⎝3 3 3 4 4 1⎠ ⎝3 3 3 4 4 1⎠
⎛1 1 1 1 1 1⎞ 1 1 1 1 1 1
⎛ ⎞
⎜0 0 0 1 1 –2 ⎟ ⎜ 0 0 0 1 1 –2⎟
⎜ ⎟ ⎜ ⎟
⎜0 0 0 –3 –3 6 ⎟ → ⎜0 0 0 0 0 0⎟ →
⎜ ⎟ ⎜ ⎟
⎜0 0 0 –1 –1 2⎟ ⎜0 0 0 0 0 0⎟
⎝0 0 0 1 1 –2 ⎠ ⎝0 0 0 0 0 0⎠
⎛1 1 1 0 0 3⎞
⎜0 0 0 1 1 –2⎟
⎜ ⎟
⎜0 0 0 0 0 0⎟
⎜ ⎟
⎜0 0 0 0 0 0⎟
⎝0 0 0 0 0 0⎠
© FUOC • P00/75004/00191 69 Matrices, vectores y sistemas de ecuaciones lineales
( x 1 , x 2 , x 3, x 4, x 5 ) = ( 3 – x 2 – x 3, x 2, x 3 , –2 – x 5 , x 5 ) =
·
= ( 3, 0, 0, 0, –2, 0 ) + x 2 ( –1, 1, 0, 0, 0 ) + x 3 ( –1 , 0, 1, 0, 0 ) + x 5 ( 0, 0, 0, –1, 1 ) .
17.
10 3 3 3 Determinada 0
11 3 3 3 Determinada 0
12 4 3 3 Indeterminada 1
13 3 2 3 Incompatible -
14 5 2 2 Indeterminada 3
18.
⎛ 2 –1 3 0 ⎞ ⎛ –1 3 7 0 ⎞
⎜ –1 3 7 0 ⎟⎟ → ⎜⎜ 2 –1 3 0 ⎟⎟ →
⎜
⎝ –1 8 24 0 ⎠ ⎝ –1 8 24 0 ⎠
⎛ –1 3 7 0 ⎞ ⎛ –1 3 7 0 ⎞
⎜ 0 5 17 0 ⎟⎟ → ⎜⎜ 0 5 17 0 ⎟⎟
⎜
⎝ 0 5 17 0 ⎠ ⎝ 0 0 0 0 ⎠
⎛ x⎞ ⎛ –16 ⎞
⎜ y⎟ = ⎜ –17 ⎟ t .
⎜ ⎟ ⎜ ⎟
⎝ z⎠ ⎝ 5⎠
19. Teniendo en cuenta que todas las ecuaciones son múltiples de la primera, el sistema se
reduce a (1 1 5 | 0), y la solución es x = – y –5 z.
20.
⎛1 1 5 0 ⎞
⎛ 1 1 5 0⎞ ⎛1 1 5 0⎞
⎜ ⎟
⎜ 4 7 –3 0 ⎟⎟ → ⎜⎜ 0 ⎟ 23
⎜ 3 –23 0 ⎟ → ⎜0 3 – 0 ⎟.
⎝ –1 –3 ⎜ 25 ⎟
2 0⎠ ⎝ 0 –2 7 0⎠ ⎝ 0 0 –--3---- 0 ⎠
21.
⎛1 1 5 –3 0 ⎞ ⎛1 1 5 –3 0⎞
⎜ ⎟ ⎜ ⎟
⎜2 –3 1 2 0 ⎟
→ ⎜⎜
0 –5 –9 8 0⎟
→
⎜4 –1 11 –4 ⎟
0 ⎟ 0 ⎟⎟
⎜0 –5 –9 8
⎜
⎝1 –4 –4 –1 0 ⎠ ⎝0 –5 –9 2 0⎠
⎛1 1 5 –3 0 ⎞ ⎛ x⎞ ⎛ –16⎞
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
0 –5 –9 8 0 ⎟ y –9⎟
→⎜ →⎜ ⎟ =⎜ s.
⎜0 0 0 0 0 ⎟ ⎜ z ⎟ ⎜ 5⎟
⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝0 0 0 –6 0 ⎠ ⎝ t⎠ ⎝ 0⎠
22. Si x 1, ..., xn son soluciones del sistema A · x = 0, querrá decir que A · x 1 = 0, ..., A · xn = 0
y, por lo tanto:
A ( t 1 x 1 + t 2 x 2 + ... + t n x n ) = A t 1 x 1 + A t 2 x 2 + ... + At n x n =
= t 1 A x 1 + t 2 Ax 2 + ... + t n Ax n = t 1 0 + t 2 0 + ... + t n 0 = 0 .
© FUOC • P00/75004/00191 70 Matrices, vectores y sistemas de ecuaciones lineales
23.
a) En primer lugar, es necesario que el determinante del sistema sea nulo (det A = 0), ya
que, si es diferente de cero, sería determinado y sólo tendría la solución trivial nula.
En segundo lugar, si det A = 0, entonces una solución es:
x 1 = A 11 t , x 2 = A 12 t , x n = A 1n t
donde A1 j es el adjunto del elemento a 1j, ya que en la primera ecuación, al sustituir las x,
tendremos:
a 11 A 11 t + a 12 A 12 t + ... + a 1n A 1n t = t det A = 0
y para las otras aparecerá la suma de productos de los elementos de una fila por los adjun-
tos de una paralela que, por la propiedad 8 de los determinantes, es cero.
Puede suceder, sin embargo, que todos estos adjuntos sean 0, lo cual nos indicaría que:
24. Para que el sistema tenga soluciones no triviales es necesario que el determinante del sis-
tema sea cero, a fin de que el rango sea inferior a 3. La condición es, de este modo:
( – 20 – λ ) 9 –5
3 2
–51 23 – λ –12 = – λ + 5 λ – 8 λ + 4 = 0.
3 –1 2 –λ
– 21 x + 9 y – 5 z = 0 ⎫
⎪
– 51 x + 22 y – 12 z = 0 ⎬ .
⎪
3x – y + z = 0⎭
o bien:
⎛ x⎞ ⎛ 2⎞
⎜ y⎟ = ⎜ 3 ⎟ t .
⎜ ⎟ ⎜ ⎟
⎝ z⎠ ⎝ –3⎠
– 22 x + 9 y – 5 z = 0 ⎫
⎪
– 51 x + 2 1y – 12 z = 0 ⎬.
⎪
3x – y = 0⎭
21 – 12 –5 1 –12 –51 21
x = t = –12 t , y = – t = –36 t, z = t = –12 t
–1 0 3 0 3 –1
o bien:
⎛ x⎞ ⎛ 1⎞
⎜ y⎟ = ⎜ 3 ⎟ t .
⎜ ⎟ ⎜ ⎟
⎝ z⎠ ⎝ 1⎠
© FUOC • P00/75004/00191 71 Matrices, vectores y sistemas de ecuaciones lineales
2 –3 –1 2
1 2 –3 4
= – 13 λ – 26,
1 –2 –1 3
λ –2 2 1
por lo tanto, este determinante se anula sólo en el caso de que λ = –2. Podemos observar
que el rango de A es 3, independientemente del valor del número real λ.
Observemos que si λ ≠ –2, el rango de la matriz ampliada será 4 y el rango de la matriz del
sistema será, como máximo, 3 (ya que la matriz del sistema sólo tiene tres columnas). Ten-
dremos, de este modo, un sistema incompatible.
Para el caso λ = –2, sustituimos este valor en el sistema y lo resolvemos siempre (es un sis-
tema “normal”: no hay parámetros, ni lambdas, ni nada). Os daréis cuenta de que es un
sistema compatible indeterminado.
λ 1 1
3 2
1 λ 1 = λ – 3 λ + 2 = (λ – 1 ) (λ + 2 ) .
1 1 λ
⎛ –2 1 1 1⎞ ⎛1 –2 1 –2⎞
⎜ 1
⎜ –2 1 –2⎟⎟ → ... → ⎜⎜ 0 –3 3 –3⎟⎟ ,
⎝ 1 1 –2 4⎠ ⎝0 0 0 4⎠
⎛ x⎞ ⎛ –3 ⁄ 7⎞
• Para λ = –5 (resolviéndolo por Gauss) obtenemos ⎜⎜ y⎟⎟ = ⎜⎜ 1⎟⎟ t.
⎝ z⎠ ⎝ –22 ⁄ 7⎠
• Para λ diferente de estos valores el sistema (¡es homogéneo!) tiene únicamente la solu-
ción trivial.
28.
a) Observad que:
• restando a la segunda columna la primera multiplicada por x,
• restando a la tercera columna la segunda multiplicada por x,
• desarrollando el determinante por la primera columna,
• aplicando la propiedad 6 de los determinantes, para extraer del determinante los esca-
lares que multiplican todos los elementos de las filas segunda y tercera,
• calculando el determinante de orden 2 que nos queda y arreglando los signos,
obtenemos:
2
1 x x 1 0 0
y –x y2 – x y
1 y y 2 = 1 y–x y2 – xy = 2 =
2 z –x z –xz
1 z z 2 1 z–x z – xz
y –x y (y – x ) 1 y
= = (y – x ) (z – x ) =
z –x z (z – x ) 1 z
= ( y – x )(z – x ) (z – y ) = (x – y ) (z – x )( y – z ) .
© FUOC • P00/75004/00191 72 Matrices, vectores y sistemas de ecuaciones lineales
b) Para el caso general, los mismos pasos que hemos seguido cuando hay tres columnas
nos proporcionan determinantes cada vez menores hasta llegar al valor del determinante
grande, que es:
⎧
⎪
⎨
⎪
⎩
29. El resultado es a 11 · a 22 · ... · a n n, es decir, el producto de los elementos de la diagonal.
30.
⎛ 1 –2 4 1 0 0⎞ ⎛1 –2 4 1 0 0⎞
⎜ 2
⎜ 5 –3 0 1 0⎟⎟ → ⎜⎜ 0 9 –11 –2 1 0⎟⎟ →
⎝–3 1 1 0 0 1⎠ ⎝0 –5 16 3 0 1⎠
⎛1 –2 4 1 0 0⎞ ⎛1 –2 4 1 0 0 ⎞
→ ⎜⎜ 0 9 –11 –2 1 0⎟⎟ → ⎜⎜ 0 1 –11 ⁄ 9 2 ⁄9 1⁄9 0 ⎟ →
⎟
⎝0 0 89 ⁄ 9 17 ⁄ 9 5⁄9 1⎠ ⎝0 0 1 17 ⁄ 89 5 ⁄ 89 9 ⁄ 89 ⎠
⎛ 2 3 1 2 –14⎞
--1
----⎜ 1 1 6 11⎟ .
89 ⎜ ⎟
⎝17 5 9⎠
31.
a)
⎛ cos ( α ) – sin ( α )⎞ – 1
= ⎛⎝
cos ( α ) sin ( α )⎞
.
⎝ sin ( α ) cos( α )⎠ – sin ( α ) cos ( α )⎠
a = 48
------ b = –36
------ c = --5
---- o bien a = –48
------ b = 36
------ c = –--5
----.
65 65 13 65 65 13
Ejercicios de autoevaluación
1.
⎛ 1⎞
a) ⎜ 4⎟ .
⎜ ⎟
⎝ 6⎠
b) Incompatible.
⎛ 11 ⁄ 5⎞ ⎛ –1 ⁄ 5 ⎞
c) ⎜ 65 ⁄ 5⎟ + t ⎜⎜ –7 ⁄ 5 ⎟⎟ .
⎜ ⎟
⎝ 0⎠ ⎝ 1⎠
⎛ 0⎞
d) ⎜⎜ 0⎟⎟ .
⎝ 0⎠
⎛ –1 ⁄ 5⎞
e) t ⎜⎜ –7 ⁄ 5⎟⎟ .
⎝ 1⎠
2.
• Si m ≠ 0 y m ≠ –1, el sistema es compatible y determinado con soluciones:
( x, y , z ) = ⎛ ---2
----m –---4
------, m 3--⎞ .
+-----
-------, ----- -------
⎝m + 1 m + 1 m + 1⎠
© FUOC • P00/75004/00191 73 Matrices, vectores y sistemas de ecuaciones lineales
( x, y , z ) = ( 0,5,0 ) + t ( 0 , –3 , 1) .
• Si m = –1 el sistema es incompatible.
3.
• Si b ≠ 0, a ≠ 1 y a ≠ –2, el determinante de la matriz del sistema es diferente de cero y,
por lo tanto, el sistema es compatible y determinado.
• Si b = 0, el sistema es incompatible independientemente del valor de a.
• Si a = 1 y b ≠ 1, el sistema es incompatible.
• Si a = 1 y b = 1, el sistema es compatible e indeterminado.
• Si a = –2 y b ≠ –2, el sistema es incompatible.
• Si a = –2 y b = –2, el sistema es compatible e indeterminado.
4.
a) Para que la condición λ 1 v1 + λ 2 v 2 + λ 3 v3 + λ 4 v 4 = 0 implique λ 1 = λ 2 = λ3 = λ 4 = 0, el
sistema homogéneo:
λ 1 + 2 λ2 – λ 3 + 4 λ4 = 0 ⎫
⎪
2 λ1 – λ 2 + 3 λ3 + 3 λ4 = 0 ⎪
⎬
λ2 – 2 λ3 = 0⎪
– λ 1 + 3 λ 2 + 2 λ 3 + 7 λ 4 = 0 ⎪⎭
⎛ 1 2 –1 4⎞
⎜ ⎟
⎜ 2 –1 3 3⎟
⎜ 0 1 –2 0⎟⎟
⎜
⎝ –1 3 2 7⎠
por lo tanto, los vectores dados serán linealmente independientes si det A ≠ 0 y linealmen-
te dependientes si det A = 0. En los apartados siguientes veremos que estamos en la segun-
da opción y, por lo tanto, los vectores dados son linealmente dependientes.
b) Para determinar el rango y obtener una base del subespacio generado por { v1 ,v2,v 3,v4 },
aplicamos el método de Gauss:
⎛ 1 2 0 –1⎞ ⎛1 2 0 –1⎞
⎜ ⎟ ⎜ ⎟
⎜ 2 –1 1 3⎟
→ ⎜0 –5 1 5⎟
→
⎜ –1 3 –2 2⎟⎟ ⎜0 5 –2 1⎟⎟
⎜ ⎜
⎝ 4 3 0 7⎠ ⎝0 –5 0 11⎠
⎛1 2 0 –1⎞ ⎛1 2 0 –1⎞
⎜ ⎟ ⎜ ⎟
0 –5 1 5⎟ ⎜0 –5 1 5⎟
→ ⎜⎜ → .
0 0 –1 6⎟ ⎜0 0 –1 6⎟
⎜ ⎟ ⎜ ⎟
⎝0 0 –1 6⎠ ⎝0 0 0 0⎠
e 1 = ( 1, 2, 0, –1 ) e 2 = ( 0 , –5 , 1 , 5) e 3 = ( 0, 0, –1, 6) .
c) Sea ahora la base canónica { u1 ,u 2,u 3,u 4}. Podemos sustituir u 1 por e1, ya que su compo-
nente según u 1 es 1 ≠ 0. Tenemos:
e 1 = u 1 + 2· u 2 – u 4, es decir: u1 = e1 – 2u 2 + u 4
w = w1 u1 + w2 u2 + w 3u 3 + w4 u4
w = w 1 ( e 1 – 2 ⋅ u 2 + u 4 ) + w 2 u 2 + w 3 u 3 +w 4 u 4 =
= w 1 e 1 + ( –2 w 1 + w 2 ) u 2 + w 3 u 3 + ( w 1 + w 4 ) u 4 .
© FUOC • P00/75004/00191 74 Matrices, vectores y sistemas de ecuaciones lineales
e 2 = 0 ⋅ e 1 + ( –2 ⋅ 0 – 5) ⋅ u 2 + 1 ⋅ u 3 + ( 0 + 5 ) ⋅ u 4 = – 5 u 2 + u 3 + 5 u 4
1 1
u 2 = –--e 2 + ---u 3 + u 4 .
5 5
1
w = w 1 e 1 + ( –2 w 1 + w 2 ) ⎛⎝ –1
---e 2 + ---u 3 + u 4⎞⎠ + w 3 u 3 + ( w 1 + w 4 )u 4 =
5 5
2 1 2 1
= w 1 e 1 + ⎛⎝ ---w 1 – ---w 2 ⎞⎠ e 2 + ⎛⎝ –---w 1 + ---w 2 + w 3⎞⎠ u 3 + ( –w 1 + w 2 + w 4 )u 4.
5 5 5 5
e 3 = –u 3 + 6 u 4
u 3 = –e 3 + 6 u 4 .
2 1 1
w = w 1 e 1 + ⎛⎝ ---w 1 –---w 2 ⎞⎠ e 2 + ⎛⎝ –2
---w 1 + ---w 2 + w 3⎞⎠ ( – e 3 + 6 u 4 ) + ( – w 1 + w 2 + w 4 ) u 4 =
5 5 5 5
2 1 2 17 11
= w 1 e 1 + ⎛⎝ -- w 1 –-- w 2⎞⎠ e 2 + ⎛⎝ --w 1 –1
---w 2 –w 3 ⎞⎠ e 3 + ⎛⎝ – ------w 1 + ------w 2 + 6 w 3 + w 4⎞⎠ u 4.
5 5 5 5 5 5
d) Aplicando la fórmula anterior podemos expresar los vectores v1,v 2,v3 ,v4 en la base
{e1,e 2,u 3,u 4 }:
2 2 2 2 17 22
v1 = ( 1, 2, 0, 1) = e 1 + ⎛⎝ --- – ---⎞⎠ e 2 + ⎛⎝ --- – -- – 0⎞⎠ e 3 + ⎛⎝ –------ + ------ + 0 – 1⎞⎠ u 4 = e 1
5 5 5 5 5 5
4 1 4 1 34 11
v2 = ( 2, –1, 1, 3 ) = 2 e1 + ⎛⎝ --- + ---⎞⎠ e2 + ⎛⎝ --- + --- – 1⎞⎠ e 3 + ⎛⎝ –------ – ------+ 6 + 3⎞⎠ u 4 = 2 e 1 + e 2 .
5 5 5 5 5 5
2 3 2 3 17 33
v 3 = ( –1, 3, –2 , 2 ) = –e 1 + ⎛⎝ – --- – ---⎞⎠ e 2 + ⎛⎝ – -- – --- + 2⎞⎠ e 3 + ⎛⎝ ------ + ------ – 12 + 2⎞⎠ u 4 = –e1 – e 2 + e 3.
5 5 5 5 5 5
v 4 = ( 4, 3, 0, 7) = 4 e 1 + ⎛ 8 ---⎞ e + ⎛ 8
--- – 3 -- + 0⎞ e 3 + ⎛ –68
--- – 3 ------+ 7⎞ u 4 = 4 e 1 + e 2 + e 3 .
------+ 33
⎝ 5 5⎠ 2 ⎝ 5 5 ⎠ ⎝ 5 5 ⎠
5.
⎛3 2 2⎞
= ⎜⎜ 1 1⎟⎟ .
–1
A 3
⎝5 3 4⎠
6.
a b⎞ –2 c c⎞ –2 b b⎞
B = ⎛⎝ C = ⎛⎝ D = ⎛⎝
2a 2 b⎠ –2 d d⎠ –4 b 2b⎠
8.
⎛3 ⎞
⎜2 ⎟
1 2 3 5⎞
a) A = ⎛ y B = ⎜ ⎟
⎝2 6 4 3⎠ ⎜5 ⎟
⎜ ⎟
⎝2 ⎠
© FUOC • P00/75004/00191 75 Matrices, vectores y sistemas de ecuaciones lineales
donde en la primera fila de A se encuentran las compras del centro de Lérida y en la se-
gunda, las del centro de Mataró, y la matriz B consta de los precios de cada uno de los pro-
ductos. Si multiplicamos:
⎛3 ⎞
1 2 3 5⎞ ⎜⎜2 ⎟⎟ 32
A ⋅ B = ⎛⎝ ⋅ = ⎛⎝ ⎞⎠
2 6 4 3 ⎠ ⎜5 ⎟ 44
⎜ ⎟
⎝2 ⎠
b) Si C = (1 1),
32
C ⋅ (A ⋅ B ) = ( 1 1) ⋅ ⎛⎝ ⎞⎠ = ( 76 ) .
44
1 2 3 5⎞
H = C ⋅ A = ( 1 1) ⋅ ⎛⎝ = (3 8 7 8 ) .
2 6 4 3⎠
⎛3 4⎞
⎜ ⎟
2 2⎟
E = ⎜ tiene los precios de cada mes en una columna.
⎜5 6⎟
⎜ ⎟
⎝2 7⎠
32 61⎞
A ⋅ E = ⎛⎝ nos da el gasto de cada centro en cada compra.
44 65⎠
e) (1 1) · (D · E ) = (76 126) es lo que han gastado entre los dos centros cada una de las
compras.
1
( ( 1 1 ) ⋅ ( D · E ) ) ⋅ ⎛ ⎞ = 7 6 + 126
⎝ 1⎠
es el gasto total.
a) El sistema será:
b) 7,5 kilogramos de plata, 20,5 kilogramos de cobre y 62 kilogramos de oro. Las propor-
ciones correspondientes son, aproximadamente, 8,33%, 22,78% y 68,89%.
c) 9%, 24,5%, 66,5%
d) Si resolvemos el sistema:
obtenemos x = –11, y = 21, z = –9. Directamente no tiene sentido considerar –11 kilos de
A, pero sí que podemos considerar que la solución consiste en tomar 21 kilos de B, extraer
11 kilos de A y 9 kilos de C. Podéis comprobar que lo que queda tiene el 20% de plata y el
80% de oro. Observad que estamos suponiendo que disponemos de un procedimiento
para extraer de una mezcla cualquiera de plata, cobre y oro ciertas cantidades de las alea-
ciones A, B y C, y esto (¡los químicos lo saben muy bien!, por cierto, ¿aún no ha llegado
el que habíamos avisado?) no siempre es posible.
e) ( 0,2; 0,3; 0,5 ).
© FUOC • P00/75004/00191 76 Matrices, vectores y sistemas de ecuaciones lineales
10. Observad que en este caso la matriz del sistema tiene determinante nulo.
a) El sistema es incompatible; por lo tanto, no hay forma de obtener una mezcla con las
características pedidas.
b) El sistema es compatible e indeterminado y, por lo tanto, tiene muchas soluciones que
se pueden expresar de esta forma:
x=5 t---– 1
--, y = 6 1
----- ----- -----– ----0t
----- ----, z = t
5 5
Por ejemplo, si n fuese 40 y m fuese 55, la solución sería (0,15; 0,3; 0,55).
• Si m = 35 y n ≠ 29 el sistema es incompatible.
• Si m = 35 y n = 29 el sistema es compatible e indeterminado con la solución siguiente:
( x, y , z) = ⎛ 0, 3
---, 2
---⎞ + t ( 1, –2 , 1 ).
⎝ 5 5⎠
Glosario
Cramer, regla de
Fórmulas cerradas, en forma de cociente de determinantes, que dan el valor de las incógnitas
en un sistema de n ecuaciones y n incógnitas compatible y determinado.
determinante
Dada una matriz cuadrada, el determinante es un número que se calcula a partir de sus ele-
mentos y que proporciona información sobre la independencia lineal de sus filas (y de sus
columnas)
Gauss, método de
Método de eliminación de incógnitas que permite reducir un sistema lineal de forma escalonada.
matriz
Tabla de números en forma de n filas y m columnas. Esto será una matriz n × m.
matriz de un sistema
Matriz formada por los coeficientes de las variables x 1, ..., x n en cada una de las ecuaciones
del sistema. Cada fila de la matriz corresponde a una ecuación.
matriz inversa
Aquella que presentan las matrices cuadradas que tienen determinante diferente de cero. Su
producto, tanto por la izquierda como por la derecha, con la matriz dada es la matriz identidad.
menor
Dada una matriz, se llama menor o determinante menor a cualquier determinante de la matriz
cuadrada que se forme con los elementos de la matriz correspondientes a k filas y k columnas,
suprimiendo las restantes. Las filas y columnas no deben coincidir necesariamente.
menor complementario
Dada una matriz cuadrada, llamamos menor complementario de un elementoa ij, al determi-
nante menor que se forma suprimiendo la fila i y la columna j.
pivote
Elemento diferente de cero que se elige para reducir a cero los restantes coeficientes en una
etapa dada del método de Gauss.
Rouché-Fröbenius, teorema de
Teorema que permite saber si un sistema lineal es compatible o no lo es. En el caso de que
sea compatible, permite saber si es determinado o indeterminado, y en el último caso, cuán-
tos grados de libertad tiene la solución general.
sistema lineal
Conjunto de ecuaciones de primer grado en m variables.
sistema compatible
Sistema de ecuaciones que tiene alguna solución. Si no tiene ninguna decimos que es incom-
patible. Por abuso de lenguaje, en el caso de un sistema homogéneo, cuando únicamente tie-
ne la solución trivial 0, decimos que es incompatible.
sistema determinado
Sistema que tiene una única solución. Si tiene más de una, decimos que es indeterminado.
sistema homogéneo
Sistema de ecuaciones sin términos independientes o, lo que es equivalente, con términos
independientes nulos.
subespacio vectorial
Subconjunto de un espacio vectorial que tiene estructura de espacio vectorial.
vector
Llamamos vectores a los elementos de un espacio vectorial. Los vectores del espacio vectorial
numérico n -dimensional están formados por una matriz de una sola columna (o una sola fi-
la) de n elementos.
Bibliografía
Anton, H. (1991). Introducción al álgebra lineal . (3.a edición). México: Ed. Limusa.
Hace un tratamiento elemental del álgebra lineal para estudiantes de primeros cursos de li-
cenciatura no específica. No requiere conocimientos de cálculo, si bien incluye ejercicios
para aquellos estudiantes que sí los tienen. No contiene todas las demostraciones, pero sí las
© FUOC • P00/75004/00191 78 Matrices, vectores y sistemas de ecuaciones lineales
que considera pedagógicas y elementales; en una parte optativa incluye otras demostraciones
no tan elementales. Es un texto bastante básico.
Birkhoff, G.; Mac Lane, S. (1985). Álgebra moderna. (3.a edición). Barcelona: Ed. Vicens-
Vives.
Es el libro de referencia general para todo el curso. Abarca todo el temario del curso y otros
temas relacionados. Se trata de un libro clásico (1941), traducido al castellano y del que se
han hecho muchas ediciones. Es un texto muy bien hecho que está a un nivel asequible para
los estudiantes que seguís este curso; contiene también muchos ejercicios propuestos.
Grossman, S.I. (1988). Álgebra lineal. (2.a edición). México: Grupo Editorial Iberoamericano.
Es un texto parecido al de Anton, pero tiene un cariz más específicamente científico. Hay
ejercicios, que están señalados, que requieren conocimientos de cálculo (análisis). Además,
contiene métodos numéricos para álgebra lineal.