Selección de Parámetros PID
Selección de Parámetros PID
Selección de Parámetros PID
Para calcular los parmetros proporcional, integral y derivativo de un controlador PID por
el mtodo de Ziegler Nichols, se requiere encontrar la respuesta del sistema en lazo
cerrado al escaln, para la ganancia proporcional crtica (Kcr), donde la respuesta presenta
oscilaciones con amplitud constante.
Para determinar el valor de Kcr, se recomiendan dos mtodos, por medio del criterio de
estabilidad de Routh Hurtwitz o por el mtodo del lugar de las races. Despus de
determinar el valor de Kcr y medir el periodo de oscilacin de la respuesta para esa
ganancia (Tcr), se determinan los valores de las ganancias del controlador PID de la
siguiente manera:
Tipo de PID Kp Ti Td
P 0.5*Kcr Inf 0
PI 0.45*Kcr Tcr/1.2 0
PID 0.6*Kcr Tcr/2 Tcr/8
Tabla 1. Configuracin de los controladores P, PI y PID.
Ejemplo:
1
Gp(s) =
s + 9s + 23s + 15
s + 9s + 23s + 15 + K = 0
A partir de los coeficientes del polinomio se obtiene el siguiente arreglo de Routh Hurwitz:
s 1 23 0
s 9 15+K 0
s 192-K 0
s 15+K
Para asegurar la estabilidad, la segunda columna del arreglo debe ser positiva, entonces se
obtienen las siguientes restricciones para el valor de K:
15 + K > 0
192 K > 0
Entonces K debe ser, K>-15 y K<192. El valor critico se encuentra cuando K=192, entonces
Kcr=192.
Por medio del lugar de las races se puede determinar un valor cercano al valor critico, se
debe buscar un valor de la ganancia K que haga que la parte real de los polos sea igual a
cero. Por ejemplo en la figura 2 se puede observar dicho punto,
Paso 2: Calcular el valor de Tcr. Se puede realizar de dos formas, una partir del lugar de las
races tomar el valor de la frecuencia para el valor de Kcr; en la figura 2 se observa que
este valor es = 4.79, para calcular el Tcr se emplea la siguiente ecuacin:
2
Tcr =
Entonces Tcr=1.31 segundos. Este tiempo tambin se puede medir, a partir de dos picos
consecutivos en la respuesta en lazo cerrado del sistema para la ganancia Kcr por medio
de la funcin ginput. En la figura 3 se puede observar cmo se mide el tiempo Tcr.
1.6
1.4
1.2
1
Amplitude
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Time (sec)
A partir de los valores Kcr y Tcr, se puede disear el controlador. A continuacin vemos las
respuestas al escaln para el sistemas en lazo cerrado, con el controlador proporcional
(ver figura 4), el controlador PI (ver figura 5) y el PID (ver figura 6). Recordar que la funcin
de transferencia de cada controlador es:
1
Cs = K 1 + + T s
T s
num = 1;
den = [1 9 23 15];
planta = tf(num,den);
Kcr = 192;
Tcr=1.31;
s=tf([1 0],1);
%Controlador P
figure(1)
planta_p = feedback(0.5*Kcr*planta,1);
step(planta_p)
title('Controlador P')
%Controlador PI
figure(2)
Kp=0.45*Kcr;
Ti=Tcr/1.2;
ftpi=Kp*(1+1/(Ti*s));
planta_pi = feedback(ftpi*planta,1);
step(planta_pi)
title('Controlador PI')
%Controlador PID
figure(3)
Kp=0.6*Kcr;
Ti=Tcr/2;
Td=Tcr/8;
ftpid=Kp*(1+1/(Ti*s)+Td*s);
planta_pi = feedback(ftpid*planta,1);
step(planta_pi)
title('Controlador PID')
Controlador P
1.4
1.2
0.8
Amplitude
0.6
0.4
0.2
0
0 2 4 6 8 10 12
Time (sec)
Figura 4. Respuesta al escaln del sistema en lazo cerrado con controlador proporcional.
Controlador PI
1.8
1.6
1.4
1.2
1
Amplitude
0.8
0.6
0.4
0.2
0
0 5 10 15 20 25
Time (sec)
Figura 5. Respuesta al escaln del sistema en lazo cerrado con controlador PI.
Controlador PID
1.5
1
Amplitude
0.5
0
0 1 2 3 4 5 6 7
Time (sec)
Figura 6. Respuesta al escaln del sistema en lazo cerrado con controlador PID.
Ejercicios:
Calcular los parmetros de los controladores P, PI y PID para las siguientes funciones de
transferencia:
1.
1
Gp(s) =
s + 3s + 11s + 12
2.
1
Gp(s) =
s + 9s + 23s + 15
3.
8
Gp(s) =
s + 7s + 14s + 8