lOMoARcPSD|3870267

Práctico - Ejercicios resueltos. Sistemas de ecuaciones

lineales
Algebra lineal (Universidad de Sevilla)

StuDocu no está patrocinado ni avalado por ningún colegio o universidad.

Descargado por Dario Calvache ([email protected])
 lOMoARcPSD|3870267

ÁLGEBRA LINEAL. EJERCICIOS RESUELTOS

Grado de Estadı́stica, curso 10 .

1. Sistemas de ecuaciones lineales

Ejercicio 1 Use eliminación Gaussiana con sustitución hacia atrás para resolver los siguientes sistemas:
 
 x1 +x2 +x3 = 1  2x1 −x2 =0
a) x1 +2x2 +2x3 = 1 , b) −x1 +2x2 −x3 =0 ,
x1 +2x2 +3x3 = 1 −x2 +x3 =1
 

 x1 +x 2 +x3 +x4 =1
4x2 −3x3 = 3


x1 +x2 +3x3 +3x4 =3
 
c) −x1 +7x2 −5x3 = 4 d)
x1 +x2 +2x3 +3x4 =3
−x1 +8x2 −6x3 = 5
 

x1 +3x2 +3x3 +3x4 =4


Solución 1 a) Formamos la matriz ampliada con la matriz de coeficientes del sistema y el vector de términos
independientes:  

1 1 1 1
A b = 1 2 2 1 .
1 2 3 1
Pivotamos sobre el elemento (1, 1), que es no nulo:
 

F2 −F1 ,F3 −F1 1 1 1 1
A b −−−−−−−−→  0 1 1 0 .
0 1 2 0

El elemento de la posición (2, 2) es no nulo, por lo que anulamos los elementos por debajo de él:
   
1 1 1 1 1 1 1 1
F3 −F2
 0 1 1 0  −−−→  0 1 1 0  .
0 1 2 0 0 0 1 0

Ya tenemos el sistema en forma triangular, que es equivalente a


 x1 +x2 +x3 = 1,
x2 +x3 = 0,
x3 = 0.


Ahora realizamos sustitución hacia atrás:

x2 = −x3 = 0, x1 = 1 − x2 − x3 = 1.

b) Procedemos como antes:  

2 −1 0 0
A b =  −1 2 −1 0  .
0 −1 1 1
El elemento en la posición (1, 1) es no nulo y pivotamos sobre él:
   
2 −1 0 0 2 −1 0 0
F2 +(1/2)F1
 −1 2 −1 0  −−−−−−→  0 1,5 −1 0  .
0 −1 1 1 0 −1 1 1

El elemento en la posición (2, 2) es no nulo y

   
2 −1 0 0 2 −1 0 0
F3 +(1/1,5)F2
 0 1,5 −1 0  −−−−−−−→  0 1,5 −1 0 .
0 −1 1 1 0 0 0,333333 1

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

El sistema equivalente que nos queda es


 2x1 −x2 = 0,
1,5x2 −x3 = 0,
0,33333x3 = 1.


Y la sustitución hacia atrás

1 1
x3 = 3, x1 = x3 = 2, x1 = x2 = 1.
1,5 2
c) Como antes, partimos de  

0 4 −3 3
A b =  −1 7 −5 4  .
−1 8 −6 5
El elemento en la posición (1, 1) es nulo, por lo que debemos hacer un intercambio de filas, como, por ejemplo,
la primera con la segunda.    
0 4 −3 3 −1 7 −5 4
F12
 −1 7 −5 4 −→ 0 4 −3 3  .
−1 8 −6 5 −1 8 −6 5
Ahora anulamos los elementos de la primera columna:
   
−1 7 −5 4 −1 7 −5 4
F3 −F1
 0 4 −3 3  −−−→  0 4 −3 3  .
−1 8 −6 5 0 1 −1 1
El elemento en la posición (2, 2) es no nulo y anulamos los términos por debajo de él.
   
−1 7 −5 4 −1 7 −5 4
F3 −1/4F2
 0 4 −3 3  −−−−−→  0 4 −3 3 .
0 1 −1 1 0 0 −0,25 0,25
Y realizamos la sustitución hacia atrás:
1
x3 = −1, x2 = (3 + 3x3 ) = 0, x1 = −(4 + 5x3 − 7x2 ) = 1.
4
d)  
1 1 1 1 1

 1 1 3 3 3 
A b =
 1
.
1 2 3 3 
1 3 3 3 4
El elemento (1, 1) es no nulo:
   
1 1 1 1 1 1 1 1 1 1
 1 1 3 3 3  F2 −F1 ,F3 −F1 ,F4 −F1  0 0 2 2 2 
 − −−−−−−−−−−−−→  .
 1 1 2 3 3   0 0 1 2 2 
1 3 3 3 4 0 2 2 2 3
El elemento de la posición (2, 2) es nulo y buscamos uno no nulo debajo de él. Se encuentra en la cuarta fila, por
lo que realizamos un intercambio de filas:
   
1 1 1 1 1 1 1 1 1 1
 0 0 2 2 2  F24  0 2 2 2 3 
 0 0 1 2 2 −
  → .
 0 0 1 2 2 
0 2 2 2 3 0 0 2 2 2
Ya hay ceros por debajo de la posición (2, 2) y bajamos a (3, 3), que es no nulo.
   
1 1 1 1 1 1 1 1 1 1
 0 2 2 2 3  F4 −2F2
 0 2 2 2 3 
  −−−−→  .
 0 0 1 2 2   0 0 1 2 2 
