Práctica de Laboratorio #5

PRÁCTICA DE LABORATORIO N° 5
A) ESTABILIDAD DE UN SISTEMA DE CONTROL DE LAZO CERRADO USANDO EL

LUGAR GEOMETRICO DE LAS RAICES
Objetivo de la Práctica
El objetivo de la práctica es obtener el Lugar de las raíces de un sistema dado por una función
de transferencia de lazo abierto y determinar si el mismo reúne de las condiciones de
estabilidad completa,
k
GH ( s ) =
¿¿
Sistema de la Práctica
El sistema que se analiza en la práctica es mostrado en la figura
Procedimiento
Se obtendrá el lugar de las raíces de un sistema cuya función de transferencia de lazo abierto
es:
k k
GH ( s ) = = 3
s( s+ 1)(s+2) s + 3 s2 +2 s
Se van a encontrar las raíces de la ecuación característica para varios valores de k (0.5, 2, 4,
6, 8), se encontrará que al aumentar la ganancia, las raíces dominantes se van acercando al
eje jω hasta que pasan al semiplano derecho del plano s.
Se va a encontrar la respuesta en el tiempo para cada valor de ganancia de lazo abierto

anteriormente dadas k (0.5, 2, 4, 6, 8) empleando Matlab. Escribirá el programa para cada caso
y simulará la respuesta al escalón.
Programa MATLAB A1
clear, close all, clc
% Lugar de las Raíces
% Se obtendrá el lugar de las Raíces de un sistema cuya función de transferencia
% de Lazo abierto es GH(s)=k/(s(s+1)(s+2))=k/(s^3+3s^2+2s)
num=[1];
den=[1 3 2 0];
% Introduzca el siguiente comando para obtener el lugar de las Raíces
rlocus(num,den)
% Introduzca la rejilla y él título de la gráfica
grid
title('Lugar de las Raíces de GH(s)=k/(s(s+1)(s+2))')
xlabel('EjeReal')
ylabel(‘Eje Imaginario')
% Se van a encontrar las Raíces de la Ecuación Característica para cuando k=0
roots(den)
% Con el siguiente comando se puede encontrar el valor de k y las raíces de la ecuación
% característica en cualquier punto del Lugar de las raíces.
[k,poles]=rlocfind(num,den)
% Éste comando se aplicará tantas veces como sea necesario.
Programa MATLAB A2
% Se van a encontrar las raíces de la ecuación característica para cuando k=0.5
%—La función de trasferencia de lazo abierto para k=0.5 es GH(s)=0.5/(s^3+3s^2+2)
num=[0.5];
den=[1 3 2 0];
% Introduzca el siguiente comando para encontrar el numerador y denominador de lazo
%—cerrado—
[numCL,denCL]=cloop(num,den)
%—Introduzca el siguiente comando para encontrar las raíces de ecuación característica
roots(denCL)
% Se va a encontrar la respuesta a una entrada escalón del sistema cuya función de
%—transferencia de lazo cerrado es C(s)/R(s)=0.5/(s^3+3s^2+2s+0.5)
%Introduzca el comando de la respuesta a un escalón step
%—Introduzca el tiempo t especificado por el usuario en el comando step
t=0:0.1:20;
step(numCL,denCL,t)
% Introduzca la rejilla y el título de la gráfica
grid
title('Respuesta en el tiempo a una entrada escalón de C(s)/R(s)=0.5/(s^3+3s^2+2s+0.5)')
xlabel('tiempo')
ylabel('salida')
Reporte.
1.- Para cada función de transferencia cuando k es (0.5, 2, 4, 6, 8) ¿dónde están ubicados los
polos de lazo cerrado del sistema en el plano complejo? Encuentre ξ y ωn para cada función de
transferencia.
2.- Obtenga la curva de respuesta en el tiempo para cada función de transferencia cuando la
entrada aplicada es un escalón unitario.
3.- Encuentre los valores de tiempo de pico (tp), tiempo de crecimiento (tr), tiempo de
establecimiento (ts 2%) y sobrepaso máximo porcentual (%Mp), para cada una de las curvas
de respuesta en el tiempo. Compare los valores anteriores con los obtenidos usando las
formulas obtenidas en clase. Si los resultados no coinciden explique la causa.
B) EFECTO DE LA ADICION DE POLOS O CEROS SOBRE EL LUGAR GEOMETRICO DE
LAS RAICES.
El objetivo de la práctica es notar que el efecto de la adición de los polos y ceros sobre el lugar
geométrico de las raíces.
Sistema de la Práctica
El sistema que se analiza en la práctica es un sistema con retroalimentación unitaria.
Procedimiento
Se obtendrá el lugar de las raíces de un sistema mostrado en la figura anterior cuya función de
transferencia de lazo abierto es:
k ( s+ 3) k (s +3)
GH ( s ) = = 3
s( s+ 1)(s+2) s + 3 s2 +2 s
Se van a encontrar las raíces de la ecuación característica para varios valores de k (0.5, 2, 4, 6,
8), se encontrará que al aumentar la ganancia, las raíces dominantes van aumentando su parte
imaginaria, jω, haciendo la respuesta más oscilatoria pero siempre estable.
Se va a encontrar la respuesta en el tiempo para cada valor de ganancia de lazo abierto

