RELACIÓN DE EJERCICIOS DEL TEMA 1: MATRICES.

Ejercicio 1:

Dada la matriz , calcula una matriz tal que es la forma normal de

Hermite (forma escalonada reducida por filas) de A. Calcula también el rango de A.

Consideramos la matriz A:

A=[0 0 -1 2 -1; 0 0 0 2 -3 ; 0 1 2 -1 0]

A = 3×5
0 0 -1 2 -1
0 0 0 2 -3
0 1 2 -1 0

Por defecto, Matlab considera variables numéricas y esto puede implicar errores de aproximación. Para que
considere la matriz anterior como una matriz simbólica, tenémos que decírselo expresamente

A=sym(A)

A =

Consideramos también la matriz identidad :

I3=eye(3)

I3 = 3×3
1 0 0
0 1 0
0 0 1

Concatenamos ambas matrices. Observa como, al operar una variable simbólica con una numérica, Matlab
considera que el resultado es simbólico.

B=[A I3]

B =

Y lo transformamos en una matriz escalonada reducida por filas:

1
 B=Hermite(B)

B =

La matriz Q que pide el ejercicio es la parte de la matriz obtenida que se corresponde con donde inicialmente
estaba la identidad.

Q=B(:,6:8)

Q =

Y la matriz de Hermite de A es:

H=B(:,1:5)

H =

Podemos comprobar que Q A = H

Q*A

ans =

El rango de la matriz es el número de filas no nulas en H, es decir, 3. Matlab tiene predefinida una función que
calcula el rango de una matriz:

rank(A)

ans = 3

2
Ejercicio 2:
Determina el rango de las siguientes matrices según el valor del parámetro a, utilizando transformaciones
elementales:

Apartado a

Para empezar definimos a como un parámetro simbólico.

syms a

Ahora definimos la primera matriz. Como una de las entradas es un símbolo, Matlab asume que la matriz es
simbólica.

A=[1 3; 4 a]

A =

Si la transformamos en escalonada reducida por filas obtenemos

Hermite(A)

ans =

Concluiríamos que la matriz tiene rango 2.

OjO: Este resultado no es correcto. El problema radica en que, al tratar con expresiones simbólicas, donde
pueden intervenir parámetros, puede que Matlab haga divisones (por cero) no permitidas.

Vamos a resolver el ejercicio paso a paso:

Primero sumamos a la fila 2 la fila 1 multiplicada por (-4)

TransElemFilaTipoIII(A,2,1,-4)

ans =

Concluimos que A tiene rango 1 si y rango 2 en otro caso.

Volvamos al error de Matlab. Si le pedimos que nos calcule el rango de A, nos dice que es 2
independientemente del valor de a.

3
 rank(A)

ans = 2

Sin embargo, si sustituimos el valor de a por y le volvemos a calcular el rango, nos dice que es 1!!!!!

A=subs(A,a,12)

A =

rank(A)

ans = 1

Importante: Por muy bueno que sea al software que utilicemos siempre es posible que obtangamos resultados
erróneos/incompletos/imprecisos. Tenemos que ser críticos con los resultados que obtengamos, y, para ello, es
necesario conocer los mecanismos utilizados.

Apartado b

Como no interviene el parámetro a, podemos hacerlo directamente.

B=sym([0 1; 1 0])

B =

Hermite(B)

ans =

Esta matriz tiene rango 2.

Apartado c

Definimos la matriz:

syms a
C=[2 a; 6 4+a; 4 6]

C =

A la fila 2 le restamos la fila 1 multiplicada por 3 ( ), y a la fila 3 le restamos la fila 1 multiplicada

por 2 ( ).

4
 C1=TransElemFilaTipoIII(C,2,1,-3);
C2=TransElemFilaTipoIII(C1,3,1,-2)

C2 =

Dividimos la fila 2 y la 3 por 2: y

C3=TransElemFilaTipoII(C2,2,1/2);
C4=TransElemFilaTipoII(C3,3,1/2)