0 0 2 2 2 0 0 0 −2 −2

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Procedemos ahora a la sustitución hacia atrás:

1 1 1
x4 = 1, x3 = 2 − 2x4 = 0, x2 = (3 − 2x4 − 2x3 ) = , x1 = 1 − x4 − x3 − x2 = − .
2 2 2
Ejercicio 2 Consideremos los tres sistemas donde los coeficientes son los mismos pero los términos a la derecha son
diferentes:
4x −8y +5z = 1 0 0
4x −7y +4z = 0 1 0
3x −4y +2z = 0 0 1
Resuelva los tres sistemas a la vez realizando eliminación Gaussiana sobre una matriz de la forma
 
A b1 b2 b3 .

Solución 2 El algoritmo de eliminación Gaussiana se basa en la matriz de coeficientes, por lo que las operaciones de
filas que se hagan para resolver un sistema son las mismas para otro en lo que varı́e sea el término independiente.
Consideremos entonces la matriz [A|B], con B la matriz formada por los términos independientes, apliquemos las
transformaciones.
 
4 −8 5 1 0 0
F2 −F1 ,F3 − 43 F1  
[A|B] →  0 1 −1 −1 1 0 
 
0 2 −7/4 −3/4 0 1
 
4 −8 5 1 0 0
F3 −2F2 
 0 1 −1 −1 1 0  .

→  
0 0 1/4 5/4 −2 1

La solución de los sistemas es

z = 5, y = (−1 + z) = 4, x = 14 (1 − 5z + 8y) = 2,
z = −8, y = (1 + z) = −7, x = 14 (0 − 5z + 8y) = −4,
z = 4, y = (0 + z) = 4, x = 14 (0 − 5z + 8y) = 3.

Ejercicio 7 Mediante la resolución de un sistema de 3 ecuaciones con 3 incógnitas, calcule los coeficientes de la
ecuación de una parábola y = α + βx + γx2 que pasa por los puntos (1, 1), (2, 2) y (3, 0).

Solución 7 Tenemos tres incógnitas α, β y γ que verifican

1 = α + β + γ,
2 = α + 2β + 4γ,
0 = α + 3β + 9γ.

Aplicamos eliminación Gaussiana para obtener

   
1 1 1 1 1 1 1 1
  F2 −F1 ,F3 −F1  
 1 2 4 2 
[A|b] =  →  0 1 3 1 
  
1 3 9 0 0 2 8 −1
 
1 1 1 1
F3 −2F2  
→  0
 1 3 1 ,
0 0 2 −3

y mediante sustitución hacia atrás nos queda

3 11
γ = − , β = 1 − 3γ = , α = 1 − β − γ = −3.
2 2

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Ejercicio 13 Resuelva el siguiente sistema mediante Gauss-Jordan:

2x1 +4x2 +2x4 = 6,

+3x2 +3x3 +x4 = 4,
2x1 +7x2 +9x3 +7x4 = 8,
6x3 +5x4 = −4.

Solución 13 Trabajamos con la matriz ampliada del sistema.

   
2 4 0 2 6 1 2 0 1 3
 0 3 3 1 4  1
2 F1
 0 3 3 1 4 
[A|b] = 
 2
 − →  
7 9 7 8   2 7 9 7 8 
0 0 6 5 −4 0 0 6 5 −4
 
1 2 0 1 3
F3 −2F1  0 3 3 1 4 
−−−−→ 
 0

3 9 5 2 
0 0 6 5 −4
 
1 2 0 1 3
 0 1 4 
1
3 F2 1 1 3 3 
−→ 
 0 3 9 5 2 
0 0 6 5 −4
1 1
 
1 0 −2 3 3
 0 1 4
F3 −3F2 ,F1 −F2 1 1 3 3

−−−−−−−−−→ 
 0

0 6 4 −2 
0 0 6 5 −4
1 1
 
1 0 −2 3 3
 0 1 4
1
6 F3 1 1 3 3

−→ 2
− 13
 
 0 0 1 
3
0 0 6 5 −4
5
− 13
 
1 0 0 3
F4 −6F3 ,F1 +2F3 ,F2 −F3  0 1 0 − 31 5
3

−−−−−−−−−−−−−−→  2
− 13

 0 0 1 
3
0 0 0 1 −2
 
1 0 0 0 3
F1 − 35 F4 ,F2 + 13 F4 ,F3 − 32 F4  0 1 0 0 1 
−−−−−−−−−−−−−−−−→   0

0 1 0 1 
0 0 0 1 −2

Por tanto, la solución del sistema es x1 = 3, x2 = 1, x3 = 1, x4 = −2.

Ejercicio 16 Use el método de Gauss-Jordan para resolver a la vez los tres sistemas

2x1 −8x2 =1 0 0
−x1 +2x2 −x3 =0 1 0
−x2 +x3 =0 0 1

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Solución 16
   
2 −8 0 1 0 0 1 1 −4 0 1/2 0 0
 −1 2 −1 0 1 0  2 F1 →  −1 2 −1 0 1 0  F2 + F1 →
0 −1 1 0 0 1 0 −1 1 0 0 1
   
1 −4 0 1/2 0 0 1 −4 0 1/2 0 0 F1 + 4F2
 0 −2 −1 1/2 1 0  − 1 F2 →  0 1 1/2 −1/4 −1/2 0  →
