Compensadores

INSTITUTO POLITÉCNICO NACIONAL

ESCUELA SUPERIOR DE INGENIERÍA MECÁNICA Y ELÉCTRICA

UNIDAD ZACATENCO

INGENIERÍA EN CONTROL Y AUTOMATIZACIÓN

TEORIA DEL CONTROL II

PROYECTO FINAL

DISEÑO DE COMPENSADORES

PROFESOR:

NORMA LOZADA C.

ALUMNO:

DÍAZ CAMACHO RICARDO DANIEL GRUPO 6AV1

CORREO: ricardodiazc@hotmail.com

OBJETIVO:

Diseñar los compensadores de la siguiente función realizando los cálculos y comparando con matlab las graficas obtenidas de la función compensada

DESCRIPCIÓN DEL PROCESO:

Realizar el diseño de los compensadores de Adelanto, Atraso y Atraso-Adelanto de la siguiente función

G(s)=1/(s(s+2))

DESARROLLO:

Obtenemos el LGR de la función de transferencia sin compensador

LGR DE LA FUNCION DE LA PLANTA

CODIGO MATLAB PARA LGR DE LA FUNCION SIN COMPENSACION

num=[1];

den=[1 2 0];

roots(den)

rlocus(num,den)

grid on;

Anexo los cálculos del LGR realizados previamente a mano.

COMPENSADOR DE ADELANTO

Elaboramos los cálculos para obtener el compensador de adelanto de la función dada anteriormente.

Y el procedimiento del mismo en matlab, lo cual nos arroja las graficas siguientes

CODIGO MATLAB PARA COMPENSADOR DE ADELANTO

%Compensador de Adelanto

num=[1];

den=conv([1 0],[1 2]);

g_s=tf(num,den)

%muestra los polos del sistema

roots(den)

p1=ans(1);

p2=ans(2);

%condiciones iniciales

tp=1;

Mp=3;

%Calculando frecuencia amortiguada (wd) y la ubicacion del polo (sp)

wd=pi/tp;

sp=complex(-tp,wd);

%Calculando factor de amortiguamiento

E=log(Mp/100);

E=-pi/E;

E=E^2+1;

E=sqrt(1/E);

%Se establece factor de amortiguamiento

E=0.8

%Se calcula la frecuencia natural

wn=wd/sqrt(1-E^2);

Ewn=wn*E;

%Se establece el polo deseado

sd=complex(-Ewn,wd)

%Se desea obtener el angulo de los vectores al polo deseado

%(teta t y alfa a, asi como la distancia del cero x y del polo y

t1=angle(sd);

t1=(t1*180)/pi

t2=angle(sd+2);

t2=(t2*180)/pi

%La sumatoria de angulos de ceros menos la sumatoria de angulos

%de polos es =-180

%El angulo alfa puede tomar valores de 87.85 hasta antes de O2 ya que si

%toma el mismo valor de O2 el polo y el cero se eliminan, por lo tanto se

%le asigna a alfa (a) valor de 89

a=89;

t3=180+a-t1-t2

%Se encuentran las distancias del polo y el cero del compensador

a1=(180*tan(a))/pi;

x=wd/a1;

x=x+Ewn

t3a=(180*tan(t3))/pi

y=wd/t3a

y=y+Ewn

y=161.55

k=[1003.6107]

%Se determinan el polo y el cero del compensador

p1c=conv(k,[1 x])

p2c=[1 y]

%Obtenemos la funcion de transferencia con el compensador y su lugar de las

%raices compensado

c_s=tf(p1c,p2c)

ftgc=g_s*c_s

numsc=p1c

densc=conv (den,p2c)

rlocus(numsc,densc)

grid on;

Anexo cálculos previamente hechos del Compensador de Adelanto.

COMPENSADOR DE RETARDO

CODIGO MATLAB PARA COMPENSADOR DE RETARDO

%Compensador de Atraso

num=[1];

...