C4 =

A la vista de esta matriz concluimos que, independientemente del valor de a, el rango siempre es 2.

Ejercico 3:
Calcula la forma normal de Hermite por filas y el rango de la siguiente matriz según los valores del parámetro a:

Apartado a

Comenzamos definiendo la matriz:

syms a
A=[1 1 -1;3 a a; 4 a 0]

A =

Restamos, a la fila 2, la 1 multiplicada por 3 y, a la fila 3, la 1 multiplicada por 4.

A=TransElemFilaTipoIII(A,2,1,-3);
A=TransElemFilaTipoIII(A,3,1,-4)

A =

5
Ahora distinguimos dos casos:

CASO 1: si tenemos que la matriz de Hermite es

Hermite(sym([1 1 -1;3 3 3; 4 3 0]))

ans =

En este caso la matriz tiene rango 3.

CASO 2: si , podemos dividir la fila 1 por .

A=TransElemFilaTipoII(A,2,1/(a-3))

A =

Ahora podemos restar a la fila 3 la 2 multiplicada por (a-4)

A=TransElemFilaTipoIII(A,3,2,-(a-4))

A =

Igualamos a cero el elemento de la posición (3,3) y resolvamos la ecuación:

solve(A(3,3)==0, a)

ans =

Concluimos que el rango de A es 2 si o . En otro caso, el rango es 3.

Comprobemos que, efectivamente, si o , el rango es 2:

subs(A,a,0)

ans =

6
 subs(A,a,5)

ans =

Apartado b

Antes de definir la matriz, volvemos a declarar la variable a como simbólica.

syms a
B=sym([2 a 1;6 3 4])

B =

Dividimos la fila 1 por 2.

B1=TransElemFilaTipoII(B,1,1/2)

B1 =

Ahora restamos a la fila 2 la fila 1 multiplicada por 6. Esto se puede hacer independientemente del valor de a.

B2=TransElemFilaTipoIII(B1,2,1,-6)

B2 =

Concluimos que la matriz tiene rango 2 independientemente del valor de a.

Ejercicio 4
Calcula la inversa de las siguientes matrices, cuando sea posible, utilizando transformaciones elementales:

7
Apartado a

Empezamos introduciendo la matriz del apartado a.

A=sym([0 1;1 0])

A =

Le realizamos operaciones elementales a A y, las mismas operaciones, a la identidad.

Id=sym(eye(2))

Id =

A1=TransElemFilaTipoI(A,1,2)

A1 =

Q=TransElemFilaTipoI(Id,1,2)

Q =

Como A ya se ha transformado en la identidad, la inversa de A es la que hemos llamado Q.

Podemos comprobar que ambos productos, AQ y QA, son la identidad.

Q*A

ans =

A*Q

ans =

Para simplificar el procedimiento, se suelen concatenar la matriz A y la identidad. De esa forma hacemos las
operaciones elementales fila a ambas matrices simultanemente. Vamos a hacerlo así en el apartado b.

Apartado b

Definimos la matriz B como matriz simbólica:

8
B=sym([1 0 -1; 0 2 3; 1 -1 2])

B =

B1=[B eye(3)]

B1 =

B1=TransElemFilaTipoIII(B1,3,1,-1)

B1 =

B1=TransElemFilaTipoII(B1,3,-1)

B1 =

B1=TransElemFilaTipoIII(B1,2,3,-2)

B1 =

B1=TransElemFilaTipoI(B1,2,3)

B1 =

B1=TransElemFilaTipoII(B1,3,1/9)

B1 =

9
 B1=TransElemFilaTipoIII(B1,1,3,1)

B1 =

B1=TransElemFilaTipoIII(B1,2,3,3)

B1 =

Por tanto, la inversa es:

B1(:,4:6)

ans =

Matlab tiene implementadas las principales funciones para matrices. Para calcular la inversa de una matriz
podemo usar la función inv() o elevar la matriz a -1.

inv(B)

