Control Digital Con Matlab
Control Digital Con Matlab
Control Digital Con Matlab
La figura de abajo muestra el tpico sistema continuo realimentado que hemos considerado hasta ahora en este tutorial. Casi la totalidad de los controladores continuos pueden implementarse usando electrnica analgica.
El controlador analgico, encerrado en el cuadrado a trazos, puede reemplazarse por un controlador digital, como se muestra abajo, el cual hace la misma tarea de control que el controlador analgico. La diferencia bsica entre estos controladores es que el sistema digital opera con seales discretas (o muestras de la seal sensada) en lugar de seales continuas.
Los diferentes tipos de seales en el esquema digital de arriba pueden representarse por las figuras siguientes.
El propsito de este Tutorial de Control Digital es mostrarle cmo trabajar con funciones discretas ya sea en la forma funcin de transferencia o en la forma espacio de estado para disear sistemas de control digital.
El reloj conectado a los conversores D/A y A/D suministra un pulso cada T segundos y cada D/A y A/D enva una seal solo cuando llega el pulso. El propsito de tener este pulso es que Hzoh(z) tiene slo
muestras u(k) para tratar y produce slo muestras como salida y(k); por lo tanto, Hzoh(z) puede ser implementado como funcin discreta. La filosofa del diseo es la siguiente. Queremos hallar una funcin discreta Hzoh(z) de modo que para una entrada constante a tramos al sistema continuo H(s), la salida muestreada del sistema continuo sea igual a la salida discreta. Suponga que la seal u(k) representa una muestra de la seal de entrada. Existen tcnicas para tomar esta muestra u(k) y mantenerla para producir una seal continua uhat(t) (por u circunflejo). El croquis abajo muestra que uhat(t) se mantuvo constante a u(k) en el intervalo kT a (k+1)T. La operacin de mantener uhat(t) constante durante el tiempo de muestreo se denomina mantenedor de orden cero.
La seal mantenida uhat(t) pasa por H2(s) y el A/D para producir la salida y(k) que ser la misma seal a tramos como si fluyera u(t) continua a travs de H(s) para producir la salida continua y(t).
Mediante Hzoh(z), podemos disear sistemas de control digital tratando solamente con funciones discretas. Nota: Existen ciertos casos donde la respuesta discreta no coincide con la respuesta analgica debido a los circuitos de retencin implementados en sistemas de control digital. Para ms informacin, vea Efecto de retardo asociado al mantenedor.
retencin de orden cero explicada arriba. El comando bsico para este comando c2dm es alguno de los siguientes.
[numDz,denDz] = c2dm (num,den,Ts,'zoh') [F,G,H,J] = c2dm (A,B,C,D,Ts,'zoh')
El tiempo de muestreo (Ts en seg/muestra) debera ser menor que 1/(30*BW), donde BW es el ancho de banda a lazo cerrado .
1. Funcin de Transferencia
suponga que tiene la siguiente funcin de transferencia continua
Asumiendo que el ancho de banda a lazo cerrado es mayor que 1 rad/seg., elegiremos el tiempo de muestreo (Ts) igual a 1/100 seg. Ahora, cree un archivo-m nuevo e ingrese los siguientes comandos.
M=1; b=10; k=20; num=[1]; den=[M b k]; Ts=1/100; [numDz,denDz]=c2dm(num,den,Ts,'zoh')
Luego de correr este archivo-m en la ventana de comandos le debera dar las siguientes matrices numDz y denDz.
numDz = 1.0e-04 * 0 denDz = 1.0000 -1.9029 0.9048 0.4837 0.4678
Nota: Las matrices numerador y denominador estarn representadas en potencias descendentes de z. Para ms informacin sobre representacin, refirase por favor a Representacin Matlab. Ahora tiene la funcin de transferencia en la forma discreta.
2. Espacio de Estado
suponga que tiene el siguiente modelo en espacio de estado continuo
Todas las constantes son las mismas que anteriormente El siguiente archivo-m convierte el espacio de estado continuo a espacio de estado discreto de arriba.
M=1; b=10; k=20; A=[0 -k/M B=[ 0; 1/M]; C=[1 0]; D=[0]; Ts=1/100; [F,G,H,J] = 1; -b/M];
c2dm (A,B,C,D,Ts,'zoh')
Cree un nuevo archivo-m y copie aquellos comandos. Corriendo este archivo-m, la ventana de comandos del Matlab le debera dar las siguiente matrices.
F = 0.9990 -0.1903 G = 0.0000 0.0095 0.9039
0.0095 H = 1 J = 0 0
Ahora tiene el modelo en espacio de estado de tiempo discreto. Nota: Para ms informacin sobre el espacio de estado discreto , refirase a Espacio de Estado Discreto .
La figura de abajo muestra el mapeo de lneas de coeficiente de amortiguamiento constante (zeta) y la frecuencia natural (Wn) del plano s al plano z mediante la expresin de abajo.
Si se dio cuenta en el plano z , el lmite de estabilidad ya no es el eje imaginario, sino el crculo unitario |z| =1. El sistema es estable cuando todos los polos se ubican dentro del crculo unitario e inestable cuando cualquier polo se ubica afuera de l. Para analizar la respuesta transitoria de ubicaciones de polo en el plano z , todava puede aplicarse las siguientes tres ecuaciones usadas en diseos de sistemas continuos .
donde
zeta = coeficiente de amortiguamiento Wn = Frecuencia natural (rad/seg.) Ts = Tiempo de establecimiento Tr = Tiempo de Subida Mp = Mximo sobrepico
Importante: La frecuencia natural (Wn) en el plano z posee las unidades de rad/muestra, pero cuando usa las ecuaciones que aparecen abajo, la Wn debe estar en unidades de rad/seg.. Suponga que tenemos la siguiente funcin de transferencia discreta
Cree un nuevo archivo-m y entre los siguientes comandos. Luego de correr este archivo-m en la ventana de comandos ,ver la figura siguiente con las lneas de coeficiente de amortiguamiento y la frecuencia natural constantes .
numDz=[1]; denDz=[1 -0.3 0.5];
pzmap(numDz,denDz)
axis([-1 1 -1 1])
zgrid
De esta figura, vemos que los polos se ubican aproximadamente a la frecuencia natural de 9pi/20T (rad/muestra) y la razn de amortiguacin de 0.25. Asumiendo que tenemos un tiempo de muestreo de 1/20 seg (que lleva a Wn = 28.2 rad/seg.) y usando las tres ecuaciones de abajo, podemos determinar que este sistema debera tener un tiempo de elevacin de 0.06 seg, el tiempo de establecimiento de 0.65 seg. y el sobrepico mximo del 45% (0.45 ms que el valor de estado estacionario). Obtengamos la respuesta al escaln y veamos si estos son correctos. Agregue los comandos siguientes al archivo-m anterior y ejectelo nuevamente en la ventana de comandos. Debera obtenerse la siguiente respuesta al escaln .
[x] = dstep (numDz,denDz,51); t = 0:0.05:2.5; stairs (t,x)
Como se puede ver del grfico, el tiempo de elevacin, el tiempo de establecimiento y el sobrepico quedan como hemos esperado . Probaremos aqu que podemos usar los lugares de los polos y las tres ecuaciones de arriba para analizar la respuesta transitoria del sistema. Para mayores anlisis con ubicaciones de polos y respuesta transitoria, vea Respuesta Transitoria.
donde G(z) es el compensador implementado en el controlador digital y Hzoh(z) es la funcin de transferencia de la planta en z . La mecnica para dibujar el lugar de las races es la misma en el plano z que en el plano s . Recordemos del Tutorial del Lugar de Races continuo, hemos usado la funcin del Matlab denominada sgrid para hallar la regin del lugar de races que proporciona la ganancia correcta (K). Para el anlisis discreto del lugar de races , usamos la funcin zgrid que posee las mismas caractersticas que sgrid. El comando zgrid(zeta, Wn) dibuja lneas de coeficiente de amortiguamiento (zeta) y frecuencia natural constante (Wn). Suponga que tenemos la siguiente funcin de transferencia discreta
y los requerimientos de tener coeficiente de amortiguamiento mayor que 0.6 y frecuencia natural mayor que 0.4 rad/muestra (esto puede hallarse a partir de los requerimientos de diseo: tiempo de muestreo (seg/muestra) y las tres ecuaciones vistas en la seccin precedente ). Los comandos siguientes dibujan el lugar de races con las lneas de coeficiente de amortiguamiento y de
frecuencia natural constantes. Cree un nuevo archivo-m y entre los siguientes comandos. Luego de correr este archivo-m le debera dar el siguiente grfico de lugar de races.
numDz=[1 -0.3]; denDz=[1 -1.6 0.7];
rlocus (numDz,denDz)
axis ([-1 1 -1 1]) zeta=0.4; Wn=0.3; zgrid (zeta,Wn)
De esta figura, ya habr notado que el sistema es estable porque todos los polos se ubican dentro del crculo unitario. Asimismo, ve dos lneas punteadas de coeficiente de amortiguamiento y de frecuencia natural constantes. La frecuencia natural es mayor que 0.3 fuera de la lnea de Wn constante, y el coeficiente de amortiguamiento es mayor que 0.4 dentro de la lnea de zeta constante. En este ejemplo, ya tenemos dibujado el lugar de races en la regin deseada. Por lo tanto, una ganancia (K) elegida en la regin deseada le debera dar la respuesta que satisfaga los requerimientos de diseo.
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies: