Análisis de Lugar de Raíces en Control
Análisis de Lugar de Raíces en Control
Tema 3
Análisis y diseño de sistemas de control
mediante el método del Lugar de las Raíces
Control y Automatización
Grado en Ingeniería de Tecnologías Industriales
Resolución de problemas
26 de marzo de 2021
Control
ControlyyAutomatización
Automatización
2
Resolución de problemas:
➢ Problema 3.2
➢ Problema 3.3
➢ Problema 3.4
➢ Problema 3.5
➢ Problema 3.6
➢ Problema 3.8
Índice ➢ Problema 3.14
➢ Problema 3.15
Resolución de ejercicios de examen:
➢ Ejercicio 3 – Examen 17/12/2019 (Grado de
Ingeniería en Robótica Software)
Control y Automatización
3
Resolución de problemas
Nombramos:
s+2 1
G s = y H s =
s+4 s+1
En MATLAB:
>> rlocus(G*H)
Control
ControlyyAutomatización
Automatización
6
Resolución de problemas
Nombramos:
1 1
G s = y H s =
s s+4 s+5
Se obtiene una fila de ceros si K=180. Por tanto, sustituyendo dicho valor en el
polinomio auxiliar (construido a partir de la fila de encima), podemos obtener el
valor de los puntos de corte.
dK
=−3s2 −18s−20=0 → s=−1,47 (s=−4,52 no pertenece al LDR)
ds
En MATLAB: >> rlocus(G*H)
Control
ControlyyAutomatización
Automatización
11
Resolución de problemas
Nombramos:
1 s+3
G s = y H s =
2
s −4 s+5
Como en los ejercicios previos, se irán dando los pasos necesarios para
realizar el bosquejo del LDR.
Paso 1: IDENTIFICACIÓN Y UBICACIÓN DE POLOS Y CEROS en lazo abierto.
Se dibujan los polos y ceros de la función de transferencia en lazo abierto.
En MATLAB: >> pzmap(G*H).
Control
ControlyyAutomatización
Automatización
12
Resolución de problemas
Control
ControlyyAutomatización
Automatización
15
Resolución de problemas
>> rlocus(G*H)
Control
ControlyyAutomatización
Automatización
16
Resolución de problemas
Nombramos:
s+2 s+3
G s = y H s =
s s+1
En MATLAB:
>> rlocus(G*H)
Control
ControlyyAutomatización
Automatización
20
Resolución de problemas
Problema 3.6. Dibuja el lugar de las raíces del sistema mostrado en la figura.
Analiza la respuesta en función de los posibles valores de K.
s+1
Se tiene: G s = y H s =1
s s+2 s2 +2s+2
Paso 1: IDENTIFICACIÓN Y UBICACIÓN DE POLOS Y CEROS en lazo abierto.
MATLAB: >> pzmap(G*H). Véase en la siguiente diapositiva.
Paso 2: NÚMERO DE RAMAS.
Se tiene: n=4 y m=1.
Por tanto, el número de ramas es: 4.
Control
ControlyyAutomatización
Automatización
21
Resolución de problemas
P s =2,77s2 +8,94=0 →
s=± −3,23=±1,79j
Paso 6: PUNTOS DE RUPTURA
s+1 O SALIDA DEL EJE REAL.
1+KG s H s =0 → 1+K =0
2
s s+2 s +2s+2 s=−1,53±0,53j y s=−0,46±0,53j
s s+2 s2 +2s+2 s4 +4s3 +6s2 +4s no pertenecen al eje real)
→ K=− =−
s+1 s+1 No hay puntos de ruptura sobre el eje real.
dK −4s3 −12s2 −12s−4 s+1 + s4 +4s3 +6s2 +4s −3s4 −12s3 −18s2 −12s−4
= = =0
ds s+1 2 s+1 2
Control
ControlyyAutomatización
Automatización
24
Resolución de problemas
Control
ControlyyAutomatización
Automatización
25
Resolución de problemas
MATLAB:
>> rlocus(G*H)
Problema 3.8. Dibuja el lugar de las raíces del sistema mostrado en la figura.
Analiza la respuesta en función de los posibles valores de K.
Se tiene: G s =
s+3
y H s =1 (realimentación negativa y unitaria)
s s+1 s2+4s+5
P s =7,10s2 +13,47=0 →
s=± −1,89=±1,37j
0.5
Imaginary Axis (seconds )
-1
153,5º
135º
0 45º
-0.5
Se selecciona el polo imaginario -2+j.
Imponemos la fórmula generalizada:
-1
90º σ θpolos − σ θceros=180º
θp +153,5º+135º+90º-45º=180º; θp = 206,5º
-1.5
-3 -2.5 -2 -1.5 -1 -0.5 0
Real Axis (seconds -1 )
El polo complejo con parte imaginaria negativa saldrá con un ángulo “simétrico”.
Control
ControlyyAutomatización
Automatización
31
Resolución de problemas
MATLAB:
>> rlocus(G*H)
Control
ControlyyAutomatización
Automatización
32
Resolución de problemas
Control
ControlyyAutomatización
Automatización
34
Resolución de problemas
4 1 2 1 dK 1 1
1+KG s H s =0 → 1+K =0 → K=− s + s =− s− =0 → s=−1
s s+2 4 2 ds 2 2
Control
ControlyyAutomatización
Automatización
36
Resolución de problemas
En MATLAB:
>> Gtd=tf([4],[1 20 0])
>> H=1
>> rlocfind(Gtd*H,-10+12.53*i)
1 1 20
e) De igual forma que en c), resulta: eSS = lim sE s = lim s =
s→0 s→0 s2 1+ 257 257
El error de velocidad se ha reducido 6,5 veces gracias a la implementación s s+20
del compensador de adelanto propuesto en el apartado d).
Control
ControlyyAutomatización
Automatización
41
Resolución de problemas
1 1 s+zc 1 zc
0,2= lim s s+zc → lim Kc =5→ =10
s→0 s2 1 s→0 s+pc s+1 s+2 pc
1+Kc
s+pc s s+1 s+2
Por simplificar, hemos considerado Kc=1. Por tanto, el error de velocidad se
reducirá 10 veces con cualquier compensador de retardo que cuente con un cero
10 veces mayor que el polo. Faltaría por considerar el requerimiento consistente en
que “se mantenga invariable la respuesta del sistema ante una entrada escalón”.
Para ello, debemos usar los fundamentos teóricos de los polos dominantes (ver
ayuda). Para mantener invariable la respuesta original, debemos introducir ceros y
polos “no dominantes”, es decir, cercanos al origen. Por ejemplo, elegimos: zc=0,05
y pc=0,005. Por tanto:
A continuación, comprobamos utilizando MATLAB como el lugar
s+0,05
Gc s = de las raíces se mantiene prácticamente invariable y la respuesta
s+0,005 temporal del sistema en régimen permanente ante una rampa se
acerca al valor ideal (reducción del error).
Control
ControlyyAutomatización
Automatización
43
Resolución de problemas
Control
ControlyyAutomatización
Automatización
44
Resolución de ejercicios de examen
Control
ControlyyAutomatización
Automatización
45
Resolución de ejercicios de examen
Control
ControlyyAutomatización
Automatización
46
Resolución de ejercicios de examen
Control
ControlyyAutomatización
Automatización
47
Resolución de ejercicios de examen
Control
ControlyyAutomatización
Automatización
48
Resolución de ejercicios de examen
Control
ControlyyAutomatización
Automatización
49
Resolución de ejercicios de examen
Control
ControlyyAutomatización
Automatización
50
Resolución de ejercicios de examen
Control
ControlyyAutomatización
Automatización