14-Simulacion Proc Ind - Respuesta en Frecuencia
14-Simulacion Proc Ind - Respuesta en Frecuencia
14-Simulacion Proc Ind - Respuesta en Frecuencia
>> t=1:0.1:10;
A=1 ; w=4;
plot(t,A*sin(w*t))
11.4.2 Criterio de estabilidad de ruta
Routh (1905) publicó una técnica analítica para determinar si las raíces de un polinomio tienen
partes reales positivas. De acuerdo con su criterio general de estabilidad, un sistema en lazo
cerrado será estable solo si todas las raíces de la ecuación característica tienen partes reales
negativas. Por tanto, al aplicar la técnica de Routh para analizar los coeficientes de la ecuación
característica, podemos determinar si el sistema en lazo cerrado es estable. Este enfoque se
conoce como el criterio de estabilidad de Routh. Solo se puede aplicar a sistemas cuyas
ecuaciones características son polinomios en s. Por lo tanto, el criterio de estabilidad de Routh
no es directamente aplicable a los sistemas que contienen retardos de tiempo, lo que limita su
aplicabilidad para procesos realistas. Si e−θs se reemplaza por una aproximación de Padé,
entonces se puede realizar un análisis de estabilidad aproximado utilizando la matriz de Routh.
Se puede realizar un análisis de estabilidad exacto de los sistemas que contienen retardos de
tiempo mediante la búsqueda directa de raíces o mediante el uso de un análisis de respuesta de
frecuencia y el criterio de estabilidad de Bode o Nyquist.
Se utilizaron técnicas de transformada de Laplace para calcular respuestas transitorias a partir de
funciones de transferencia. Este capítulo se centra en una forma alternativa de analizar sistemas
dinámicos mediante el análisis de respuesta de frecuencia. Los conceptos y técnicas de respuesta
de frecuencia juegan un papel importante en el análisis de estabilidad, el diseño del sistema de
control y la evaluación de la robustez.
Históricamente, las técnicas de respuesta de frecuencia proporcionaron el marco conceptual
para las primeras teorías de control y aplicaciones importantes en el campo de las
comunicaciones.
Se presenta un procedimiento para calcular las características de respuesta en frecuencia a partir
del modelo de función de transferencia de cualquier proceso lineal.
Dos conceptos, los criterios de estabilidad de Bode y Nyquist, son aplicables para sistemas de
control de retroalimentación y análisis de estabilidad.
A continuación, presentamos dos métricas útiles para la estabilidad relativa, a saber, los
márgenes de ganancia y de fase. Estas métricas indican qué tan cerca está un sistema de control
de la inestabilidad. Un tema relacionado es la robustez, es decir, la sensibilidad del desempeño
del sistema de control a las variaciones del proceso ya la incertidumbre en el modelo del
proceso.
14.3 DIAGRAMAS DE BODE
El diagrama de Bode (o gráfico de Bode) proporciona una visualización conveniente de las
características de respuesta de frecuencia en las que AR y ϕ se representan en función de ω.
Por lo general, ω se expresa en unidades de radianes/tiempo para simplificar los cálculos de la
tangente inversa (p. ej., Ec. 14-18b) donde los argumentos deben ser adimensionales, es decir,
en radianes. Ocasionalmente, se utiliza una frecuencia cíclica, ω/2π, con unidades de
ciclos/tiempo.
El ángulo de fase ϕ es normalmente expresado en grados en lugar de radianes.
Por razones que se harán evidentes en el siguiente desarrollo, el diagrama de Bode consta de: (1)
una gráfica logarítmica de AR frente a ω y (2) una gráfica semilogarítmica de ϕ frente a ω.
Estos gráficos son particularmente útiles para el análisis rápido de las características de respuesta
y la estabilidad de los sistemas de circuito cerrado.
Procesos de Primer Orden
Procesos de Segundo Orden
%Make Bode plot for G = 5 (0.5s + 1)e^–0.5s/(20s + 1)(4s + 1)
close all
gain = 5;
tdead = 0.5;
num = [0.5 1];
den = [80 24 1];
G = tf (gain*num, den) %Define the transfer function
points = 500; %Define the number of points
ww = logspace (-2, 2, points); %Frequencies to be evaluated
[mag, phase, ww] = bode (G,ww); % Generate numerical
%values for Bode plot
AR = zeros (points, 1); % Preallocate vectors for Amplitude
%Ratio and Phase Angle
PA = zeros (points, 1);
for i = 1 : points
AR(i) = mag (1,1,i)/gain; %Normalized AR
PA(i) = phase (1,1,i) - ((180/pi) *tdead*ww(i));
end
figure
subplot (2,1,1)
loglog(ww, AR)
axis ([0.01 100 0.001 1])
title ('Frequency Response of a SOPTD with Zero')
ylabel('AR/K')
subplot (2,1,2)
semilogx(ww,PA)
axis ([0.01 100 -270 0])
ylabel('Phase Angle (degrees)')
xlabel('Frequency (rad/time)')
Interpretación: https://www.youtube.com/watch?v=SIS3cw5hj8g
>> H = tf([1 0.1 7.5],[1 0.12 9 0 0]);
bode(H)
El diagrama de Bode de ganancias es una recta con pendiente -20, dB/década. Pasa por 0 dB en
la frecuencia de 1 rad/s.
Derivador: G(s) = s
El diagrama de Bode de ganancias es una recta con pendiente 20 dB/década. La recta pasa por
0dB en la frecuencia de 1 rad/s
Polo simple estable
Polo simple estable
>>clear tf
>>H = tf([0 0.3],[1 0.3]);
bode(H)
14.4 Características de las respuestas de los controladores FeedBack.
Para usar el análisis de respuesta de frecuencia para diseñar sistemas de control, se deben
conocer las características relacionadas con la frecuencia de los controladores de
retroalimentación para las formas más utilizadas del controlador PID (Cap 8).
En las siguientes derivaciones, generalmente asumimos que el controlador es acción inversa
(Kc > 0). Si un controlador es de acción directa (Kc < 0), el gráfico AR no cambia porque |Kc| se
utiliza para calcular la magnitud. Sin embargo, el ángulo de fase se desplaza −180∘ cuando Kc es
negativo. Por ejemplo, un controlador proporcional de acción directa (Kc < 0) tiene un ángulo de
fase constante de −180∘.
Es posible usar el valor absoluto de Kc para calcular ϕ al diseñar sistemas de control de lazo
cerrado, porque las consideraciones de estabilidad (Cap 11) requieren que Kc < 0 solo cuando
KvKpKm < 0.
Esta elección garantiza que el la ganancia de lazo abierto (KOL = KcKvKpKm) siempre será positiva. El
uso de esta convención produce convenientemente ϕ = 0∘ para cualquier controlador
proporcional y, en general, elimina la necesidad de considerar la contribución de cambio de fase
−180∘ de la ganancia negativa del controlador.
Ejercicio: Una máquina de producción industrial tiene la siguiente función de transferencia.
𝑏1 𝑠 + 𝑏0
𝐺 𝑠 = 2
𝑠 + 𝑎1 𝑠 + 𝑎0
Use el diagrama DiagramaBodeMaquina.m
ft =
0.5 s + 100
-------------------
0.5 s^2 + 5 s + 323
Ejercicio: Un sistema de trituración de minerales tiene la siguiente función de transferencia.
𝑏1 𝑠 + 𝑏0
𝐺 𝑠 =
𝑎4 𝑠 4 + 𝑎3 𝑠 3 + 𝑎2 𝑠 2 + 𝑎1 𝑠 + 𝑎0
Use el diagrama DiagramaBodeTrituradora2dof.m para analizar su comportamiento y estabilidad
ft =
2.902e08 s + 4.559e09
--------------------------------------------------------------
7312 s^4 + 9.865e04 s^3 + 3.042e07 s^2 + 2.902e08 s + 4.559e09