2
0 −1 1 0 0 1 0 −1 1 0 0 1 F3 + F2
   
1 0 2 −1/2 −2 0 1 0 2 −1/2 −2 0
F − 2F
 0 1 1/2 −1/4 −1/2 0  →  0 1 1/2 −1/4 −1/2 0  1 1 3 →
2 F2 − 2 F3
0 0 3/2 −1/4 −1/2 1 3 F3 0 0 1 −1/6 −1/3 2/3
 
1 0 0 −1/6 −4/3 −4/3
 0 1 0 −1/6 −1/3 −1/3 
0 0 1 −1/6 −1/3 2/3
La solución de los sistemas es
     
−1/6 −4/3 −4/3
x(1) =  −1/6  , x(2) =  −1/3  , x(3) =  −1/3  .
−1/6 −1/3 2/3
Ejercicio 19 Calcule las soluciones de la ecuación ax2 + bx + c = 0, con aritmética de 8 dı́gitos significativos, para
a = 1, b = −105 , c = 1. Compare el resultado con la solución exacta.
Solución 19 Las raı́ces exactas, redondeadas a 11 decimales, son
x1 = 99999,99990, x2 = 0,000010000000001.
Vamos a realizar los cálculos con la fórmula tradicional
1 p 1 p
x1 = (−b + b2 − 4ac), x2 = (−b − b2 − 4ac).
2a 2a
En primer lugar, calculemos el discriminante ∆ = b2 − 4ac:
fl √
a · c = 1, b2 = 1010 , ∆ = 1010 − 4 = 999999996 → 1010 , ∆ = 105 .
Entonces
1 1
x1 = (105 + 105 ) = 105 (aceptable), x2 = (105 − 105 ) = 0(no aceptable).
2 2
De nuevo, el cálculo de x2 nos ha llevado a una cancelación catastrófica. Existen varias alternativas para calcular
las raı́ces de una ecuación cuadrática que evitan esta cancelación. Una de ellas es usar el signo del coeficiente b para
determinar cuál de las fórmulas causa la menor cancelación y de esta forma obtener una de las raı́ces. En concreto,
1 p
x1 = − (b + signo(b) b2 − 4ac).
2a
La otra raı́z se calcula como
c
x2 = .
ax1
Para este problema,
x1 = 105 , x2 = 1/105 = 10−5 ,
y ambas respuestas son aceptables.
Ejercicio 21 Vamos a resolver el sistema
+ 21 x2 + 13 x3 + 41 x4


 x1 = 1,

 1 x1 + 13 x2 + 14 x3 + 51 x4

= 1,

2
1


 3 x1 + 14 x2 + 15 x3 + 61 x4 = 1,

1
+ 15 x2 + 16 x3 + 71 x4

4 x1 = 1.


Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

con la precisión simple de Matlab. La matriz de coeficientes del sistema y el término independiente se pueden generar
con los comandos
>> A = single(hilb(4))
>> b = single(ones(4,1))

1. Aplique la eliminación gaussiana para resolver el sistema anterior con Matlab.

2. La solución con aritmética exacta es x1 = −4, x2 = 60, x3 = −180, x4 = 140. ¿Qué ha ocurrido?


Solución 21 Partimos de la matriz ampliada A b , escribimos la secuencia de comandos para obtener la matriz
en forma triangular y aplicamos sustitución hacia atrás.
>> A = single(hilb(4));
>> b = single(ones(4,1));
>> B1 = [A,b];
>> B1(2,:) = B1(2,:) - B1(2,1)/B1(1,1) * B1(1,:);
>> B1(3,:) = B1(3,:) - B1(3,1)/B1(1,1) * B1(1,:);
>> B1(4,:) = B1(4,:) - B1(4,1)/B1(1,1) * B1(1,:);
>> B2 = B1
B2 =

1.0000 0.5000 0.3333 0.2500 1.0000

0 0.0833 0.0833 0.0750 0.5000
0 0.0833 0.0889 0.0833 0.6667
0 0.0750 0.0833 0.0804 0.7500

>> B2(3,:) = B2(3,:) - B2(3,2)/B2(2,2)*B2(2,:);

>> B2(4,:) = B2(4,:) - B2(4,2)/B2(2,2)*B2(2,:);
>> B3 = B2

B3 =

1.0000 0.5000 0.3333 0.2500 1.0000

0 0.0833 0.0833 0.0750 0.5000
0 0 0.0056 0.0083 0.1667
0 0 0.0083 0.0129 0.3000

>> B3(4,:) = B3(4,:) - B3(4,3)/B3(3,3) * B3(3,:);

>> B3

B3 =

1.0000 0.5000 0.3333 0.2500 1.0000

0 0.0833 0.0833 0.0750 0.5000
0 0 0.0056 0.0083 0.1667
0 0 0 0.0004 0.0500

>> x4 = B3(4,5)/B3(4,4)

x4 =

140.0002

>> x3 = 1/B3(3,3) * (B3(3,5) - x4*B3(3,4))

x3 =

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

-180.0002

>> x2 = 1/B3(2,2) * (B3(2,5) - x4*B3(2,4) - x3*B3(2,3))

x2 =

60

>> x1 = 1/B3(1,1) * (B3(1,5) - x4*B3(1,4) - x3*B3(1,3) - x2*B3(1,2))

x1 =