ans =

B^(-1)

ans =

10
Apartado c

Lo vamos a hacer ya directamente usando la función Hermite que hemos programado.

C=sym([1 2 3 4; 0 2 3 4; 0 0 3 4; 0 0 0 4])

C =

C1=[C eye(4)]

C1 =

C1=Hermite(C1)

C1 =

Matlab trae predefinda la función rref que hace los mismos cálculos que la funcion Hermite que hemos
programado nosostros. Veámoslo:

rref(C1)

ans =

11
La inversa es:

C1(:,5:8)

ans =

Apartado d

D=sym([1 3 -4;1 5 -1;3 13 -6])

D =

D1=[D eye(3)]

D1 =

D1=Hermite(D1)

D1 =

OjO: la matriz D no tiene inversa porque tiene rango 2. Veamos qué devuelve Matlab si le pedimos que calcule
la inversa de D con los dos métodos mencionados anteriormente:

D^(-1)

12
 ans =

inv(D)

ans =

Ejercicio 5
Para una matriz diagonal de orden n, deduce:

1. la expresión de su potencia n-sima;

2. el valor de su determinante;
3. la condición necesaria y suficiente para que sea inversible; en tal caso, la expresión de su inversa.

Esta pregunta es teórica. Analizando el algortimo de la multiplicación de matrices, concluimos que

Finalmente, si el determinate no es cero, la inversa es

Ponemos algunos ejemplos:

syms a b c d
D=sym(diag([ a b c d]))

D =

13
La quinta potencia de D es:

D^5

ans =

Su determinante es:

DeterminantePorGauss(D)

ans =

Y su inversa es:

D^(-1)

ans =

Ejercicio 6
Discute y resuelve, cuando sea posible, utilizando el método de Gauss y el de Gauss-Jordan, los siguientes
sistemas de ecuaciones:

Definimos la matriz de coeficientes y la matriz ampliada:

A=sym([3 2 1; 5 3 3; 1 1 -1]);
B=[1;2;1];
Aamp=[A B]

Aamp =

14
Aplicando Gauss obtenemos la matriz

A2=TransElemFilaTipoI(Aamp,1,3)

A2 =

A2=TransElemFilaTipoIII(A2,2,1,-5);
A2=TransElemFilaTipoIII(A2,3,1,-3)

A2 =

A2=TransElemFilaTipoI(A2,2,3);
A2=TransElemFilaTipoII(A2,2,-1)

A2 =

A2=TransElemFilaTipoIII(A2,3,2,2)

A2 =

Y esto quiere decir que el sistema es incompatible porque la última ecuación se ha transformado en .

Ya no tiene sentido aplicar Gauss-Jordan.

Apartado b

Comenzamos definiendo las matrices implicadas:

A=[1 -3 -2 4; 3 -8 -3 8; 2 -3 5 -4];
B=sym([5;18;19]);
Aamp=[A B]

Aamp =

15
 A2=TransElemFilaTipoIII(Aamp,2,1,-3);
A2=TransElemFilaTipoIII(A2,3,1,-2)

A2 =

A2=TransElemFilaTipoIII(A2,3,2,-3)

A2 =

El método de Gauss transforma ahora esta matriz en ecuaciones para despejar los valores desde abajo hacia
arriba.

Viendo el rango de las matrices ya sabemos que el sistema es compatible indeterminado.

Además, el sistema tiene dos grados de libertad, es decir, el conjunto de las soluciones se puede expresar
usando dos parámetros.

Lo más sencillo es asignar parámetros a las incógnitas cuyas columnas no tinen pivote y despejar las otras:

Asignamos y .

De obtenemos que

y de obtenemos que .

Terminamos ahora con los últimos pasos del método de Gauss-Jordan:

A2=TransElemFilaTipoIII(A2,1,2,3)

A2 =

Y, de aquí, obtenemos directamente la solución:

; ; y .

Apartado c

16
En este, y los siguientes apartados, lo vamos a hacer directamente. Definimos la matriz ampliada

