Control Moderno - Ing. Electrónica Ejercicio Resuelto 2

Control Moderno - Ing. Electrónica
Ejercicio Resuelto 2: Transformaciones
1.
Introducción
Las transformaciones nos ayudan a percibir un problema de otra forma, sin que ello implique en
absoluto una modificación en las caracterı́sticas del mismo. Si a un sistema se le aplica una transformación lineal, éste conservará todas sus cualidades, como son: los autovalores, observabilidad,
controlabilidad, etc. El cambio de base del sistema de coordenadas de referencia que implica una
transformación es equivalente a una multiplicación de matrices: x = P z , donde existe una matriz
que relaciona ambas bases entre sı́. Esta matriz P es de suma importancia ya que es la que nos
permitirá pasar de un sistema a otro para poder aprovechar las ventajas de cálculo y retornar al sistema original con los resultados finales. Veamos que sucede con un modelo de estados al aplicar una
transformación lineal x = P z, donde x es el vector de estados original y z el vector transformado.
Sea el modelo de estados
ẋ = Ax + Bu, y = Cx.
(1)
Si aplicamos la trasformación x = P z y reemplazamos en (1), obtendremos
P ż = AP z + Bu,
y = CP z.
(2)
Premultiplicando por P −1
ż = P −1 AP z + P −1 Bu,
y = CP z
(3)
y reordenando, obtendremos un nuevo modelo de estados
ż = Az z + Bz u,
y = Cz z
(4)
donde Az = P −1 AP , Bz = P −1 B y Cz = CP .
2.
Aplicaciones
Estas transformaciones son útiles para expresar el modelo original en otras variables de estados
donde ciertas caracterı́sticas propias sean más fácilmente observadas. Tal vez las más utilizadas son
aquellas que nos permiten obtener un modelo de estados en forma diagonal (matriz A diagonal),
en forma de Jordan o en forma canónica controlable.
2.1. Forma diagonal
Consideremos el siguiente sistema
ẋ =
"
−14 −4
4 −4
#
x+
"
1
1
#
u,
y=
h
1 1
i
x.
(5)
Buscamos un modelo de estados alternativo cuya matriz A sea diagonal. Si el sistema presenta
autovalores reales y distintos (es este caso) la nueva matriz Ad será1
Ad =
1
"
λ1 0
0 λ2
#
,
¿Por qué la matriz Ad tiene esta estructura, con los autovalores en la diagonal?
1
(6)
donde λ1 , λ2 son los autovalores de A.
Para obtener el modelo diagonal debemos hallar la matriz de transformación P , tal que
P
−1
AP = Ad =
"
λ1 0
0 λ2
#
.
(7)
Esta condición puede reescribirse como
AP = P
"
#
λ1 0
0 λ2
(8)
y, si consideramos P compuesta por dos vectores columna v1 y v2 , es decir, P =
entonces
"
#
i h
i λ
h
0
1
A v1 v2 = v1 v2
0 λ2
h
i
v1 v2 ,
(9)
Luego, para hallar los vectores vi debemos resolver dos sistemas de ecuaciones de la forma
Avi = λi vi ,
(10)
donde vi es el autovector asociado al autovector λi de la matriz A.
Apliquemos estas ideas al ejemplo. Primero debemos encontrar los autovalores del sistema, esto es,
debemos resolver
|λI2×2 − A| = 0,
(λ + 14)
−4
4
(λ + 4)
de donde obtenemos λ1 = −6 y λ2 = −12.
= λ2 + 18λ + 72 = 0,
Una vez hallados los autovalores, debemos plantear la ecuación (10) para hallar los autovectores
asociados a cada uno de ellos.
Para λ1 = −6, la (10) corresponde a
reemplazando A y v1 =
h
v11 v21
"
iT
−14 −4
4 −4
Av1 = −6v1 ,
#"
v11
v21
#
=
"
−6v11
−6v21
→
(
#
.
Expandiendo,
(
−14v11 − 4v21 = −6v11
4v11 − 4v21
= −6v21
v11 = −v21 /2
v21 =
v21
Adoptando, por ejemplo, v21 = 1 entonces v11 = 1/2. El vector resultante puede dividirse por su
h
√
√ iT
norma para normalizarlo, con lo cual, v1 = 2/ 5 −1/ 5 .
Para λ2 = −12, debemos resolver Av2 = −12v2 , es decir,
"
donde v2 =
h
eligiendo v22
v12 v22
−14 −4
4 −4
#"
v12
v22
#
=
"
−12v12
−12v22
#
,
iT
. El sistema de ecuaciones anterior es equivalente a v12 = −2v22 . Luego,
h
√
√ iT
= 1 y normalizando v2 = 1/ 5 −2/ 5
.
Control Moderno - Ing. Electrónica: Ejercicio Resuelto 2: Transformaciones
2
Finalmente, la matriz de transformación P resulta
"
√
√ #
2/√5
1/√5
P =
−1/ 5 −2/ 5
y el nuevo modelo de estados
ż = Ad z + Bd u,
y = Cd z,
donde
Ad = P
−1
AP =
"
−6
0
0 −12
#
,
Bd = P
−1
B=
"
2
−1
#
,
Cd = CP =
h
1/2 1
i
.
2.2. Forma de Jordan
Consideremos el siguiente modelos de estados




−1 −2
0
0




ẋ =  0 −1
0  x +  1  u,
3 −3 −5
1
y=
h
1 0 0
i
.
(11)
Al resolver la ecuación
|A − λI3×3 | = 0
(12)
obtenemos dos autovalores distintos λ1 = −5 y λ2 = −1, este último de multiplicidad 2. Cuando
A tiene autovalores reales pero no todos son distintos, no siempre tendremos los suficientes autovectores linealmente independientes para expresar la matriz A en forma diagonal. Si no tenemos
autovectores linealmente independientes la matriz de transformación P será singular. En este caso,
podemos hallar los denominados autovectores generalizados y expresar la matriz A en la forma de
Jordan, esto es


λ1 0 0


Az =  0 λ2 1  .
0 0 λ2
Observemos que Az estará formada por bloques, los bloques correspondientes a los autovalores
de multiplicidad 1 serán diagonales y los correspondientes a autovalores de multiplicidad mayor
tendrán la forma
"
#
λi 1
.
0 λi
Comencemos con el autovector asociado a λ1 = −5. Como este tiene multiplicidad 1 podemos hallar
un autovector asociado a λ1 usando (10), es decir,
(A + 5I3×3 )v1 = 0,
definiendo v1 =
h
v11 v21 v31
iT
(13)
, la (13) es equivalente a


 4v11 − 2v21
4v
21

 3v − 3v
11
21
de donde se verifica que el vector v1 =
h
0 0 1
iT
= 0,
= 0,
= 0,
(14)
es solución de (13).
Control Moderno - Ing. Electrónica: Ejercicio Resuelto 2: Transformaciones
3
El problema surge en los autovectores asociados al autovalor doble λ2 = −1. Si buscamos las
restantes columnas de P a partir de (10) sólo obtendremos vectores linealmente dependientes. Para
salvar el inconveniente reemplazamos (10) por las siguientes ecuaciones
(A − λ2 I3×3 )v2 = 0,
(15)
(A − λ2 I3×3 )v3 = v2 .
(16)
Reemplazando en (15) los valores del ejemplo considerado y v2 =




h
v12 v22 v32
iT
tenemos

0 −2
0
v12
0


 

0
0   v22  =  0  ,
 0
3 −3 −4
0
v32
(17)
que simplificando resulta en dos ecuaciones
(
−2v22 = 0,
3v12 − 3v22 − 4v32 = 0,
(18)
que son equivalente a 3v12 = 4v32 . Luego, adoptando v32 = 3 y normalizando, tenemos
v2 =
h
4/5 0 3/5
iT
.
Obtenido el vector v2 , reemplazamos en (16) y resolvemos el sistema de ecuaciones resultante

Simplificando,




0 −2
0
v13
4


 

0
0   v23  =  0  .
 0
v33
3 −3 −4
3
(
−2v23 = 4,
3v13 − 3v23 − 4v33 = 3.
(19)
(20)
Luego, esto implica que v23 = −2 y que 3v13 − 4v33 = −3. Adoptando v13 = 3 y normalizando,
h
√
√
√ iT
v2 = 3/ 22 −2/ 22 3/ 22 .
Finalmente, a partir de los tres autovectores obtenemos la matriz de transformación

√ 
0 4/5
3/√22


P = 0
0 −2/√22  .
1 3/5
3/ 22
Con esta transformación, la matriz Az tiene la forma Jordan correspondiente


−5
0
0


1 .
Az =  0 −1
0
0 −1
Control Moderno - Ing. Electrónica: Ejercicio Resuelto 2: Transformaciones
4