-4.0000
La diferencia con la solución exacta viene determinado por los errores acumulados en el cálculo, pero en último
término es porque este es un ejemplo de matriz mal condicionada. Si, por ejemplo, se cambia el valor A(4,4) = 0.1429,
la nueva solución es x1 = −3,25, x2 = 51,003, x3 = −157,5009, x4 = 125,007. ası́, pequeños cambios en la matriz de
coeficientes producen grandes cambios en la solución.
Ejercicio 22 Consideremos el sistema
10−3 x −y =1
x +y =0
1. Calcule una solución con aritmética de 3 dı́gitos y sin pivoteo.
2. Calcule una solución con aritmética de 3 dı́gitos con pivoteo escalado.
3. Calcule una solución exacta.
4. Compare los dos primeros apartados con la solución exacta redondeada a tres dı́gitos.
Solución 22 1. El multiplicador es m = fl( 101−3 ) = 103 y las operaciones son

fl(m · 10−3 ) = 1,0, fl(−1 · m) = −103 , fl(1 · m) = 103 ,

fl(1 − 1,0) = 0, fl(1 + 103 ) = 103 , fl(0 − 103 ) = −103 .
La matriz ampliada se transforma en  
10−3 −1 1
,
0 103 −103
y la solución calculada es y = −1, x = 0.
2. Con pivoteo escalado: Tenemos la matriz  
10−3 −1 1
,
1 1 0
formamos el vector auxiliar    
10−3 10−3
s1 = 1
1 = ,
1
1
obtenido por la división de cada elemento de la primera columna por el elemento de mayor módulo de su fila
correspondiente, en valor absoluto. Vemos que la segunda entrada es mayor que la primera, por lo que efectuamos
el intercambio de las dos primeras filas:
   
10−3 −1 1 F12 1 1 0
−
→
1 1 0 10−3 −1 1
Procedemos a anular el coeficiente de la posición (2, 1):
   
1 1 0 F2 − 101 F1 1 1 0
−3

−−−−−−→
10−3 −1 1 0 fl(−1, 001) = −1 1
1 1
y resolvemos y = −1 = −1, x = 1 = 1.

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

3. La solución exacta del sistema es

1000 1000
x= ≈ 0,999000, y = − ≈ −0,999000.
1001 1001

4. La solución exacta redondeada a tres dı́gitos es xc = 0,999, yc = −0,999. El segundo método se aproxima mucho
más que el primero.

Ejercicio 23 Consideremos el sistema

x − y = −2,
−9x + 10y = 12.
1. Con aritmética de tres dı́gitos, resuelva el sistema con pivoteo parcial sin escalado.
2. Con aritmética de tres dı́gitos, resuelva el sistema con pivoteo escalado.

Solución 23 Consideremos la matriz ampliada, sobre la que vamos a efectuar las operaciones:
" #
1 −1 −2
[A|b] = .
−9 10 12

1. Con pivoteo parcial sin escalado. La primera operación es el intercambio de las filas 1 y 2, pues el elemento de
mayor módulo se encuentra en la posición (2, 1).
" #
F12 −9 10 12
[A|b] → .
1 −1 −2

fl
El multiplicador es α = −1/9 → −0,111. Entonces

fl

αF1 = 1 −10/9 −12/9 → 1 −1,11 −1,33 ,

y


F2 − αF1 = 0 −1 + 1,11 −2 + 1,33 = 0 0,11 −0,67 .
El sistema queda en la forma
−9x +10y = 12,
,
0,11y = −0,67,
de donde y = − fl(0,67/0,11) = −6,09, −9x = 12 − 10y = 72,9, x = −8,1.
2. Con pivoteo escalado. El vector  
1/1
a1 =
9/10
nos indica que no hay que hacer intercambio de filas. Hacemos la transformación F2 + 9F1 y obtenemos
 
1 −1 −2
,
0 1 −6

de donde y = −6, x = −8.

Ejercicio 38 Mediante la técnica de pivoteo parcial escalado, resuelva simultáneamente los sistemas Ax = b1 , Ax =
b2 , donde " # " # " #
2 1000 4000 4000
A= , b1 = , b2 = .
6 −1 2 3
Solución 38 Formamos la matriz ampliada
" #

2 1000 4000 4000
A b1 b2 .
6 −1 2 3

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

para la resolución simultánea, aplicamos el pivoteo escalado. Formamos el vector auxiliar

 2   
1000 2 · 10−3
s1 = 6 = ,
6 1

obtenido por la división de cada elemento de la primera columna por el elemento de mayor módulo de su fila co-
rrespondiente, en valor absoluto. Vemos que la segunda entrada es mayor que la primera, por lo que efectuamos el
intercambio de las dos primeras filas:
" # " #
2 1000 4000 4000 F12 6 −1 2 3
−
→ .
6 −1 2 3 2 1000 4000 4000

Procedemos a anular el coeficiente de la posición (2, 1):

" # " #
6 −1 2 3 F2 − 62 F1 6 −1 2 3
−−−−→ .
3001 11998
2 1000 4000 4000 0 3 3 3999

Ya tenemos la matriz en forma triangular y resolvemos los sistemas por sustitución hacia atrás. En el primer sistema,
11998 3 11998 1 3000
x2 = · = ≈ 3,998001, x1 = (2 + x2 ) = ≈ 0,999667.
3 3001 3001 6 3001
En el segundo,
3 11997 1 3500
y2 = 3999 · = ≈ 3,997667, y1 = (3 + x2 ) = ≈ 1,166278.
3001 3001 6 3001

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Sistemas rectangulares y formas escalonadas