Aamp=sym([1 1 1 -1 -1; 3 -2 5 0 10; 0 3 0 5 7; 1 0 1 -1 0])

Aamp =

EscalonaPorFilas(Aamp)

ans =

Observamos que el sistema es compatible determinado. Despejamos las incógnitas de abajo hacia arriba:

Si utilizamos el método de Gauss-Jordan, obtenemos la solución directamente:

Hermite(Aamp)

ans =

Apartado d

17
Este sistema es lineal homogeneo y, por tanto, compatible y, dado que hay más incógnitas que ecuaciones,
será indeterminado.

Aamp=sym([1 -1 1 -1 1 0; 1 1 2 2 -1 0;-1 5 1 7 -5 0])

Aamp =

EscalonaPorFilas(Aamp)

ans =

Efectivamente, el sistema es compatible indeterminado y tiene 3 grados de libertad. Despejamos de abajo

hacia arriba.

, ,

para todo

Hermite(Aamp)

ans =

El método de Gauss-Jordan nos da la solución directamente:

, ,

18
 ,

para todo

Apartado e

Aamp=sym([1 1 -2 4 5; 2 2 -3 1 3; 3 3 -4 -2 1])

Aamp =

EscalonaPorFilas(Aamp)

ans =

Es un sistema compatible indeterminado con dos grados de libertad. Las incognitas que podemos despejar
directamente son

para todo

Lo resolvemos por Gauss-Jordan:

Hermite(Aamp)

ans =

Y obtenemos

19
para todo

Apartado f

Directamente podemos decir que es un sistema comopetible indeterminado. Lo resolvemos por Gauss-Jordan:

Aamp=sym([1 2 -3 1 0;1 -3 1 -2 0; 2 1 -3 5 0])

Aamp =

Hermite(Aamp)

ans =

Tiene un grado de libertad y la solución es

para todo

Apartado g

Lo resolvemos primero por Gauss

Aamp=sym([1 1 1 6; 2 -1 2 3; 3 2 -3 3])

Aamp =

20
 EscalonaPorFilas(Aamp)

ans =

Es un sistema compatible determinado y la solución es:

Si lo resolvemos por Gauss-Jordan obtenemos la solución directamente:

Hermite(Aamp)

ans =

Solución: , , .

Ejercicio 7
Discute los siguientes sistemas de ecuaciones lineales según los valores de los parámetros a y b:

Apartado a

Declaramos a como simbólica y definimos la matriz ampliada:

syms a
Aamp=[(1+a) 1 1 1; 1 (1+a) 1 a; 1 1 (1+a) a^2]

Aamp =

Como ya hemos docho anteriormente, tenemos que aplicar Gauss paso a paso porque Matlab, trabajando en
modo simbólico, podría hacer divisiones por cero. Comenzamos intercambiando la fila 1 con la 3:

A2=TransElemFilaTipoI(Aamp,1,3)

21
 A2 =

Ahora restamos a la fila 2 la fila 1 y a la fila 3 la 1 multiplicada por . Estas operaciones se pueden hacer
indpendientemente del valor de a.

A2=TransElemFilaTipoIII(A2,2,1,-1);
A2=TransElemFilaTipoIII(A2,3,1,-(a+1))

A2 =

Sumamos ahora a la fila 3 la fila 2:

A2=TransElemFilaTipoIII(A2,3,2,1)

A2 =

Para ver mejor la matriz, vamos a pedir a Matlab que simplifique las expresiones:

simplify(A2)

ans =

Observamos que, si y , el sistema es compatible determinado. Veamos que sucede cuando a es 0

o -3.

Empezamos sustituyendo a por 0:

A3=subs(A2,a,0)

A3 =

En este caso nos sale un sistema incompatible. Veamos qué sucede si a es -3:

A3=subs(A2,a,-3)

22
 A3 =

También nos sale que es un sistema incompatible.

