Solucionario 3era Unidad
Solucionario 3era Unidad
Solucionario 3era Unidad
FACULTAD DE INGENIERIA
ESC. PROF. ING. ELECTRÓNICA
Problema Nº 01
Considere el sistema con realimentación unitaria y F.T. directa:
β>1
Entonces la F.T. de lazo abierto del sistema compensado es Gc(s)G(s).
Pero sabemos que:
UNIVERSIDAD PRIVADA DE TACNA
FACULTAD DE INGENIERIA
ESC. PROF. ING. ELECTRÓNICA
Ts 1 K
Kv lim sGc( s )G ( s ) lim s 4.5
s 0 s 0 Ts 1 s( s 1)(0.5s 1)
» K=4.5;
» p1=[1 0];
» p2=[1 1];
» p=conv(p1,p2);
» p3=[0.5 1];
» den=conv(p,p3);
» num=K;
» bode(num,den)
UNIVERSIDAD PRIVADA DE TACNA
FACULTAD DE INGENIERIA
ESC. PROF. ING. ELECTRÓNICA
» margin(num,den)
Bode Diagrams
Gm=-3.5 dB (Wcg=1.4); Pm=-10.3 deg. (Wcp=1.7)
100
50
Phase (deg); Magnitude (dB)
-50
200
-200
-400
-2 -1 0 1
10 10 10 10
Frequency (rad/sec)
Se caracteriza por:
Un par polo-cero muy próximos al origen.
Muy próximos entre sí.
El polo más cerca del origen.
s 0.05
= 0.45
s 0.005
( s 0.05) 1
0.45
G (s)G(s)=
c
( s 0.005) s ( s 1)(0.5s 1)
Programa:
K=0.45;
» z=[1 0.05];
» p1=[1 0.005];
UNIVERSIDAD PRIVADA DE TACNA
FACULTAD DE INGENIERIA
ESC. PROF. ING. ELECTRÓNICA
» p2=[1 1];
» p3=[0.5 1];
» p4=[1 0];
» num=K*z;
» p=conv(p1,p2);
» p5=conv(p3,p4);
» den=conv(p,p5);
» bode(num,den);
» margin(num,den)
Bode Diagrams
Gm=15.9 dB (Wcg=1.4); Pm=49.7 deg. (Wcp=0.4)
100
50
Phase (deg); Magnitude (dB)
-50
-100
0
-100
-200
-300
-4 -3 -2 -1 0 1
10 10 10 10 10 10
Frequency (rad/sec)
UNIVERSIDAD PRIVADA DE TACNA
FACULTAD DE INGENIERIA
ESC. PROF. ING. ELECTRÓNICA
Problema Nº 02
SOLUCION Nº02
El criterio de Cauchy
El criterio de Cauchy (del análisis complejo) establece que cuando se recorre un
camino cerrado en el plano complejo, y se lo mapea através de una función compleja
G(s), la cantidad de veces que el gráfico de G(s) rodea el origen es igual a la cantidad de
ceros de G(s) menos la cantidad de polos de G(s) encerrados por el camino de
frecuencias . Los rodeos al origen se cuentan como positivos si están en la misma
dirección que el camino cerrado original o negativos si están en la dirección contraria.
G(s)
---------
1 + G(s)
Si 1+ G(s) rodea el origen, entonces G(s) encerrará el punto -1. Como nos interesa la
estabilidad a lazo cerrado, queremos saber si habrá algún polo a lazo cerrado (ceros de
1 + G(s)) en el semiplano derecho. Más adelante se verán mayores detalles sobre su
determinación. Enonces:
Recuerde que del criterio de Cauchy el número N de veces que el gráfico de G(s)H(s)
rodea -1 es igual al número Z de ceros de 1 + G(s)H(s) encerrado por el camino de
frecuencias menos el número P de polos de 1 + G(s)H(s) encerrado por el camino de
frecuencias (N = Z - P). Teniendo especial cuidado en las funciones de transferencia a
lazo cerrado y abierto, así como en numeradores y denominadores, debe
autoconvencerse que:
Otra forma de verlo es imaginarse sentado sobre el punto -1 y seguir el diagrama del
principio al final. Ahora pregúntese: Cuántas veces giré la cabeza 360 grados
completos? De nuevo, si el movimiento fué a-reloj, N es positivo, y si el movimiento
fué contra-reloj, N es negativo.
K(s^2 + 10 s + 24)
H(s)= 1 ; G(s) = ---------------
s^2 - 8 s + 15
Este sistema tiene una ganancia K que puede variarse para modificar la respuesta del
sistema a lazo cerrado. Sin embargo, veremos que podemos variar esta ganancia sólo
entre ciertos límites, dado que debemos asegurar que nuestro sistema a lazo cerrado
será estable. Esto es lo que buscaremos: el rango de ganancias que estabiliza a este
sistema a lazo cerrado.
roots([1 -8 15])
ans =
5
3
Los polos de la función de transferencia a lazo abierto son ambos positivos. Por lo
tanto, necesitamos dos rodeos contra-reloj (N = -2) del Diagrama de Nyquist para tener
un sistema estable a lazo cerrado (Z = P + N). Si la cantidad de rodeos es menor que
dos o los rodeos no son contra-reloj, nuestro sistema será inestable.
hay dos rodeos en contra-reloj a -1. Por lo tanto, el sistema es estable para una
ganancia de 1. Ahora veamos cómo se comporta el sistema si incrementáramos la
ganancia a 20:
El diagrama se expandió. Por lo tanto, sabemos que el sistema será estable sin
importar qué tanto incrementamos la ganancia. Sin embargo, si decrementáramos la
ganancia, el diagrama se contraerá y el sistema puede volverse inestable. Veamos qué
sucede para una ganancia de 0.5:
El sistema ahora es inestable. Por prueba y error encontramos que este sistema se
inestabilizará con ganancias menores que 0.80. Podemos verificar nuestras respuestas
ya sea agrandando con zoom los diagramas de Nyquist; como mirando las respuestas
al escalón a lazo cerrado para ganancias de 0.79, 0.80, and 0.81.
Si tiene problemas para contar los rodeos de Nyquist, le sugerimos usar nyquist1. El
número de rodeos en contra-reloj se verá en su pantalla (recuerde que este número
representa realmente N negativos, e.d.. Si nyquist1 muestra 2, es N = -2) así como la
cantidad de POLOS REALES positivos a lazo abierto y cerrado.
Problema Nº 03
Se desea compensar en adelanto un sistema con realimentación unitaria a la función
de transferencia F.T. teniendo como planta G(s) como se indica.
ξ = 0.5 y ωn = 4.
Hacer uso del modelo de Lugar geométrico de las raíces.
NO SE HA SOLUCIONADO EL Nº03
Problema Nº 04
Considere el sistema con realimentación unitaria y F.T. directa:
» num=K;
» bode(num,den);
» margin(num,den)
Bode Diagrams
Gm=43.7 dB (Wcg=25.1); Pm=51.8 deg. (Wcp=1.6)
50
Phase (deg); Magnitude (dB)
-50
-50
-100
-150
-200
-2 -1 0 1 2
10 10 10 10 10
Frequency (rad/sec)
Se debe hallar el error estático de velocidad para este sistema que es tipo 1, la
entrada debe ser una rampa y por lo tanto tenemos que:
Ts 1 K
Kv lim sGc( s )G ( s ) lim s 4.5
s 0 s 0 Ts 1 s( s 1)(0.5s 1)
» K=4.5;
» p1=[1 0];
» p2=[1 1];
» p=conv(p1,p2);
» p3=[0.5 1];
UNIVERSIDAD PRIVADA DE TACNA
FACULTAD DE INGENIERIA
ESC. PROF. ING. ELECTRÓNICA
» den=conv(p,p3);
» num=K;
» bode(num,den)
» margin(num,den)
Bode Diagrams
Gm=-3.5 dB (Wcg=1.4); Pm=-10.3 deg. (Wcp=1.7)
100
50
Phase (deg); Magnitude (dB)
-50
200
-200
-400
-2 -1 0 1
10 10 10 10
Frequency (rad/sec)
Se caracteriza por:
Un par polo-cero muy próximos al origen.
Muy próximos entre sí.
El polo mas cerca del origen.
s 0.05
0.45
=
s 0.005
( s 0.05) 1
0.45
G (s)G(s)=
c
( s 0.005) s ( s 1)(0.5s 1)
Programa:
UNIVERSIDAD PRIVADA DE TACNA
FACULTAD DE INGENIERIA
ESC. PROF. ING. ELECTRÓNICA
K=0.45;
» z=[1 0.05];
» p1=[1 0.005];
» p2=[1 1];
» p3=[0.5 1];
» p4=[1 0];
» num=K*z;
» p=conv(p1,p2);
» p5=conv(p3,p4);
» den=conv(p,p5);
» bode(num,den);
» margin(num,den)
Bode Diagrams
Gm=15.9 dB (Wcg=1.4); Pm=49.7 deg. (Wcp=0.4)
100
50
Phase (deg); Magnitude (dB)
-50
-100
0
-100
-200
-300
-4 -3 -2 -1 0 1
10 10 10 10 10 10
Frequency (rad/sec)