Ejercicio 39 Supongamos que A es una matriz m × n. Explique brevemente si las siguientes afirmaciones son verda-
deras o falsas:
1. rango(A) ≤ mı́n{m, n}.
2. rango(A) < m si una fila de A es nula.

3. rango(A) < m si una fila de A es múltiplo de otra.

4. rango(A) < m si una fila es combinación lineal de las otras.
5. rango(A) < n si una columna es nula.

Solución 39 1. Verdadero. Como el rango es igual al número de filas no nulas en una forma escalonada, rango(A) ≤
m. Como también es el número de columnas básicas, rango(A) ≤ n.
2. Verdadero. Si una fila de A es nula, en la forma reducida quedará al final y el rango será menor que m.
3. Verdadero. Si como primera operación elemental anulamos la fila que es múltiplo de otra, estamos en el caso
anterior y rango(A) < m.
4. Verdadero. Como primera operación elemental, restamos a dicha fila la combinación lineal de las restantes y
anulamos la fila. Entonces rango(A) < m.
5. Verdadero. Las operaciones elementales no alteran el carácter nulo de la columna. Entonces no será columna
básica, por lo que rango(A) < n.

Ejercicio 42 Calcule la forma escalonada reducida por filas de cada una de las siguientes matrices y exprese las
columnas no básicas en función de las básicas.
 
  2 1 1 3 0 4 1
1 2 3  4
  2 4 4 1 5 5 
1 2 3 3  2 6 8   
 , c)  2 1 3 1 0 4 3
   
a)  2 4 6 9  , b)  2 6 0  6
.
  3 4 8 1 9 5 
2 6 7 6  1 2 5  
 0

0 3 −3 0 0 3 
3 8 6
8 4 2 14 1 13 3

Solución 42 a)      
1 2 3 3 1 2 3 3 1 2 3 3
     
 2 4 6 9 
A=  0 0 0 3   0 2 1 0 
→  →  
 2 6 7 6  0 2 1 0   0 0 0 3 
1 2 3 3 1 0 2 3 1 0 2 3
     
 0 1 1/2 0  →  0 1 1/2 0 
→
 0 1 1/2 0 
    
 0 0 0 3  0 0 0 3 0 0 0 1
1 0 2 0
 
 0 1 1/2 0  .
 
0 0 0 1
Por tanto, el rango de la matriz es igual a 3, las columnas básicas son A∗1 , A∗2 , A∗4 . La columna no básica se
expresa como
1
A∗3 = 2A∗1 + A∗2 .
2

10

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

b)      
1 2 3 1 2 3 1 2 3

 2 6 8 


 0 2 2 


 0 1 1 


 2 6 0 →


 0 2 −6 →


 0 2 −6 →

 1 2 5   0 0 2   0 0 2 
 3 8 6   0 2 −3   0 2 −3 
1 0 1 1 0 1 1 0 1

 0 1 1 


 0 1 1 


 0 1 1 


 0 0 −8 →

 0 0 1 →


 0 0 1 →

 0 0 2   0 0 2   0 0 2 
 0 0 −5 0 0 −5 0 0 −5
1 0 0

 0 1 0 


 0 0 1 

 0 0 0 
0 0 0
Todas las columnas son básicas.
c) Escribimos las transformaciones principales para llegar a la forma escalonada reducida por filas.
   
2 1 1 3 0 4 1 1 1/2 1/2 3/2 0 2 1/2
 4 2 4 4 1 5 5   4 2 4 4 1 5 5 
   
 2 1 3 1 0 4 3   2 1 3 1 0 4 3 
A=  →  
 6 3 4 8 1 9 5   6 3 4 8 1 9 5 
  
 0 0 3 −3 0 0 3   0 0 3 −3 0 0 3 
 8 4 2 14 1 13 3   8 4 2 14 1 13 3 
1 1/2 1/2 3/2 0 2 1/2 1 1/2 1/2 3/2 0 2 1/2
 0 0 2 −2 1 −3 3   0 0 1 −1 1/2 −3/2 3/2 
   
 0 0 2 −2 0 0 2   0 0 2 −2 0 0 2 
 →  →
 0 0 1 −1 1 −3 2   0 0 1 −1 1 −3 2 
   
 0 0 3 −3 0 0 3   0 0 3 −3 0 0 3 
 0 0 −2 2 1 −3 −1   0 0 −2 2 1 −3 −1 
1 1/2 0 2 −1/4 11/4 −1/4 1 1/2 0 2 −1/4 11/4 −1/4
 0 0 1 −1 1/2 −3/2 3/2   0 0 1 −1 1/2 −3/2 3/2 
   
 0 0 0 0 −1 3 −1   0 0 0 0 1 −3 1 
 →  →
 0 0 0 0 1/2 −3/2 1/2   0 0 0 0 1/2 −3/2 1/2 
   
 0 0 0 0 −3/2 9/2 −3/2   0 0 0 0 −3/2 9/2 −3/2 
 0 0 0 0 2 −6
 2 0 0 0 0 2 −6 2
1 1/2 0 2 0 2 0
 0 0 1 −1 0 0 1 
 
 0 0 0 0 1 −3 1 
 .
 0 0 0 0 0 0 0 
 
 0 0 0 0 0 0 0 