Solución: Si el sistema es incompatible. En otro caso, el sistema es comatible determinado.

Apartado b

Definimos la matriz ampliada:

syms a
Aamp=[1 2 -1 a;1 -a 1 2]

Aamp =

A la fila 2 le restamos la fila 1:

A2=TransElemFilaTipoIII(Aamp,2,1,-1)

A2 =

Y observamos que, independientemente del valor de a, la matriz de coeficientes tien rango 2 y, por tanto, la
matriz ampliada también.

Solución: El sistema es compatible indeterminado para todo .

Apartado c

Declaramos las variables simbólicas a y b, y definimos la matriz ampliada:

syms a b
Aamp=[a*(1-b) (b-1) 0 1-b; a*b 1 1 b; a b 1 1]

Aamp =

23
Intercambiamos las filas 1 y 3:

A2=TransElemFilaTipoI(Aamp,1,3)

A2 =

A la fila 2 la restamos al 1 multiplicada por b y a la 3 le sumamos la 1 multiplicada por :

A2=TransElemFilaTipoIII(A2,2,1,-b);
A2=TransElemFilaTipoIII(A2,3,1,b-1)

A2 =

Simplificamos las expresiones:

A2=simplify(A2)

A2 =

A la fila 3 le sumamos la fila 2. Observa que las operaciones que estamos realizando se pueden hacer
independientemente del valor de a y de b.

A2=TransElemFilaTipoIII(A2,3,2,1)

A2 =

Está claro que, si y , el sistema es compatible indeterminado. Determinado nunca podría

ser porque el rango de la matriz de coeficientes es siempre menor que 3. Tampoco puede ser incompatible
porque el rango de la matriz de coeficientes no puede ser mayor que el da la matriz ampliada.

Solución: el sistema es compatible indetermiado para todo .

Comprobemos que es siempre compatible indeterminado en alguno de los casos críticos:

A3=subs(A2,b,-1)

A3 =

24
 A3=subs(A2,b,1)

A3 =

A3=subs(A2,b,0)

A3 =

Podemos ir un poco más allá y decir que el sistema es compatible indeterminado con un grado de libertad
excepto cuando que tendría dos grados de libertad.

Apartado d

Comenzamos definiendo la matriz ampliada:

syms a b
Aamp=[a 1 1 1 1;1 a 1 1 b; 1 1 a 1 b^2; 1 1 1 a b^3]

Aamp =

A2=TransElemFilaTipoI(Aamp,1,4)

A2 =

A2=TransElemFilaTipoIII(A2,2,1,-1);
A2=TransElemFilaTipoIII(A2,3,1,-1);

25
 A2=TransElemFilaTipoIII(A2,4,1,-a)

A2 =

A2=TransElemFilaTipoIII(A2,4,2,1)

A2 =

A2=TransElemFilaTipoIII(A2,4,3,1)

A2 =

solve(A2(4,4)==0,a)

ans =

A3=subs(A2,a,1)

A3 =

factor(A3(4,5))

ans =

Primera conclusión: si y entonces el sistema es compatible indeterminado con tres grados de

libertad.

subs(A3,b,1)

ans =

26
Si y , el sistema es incompatible. Por ejemplo:

subs(A3,b,0)

ans =

Veamos qué sucede si

A3=subs(A2,a,-3)

A3 =

assume(b,'real')
solve(A3(4,5)==0,b)

ans =

Tenemos que, si y , el sistema es compatible indeterminado con un grado de libertad; mientras

que, si y , el sistema es incompatible.

Finalmente, si , el sistema es compatible determinado.

Solución:

• Si , el sistema es compatible determinado.

• Si y entonces el sistema es compatible indeterminado con tres grados de libertad.
• Si y , el sistema es incompatible.
• Si y , el sistema es compatible indeterminado con un grado de libertad
• Si y , el sistema es incompatible.

Ejercicio 8
Usando Gauss, clasifica el siguiente sistema en función de

27
Declaramos k como variable simbólica y definimos la matriz ampliada:

