UNIVERSIDAD AUTÓNOMA DE NUEVO LEÓN
FACULTAD DE INGENIERÍA MECÁNICA Y ELÉCTRICA
LABORATORIO DE CONTROL MODERNO
PRÁCTICA 2
Valores propios, vectores propios y matriz inversa
Objetivo
Calcular los valores y vectores propios, así como la matriz inversa de una
matriz.
Valores propios
Los valores propios o eigenvalores de una matriz se obtienen a partir del
polinomio característico det(sI-A)=0 donde s son los valores propios de la
matriz A. Una forma de obtener los valores propios es programar la ecuación
previa esto es.
Ejemplo: Obtenga los valores propios de la siguiente matriz
1 −3 1
𝐴 = −2 −2 1
0 1 2
s = %s;
I = eye(3,3);
A = [1 -3 1; -2 -2 1; 0 1 2];
pc = det(s*I-A);
pc =
19 -11s -1s² +1s³
valp = roots(pc)
con esto resulta que los valores propios de A son:
valp =
-3.5707506 + 0.i
2.2853753 + 0.3131613i
2.2853753 - 0.3131613i
Una forma alternativa de calcular el polinomio característico de una matriz se
utiliza el comando poly cuya sintaxis es:
p = poly(vec, “vname”)
donde vec es un vector o matriz cuadrada de quien se requiere calcular el
polinomio y vname es la variable simbólica que aparecerá como variable
independiente en el polinomio.
pc = poly(A,"x")
pc =
19 -11x -1x² +x³
En forma directa los valores propios se pueden calcular utilizando el comando
el comando spec cuya sintaxis es:
evals=spec(A)
donde A es una matriz cuadrada con elementos reales o complejos.
Ejemplo: Obtenga los valores propios de la siguiente matriz
1 −3 1
𝐴 = −2 −2 1
0 1 2
A =[1 -3 1;-2 -2 1; 0 1 2];
valp = spec(A)
valp =
-3.5707506 + 0.i
2.2853753 + 0.3131613i
2.2853753 - 0.3131613i
Valores y vectores propios
El comando spec permite calcular los valores y vectores propios de una
matriz y su sintaxis es
[vecp valp] = spec(A)
donde A es una matriz cuadrada con elementos reales o complejos, vecp es el
nombre de la matriz que contiene los vectores propios calculados y valp es el
nombre de la matriz que contiene los valores propios calculados.
Ejemplo: Obtenga los valores y vectores propios de la siguiente matriz
1 −3 1
𝐴 = −2 −2 1
0 1 2
A = [1 -3 1; -2 -2 1; 0 1 2];
[vecp valp] = spec([A])
vecp =
0.564946 + 0.i 0.0479473 + 0.5496183i 0.0479473 - 0.5496183i
0.8121465 + 0.i -0.2191553 - 0.2404937i -0.2191553 + 0.2404937i
-0.1457876 + 0.i -0.7679548 + 0.i -0.7679548 + 0.i
valp =
-3.5707506 + 0.i 0. + 0.i 0. + 0.i
0. + 0.i 2.2853753 + 0.3131613i 0. + 0.i
0. + 0.i 0. + 0.i 2.2853753 - 0.3131613i
Matriz inversa
El comando inv permite calcular la inversa de una matriz cuadrada y su
sintaxis es:
Ai = inv(A)
donde A es una matriz cuadrada.
Ejemplo. Obtenga la matriz inversa de
1 −3 1
𝐴 = −2 −2 1
0 1 2
Ai = inv(A)
Ai =
0.2631579 -0.3684211 0.0526316
-0.2105263 -0.1052632 0.1578947
0.1052632 0.0526316 0.4210526
Reporte
Resuelva los siguientes ejercicios.
1.- Para las siguientes matrices obtenga el polinomio característico y sus valores
propios. Utilice las dos formas revisadas en la práctica para obtener los valores
propios.
3 0 −2 5 −1 2
𝐴= 1 1 4 𝐵 = −2 2 3
2 1 −3 1 3 1
2.- Para las matrices del ejercicio 1 obtenga los valores y vectores propios.
Verifique que los valores propios correspondan a los obtenidos en el ejercicio
1.
3.- Obtenga la matriz inversa de las matrices del ejercicio 1.
4.- Defina sus conclusiones.