0 0 0 0 0 0 0
El rango vale 3, las columnas básicas son {A∗1 , A∗3 , A∗5 } y las expresiones de las columnas no básicas son
1
A∗2 = A∗1 ,
2
A∗4 = 2A∗1 − A∗3 ,
A∗6 = 2A∗1 − 3A∗5 ,
A∗7 = A∗3 + A∗5 .
Ejercicio 48 Pruebe que el siguiente sistema es compatible y calcule todas sus soluciones:
x1 + x2 + 2x3 − x4 = 2,
3x2 − x3 + 4x4 = 3,
x1 + 2x2 − 3x3 + 5x4 = 3,
x1 + x2 − 5x3 + 6x4 = 2.

11

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

De entre ellas, seleccione las que se encuentran sobre la esfera x21 + x22 + x23 + x24 = a, con a > 0.

Solución 48 Procedemos a calcular la forma escalonada reducida por filas de la matriz ampliada asociada al sistema.
 
1 0 0 0 1
 
rref  0 1 0 1 1 

[A|b] →  .

 0 0 1 −1 0 
 
0 0 0 0 0

Como no hay pivote en la última columna, el sistema es compatible. Las variables básicas son x1 , x2 , x3 y la variable
libre es x4 . Por tanto, el conjunto de soluciones del sistema es de la forma

x1 = 1
x2 = 1 − t,
x3 = t,
x4 = t,
t ∈ R.

Para calcular las que estén en la esfera, sustituimos en la ecuación cuadrática para expresarla en función del parámetro
t. Nos queda
1 + (1 − t)2 + t2 + t2 = a ⇒ (2 − a) − 2t + 3t2 = 0.
El discriminante es ∆ = 4(3a − 5), por lo que, para tener solución, debe ser a ≥ 53 . En tal caso, las soluciones son
√ √
t1 = 13 (1 + 3a − 5), t2 = 13 (1 − 3a − 5) y los puntos buscados son pi = (1, 1 − ti , ti ti )t , i = 1, 2.

Ejercicio 49 Si A es una matriz de orden m × n y rango(A) = m, explique por qué el sistema [A|b] es compatible
para cualquier b ∈ Km .

Solución 49 El rango de la matriz [A|b] es igual que el de A, puesto que no puede ser mayor que su número de filas,
que es m. Por tanto, el sistema es compatible: rango(A) = m = rango([A|b]), es decir, la última columna de [A|b] no
es básica.

Ejercicio 50 Consideremos dos sistemas compatibles, cuyas matrices ampliadas son de la forma [A|b] y [A|c], esto
es, tienen la misma matriz de coeficientes. ¿Es el sistema [A|b + c] compatible? Explique por qué.

Solución 50 Si los sistemas son compatibles, entonces la columna b es combinación lineal de las columnas de A y la
columna c también, esto es, podemos escribir
n
X n
X
b= βj A∗j , c = γj A∗j , para ciertos βj , γj .
j=1 j=1

Pero entonces b + c = nj=1 (βj + γj )A∗j , es decir, b + c se expresa como combinación lineal de las columnas de A, lo
P
que implica que el sistema [A|b + c] es compatible.

Ejercicio 51 Supongamos que [A|b] se reduce a una forma escalonada por filas [B|c]. ¿Es B una forma escalonada
por filas de A?

Solución 51 De forma independiente a que la columna c contenga un pivote, los restantes se encuentran en B y por
tanto B tiene forma escalonada por filas. Por ello, las transformaciones para llegar a [B|c] son las mismas que hay que
aplicar a A para producir una forma escalonada. Ası́, B es una forma escalonada obtenga a partir de A.

Ejercicio 60 Calcule la solución general del sistema


 x3 +x4 = 0,
−2x1 −4x2 +x3 = 0,
3x1 +6x2 −x3 +x4 = 0.


12

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Solución 60 Calculamos la forma escalonada reducida por filas de la matriz de coeficientes y a partir de dicha
expresión obtenemos todas las soluciones:
     
0 0 1 1 −2 −4 1 0 1 2 −1/2 0
     
 −2 −4 1 0  →  0 0 1 
1→  0 0 1 1 

   
 3 6 −1 1  3 6 −1 1  3 6 −1 1
1 2 −1/2 0 1 2 0 1/2 1 2 0 1/2
     
 0 0 1 1   0 0 1 1   0 0 1 1
→ →

   
 0 0 1/2  1 0 0 0 1/2 0 0 0 1
1 2 0 0
 
 0 0 1 0 .
 
0 0 0 1

El sistema queda en forma equivalente como


 x1 = −2x2 ,

 x1 +2x2 = 0,

x2 = x2 ,

x3 = 0, ⇒
x3 = 0,
x4 = 0.
 

x4 = 0.


Por tanto, el conjunto de soluciones es

   
−2x2 −2
 x2 
 = x2  1  .
 
u=
 0   0 
0 0

Ejercicio 61 Calcule b − Ay y b − Az para

       
0,217 0,780 0,563 0,341 0,999
b= ,A = ,y = ,z = .
0,254 0,913 0,659 −0,087 −1,0

¿Es mejor y que z como solución de Ax = b? Considere los siguientes puntos de vista:
1. Compare los residuos b − Ay y b − Az.
2. Compare y y z con la solución exacta x = (1, −1)t .

Solución 61 Lo que se pide es la comparación de los vectores de residuos para cada una de las soluciones calculadas
y y z. Es fácil ver que
" # " #
0,00000099999999997324452 0,000779999999999919648
ry = b − Ay = , rz = b − Az = .
−5,55111512312578272 × 10−17 0,000912999999999997148