syms k
Aamp=[k 1 1 1;1 k 1 1; 1 1 k 1]

Aamp =

A2=TransElemFilaTipoI(Aamp,1,3)

A2 =

A2=TransElemFilaTipoIII(A2,2,1,-1);
A2=TransElemFilaTipoIII(A2,3,1,-k)

A2 =

A2=TransElemFilaTipoIII(A2,3,2,1)

A2 =

assume(k,'real')
solve(A2(3,3)==0,k)

ans =

Distinguimos entonces tres casos. Si entonces el sistema es compatible indeterminado con dos grados de
libertad.

subs(A2,k,1)

ans =

28
Si entonces el sistema es incompatible.

subs(A2,k,-2)

ans =

En otro caso, el sistema es compatible determinado.

Ejercicio 9
Determina si hay valores de a y b que hagan que la siguiente ecuación matricial tenga solución,
siendo

Para despejar la X, vamos a pasar todos los términos que la tengan a la parte derecha de la ecuación.
. Ahora necesitamos incluir a la matriz identidad, , para poder sacar factor común
la X. De esta forma, la eqcuación queda como sigue: . Definimos estas matrices, no sin antes
declarar las variables simbólicas:

syms a b
A=[0 0 2;0 a 0;-2 0 0]

A =

A2=A-2*eye(3)

A2 =

B=[1 b 1; 3 5 a; -2 2 7]

B =

29
 A3=[A2 B]

A3 =

Realizamos operaciones elementales fila y analizamos las situaciones posibles.

A3=TransElemFilaTipoIII(A3,3,1,-1)

A3 =

Si el sistema es incompatible y, si , el sistema tiene solución porque la matriz es regular. La

solución será

X=inv(A2)*B

X =

Esta opreción viene predifinida en Matlab. En concreto, para matrices, Matlab define dos divisiones que son
y . La que nos interesa es la siguiente:

X=A2\B

X =

Ejercicio 10

30
Calcula los siguientes determinantes

Apartado a

Existen distintos métodos para calcular el determinante de una matriz. De todos ellos, como ya se ha dicho en
clase, el más eficente es el que utiliza transformaciones elementales fila para escalonar la matriz.

OjO: hay que tener en cuanta que las transfoamciones tipo I cambian el signo al determinante y que las tipo II
multiplican el determinante por un factor.

Aunque ya hemos implementado una función en Matlab para calcular el determinante

(DeterminatePorGauss) y el propio Matlab lo trae implementado (det), lo vamos a hacer paso a paso:

A=sym([1 -1 2 3; 4 0 2 5; -1 2 3 7; 5 1 4 0])

A =

A1=TransElemFilaTipoIII(A,2,1,-4);
A2=TransElemFilaTipoIII(A1,3,1,1);
A3=TransElemFilaTipoIII(A2,4,1,-5)

A3 =

Como solo hemos transformaciones elementales tipo III, el determinante no ha cambiado, es decir,

A4=TransElemFilaTipoI(A3,2,3)

A4 =

En este caso sí ha cambiado el determinante. Ahora tenemos que

31
 A5=TransElemFilaTipoIII(A4,3,2,-4);
A6=TransElemFilaTipoIII(A5,4,2,-6)

A6 =

Y tenemos que . Para que resulten más sencillos los siguientes pasos, vamos a
dividir la fila 3 por . Con ello, estamos dividiendo también por el determinante y tendremos que tenerlo
en cuenta.

A7=TransElemFilaTipoII(A6,3,-1/26)

A7 =

Ahora tenemos que

A8=TransElemFilaTipoIII(A7,4,3,36)

A8 =

Y, como esta matriz ya es escalonada por filas, su determinante es el producto de los elementos de la

diagonal. Por tanto,

Vamos a comprobar que es correcto:

det(A)

ans =

Apartado b

Definimos las matrices y lo hacemos de la misma forma:

A=sym([1 -1 2; 3 4 2;-2 3 4])

