En esta sección, nuestro sistema lineal será una matriz de n×n de números complejos. Algunos conceptos de este modulo están basado en los conceptos básicos de álgebra lineal.
Sea
A
A una matriz de n×n donde
A
A es un operador lineal en los vectores de
Cn
n
.
Ax=b
A
x
b
(1)
donde
x
x y
b
b
son vectores de n×1 (
figura 1).
- Definition 1: Eigenvector
Un eigenvector de
A
A es un vector
v∈Cn
v
n
tal que
Av=λv
A
v
λ
v
(2)
donde
λ
λ es llamado el
eigenvalor correspondiente.
A
A solo cambia la
longitud de
v
v, no su dirección.
A través de las siguientes figura 2 y figura 3,
veamos las diferencias de la ecuación 1 y de la ecuación 2.
Si
v
v es un eigenvector de
A
A, entonces solo su longitud cambia.
Véase figura 3 y note que la longitud de nuestro vector esta simplemente escalada por una variable
λλ, llamada eigenvalor:
Cuando tratamos con una matriz
A A, los eigenvectores son los vectores posibles más
simples para trabajar.
Por inspección y entendimiento de eigenvectores, encuentre los dos eigenvectores
v
1
v
1
y
v
2
v
2
,
de
A=(
30
0-1
)
A
3
0
0
-1
También ¿cuáles son los eigenvalores correspondientes,
λ
1
λ
1
y
λ
2
λ
2
? No se preocupe si tiene problemas viendo estos valores de la información dada hasta ahora, veremos otras maneras mas rigurosas de encontrar estos valores.
Los eigenvectores que debió encontrar son:
v
1
=10
v
1
1
0
v
2
=01
v
2
0
1
Y los eigenvalores correspondientes son:
λ
1
=3
λ
1
3
λ
2
=-1
λ
2
-1
Muestre que estos dos vectores,
v
1
=11
v
1
1
1
v
2
=1-1
v
2
1
-1
son eigenvectores de
A
A, donde
A=(
3-1
-13
)
A
3
-1
-1
3
. También encuentre los eigenvalores correspondientes.
Para poder probar que estos dos vectores son eigenvectores, mostraremos que estas afirmaciones cumplen con los requisitos que indica la definición.
A
v
1
=(
3-1
-13
)11=22
A
v
1
3
-1
-1
3
1
1
2
2
A
v
2
=(
3-1
-13
)1-1=4-4
A
v
2
3
-1
-1
3
1
-1
4
-4
Este resultado nos muestra que
A
A
solo escala los dos vectores (es decir
cambia sus longitudes) y esto prueba que la ecuación 2 es cierta para los siguientes dos eigenvalores que se le pidió que encontrara:
λ
1
=2
λ
1
2
λ
2
=4
λ
2
4
.
Si quiere convencerse más, entonces también se pueden graficar los vectores y su producto correspondiente con
AA para ver los resultados como una versión escalada de los vectores originales
v
1
v
1
y
v
2
v
2
.
En los ejemplos anteriores, confiamos en su entendimiento de la definición y de algunas observaciones para encontrar y probar los valores de los eigenvectores y eigenvalores. Sin embrago como se puede dar cuenta, encontrar estos valores no siempre es fácil. A continuación veremos un método matemático para calcular eigenvalores y eigenvectores de una matriz.
Encontrar
λ∈C
λ
tal que
v≠0
v
0
,
donde
0
0 es el “vector cero”. Empezaremos con la ecuación 2, trabajemos de la siguiente manera mientras encontramos una manera explicita de calcular
λ
λ.
Av=λv
A
v
λ
v
Av−λv=0
A
v
λ
v
0
(A−λI)v=0
A
λ
I
v
0
En el paso previo, usamos el hecho de que
λv=λIv
λ
v
λ
I
v
donde II es la matriz identidad.
I=(
10…0
01…0
00⋱⋮
0……1
)
I
1
0
…
0
0
1
…
0
0
0
⋱
⋮
0
…
…
1
Por lo tanto,
A−λI
A
λ
I
es justo una matriz nueva.
Dada la siguiente matriz,
AA, entonces podemos encontrar nuestra nueva matriz,
A−λI
A
λ
I
.
A=(
a1,1a1,2
a2,1a2,2
)
A
a
1
1
a
1
2
a
2
1
a
2
2
A−λI=(
a1,1−λa1,2
a2,1a2,2−λ
)
A
λ
I
a
1
1
λ
a
1
2
a
2
1
a
2
2
λ
Si
(A−λI)v=0
A
λ
I
v
0
para algún
v≠0
v
0
, entonces
A−λI
A
λ
I
es no invertible. Esto quiere decir:
det(A−λI)=0
A
λ
I
0
este determinante (el mostrado arriba) se vuelve una expresión polinomial (de grado nn). Véase el siguiente ejemplo para entender mejor.
Empezando con la matriz AA
(mostrada a continuación), encontremos la expresión polinomial, donde nuestros eigenvalores serán variables dependientes.
A=(
3-1
-13
)
A
3
-1
-1
3
A−λI=(
3−λ-1
-13−λ
)
A
λ
I
3
λ
-1
-1
3
λ
det(A−λI)=3−λ2−-12=λ2−6λ+8
A
λ
I
3
λ
2
-1
2
λ
2
6
λ
8
λ=24
λ
2
4
Empezando con la matriz AA
(mostrada a continuación),encontremos la expresión polinomial, donde nuestros eigenvalores serán variables dependientes.
A=(
a1,1a1,2
a2,1a2,2
)
A
a
1
1
a
1
2
a
2
1
a
2
2
A−λI=(
a1,1−λa1,2
a2,1a2,2−λ
)
A
λ
I
a
1
1
λ
a
1
2
a
2
1
a
2
2
λ
det(A−λI)=λ2−(a1,1+a2,2)λ−a2,1a1,2+a1,1a2,2
A
λ
I
λ
2
a
1
1
a
2
2
λ
a
2
1
a
1
2
a
1
1
a
2
2
Si no lo han notado, calcular los eigenvalores es equivalente a calcular las raíces de
det(A−λI)=
c
n
λn+
c
n
−
1
λn−1+…+
c
1
λ+
c
0
=0
A
λ
I
c
n
λ
n
c
n
−
1
λ
n
1
…
c
1
λ
c
0
0
Por lo tanto usando unos pequeños cálculos para resolver las raíces de nuestro polinomio, podemos encontrar los eigenvalores de la matriz.
Dado un eigenvalor,
λ
i
λ
i
, el eigenvector asociado esta dado por
Av=
λ
i
v
A
v
λ
i
v
A
v
1
⋮
v
n
=
λ
1
v
1
⋮
λ
n
v
n
A
v
1
⋮
v
n
λ
1
v
1
⋮
λ
n
v
n
conjunto de nn ecuaciones con
nn incognitas. Simplemente se resuelven las solve the nn
ecuaciones para encontrar los eigenvectores.
El decir que los eigenvectores de AA,
v
1
v
2
…
v
n
v
1
v
2
…
v
n
,
generan el subespacio
Cn
n
, significa que
v
1
v
2
…
v
n
v
1
v
2
…
v
n
son linealmente independientes y que podemos escribir cualquier
x∈Cn
x
n
como
x=
α
1
v
1
+
α
2
v
2
+…+
α
n
v
n
x
α
1
v
1
α
2
v
2
…
α
n
v
n
(3)
donde
α
1
α
2
…
α
n
∈C
α
1
α
2
…
α
n
Todo lo que estamos haciendo es reescribir
xx en términos de los eigenvectores de
AA. Entonces,
Ax=A(
α
1
v
1
+
α
2
v
2
+…+
α
n
v
n
)
A
x
A
α
1
v
1
α
2
v
2
…
α
n
v
n
Ax=
α
1
A
v
1
+
α
2
A
v
2
+…+
α
n
A
v
n
A
x
α
1
A
v
1
α
2
A
v
2
…
α
n
A
v
n
Ax=
α
1
λ
1
v
1
+
α
2
λ
2
v
2
+…+
α
n
λ
n
v
n
=b
A
x
α
1
λ
1
v
1
α
2
λ
2
v
2
…
α
n
λ
n
v
n
b
por lo tanto podemos escribir,
x=∑i
α
i
v
i
x
i
α
i
v
i
Y esto nos lleva a la siguiente representación del sistema:
donde en la figura 4 tenemos,
b=∑i
α
i
λ
i
v
i
b
i
α
i
λ
i
v
i
Descomponiendo nuestro vector, xx, en una combinación de eigenvectores, la solución de
Ax
A
x
esta dada en piezas “fáciles de digerir".
Para la siguiente matriz, AA y
vector xx, resuélvase por sus productos. Trate de resolverlos por los dos diferentes métodos: directamente y usando eigenvectores.
A=(
3-1
-13
)
A
3
-1
-1
3
x=53
x
5
3
Método Directo (usese la multiplicación básica de matrices)
Ax=(
3-1
-13
)53=124
A
x
3
-1
-1
3
5
3
12
4
Eigenvectores (use los eigenvectores
y eigenvalores que se encotraron anteriormente para esta misma matriz)
v
1
=11
v
1
1
1
v
2
=1-1
v
2
1
-1
λ
1
=2
λ
1
2
λ
2
=4
λ
2
4
Como se muestra en la ecuación 3, queremos representar xx como la suma de sus eigenvectores escalados. Para este caso tenemos:
x=4
v
1
+
v
2
x
4
v
1
v
2
x=53=411+1-1
x
5
3
4
1
1
1
-1
Ax=A(4
v
1
+
v
2
)=
λ
i
(4
v
1
+
v
2
)
A
x
A
4
v
1
v
2
λ
i
4
v
1
v
2
Por lo tanto, tenemos
Ax=4×211+41-1=124
A
x
4
2
1
1
4
1
-1
12
4
Nótese que el método usando eigenvectores
no requiere multiplicación de matrices. . Esto puede parecer mas complicado hasta ahora, pero, imagine que AA es de dimensiones muy grandes.
"Señales y Sistemas is a Spanish translation of Dr. Rich Baraniuk's collection Signals and Systems (col10064). The translation was coordinated by an an assistant electrical engineering professor […]"