Si comparamos estos residuos, vemos que las entradas de ry son bastante más pequeñas que las de rz , lo que nos
puede llevar a pensar que la solución y es mejor que la solución z.
La solución exacta del sistema es x = (1, −1)t y la distancia de y a x es mayor que la de z a x. ¿Por qué los
residuos nos dicen una cosa y las soluciones calculadas otra? Debido al mal condicionamiento de la matriz.

Ejercicio 62 Calcule el rango de la matriz

a2
 
1 1 2
A= 1 1−a 2 0 
2 2−a 6−a 4

en función del parámetro a.

13

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Solución 62 Efectuaremos las operaciones elementales por fila sobre la matriz A para llevarla a su forma escalonada
reducida por filas. Cuando un pivote dependa del valor del parámetro a, tendremos que distinguir casos para obtener
los diferentes valores del rango de A.

F2 − F1 a2
 
1 1 2
F3 − 2F1  
A →  0 −a
A1 =  0 −a2 .

0 −a 2 − a 4 − 2 a2

En la segunda columna, por debajo de la posición (2, 2), dependemos del valor de a para seleccionar un pivote. Si
a 6= 0, podemos escoger la posición (2, 2) o la (3, 2) como posición pivotal.

Caso 1. a 6= 0. Continuamos con el proceso:

− a1 F2
F1 − F2 a2 − a
 
1 0 2
F3 + aF1  
A1 →  0
A2 =  1 0 a .

0 0 2−a 4 − a2

En la posición (3, 3) tenemos ahora una entrada que se anula si 2 − a = 0. Procedemos entonces a considerar

• Caso 1.1. 2 − a 6= 0. Es inmediato que entonces tenemos 3 pivotes y rango(A) = 3.

• Caso 1.2. 2 − a = 0. Entonces

a2 − a
   
1 0 2 1 0 2 2
   
 0 1 0 a
A2 =  = 0 1 0 2 
  ,
2
0 0 2−a 4−a 0 0 0 0

y rango(A) = 2.

Caso 2. a = 0. La matriz queda como  

1 1 2 0
 
 0 0 0 0 
 ,
0 0 2 4
cuya forma reducida por filas es igual a  
1 1 0 −4
 
 0 0 1 2 
 .
0 0 0 0
Entonces, rango(A) = 2.

En resumen,
si a 6= 0, 2, entonces rango(A) = 3;
si a = 0, entonces rango(A) = 2;
si a = 2, entonces rango(A) = 2.

Ejercicio 63 Discuta el rango de la matriz  

a 1 1 2
 1 a 1 2 .
1 1 a 2
en función de los valores del parámetro a.

14

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Solución 63 Vamos a calcular una forma escalonada para determinar su rango. Partimos de la matriz dada, que
llamaremos C. Observamos que en la posición (1, 1) tenemos el parámetro a y por debajo hay elementos que no
dependen de él. Para evitar el caso a = 0, efectuamos un intercambio de filas para colocar en la posición pivotal un
elemento que podamos asegurar que no es nulo y realizamos las operaciones por filas para anular los términos por
debajo.

F12
 
1 a 1 2
F2 − aF1 , F3 − F1  2

C −−−−−−−−−−−−−−−→   0 1−a 1−a 2−2a 
.
0 1−a a−1 0

Encontramos ahora en la posición (2, 2) el término 1 − a2 , que se anula en a = ±1 y en la posición (2, 3) el elemento
1 − a, que se anula para a = 1. Para evitar, de nuevo, un caso, efectuamos la permutación entre la segunda y tercera
filas.
 
1 a 1 2
F23  
−→   0 1−a a−1 0  = C1 .

0 1 − a2 1 − a 2 − 2 a

Tenemos ahora que distinguir dos casos:

Caso 1. a 6= 1. En este caso podemos usar 1 − a como pivote y anular la entrada de la posición (2, 3).
 
1 a 1 2
F3 −(1+a)F2  
−−−−−−−→   0 1−a a−1 0 .

0 0 2 − a − a2 2 − 2 a

El elemento en la posición (3, 3) es igual a 2 − a − a2 = −(a + 2)(a − 1). Por tanto, dentro del caso en el que
estamos, solamente se puede anular para a = −2 (recordemos que estamos suponiendo que a 6= 1). Estudiamos
entonces las dos posibilidades:
• a 6= −2. En este caso, el elemento (3, 3) es no nulo y el rango de la matriz vale 3. No nos importa lo que
ocurra en la posición (3, 4).
• a = −2. El elemento (3, 3) se anula, pero con este valor 2 − 2a = 6 6= 0, por lo que tenemos un pivote en la
posición (3, 4). Ası́, el rango vale 3.
Caso 2. a = 1. Con este valor, la matriz C1 es
 
1 1 1 2
 0 0 0 0 .
0 0 0 0

Es una forma escalonada, con una sola fila no nula, por lo que el rango es igual a 1.

Ejercicio 64 Discuta el sistema

2x + 3y + 3z = 0,
4x + 8y + 12z = −4,
6x + 2y + αz = 4,
en función de los valores de α. Calcule las soluciones en función del parámetro.

Solución 64 Si efectuamos reducciones a la matriz ampliada,

 
2 0 −6 6
[A|b] →  0 2 6 −4  .
0 0 α + 12 −10