32
 A =

A1=TransElemFilaTipoIII(A,2,1,-3);
A2=TransElemFilaTipoIII(A1,3,1,2)

A2 =

Tenemos que

A3=TransElemFilaTipoI(A2,2,3)

A3 =

Ahora tenemos que

A4=TransElemFilaTipoIII(A3,3,2,-7)

A4 =

Concluimos que . Lo comprobamos:

det(A)

ans =

Apartado c

En este caso tenemos que declarar antes las variables a y b como simbólicas:

syms a b
A=[a^2 a*b b^2; 2*a a+b 2*b; 1 1 1]

A =

33
 A1=TransElemFilaTipoI(A,1,3)

A1 =

Con lo que

A2=TransElemFilaTipoIII(A1,2,1,-2*a);
A3=TransElemFilaTipoIII(A2,3,1,-a^2)

A3 =

Ahora tenemos que

Como

simplify(A3(3,2))

ans =

A4=TransElemFilaTipoIII(A3,3,2,-a)

A4 =

Matlab no ha simplificado las expresiones que hay en la matriz. Vamos a forzarle a que lo haga:

A4=simplify(A4)

A4 =

Por tanto, el determinannte es y esto es cierto para todo

. Como en los casos anteriores, vamos a comprobarlo.

det(A)

ans =

simplify(det(A))

34
 ans =

Apartado d

Para empezar, vamos a ver que sucede con una matriz que siga este esquema:

A=sym([0 1 1 1; 1 0 1 1; 1 1 0 1; 1 1 1 0])

A =

A1=TransElemFilaTipoI(A,1,4)

A1 =

Con esto tenemos que

A2=TransElemFilaTipoIII(A1,2,1,-1);
A3=TransElemFilaTipoIII(A2,3,1,-1)

A3 =

Como las transformaciones tipo III no cambian el determinante,

A4=TransElemFilaTipoIII(A3,4,2,1)

A4 =

A5=TransElemFilaTipoIII(A4,4,3,1)

A5 =

35
Por tanto, . Podemos repetir estos mismos pasos en una matriz

A partir de aquí solo vamos a usar trasnformaciones elementales tipo III y, por tanto, el determinate no cambia.
Empezamos por el cambio para todo .

Y ahora, para todo , realizamos la siguiente operación:

Ejercicio 11
Calcula el determinante

a) Desarrollando por la cuarta fila.

b) Desarrollando por la fila o columna para la que sea necesario calcular menos adjuntos.

c) Desarrollando por la segunda columna realizando antes operaciones elementales de forma que solamente
sea necesario calcular un adjunto.

36
Apartado a

y el resultado es

-5*det([6 2 0; 3 -1 4; 0 1 2])-3*det([1 2 0; 0 -1 4;2 1 2])-4*det([1 6 0;0 3 4; 2 0 2])

ans = -6

Apartado b

La columna cuatro es la que mas ceros tiene y la vamos a usar para calcular el determinante:

y el resultado es

4*det([1 6 2; 2 0 1;5 -3 4])-2*det([1 6 2;0 3 -1;5 -3 4])

ans = -6.0000

Apartado c

y el resultado también es

3*det([1 4 -8; 2 1 2; 5 3 4])

ans = -6.0000

Ejercicio 12
Se dice que una matriz cuadrada A es nilpotente si existe tal que . Calcula el determinante de
una matriz nilpotente.

Como para todo par de matrices cuadradas, tenemos que implica que y,
por tanto, porque es un cuerpo y en los cuerpos no hay divisores de cero.

Ejercicio 13

37
Se dice que una matriz es idempotente si . Calcula los posibles valores de su determinante.

Razonando igual que antes, si , debe cumplir . Resolvemos esta ecuación: equivale a
y, por no existir divisores de cero, esto solo es posible si o . Solución: el determinate
debe ser 0 o 1.

Ejercicio 14
Sabiendo que A y B son matrices cuadradas del mismo orden tales que y , calcula:

38