Tarea1 LuisÁlvarez
Tarea1 LuisÁlvarez
Tarea1 LuisÁlvarez
MAESTRÍA EN ELECTRICIDAD
COHORTE 2
PRÁCTICA 1
AUTOR
MATERIA
DOCENTE
2021
2
Al remplazar los valores de 𝑇𝐴 , 𝑇𝐸 y 𝑇′𝐷𝑂 y efectuando los cálculos en (4), se tiene lo siguiente:
0.0052𝑠 3 + 0.1224𝑠 2 + 0.66𝑠 + (𝐾 + 1) = 0 (5)
Aplicando el criterio de Routh – Hurwitz para (5), se obtiene lo siguiente:
𝑠3 0.0052 0.66 Cálculo intermedio para determinar valor de 𝑠1 :
0.1224 ∗ 0.66 − 0.0052 ∗ (𝐾 + 1)
𝑠2 0.1224 K+1 = 0.66 − 0.0425 ∗ (𝐾 + 1)
0.1224
Para que la respuesta del sistema en lazo cerrado sea estable, el coeficiente obtenido para 𝑠1
que contiene al parámetro K debe ser mayor a 0, pero para que sea marginalmente estable basta
con igualarlo a 0. A continuación, se muestra lo anterior:
0.6175 − 0.0425𝐾 = 0
0.6175
𝐾= = 14.5294
0.0425
Para obtener una respuesta marginalmente estable para este caso, K debe ser igual a 14.5294.
Con lo anterior, 𝐻(𝑠) queda expresado de la siguiente manera:
𝐾
𝐻(𝑠) =
0.0052𝑠 + 0.1224𝑠 2 + 0.66𝑠 + (𝐾 + 1)
3
14.5294 (6)
𝐻(𝑠) =
0.0052𝑠 3 + 0.1224𝑠 2 + 0.66𝑠 + 15.5294
Utilizando el software de Matlab para corroborar los resultados obtenidos, se muestra la Figura
1. que contiene los comandos empleados y las gráficas de ubicación de las raíces en el plano
complejo y la respuesta del sistema ante una señal de impulso en la entrada en t = 0 s.
a)
b) c)
Figura 1. Análisis de respuesta marginalmente estable en lazo cerrado del sistema de control AVR propuesto utilizando
Matlab para una ganancia K = 14.5294. a) Ventada de comandos de Matlab. b) Ubicación de los polos en el plano complejo.
c) Respuesta del sistema a una señal de impulso en la entrada en t = 0 s.
4
Como se puede observar en la Figura 1., los polos de 𝐻(𝑠) se ubican en el semiplano complejo
izquierdo, sin embargo, dos de ellos se encuentran sobre el eje imaginario que producen una
respuesta inestable para una señal de impulso como entrada, pero transcurridos alrededor de
6000 s., el sistema retorna a un punto estable.
Por lo anterior, se disminuye la ganancia K tomando como límite el valor determinado
anteriormente para una respuesta marginalmente estable con el propósito de evidenciar una
respuesta estable y la ubicación de los polos de 𝐻(𝑠) en el semiplano complejo izquierdo y
ninguno de ellos sobre el eje imaginario. Se debe tomar en cuenta que K puede disminuir hasta
un valor de -1 para que los coeficientes de la columna pivote, definida por el método de Routh
– Hurwitz, no presenten variación de signos entre ellos. En la Figura 2. se despliega los
resultados obtenidos variando el parámetro K a 10.0 en el entorno de Matlab.
a)
b) c)
Figura 2. Análisis de respuesta estable en lazo cerrado del sistema de control AVR propuesto utilizando Matlab para una
ganancia K = 10.0. a) Ventada de comandos de Matlab. b) Ubicación de los polos en el plano complejo. c) Respuesta del
sistema a una señal de impulso en la entrada en t = 0 s.
Con lo anterior, se comprueba que para valores de K menores a 14.5924 la respuesta de 𝐻(𝑠)
es estable porque sus polos se ubican en el semiplano complejo izquierdo y ninguno sobre el
eje imaginario.
5
2. Si se aplica una compensación, en donde TC = 0.36 s. Analizar en Matlab las gráficas de los
polos y su respuesta en el tiempo.
Para este caso, 𝐺𝑐 (𝑠) queda expresado de acuerdo a lo siguiente:
𝐾(1 + 𝑠𝑇𝐶 ) (7)
𝐺𝑐 (𝑠) =
(1 + 𝑠𝑇𝐴 )(1 + 𝑠𝑇𝐸 )(1 + 𝑠𝑇′𝐷𝑂 )
Considerando que, 𝑇𝐶 = 𝑇𝐸 , (7) queda definido como:
𝐾 (8)
𝐺𝑐 (𝑠) =
(1 + 𝑠𝑇𝐴 )(1 + 𝑠𝑇′𝐷𝑂 )
y, remplazando (8) en (3) se tiene:
𝐾
1+ =0
(1 + 𝑠𝑇𝐴 )(1 + 𝑠𝑇′𝐷𝑂 )
(1 + 𝑠𝑇𝐴 )(1 + 𝑠𝑇 ′ 𝐷𝑂 ) + 𝐾 = 0
𝑠 2 𝑇𝐴 𝑇′𝐷𝑂 + 𝑠(𝑇𝐴 + 𝑇′𝐷𝑂 ) + (𝐾 + 1) = 0 (9)
Al remplazar los valores de 𝑇𝐴 y 𝑇′𝐷𝑂 y efectuando los cálculos en (9), se obtiene lo siguiente:
0.0144𝑠 2 + 0.3𝑠 + (𝐾 + 1) = 0 (10)
Aplicando el criterio de Routh – Hurwitz para (10), se tiene lo siguiente:
𝑠2 0.0144 K+1 Cálculo intermedio para determinar valor de 𝑠 0 :
0.3 ∗ (𝐾 + 1)
𝑠1 0.3 0 =𝐾+1
0.3
𝑠0 K+1
Al igual que en el numeral 1., para que la respuesta del sistema en lazo cerrado sea estable, el
coeficiente obtenido para 𝑠 0 que contiene al parámetro K debe ser mayor a 0, pero para que sea
marginalmente estable basta con igualarlo a 0. A continuación, se muestra lo anterior:
𝐾+1=0
𝐾 = −1
Para obtener una respuesta marginalmente estable para este caso, K debe ser igual a -1. Con lo
anterior, 𝐻(𝑠) queda expresado de la siguiente manera:
𝐾
𝐻(𝑠) = 2
0.0144𝑠 + 0.3𝑠 + (𝐾 + 1)
1 (6)
𝐻(𝑠) = −
0.0144𝑠 2 + 0.3𝑠
Utilizando el software de Matlab para corroborar los resultados obtenidos, se muestra la Figura
3. que contiene los comandos empleados y las gráficas de ubicación de las raíces en el plano
complejo y la respuesta del sistema ante una señal de impulso en la entrada en t = 0 s. Como se
puede observar, se obtiene un polo situado en el origen que es el que influye en la respuesta
marginalmente estable del sistema. El otro polo se encuentra sobre el eje real, pero en el
6
semiplano complejo izquierdo, por lo cual, no influye en la respuesta del sistema como el
anterior.
a)
c)
b)
Figura 3. Análisis de respuesta marginalmente estable en lazo cerrado del sistema de control AVR propuesto utilizando
Matlab para una ganancia K = -1 y un compensador. a) Ventada de comandos de Matlab. b) Ubicación de los polos en el
plano complejo. c) Respuesta del sistema a una señal de impulso en la entrada en t = 0 s.
Como en el caso del numeral 1., se procede a variar el valor de la ganancia K tomando como
límite el valor determinado anteriormente para una respuesta marginalmente estable con el
propósito de evidenciar una respuesta estable, y, la ubicación de los polos de 𝐻(𝑠) en el
semiplano complejo izquierdo con ninguno de ellos sobre el eje imaginario. En la Figura 2.
Para una respuesta estable, se debe considerar que los coeficientes de la columna pivote,
definida por el método de Routh – Hurwitz, no presenten variación de signos entre ellos, por
ello K debe ser siempre mayor a -1. En la Figura 4. se despliegan los resultados obtenidos
variando el parámetro K a 10.0 en el entorno de Matlab.
7
a)
b) c)
Figura 4. Análisis de respuesta estable en lazo cerrado del sistema de control AVR propuesto utilizando Matlab para una
ganancia K = 10.0 y un compensador. a) Ventada de comandos de Matlab. b) Ubicación de los polos en el plano complejo. c)
Respuesta del sistema a una señal de impulso en la entrada en t = 0 s.