15

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Si α + 12 = 0, el sistema es incompatible. Si α + 12 6= 0, la solución del sistema es

α+2
x = 3 ,
α + 12
α−3
y = −2 ,
α + 12
1
z = −10 .
α + 12
Ejercicio 65 Determine el carácter incompatible, con solución única o con infinitas soluciones, de los siguientes
sistemas, en función de los valores de c:
 
  x1 +2x2 −x3 = c,  cx1 +x2 +x3 = 2,
x2 +cx3 = 0,
(a) (b) x1 +3x2 +x3 = 1, (c) x1 +cx2 x3 = 2,
x1 −cx2 = 1.
3x1 +7x2 −x3 = 4. x1 +x2 +cx3 = 2.
 

Solución 65 En todos los casos vamos a estudiar la forma escalonada reducida por filas de la matriz ampliada del
sistema. Dicha matriz dependerá del parámetro c y según su configuración de pivotes clasificaremos el sistema.
(a) El sistema a) tiene como matriz ampliada
 
0 1 c 0
.
1 −c 0 1
Lo primero que debemos hacer es intercambiar la primera y segunda fila para tener un elemento no nulo en la
posición (1, 1) y nos queda la matriz " #
1 −c 0 1
.
0 1 c 0
El pivote ya es igual a 1 y además no es necesario hacer ceros por debajo de esa posición, con lo que pasamos a la
segunda fila. Encontramos en la posición (2, 2) un pivote igual a 1. Con esto ya tenemos garantı́a que el sistema
es compatible, pues no hay pivote en la columna correspondiente al término independiente y las variables básicas
son x1 y x2 , mientras que la libre es x3 . En consecuencia, el sistema es compatible indeterminado. Por si hay
alguna duda, la forma escalonada reducida por filas es
1 0 c2 1
" #
.
0 1 c 0
Observemos que no hay que distinguir el caso c = 0, pues no afecta a los pivotes ni a la variable libre.
(b) El sistema b) tiene como matriz ampliada
 
1 2 −1 c
 
 1 3
[A|b] =  1 
1.
3 7 −1 4
Procedemos a calcular su forma escalonada reducida por filas.

F2 − F1
 
1 2 −1 c
F3 − 3F1  
[A|b] −−−−−−−−→  0 1 2 1−c 
 
0 1 2 4−3c
F3 − F2
 
1 0 −5 3 c − 2
F1 − 2F2  
 0 1
−−−−−−−−→  2 1−c .
0 0 0 3−2c
Si 3 − 2c = 0, es decir, c = 3/2, no hay pivote en la columna correspondiente al término independiente. El sistema
es compatible, con variable libre x3 . Entonces es compatible indeterminado.
Si c 6= 3/2, hay un pivote en la última columna y el sistema es incompatible.

16

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

(c) La matriz ampliada del sistema c) es  

c 1 1 2
 
 1 c 1 2 
 .
1 1 c 2
Como existen elementos no nulos en la primera columna que no dependen del valor de c, efectuamos un inter-
cambio de las filas 1 y 2 para evitar el caso c = 0.
 
1 c 1 2
F12 
 
−→ c 1 1 2 
.
1 1 c 2
Ahora anulamos los elementos por debajo del pivote:
F2 − cF1
 
1 c 1 2
F3 − F1 
−−−−−−−→  0 1 − c2

 1−c 2−2c 
.
0 1−c c−1 0

Como pivote en la posición (2, 2) tenemos el elemento 1 − c2 , que se anula para c = 1 o c = −1. Por debajo
encontramos otro elemento que depende de c, pero que solamente se anula para c = 1. Efectuamos entonces un
intercambio entre la segunda y tercera filas para simplificar los casos.
 
1 c 1 2
F23  
−
→  0 1 − c −1 + c 0 . (1)

2
0 1−c 1−c 2−2c
Caso 1. c 6= 1. Dividimos la segunda fila por 1 − c y anulamos las posiciones por debajo y encima del pivote.
 
1 c 1 2
1
1−c F2
 
−−−→  0 1 −1+c
0 
 1−c 
2
0 1−c 1−c 2−2c
2
F3 − (1 − c )F2
 
1 0 c+1 2
F1 − cF2  
−−−−−−−−−−−−−→   0 1 −1 0 .

0 0 −c2 − c + 2 2 − 2 c

En la posición (3, 3) encontramos el elemento −c2 − c + 2, que se anula para c = 1 y para c = −2. La primera
posibilidad no puede ser, pues estamos suponiendo que c 6= 1. Entonces
Si c = −2, entonces la última matriz queda de la forma
 
1 0 −1 2
 0 1 −1 0  ,
0 0 0 6
por lo que existe un pivote en la columna del término independiente. En este caso, el sistema es incompatible.
Si c 6= −2, la posición (3, 3) contiene un pivote y el sistema es compatible determinado.
Caso 2. c = 1. La matriz que aparece en (1) toma la forma
 
1 1 1 2
 0 0 0 0 ,
0 0 0 0
que ya está en forma escalonada. Hay una variable básica (x1 ) y dos variables libres (x2 , x3 ). El sistema es
compatible indeterminado.
En resumen, podemos decir que

17

Descargado por Dario Calvache ([email protected])

 lOMoARcPSD|3870267

Si c = 1, el sistema es compatible indeterminado.

Si c = −2, el sistema es incompatible.
Si c 6= 1, −2, el sistema es compatible determinado.

18

Descargado por Dario Calvache ([email protected])