anteriormente dadas k (0.5, 2, 4, 6, 8). Para ello se tiene que encontrar la función de
transferencia de lazo cerrado y con ella encontrar la respuesta en el tiempo.
Programa MATLAB B1
% Lugar de las Raíces
% Se obtendrá el lugar de las Raíces de un sistema cuya función de transferencia
% de Lazo abierto es GH(s)=k(s+3)/(s(s+l)(s+2))=k(s+3)/(s^3+3s^2+2s)
num=[l 3];
den=[l 3 2 0];
% Introduzca el siguiente comando para obtener el lugar de las Raíces
rlocus(num,den)
% Introduzca la rejilla y el título de la gráfica
grid
title('Lugar de las Raíces de GH(s)=k(s+3)/(s(s+1)(s+2))')
xlabel('EjeReal')
ylabel(‘Eje Imaginario')
% Se van a encontrar las Raíces de la Ecuación Característica para cuando k=0
roots(den)
%—Con el siguiente comando se puede encontrar el valor de k y las raíces de la ecuación
%—característica en cualquier punto del Lugar de las raíces.
[k,poles]=rlocfind(num,den)
%—Este comando se aplicará tantas veces como sea necesario.
Programa MATLAB B2
% Se van a encontrar las raíces de la ecuación característica para cuando k=0.5
%-La función de transferencia de lazo abierto es GH(s)=(0.5s+1.5)/(s^3+3s^2+2)
num=[0.5 1.5];
den=[l 3 2 0];
%-—Introduzca el siguiente comando para encontrar el numerador y denominador de lazo
%—cerrado—
[numCL,denCL]=cloop(num,den)
%—Introduzca el siguiente comando para encontrar las raíces de ecuación característica
roots(denCL)
% Se va a encontrar la respuesta a una entrada escalón del sistema cuya función de
%—transferencia de lazo cerrado es C(s)/R(s)=0.5(s+3)/(s^3+3s^2+2.5s+1.5)
%—Introduzca el comando de la respuesta a un escalón step
%—Introduzca el tiempo t especificado por el usuario en el comando step
t=0:0.1:20;
step(numCL,denCL,t)
%— Introduzca la rejilla y el título de la gráfica
grid
title('Respuesta en el tiempo a una entrada escalón de C(s)/R(s)=0.5(s+3)/(s A3+3sA2+2.5s+1.5)')
xlabel('tiempo')
ylabel('salida')
Reporte.
1. Para cada función de transferencia ¿dónde están ubicados los polos de lazo cerrado del
sistema en el plano complejo? Encuentre ξ y ωn
2. ¿Qué efecto tuvo el agregar un cero en la posición de los polos de lazo cerrado para los
diferentes valores de ganancia?
3. Obtenga la curva de respuesta en el tiempo para cada función de transferencia cuando la
entrada aplicada es un escalón.
4. Encuentre los valores de tiempo de pico (tp), tiempo de crecimiento (tr), tiempo de
establecimiento (ts para el 2%) y sobrepaso máximo porcentual (Mp), para cada una de las
curvas de respuesta en el tiempo. Compare los valores anteriores con los obtenidos usando las
formulas obtenidas en clase. Si los resultados no coinciden explique la causa.
C) RESPUESTA A LA FRECUENCIA DE UN SISTEMA DE CONTROL (DIAGRAMA DE
BODE)
El objetivo de la práctica es encontrar el diagrama de Bode con MATLAB, los cuales son
utilizados para analizar y diseñar sistemas de control, aquí será usado para encontrar la
estabilidad de lazo cerrado. Para ello se utiliza la función de transferencia de lazo abierto:
k
GH ( s ) =
¿¿
Sistema de la práctica.
El sistema que se analiza en la práctica es un sistema con retroalimentación unitaria
Se encontrarán los diagramas de Bode para cuando k=20 y para cuando k=200
Procedimiento.
Se obtendrá el diagrama de Bode del sistema cuya función de transferencia de lazo abierto es:
GH(s)=k/(s(s+1)(s+10))
Se encontrarán los diagramas de Bode para cuando k=20 y para cuando k=200,
Con los diagramas de Bode se obtendrá el Margen de fase (MF) y frecuencia de cruce de
ganancia (ωg), el Margen de ganancia (MG) y la frecuencia de cruce de fase (ωp), obteniendo
con estos valores la estabilidad de lazo cerrado del sistema.
También se encontrará la respuesta en el tiempo para los valores de ganancia de lazo abierto
como una forma de comprobar la estabilidad del sistema.
Programa MATLAB C1
% Se va a encontrar el Diagrama de Bode del sistema cuya GH(s)= k/(s(s+1)(s+10))
% o GH(s)= k/(s^3+11s^2+10s) para cuando k=20
num=[20];
den=[1 11 10 0];
% Para encontrar el Diagrama de Bode usamos los comandos
bode(num,den)
title('Diagrama de Bode de un sistema cuya GH(s)=20/(s^3+11s^2+10s) ')
% Se va a encontrar la Estabilidad del sistema, encontrando el MF y ωg y MG y ωp
% para ello se utiliza el siguiente comando.
margin(num,den)
Programa MATLAB C2
% se ya encontrar la respuesta en tiempo para una entrada escalón, del sistema
% cuya GH(s)= k/(s(s+1)(s+10)) = k/(s^3+11s^2+10s) para cuando H(s)=l y k=20
num=[20];
den=[1 11 10 0];
% para encontrar la función de transferencia de lazo cerrado usamos los comandos
[numCL,denCL]= cloop(num,den)
% Para encontrar la Respuesta en el Tiempo usamos los comandos
step(numCL, denCL)
grid
title('Respuesta a un Escalón de C(s)/R(s)=20/(s^3+11s^2+10s+20)')
xlabel('Tiempo')
ylabel('Salida')
Reporte.
1. Dibuje los diagramas de Bode para cuando la ganancia k= 20 y k=200.
2. Obtenga MF, ωg, MG y ωp para cada una de las funciones de transferencia, para cuando
k=20 y k=200, con estos valores obtenga la estabilidad del sistema.
3. ¿Cuál es la diferencia de los diagramas de Bode para cuando k=20 y k=200, en las gráficas
de magnitud y de ángulo?
4. Obtenga las gráficas de respuesta en el tiempo de lazo cerrado para cuando la entrada
aplicada es un escalón, para cuando la k=20 y k=200.
5. Defina los siguientes comandos: num, den, cloop, step, grid, title, xlabel, ylabel, margin, grid.
6. ¿Si el sistema es inestable como se encuentra el MF y el MG del sistema?
D) RESPUESTA A LA FRECUENCIA DE UN SISTEMA DE LAZO CERRADO.
El objetivo de la práctica es encontrar la respuesta en frecuencia de lazo cerrado con MATLAB,
la cual será usada para encontrar los valores de pico de resonancia (Mr) y frecuencia de
resonancia (wr), y ancho de banda (Bw).
Sistema de la Práctica.
El sistema que se analiza en la práctica es un sistema con retroalimentación unitaria cuya
función de transferencia de lazo abierto es: GH(s)=k/(s(s+1)(s+10)), el diagrama de bloques del
sistema aparece mostrado en la gráfica.
Se encontrarán los diagramas de Bode de lazo cerrado para cuando k=20 y para cuando
k=200.
Procedimiento.
Se obtendrá el diagrama de Bode del sistema cuya función de transferencia de lazo cerrado es:
C(s)/R(s) =k/((s(s+1)(s+10)+k) =k/(s3+11s2+10s+k).
Se encontrarán los diagramas de Bode para cuando k=20 y para cuando k=200. También se
calculará el pico de resonancia, la frecuencia de resonancia y el ancho de Banda.
Programa MATLAB D1
% Se va a encontrar el Diagrama de Bode del sistema cuya:
% C(s)/R(s)=k/(s(s+l)(s+10)+k) o C(s)/R(s)= 20/(s^3+11s^2+10s+20) para cuando
% k=20
% Para encontrar el Diagrama de Bode usamos los comandos
num=[20];
den=[1 11 10 20];
bode(num,den)
grid
title('Diagrama de Bode de un sistema cuya C(s)/R(s)=20/(s^3+11s^2+10s+20) ')
% Las órdenes de MATLAB que se puede utilizar para encontrar, el pico de resonancia
% y la frecuencia de resonancia son las siguientes:
w=logspace(-1,1);
[mag,phase,w]=bode(num,den,w); % o [mag,phase,w]=bode(sys,w)
[Mp,k]=max(mag)
Mr=20*log10(Mp);
wr=w(k)
% El ancho de banda se obtiene introduciendo las siguientes líneas en el programa
n=1;
while 20*log10(mag(n))>=-3;n=n+1
end
BW=w(n)
Programa MATLAB D2
% Se va encontrar la respuesta en tiempo para una entrada escalón, del sistema
% cuya C(s)/R(s)= k/(s(s+l)(s+10)+k) o C(s)/R(s)= 20/(s^3+11s^2+10s+20) para
% cuando k=20
num=[20];
den=[1 11 10 20];
% para encontrar la Respuesta en el Tiempo usamos los comandos
step(num,den)
grid
title('Respuesta a un Escalón de C(s)/R(s)=20/(s^3+11s^2+10s+20)')
xlabel('Tiempo')
ylabel('Salida')
Reporte.
1. Dibuje los diagramas de Bode de lazo cerrado para cuando la ganancia es: k= 20 y k=200.
2. Obtenga el pico de resonancia, Mr, la frecuencia de resonanica wr, y el ancho de banda BW
para cada una de las funciones de transferencia de lazo cerrado, para cuando k=20 y k=200.
3. ¿Cuál es la diferencia de los diagramas de Bode para cuando k=20 y k=200 en las gráficas
de magnitud y en las gráficas de ángulo?
4. Obtenga las gráficas de respuesta en el tiempo de lazo cerrado para cuando la entrada
aplicada es un escalón unitario, para cuando la k=20 y k=200.
5. Compare los valores de Mr y de (1+Mp) ¿por qué se hace esta comparación?
TAREA CORRESPONDIENTE A UN EXAMEN PROCEDIMENTAL.

1. Presentar un informe conteniendo las respuestas o reportes de cada uno de los temas
plateados.
2. En el informe incluir los códigos de los programas elaborados y los resultados tanto gráficos
como textuales.
3. Presentarlo en un folder guardando la limpieza y las reglas gramaticales. Incluya si así lo
estima, un breve sustento teórico de cada Tema planteado.
4. La presentación del informe será el día y hora antes del inicio del segundo examen parcial.
No hay prorrogas.

Práctica de Laboratorio #5

Cargado por

Copyright:

Formatos disponibles

Práctica de Laboratorio #5

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Práctica de Laboratorio #5

Cargado por

Copyright:

Formatos disponibles

PRÁCTICA DE LABORATORIO N° 5

A) ESTABILIDAD DE UN SISTEMA DE CONTROL DE LAZO CERRADO USANDO EL

Se va a encontrar la respuesta en el tiempo para cada valor de ganancia de lazo abierto

Se va a encontrar la respuesta en el tiempo para cada valor de ganancia de lazo abierto

TAREA CORRESPONDIENTE A UN EXAMEN PROCEDIMENTAL.

También podría gustarte