Sistemas de Control I

Prácticas de Sistemas
de Control
Roberto Martínez Montejano

Unidad académica Multidisciplinaria Zona Media
Isaac Campos Cantón

Facultad de Ciencias
Raúl Castillo Meraz

Coordinación Académica Región Huasteca Sur
Tecnología y ciencias aplicadas

.
© Roberto Martínez Montejano, Isaac Campos Cantón, Raúl Castillo Meraz.

D. R. © Universidad Autónoma de San Luis Potosí 2019.
Edición a cargo Roberto Martínez Montejano, Isaac Campos Cantón.
Diseño editorial: M. C. Juan Alberto Vértiz Hernández, Dr. Isaac Campos Cantón.
ISBN: 978-607-535-122-3 (e-Libro)
San Luis Potosí, México.
Todos los derechos reservados. Esta obra no puede ser reproducida en todo o en parte, ni
registrada o transmitida por un sistema de recuperación, en ninguna forma y medio, sea
mecánico, fotoquímico, electrónico, magnético, electroóptico por fotocopia o cualquier otro,
sin permiso previo del titular de los derechos patrimoniales.
Índice general
1. Introducción a Matlab 3
1.1. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4. Reporte de resultados . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.5. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2. Uso básico de Simulink 7

2.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3. Modelado matemático de sistemas fı́sicos:parte 1 15

3.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.1.1. Sistemas mecánicos traslacionales. . . . . . . . . . . . . . . . . . 15
3.1.2. Sistemas mecánicos rotacionales. . . . . . . . . . . . . . . . . . 17
3.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
3.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
4. Modelado matemático de sistemas fı́sicos: parte II 21

4.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
4.1.1. Sistemas eléctricos. . . . . . . . . . . . . . . . . . . . . . . . . . 21
4.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
i
ii
5. Diseño de controladores analógicos por compensación de fase. 27
5.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
5.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
5.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
5.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
5.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
6. Identificación de procesos de primer orden 31

6.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
6.1.1. Método gráfico para identificar un sistema de primer orden . . . 32
6.1.2. Segundo método gráfico para identificación de parámetros de pro-
cesos de primer orden . . . . . . . . . . . . . . . . . . . . . . . . 33
6.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
6.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
6.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
6.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
7. Sintonización de controladores PID por método experimental parte

I. 37
7.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
7.1.1. Método de Ziegler and Nichols. . . . . . . . . . . . . . . . . . . 39
7.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
7.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
7.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
7.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
8. Sintonización de controladores PID por método experimental parte

II. 43
8.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
8.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
8.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
8.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
8.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
9. Sintonización de controladores PID forma analı́tica parte I. 47

9.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
9.1.1. Ejemplo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
9.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
9.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
iii
9.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
9.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
10.Sintonización de controladores PID forma analı́tica parte II. 53

10.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
10.1.1. Ejemplo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
10.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
10.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
10.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
10.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
11.Creación de planta analógica de segundo orden. 59

11.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
11.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
11.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
11.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
11.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
12.Implementación de un PID analógico. 65

12.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
12.1.1. Etapa proporcional. . . . . . . . . . . . . . . . . . . . . . . . . . 65
12.1.2. Etapa integral. . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
12.1.3. Etapa derivadora. . . . . . . . . . . . . . . . . . . . . . . . . . . 66
12.2. Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
12.3. Material . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
12.4. Procedimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
12.6. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
Introducción
Los sistemas de control son de mucha utilidad en el campo de la investigación y en

el industrial. El conocimiento de esta disciplina ahorra mucho en tiempo y en capi-
tal monetario a los ingenieros e investigadores. El entrenamiento para comprender las
técnicas empleadas en los sistemas de control debe realizarce a través de experimentos
en un ambiente controlado, para ello, el laboratorio es el lugar propicio. Ası́, a través
de prácticas diseñadas y con un objetivo a cubrir se generará la experiencia; y como
concecuencia, cuando se presente un desafı́o en el campo de trabajo, este se abordará
de la mejor forma esperando obtener los resulatdos satisfactorios.
Estas prácticas de control abarca lo que se ha denominado teorı́a de control clásico, se

han formulado 12 prácticas que a continuación se describen:
Práctica 1 Se conocerá el software Matlab a través de aplicar instrucciones básicas y
el concepto de sistema de primer orden.
Práctica 2 Se desarrollarán bloques básicos para la construcción de modelos en el
entorno de programación visual Simulink de Matlab, además de ejercitarse con
sistemas de segundo orden.
Práctica 3 Se obtedrá el modelado matemático de sistemas mecánicos rotacionales
y traslacionales, su función de transferencia y su espacio de estados. Además
de obtener las respuestas al impulso y al escalón de dichos sistemas usando el
software Matlab.
Práctica 4 Se estudiará el modelado matemático de sistemas eléctricos obteniendo la
función de transferencia y su espacio de estados. Además de obtener la respuesta
al impulso y al escalón unitario de dichos sistemas a través del software Matlab.
Práctica 5 Se diseñará un controlador por compensación de fase utilizando un com-
pensador por adelanto de fase. Se obtendrán las caracterı́sticas de la planta en
lazo abierto con base a las especificaciones dadas y se realizará el cálculo para las
ganancias del compensador.
Práctica 6 Se identificarán los parámetros básicos en un sistema de primer orden con
retardo.
Práctica 7 Se diseñará un controlador por método gráfico de Ziegler-Nichols para
sintonizar los parámetros de control en sistemas de primer orden.
Práctica 8 Se diseñará un controlador por el método última ganancia para sintonizar
los parámetros de control.
1
2
Práctica 9 Se diseñará un controlador usando el criterio de Routh-Hurwitz para sin-
tonizar los parámetros de control de forma analı́tica.
Práctica 10 Se diseñará un controlador por la modificación de polos a través del
análisis matemático y se realizará la sintonización de sus parámetros de control.
Práctica 11 Se construirá una planta de segundo orden utilizando amplificadores ope-
racionales, se identificarán sus parámetros y se variará la constante de amorti-
guamiento.
Práctica 12 Se implementará un controlador PID utilizando amplificadores operacio-
nales y su sintonización para un sistema de segundo orden.
De esta forma se aborda un conjunto de prácticas para la teorı́a de los sistemas de
control clásico. A continuación se describe de forma detallada cada práctica.
Práctica 1
Introducción a Matlab
El modelado matemático de un sistema dinámico, es el conjunto de ecuaciones que

representan la evolución del mismo con bastante precisión o en el grado más cercano
a ello. Dicho modelo matemático no es único para un sistema determinado, ya que
puede representarse de diversas formas, por lo que se pueden tener diferentes modelos
matemáticos dependiendo de la perspectiva a analizar. La dinámica de los sistemas ge-
neralmente se representa a través de ecuaciones diferenciales, que se obtienen mediante
la aplicación de leyes fı́sicas como leyes de Newton o leyes de Kirchhoff, por mencionar
algunas.
Un sistema que se denomina lineal e invariante en el tiempo, es aquel en el que se

puede aplicar el principio de superposición y además tiene coeficientes constantes o
son funciones de una variable independiente.
En un sistema lineal invariante en el tiempo, una función de transferencia se define

como el cociente entre la transformada de Laplace de la señal de salida y la transfor-
mada de Laplace de la señal de entrada. El sistema viene dado por la siguiente ecuación
diferencial:
an y n (t)+an−1 y n−1 (t)+..+a1 ẏ(t)+a0 y(t) = bm um (t)+bm−1 um−1 (t)+..+b1 u̇(t)+b0 u(t)
(1.1)
En donde u(t) representa la entrada y y(t) la salida. La función de transferencia que

relaciona la salida con la entrada de la ecuación diferencial (1.1) es:
Y (s) bm sm + bm−1 sm−1 + .. + b1 s + b0

G(s) = = (1.2)
U (s) an sn + an−1 ss−1 + .. + a1 s + a0
A partir de este concepto, es posible representar la dinámica de un sistema mediante

ecuaciones algebraicas en el dominio s, siendo la potencia más alta en el denominador
la que originará el orden del sistema, esto es, si la función es de orden n, el sistema
3
4
será de n-ésimo orden.
Para poder conocer la respuesta del sistema, se utilizan diversas entradas prototipo,
dentro de las cuales, las más comunes son el impulso y el escalón. La salida del sistema
con esas entradas, arrojará información acerca de la estabilidad, rapidez del sistema,
sobretiro, tiempo de asentamiento, etc. El análisis de estas respuestas, se conoce como
análisis de respuesta transitoria y estacionaria. El software Matlab se puede utilizar
para conocer y analizar la respuesta del sistema dinámico modelado matemáticamen-
te, en donde se puede encontrar el conjunto de instrucciones que se detalla en la tabla 1.
1.1. Objetivo
Conocer y aplicar las intrucciones básicas de Mtalab para obtener la respuesta de un
proceso ante diversas funciones de entrada, además de estudiar los parámetros de un
sistema de primer orden.
1.2. Material
1. Materiales utilizados
Lápiz y papel
2. Equipo utilizado
Equipo de cómputo con software Matlab
1.3. Procedimiento
1. Capture las siguientes funciones de transferencia en MATLAB:
2s + 25
G1 (s) = (1.3)
12s2 + 115s + 245
315
G2 (s) = (1.4)
17s2
+ 115s + 350
2. Obtenga la expansión en fracciones parciales de las funciones G1 y G2 .
3. Dado que G(s)=Y(s)/R(s), suponga que r(t)=δ(t) (δ(t) = impulso). Obtenga la
respuesta en el tiempo y(t) si G(s)=G1 (s)*G2 (s).
Utilizando la transformada inversa de Laplace.
Utilizando la función impulse.
4. Obtenga la respuesta en el tiempo a una entrada escalón unitario a la planta
G(s)=G1 (s)*G2 (s):
Utilizando la transformada inversa de Laplace.
5
Tabla 1.1: Ejemplos de comandos en Matlab

Sintaxis Descripción Ejemplo
B=[b1 b2...bn] Asigna un vector 1×n (con B1=[2 30]; A1=[1 25 174
elementos b1, b2,...,bn) a la 360];
variable B. El polinomio es
un vector de coeficientes de
potencias descendentes.
tf(B,A) Define una función de trans- G1=tf(B1,A1)
ferencia
conv(p,q) Multiplica los polinomios MB2=conv([2],[-1.5 1];
representados por los vecto-
res p y q.
[R,P,K]=residue(B,A) Calcula los residuos, polos y [R1,P1,K1]=residue(B1,A1);
término directo de una ex-
pansión en fracciones par-
ciales de la relación de po-
linomios B(s)/A(s)
syms x Construye objetos simbóli- syms s;
cos (x es una variable a la
que no se asigna un valor en
particular).
ilaplace(arg) Obtiene la transformada in- ilaplace((3*s+5)/(s∧ 2+2));
versa de Laplace de la fun-
ción que se encuentra en el
argumento.
impulse(g) Obtiene la gráfica de la res- impulse(G1);
puesta a una entrada impul-
so de un sistema dinámico,
donde g es una función de
transferencia.
[Y,t]=step(g) Obtiene la respuesta ante [Y1,t1]=step(G1);
una entrada escalón de una
función de transferencia g.
Y es el vector de la respues-
ta y t el del tiempo.
Grid Cuadricula el área de la Grid
gráfica.
xlabel(’texto’) Coloca el nombre del eje x. xlabel(’tiempo(s)’)
ylabel(’texto’) Coloca el nombre del eje y. ylabel(’Respuesta’)
title(’texto’) Coloca el tı́tulo de la gráfi- title(’Respuesta G1 al im-
ca. pulso’)
hold on Retiene la curva para añadir hold on
otra sobre la misma gráfica.
6
Utilizando la función step.
5. Añadir un retraso de 2.5 segundos a la función de transferencia G(s). Obtenga la
respuesta al escalón unitario. Observación: use la inversa de Laplase.
6. Defina cuatro sistemas de primer orden (usted deberá de otorgar los valores ne-
cesarios a k4 , k5 , τ4 , τ6 , θ7 ). Obtenga la respuesta al escalón unitario de los cuatro
sistemas en una misma gráfica y concluya el efecto de los parámetros.
k4
G(3) = (1.5)
τ4 s + 1
k5
G(3) = (1.6)
τ4 s + 1
k5
G(3) = (1.7)
τ6 s + 1
k5 e−θ7 s
G(3) = (1.8)
τ6 s + 1
1.4. Reporte de resultados

En la fecha de entraga pactada, el alumno deberá mostrar la simulación realizada,
además de los siguientes puntos:
Se deberá de describir todo el procedimiento utilizado en los puntos mencionados
en la sección anterior.
En el punto 2, se debe incluir en el reporte la expansión en fracciones parcia-
les no como una captura de pantalla del resultado arrojado por Matlab, sino
reescribiéndolo en el editor de ecuaciones de Word.
De los puntos 3, 4, 5 y 6 se debe de incluir el procedimiento, las ecuaciones
obtenidas y las gráficas.
Se deberá entregar un reporte que contenga lo anteriormente solicitado, con las siguien-
tes secciones:
1. Portada
2. Introducción (máximo una cuartilla)
3. Desarrollo
4. Resultados
5. Discusión y análisis de los resultados
6. Referencias
1.5. Conclusiones
Redacte de forma concisa las conclusiones acerca de las actividades realizadas y del
alcance de los objetivos planteados en esta práctica.
Práctica 2
Uso básico de Simulink
2.1. Introducción
Un diagrama de bloques de un sistema, es una representación gráfica de las funciones
que tiene cada componente y flujo de señales, mostrando las relaciones existentes entre
los diversos componentes. Un diagrama de bloques tiene la ventaja de mostrar el flujo
de las señales de una forma directa. En estos diagramas, las variables se enlazan unas
con otras mediante bloques funcionales, que incluyen un sı́mbolo para la representación
matemática que realiza ese bloque sobre la señal de entrada, para producir una cierta
salida.
Las funciones de transferencia de los bloques se escribe dentro de estos, interconec-

tandose entre sı́ mediante flechas, que indican la dirección del flujo de las señales. Un
diagrama de bloques muestra una propiedad explı́cita unilateral. En la figura 2.1 se
muestra el esquema de un diagrama de bloques.
Figura 2.1: Bloque donde se especifica la función de transferencia.
Un diagrama de bloques contiene información relacionada con el comportamiento

dinámico, pero no incluye la información de la construcción fı́sica del sistema, es por
ello, que un mismo sistema puede ser representado por diferentes diagramas de bloques.
7
8
Para poder simular los diagramas de bloques, se puede utilizar el entorno Simulink
de Matlab, el cual es un entorno de programación visual en bloques de alto nivel. Se
utiliza para realizar pruebas de forma anticipada y eficiente, que incluye diversos en-
tornos como comunicaciones inalámbricas, electrónica de potencia, sistemas de control,
procesamiento de señales, robótica, procesamiento de imágenes, entre otros.
Dentro de los bloques básicos que se pueden encontrar en Simulink, se detallan algunos
en la Tabla 2.1.
Se pueden modificar los parámetros de la simulación en Simulink, para ajustarlo a

las condiciones requeridas, como se muestra en la Tabla 2.2.
En la figura 2.2 se muestran los parámetros descritos en la Tabla 2.2 como apare-
cen en Simulink.
Figura 2.2: Configuración de parámetros de Simulink.

9
Tabla 2.1: Ejemplos de bloques en Simulink

Bloque Descripción Librerı́a
Genera una señal de tipo escalón, en Sources

la que puede configurarse el tiempo de
aplicación, valor inicial y valor final.
Define una función de transferencia Continuous

donde el numerador y denominador
están dados por dos vectores de coefi-
cientes de potencias descendentes de s.
Simula el tiempo muerto del proceso. Continuous
Multiplexa las señales de entrada en Signal Routing

una señal de salida, donde se puede
configurar el número de señales de en-
trada
Genera una gráfica de la simulación con Sinks

la señal de entrada que puede contener
información de una o varias variables.
Permite transferir los resultados de la Sinks

simulación al espacio de trabajo de
Matlab. Puede configurarse el nombre
de la variable y su tipo.
Suma señales. El número y signo de las Math Opera-

señales a sumar puede configurarse. tions
10
Tabla 2.2: Parámetros de simulación en Simulink

Parámetro Descripción Valor tı́pico
Start time Tiempo de inicio de la simulación. Cero
Stop time Tiempo final de la simulación. Ajustable pa-
ra observar el
transitorio
Solver Se especifica si el tamaño de paso para Variable
options- la simulación es fijo o variable.
Type
Solver op- Método numérico que el simulador uti- Ode45
tions lizará para la solución de ecuaciones di-
ferenciales ordinarias
Finalmente, en la tabla 2.3 se describen los comandos que se pueden utilizar al mo-
mento de exportar una señal mediante simout a ventana de comandos de Matlab.
Tabla 2.3: Instrucciones de Matlab

Sintaxis Descripción Ejemplo
plot(x,y) Grafica la curva y(x). plot(tout,simout);
Arreglo(:,num) Se hace referencia a los ele- simout(:,1);
mentos de todas las filas de
la columna num.
plot(x,y,’g’) Grafica la curva y(x) en co- plot(tout,simout(:,1),’g’);
lor verde.
hold on Permite añadir curvas al hold on
mismo espacio de graficado
gtext(’letrero’) Coloca letrero donde se po- gtext(’G1’)
sicione el cursor sobre la
gráfica
hold off Deja de añadir curvas sobre hold off
la misma gráfica.
grid Cuadricula el área de grafi- grid
cado.
title(’tı́tulo’) Coloca tı́tulo especificado title(’Diferentes tipos de sis-
en la gráfica. temas de segundo orden’)
11
2.2. Objetivo
Conocer y aplicar bloques básicos para la construcción de modelos en Simulink además
de distinguir distintos tipos de sistemas de segundo orden.
2.3. Material
Lápiz y papel
2. Equipo utilizado
2.4. Procedimiento
1. Construya un diagrama de bloques en Simulink para la implementación de un
lazo de control; formado por un proceso y controlador, dados por Gp (s), Gd (s) y
Gc (s), como se muestra en la figura 2.3.
Figura 2.3: Diagrama de bloques de lazo de control retroalimentado.
En donde:
k 2
Gp (s) = = (2.1)
τs + 1 7s + 1
kd −0.5
Gd (s) = = (2.2)
τd s + 1 10s + 1
1.75s + 1.75
Gc (s) = (2.3)
s
La figura 2.4 muestra la implementación en Simulink del sistema de control mos-

trado previamente.
12
Figura 2.4: Implementación de lazo de control en Simulink.
2. Simule un cambio de referencia R(t) de 0 a 1 iniciando en t=1, además de la

presencia de una perturbación D(t) de 0 a 1 en t=5. El tiempo total de simulación
del sistema de control deberá ser de 30.
Grafique los resultados de simulación en Matlab
3. Construya el modelo descrito en la figura 2.5, en Simulink.
Figura 2.5: Diagramas de Simulink para comparar sistemas de segundo orden.
4. Utilice el diagrama de Simulink descrito en el punto anterior para comparar la

respuesta de los tres tipos de sistemas de segundo orden ante una entrada escalón.
13
5. Calcule el discriminante e identifique cada respuesta (subamortiguada, sobre-
amortiguada y crı́ticamente amortiguada) e identifique cada sistema del inciso
anterior en la gráfica del ”scope”de Simulink.
6. Grafique los resultados de la simulación en Matlab, utilizando el siguiente código:

plot(tout,simout(:,1),tout,simout(:,2),tout,simout(:,3))
grid
hold
xlabel(’tiempo’)
ylabel(’Respuesta transitoria’)
gtext(’Y 1’,’fontzise’,12)
7. El diagrama de simulink de la figura 2.6, muestra la manera de incluir un tiem-

po muerto a la respuesta del sistema. Compare la respuesta de los sistemas con
tiempo muerto y sin tiempo muerto.
Figura 2.6: Sistema en Simulink con tiempo muerto.
8. Construya las siguientes funciones de transferencia y defina los parámetros. Mues-

tre el efecto de cada parámetro.
ke−θ1 s
G1 (s) = (2.4)
(τ 2 )s2 + 2(τ )(ζ1 )s + 1
ke−θ1 s
G1 (s) = 2 2 (2.5)
(τ )s + 2(τ )(ζ2 )s + 1
14
ke−θ2 s
G1 (s) = (2.6)
(τ 2 )s2 + 2(τ )(ζ1 )s + 1

Se deberá de describir todo el procedimiento utilizado en los puntos mencionados
en la sección anterior.
Se deben incluir los cálculos realizados en el punto 5 para la obtención del dis-
criminente.
Se deberá entregar un reporte que contenga lo solicitado anteriormente, con las siguien-
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
2.6. Conclusiones
Práctica 3
Modelado matemático de sistemas

fı́sicos:parte 1
3.1. Introducción
Los modelos matemáticos pueden adoptar diferentes formas dependiendo del sistema
que se trate y dependiendo de las circunstancias especı́ficas del mismo, un modelo pue-
de ser más conveniente que otro, por ejemplo en control óptimo, es provechoso utilizar
representaciones en espacio de estados, en cambio para el análisis de la respuesta transi-
toria, es mejor utilizar la función de transferencia. Teniendo dicho modelo matemático
se pueden utilizar diversos recursos como algún software para analizarlo y sintetizarlo.
Dependiendo del sistema fı́sico a tratar, serán las leyes fı́sicas aplicadas al mismo, te-
niendo sistemas eléctricos, mecánicos, electrónicos, hidráulicos, neumáticos, entre otros.
3.1.1. Sistemas mecánicos traslacionales.

Los sistemas mecánicos son aquellos conformados por tres elementos básicos: masa,
resorte y amortiguador. La ley base sobre la cual se rigen es la primera ley de Newton,
que relaciona las fuerzas aplicadas sobre un objeto respecto a su desplazamiento:
X d2 x
Fi = m (3.1)
i
dt2
Donde Fi son las fuerzas vectoriales que actúan sobre el objeto, x es el valor de des-
plazamiento y m la masa del objeto.
En la tabla 3.1, se muestran las ecuaciones de dichos componentes:
En los sistemas mecánicos, la masa y el resorte, son elementos que almacenan energı́a.
15
16
Tabla 3.1: Ecuaciones de elementos de sistemas mecánicos traslacionales.

Bloque funcional Ecuación descriptiva Energı́a almacena-
da/potencia disipada
2
Resorte F = kx E = 12 Fk
2
Masa F = m ddt2x E = 12 mv 2
dx
Amortiguador F = c dt P = cv 2
Por su parte, el amortiguador es un elemento meramente disipativo.
Ejemplo: considere una masa que se mueve sobre un plano en la cual existe una ac-
ción de fuerza externa y fricción viscosa por la superficie de desplazamiento (b es el
coeficiente de fricción) como la que se muestra en la figura 1. Se deberá de obtener
la función de transferencia y la representación en espacio de estados, donde x(t) es el
desplazamiento y la fuerza externa es F (t).
Figura 3.1: Figura de ejemplo 1.
Aplicando la ley de Newton se obtiene (observación: dx/dt = ẋ)
F (t) − bẋ(t) = mẍ(t) (3.2)
Aplicando la transformada de Laplace con condiciones iniciales cero, se obtiene:
X(s) 1/m
G(s) = = (3.3)
F (s) s(s + b/m)
Para la representación en espacio de estados, se considera la fuerza aplicada como la

entrada u(t) = F (t) y la salida como el desplazamiento y(t) = x(t), y se definen dos
estados:
17
x1 = x(t) (3.4)
x2 (t) = ẋ(t) (3.5)
Al momento de realizar la sustitución se tiene:
ẋ1 = x2 (3.6)
b 1
ẋ2 = − x2 + u(t) (3.7)
m m
Con lo que se puede representar el modelo en el espacio de estados de forma matricial

como sigue:

ẋ1 0 1 x1 0
= + u(t) (3.8)
ẋ2 0 −b/m x2 1/m

x1
y(t) = 1 0 + 0 u(t) (3.9)
x2
3.1.2. Sistemas mecánicos rotacionales.

Los sistemas mecánicos rotacionales están conformados por tres elementos básicos: mo-
mento de inercia, resorte de torsión y amortiguador rotacional. La forma de obtener su
modelado matemático, es muy parecida a la de los sistemas tralacionales, teniendo las
ecuaciones que se describen en la tabla 3.2.
Tabla 3.2: Ecuaciones de elementos de sistemas mecánicos rotacionales

2
Resorte τ = kθ E = 12 τk
2
Masa τ = I ddt2θ E = 12 Iω 2
Amortiguador τ = c dθ
dt
P = cω 2
Las leyes fı́sicas que se aplican en estos sistemas son las mismas leyes de Newton pero
en versión de movimiento rotacional, teniendo que θ es el desplazamiento, I es el mo-
mento de inercia, τ es el par o torque y ω es la velocidad angular. La ecuación principal
es la que a continuación se muestra:
X X
τi = Fi l = I θ̈ (3.10)
i i
donde F es la fuerza lineal y l la distancia.

18
3.2. Objetivo
Obtener el modelado matemático de sistemas mecánicos rotacionales y traslacionales,
en función de transferencia y espacio de estados. Además de obtener las respuestas al
impulso y escalón de dichos sistemas en Matlab.
3.3. Material
Lápiz y papel
2. Equipo utilizado
3.4. Procedimiento
1. Obtener el modelo matemático que describe el comportamiento del siguiente sis-
tema descrito en la figura 3.2, r(t) representa la fuerza de entrada, y2 y y1 los
desplazamientos de las masas.
Donde se deben de utilizar las relaciones k2 (y2 − y1 ), b(ẏ2 − ẏ1 ) y k1 y1 . El modelo

matemático deberá quedar de la forma:
m1 ÿ1 + bẏ1 + (k2 + k1 )y1 = bẏ2 + k2 y2 + r(t) (3.11)

m2 ÿ2 + bẏ2 + k2 y2 = bẏ1 + k2 y1 (3.12)
2. Obtener la función de transferencia Y1 (s)/R(s) del sistema anterior.
3. Obtener la respuesta al impulso y la respuesta escalón del sistema en Simulink

de Matlab. Considerar m1 = m2 = 1, k1 = k2 = 2 y b = 3.
4. Obtener el modelo matemático del sistema mecánico rotacional descrito en la

figura 3.3.
El modelo deberá quedar de la forma:
τ = kθ + bθ̇ + I θ̈ (3.13)
5. Desarrolle la función de transferencia y el espacio de estados del sistema rotacio-
nal.
6. Obtener la respuesta al impulso y respuesta al escalón del sistema rotacional.
Según el tipo de respuesta clasifique el sistema (subamortiguado, sobreamorti-
guado, crı́ticamente amortiguado), usando k = b = I = 1. Obtener de forma
19
Figura 3.2: Sistema de dos masas, dos resosrtes y un amortiguador.
Figura 3.3: Sistema rotacional masa resorte amortiguador.

20
analı́tica el tipo de sistema utilizando el discriminante de la función de transfe-
rencia.

Se deberá de describir todo el procedimiento paso por paso para obtener el modelo
matemático, la función de transferencia y el espacio de estados de ambos sistemas.
Detallar el procedimiento utilizado en Simulink para obtener las respuestas co-
rrespondientes.
Las gráficas deberán guardarse en JPEG y no como captura de pantalla, para
esto hacer uso de simout.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
3.6. Conclusiones
Práctica 4
Modelado matemático de sistemas

fı́sicos: parte II
4.1. Introducción
4.1.1. Sistemas eléctricos.
Los sistemas eléctricos, al igual que los sistemas mecánicos rotacionales y traslaciona-
les, están compuestos de tres elementos básicos: resistencia, capacitor e inductor. Las
ecuaciones de dichos componentes se presentan en la tabla 4.1.
Tabla 4.1: Ecuaciones de sistemas eléctricos

di 1
E = 12 Li2
R
Inductor v = L dtR , i = L vdt
Capacitor v = C1 idt, i = C dvdt
E = 12 Cv 2
2
Resistor v = Ri, i = Rv P = vR
En donde el inductor y el capacitor, son elementos que almacenan energı́a, el capaci-

tor en forma de voltaje y el inductor en forma de corriente eléctrica; en tanto que la
resistencia es un elemento solo disipativo.
Las leyes fundamentales que se aplican a los sistemas eléctricos, son las leyes de Kir-
chhoff:
Corriente: la suma de corrientes que entran en un nodo, es igual a la suma de

corrientes que salen de él. La ecuación que se plica es:
n
X m
X
ijent = ijsal (4.1)
j=1 j=1
21
22
Voltaje: la suma de las diferencias de potencial sobre una malla es cero. Esto es:
n
X
vj = 0 (4.2)
j=1
Ejemplo: Encontrar la ecuación diferencial que gobierna al circuito eléctrico conside-

rando condiciones iniciales nulas, la variable manipulada será v(t) y la variable contro-
lada i(t). El circuito es un RLC en serie como el que se muestra en la figura 4.1.
Figura 4.1: Circuito RLC en serie.
Aplicando la ley de voltaje de Kirchhoff se obtiene:
v = vL + vR + vC (4.3)
Sustituyendo las respectivas ecuaciones para inductor, resistor y capacitor se obtiene:

Z
di 1
v = L + iR + idt (4.4)
dt C
Mediante la transformada de Laplace y relacionando corriente y voltaje, entonces la

función de transferencia es de la siguiente forma:
I(s) s/L
G(s) = = 2 (4.5)
V (s) s + sR/L + 1/LC
Para el espacio de estados se considera el voltaje de la fuente como la entrada u(t) = v(t)
y la salida la corriente en el capacitor i(t), entonces los estados estarán dados por:
Z
x1 = i(t)dt (4.6)
23
x2 = i(t) (4.7)
Al sustituirlos se obtiene:
ẋ1 = x2 (4.8)
1 R 1
ẋ2 = − x1 − x2 + u(t) (4.9)
LC L L
Quedando el espacio de estados como a continuación se muestra:

ẋ1 0 1 x1 0
= 1 R + 1 u(t) (4.10)
ẋ2 − LC L
x2 L

x1
y(t) = 1 0 + 0 u(t) (4.11)
x2
4.2. Objetivo
Obtener el modelado matemático de sistemas eléctricos, en función de transferencia y
espacio de estados. Además de obtener las respuestas al impulso y escalón de dichos
sistemas en Matlab.
4.3. Material
Lápiz y papel
2. Equipo utilizado
24
4.4. Procedimiento
1. Del modelo descrito en la figura 4.1, simule la función de transferencia en Simu-
link de Matlab con una entrada escalón para los siguientes valores:
R=L=C=1
R=2, L=C=1
R=3, L=C=1
R=0, L=C=1
2. Clasifique el sistema según el tipo de respuesta en amortiguado, subamortiguado,
crı́ticamente amortiguado, etc. Compruebe además de forma analı́tica analizando
el discriminante de la función de transferencia.
3. Obtenga el modelo matemático del circuito RLC en paralelo de la figura 4.2.
Figura 4.2: Circuito RLC en paralelo.
4. Obtenga la función de transferencia y el espacio de estados del sistema RLC en

paralelo, tome el voltaje en la resistencia como la salida.
5. Con la función de transferencia obtenga la salida del sistema con una entrada
escalón, cuando los valores de los componentes son los siguientes:
R=L=C=1
R=2, L=C=1
R=3, L=C=1
R=0, L=C=1
6. Clasifique el sistema según el tipo de respuesta en amortiguado, subamortiguado,
crı́ticamente amortiguado, etc. Compruebe además de forma analı́tica analizando
el discriminante de la función de transferencia.
25
Se deberá de describir todo el procedimiento paso por paso para obtener el modelo
matemático, la función de transferencia y el espacio de estados de ambos sistemas.
Detallar el procedimiento utilizado en Simulink para obtener las respuestas co-
rrespondientes, además de los discriminantes para clasificar los sistemas.
Las gráficas deberán guardarse en JPEG y no como captura de pantalla, para
esto hacer uso de simout.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
4.6. Conclusiones
26
Práctica 5
Diseño de controladores analógicos

por compensación de fase.
5.1. Introducción
Con el término de respuesta en frecuencia se hace referencia a la respuesta de un sistema
en estado estacionario a una entrada senoidal. Los métodos de respuesta de frecuencia
varı́an la frecuencia de la señal de entrada en un cierto rango para poder analizar la
respuesta de salida a estas variaciones de frecuencia.
Los métodos de respuesta en frecuencia se desarrollaron por Nyquist, Bode y Nichols

entre los años 1930-1940, los cuales son los más potentes en teorı́a de control clásica,
aunque también son indispensables en control robusto.
Para la presentación de las caracterı́sticas de la respuesta en frecuencia de forma gráfica

se puede utilizar un diagrama de Bode o diagrama logarı́timico. El diagrama de Bode
está conformado por dos gráficas: una de ellas la gráfica del logaritmo de la magnitud
de la función de transferencia y la otra es la gráfica del ángulo de fase; siendo ambas
dibujadas contra la frecuencia en escala logarı́tmica.
La ventaja principal de utilizar el diagrama de Bode reside en que la multiplicación

de magnitudes se convierte en suma. Además, cuenta con un método sencillo para di-
bujar una curva aproximada de magnitud logarı́tmica, que se basa en aproximaciones
asintóticas, la cual es suficiente si se desea obtener información de forma general de la
respuesta en frecuencia.
Cuando se utiliza el método de respuesta en frecuencia para controlar un sistema,

el comportamiento de la respuesta se especifica en términos de margen de fase, margen
de ganancia, magnitud de pico de resnonancia, ancho de banda, entre otros. Después
de diseñar el lazo abierto, se determinan los polos y ceros del sistema en lazo cerrado;
por lo que se deben verificar las caracterı́sticas de la respuesta en un diagrama de Bode.
27
28
Existen ocaciones en las que es necesario compensar un sistema en fase, teniendo ade-
lanto y atraso de la fase. La compensanción en adelanto produce un mejoramiento en
la respuesta transitoria y un cambio en la precisión en estado estable, aunque puede
acentuar los efectos del ruido de alta frecuencia. Por su parte la compensación en atraso
produce un mejoramiento notable en la precisión en estado estable aunque aumenta el
tiempo de la respuesta transitoria; además suprime los efectos de las señales de ruido
a altas frecuencias.
Existe otro tipo de compensación que combina adelanto y atraso, teniendo las ca-
racterı́sticas de ambos sistemas, aunque eleva el orden del sistema en 2, con lo que se
vuelve más complejo y difı́cil de controlar el comportamiento de la respuesta transitoria.
5.2. Objetivo
Diseñar un controlador por compensación de fase utilizando un compensador por ade-
lanto de fase, obteniendo las caracterı́sticas de la planta en lazo abierto y con base a
las especificaciones dadas, realizar el cálculo para las ganancias del compensador.
5.3. Material
Lápiz y papel
2. Equipo utilizado
5.4. Procedimiento
1. Del siguiente sistema, se sabe que se tiene un margen de fase entre 10◦ y 30◦ y
un ancho de banda menor a 0.01 rad/s. Utilce diagramas de bode para verificar
la información anterior. La función a analizar viene dada por:
1
F (s) = (5.1)
s3 + 20000s2+ 50s + 0.05
Notas:
Si el margen de fase es igual o superior a 55◦ , se asegura estabilidad.
Si el ancho de banda es muy pequeño, el sistema es lento.
2. Realice la gráfica de la respuesta a lazo abierto del sistema con una rampa como
entrada. La pendiente de la rampa será de 0.25
29
Determinar el error en estado estable.
Describa la respuesta del sistema.
3. Realizar un compensador por adelanto de fase cuyas principales caracterı́sticas
son:
Aumento del ancho de banda
El margen de fase puede mantenerse o ser aumentado
La función de transferencia será la siguiente:
Ts + 1
G = Kc (5.2)
αT s + 1
4. Calcule los valores del cero, polo y ganancia de la función de transferencia de este
compensador si se requiere un ancho de banda de 0.1 rad/s y un margen de fase
de 70◦ .
Investigue el tema e incluya las fórmulas utilizadas, ası́ como el procedimien-

to completo.
5. Realizar un programa en Matlab en archivo .m con las siguientes caracterı́sticas:
Recibirrá los datos del coeficiente del numerador y denominador del sistema
a compensar
Recibirá el margen de fase deseado
Recibirá el ancho de banda deseado
6. El programa deberá de devolver lo siguiente:
Función de transferencia del compensador

3 gráficas de Bode (sistema original, compensador y sistema en conjunto)
Margen de fase y ancho de banda resultantes
7. Teniendo el sistema con el compensador diseñado realice lo siguiente:
Simule el controlador a lazo abierto con entrada de rampa con pendiente de
0.25
Simular de nuevo en lazo cerrado con retroalimentación negativa y la misma
entrada, describa qué es lo que sucede.

Incluya en el reporte todas las fórmulas utilizadas y el procedimiento completo
Incluya los diagramas de Bode de todos los sistemas en el reporte
Icluir los resultados de los sistemas con y sin control
tes secciones:
30
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
5.6. Conclusiones
Práctica 6
Identificación de procesos de primer

orden
6.1. Introducción
El control de un proceso implica el uso de una serie de sensores, dependiendo del tipo
de variable a controlar, de tal manera que se puede obtener información en tiempo
real del proceso. De esta forma, se puede obtener un valor de error entre la señal de
referencia y la variable de entrada.
Actualmente, existen dispositivos para controlar procesos industriales como control

de velocidad de un motor, control de temperatura, control de humedad, por mencionar
algunos. Estos dispositivos disponen de un autoajuste para controlar los parámetros
del sistema de forma automática. Pero, algunas veces no se obtiene un correcto funcio-
namiento del proceso, por lo que es necesario diseñar el controlador de manera analı́tica.
En teorı́a de control se puede denominar como identificación de sistemas, al proceso

que nos permite obtener el valor de los parámetros de un modelo matemático a partir
de métodos experimentales. De tal forma que, se pueden encontrar métodos gráficos co-
mo el de la respuesta transitoria y métodos estadı́sticos como el de mı́nimos cuadrados.
En un sistema de primer orden con tiempo muerto (retardo), la función de transfe-

rencia se expresa de la siguiente forma:.
k
G(s) = e−t0 s (6.1)
τs + 1
31
32
6.1.1. Método gráfico para identificar un sistema de primer
orden
Se toma en consideración la figura 6.1 para explicar cómo obtener el valor de los paráme-
tros de un modelo de primer orden, considerando que la gráfica real corresponde a la
respuesta de un proceso en lazo abierto.
Figura 6.1: Obtención de parámetros de un modelo de primer orden.
El proceso consiste en lo siguiente:
∆y
El valor de k se obtiene del cociente de ∆u . Donde ∆y corresponde a la amplitud
de la respuesta del sistema y ∆u corresponde al valor del escalón de entrada.
Para el parámetro t0 se realiza lo siguiente: se traza una lı́nea tangente a la curva

de respuesta del proceso en el punto en que la velocidad de cambio es máxima
(punto de inflexión, en los casos en el que el sistema no sea de primer orden).
El valor se obtiene al medir el tiempo que tarda la respuesta desde el punto en
donde se aplica el escalón hasta el punto en donde se intercepta la lı́nea tangente
y el eje horizontal tiempo.
El valor de τ se obtiene midiendo el tiempo en donde el cambio en la respuesta

alcanza el 63.2 % de su valor de cambio total. Se traza una lı́nea horizontal en
dicho punto hasta que cruce con la curva y ahı́ se obtiene el tiempo t0 + τ . Sola-
mente se quita el valor de t0 y se obtiene el de τ .
33
6.1.2. Segundo método gráfico para identificación de paráme-
tros de procesos de primer orden
Un segundo método fue propuesto por Cecil L. Smith, ya que se puede tener una me-
nor precisión en el trazo de la tangente. Smith propone que los valores de t0 y τ se
seleccionen de tal manera que la respuesta del modelo coincida con la real en la región
de mayor cambio. Por lo que recomienda que los puntos sean t0 + τ /3 y t0 + τ . En la
figura 6.2 se pueden apreciar estos puntos como t1 y t2 respectivamente.
Figura 6.2: Segundo método para la obtención de parámetros de un modelo de primer

orden.
Ası́ los valores de t0 y de τ se obtienen mediante:
3
τ = (t2 − t1 ) (6.2)
2
t0 = t2 − τ (6.3)
Donde:
t1 = tiempo en el cual ∆y = 0.283 ∗ ∆y (6.4)
t2 = tiempo en el cual ∆y = 0.683 ∗ ∆y (6.5)
6.2. Objetivo
Diseñar un controlador por compensación de fase utilizando un compensador por ade-
lanto de fase, obteniendo las caracterı́sticas de la planta en lazo abierto y con base a
las especificaciones dadas, realizar el cálculo para las ganancias del compensador.
34
6.3. Material
Lápiz y papel
2. Equipo utilizado
6.4. Procedimiento
1. Generar una serie de datos correspondientes a un proceso de primer orden con
tiempo muerto (proponer una función de transferencia).
2. Obtener la respuesta del sistema utilizando un escalón unitario.

Identifique los parámetros de los datos generados usando las dos técnicas del
método gráfico:
Método 1: gráficas y procedimiento completo.

Método 2: gráficas y procedimiento completo.
3. Realizar tabla comparativa de los parámetros obtenidos contra los parámetros
del modelo original.
4. Grafique la función de tranferencia original en la misma gráfica con las funciones

de transferencia obtenidas mediante los métodos gráficos, con una entrada de
escalón unitario.

Incluir todas las gráficas obtenidas y procedimiento completo para la identifica-
ción de las funciones de transferencia.
Incluir una gráfica de la respuesta de todas las funciones de transferencia al

escalón unitario. Diferenciar cada respuesta con un color y etiqueta.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
35
6.6. Conclusiones
36
Práctica 7
Sintonización de controladores PID

por método experimental parte I.
7.1. Introducción
El elemento más importante de un lazo de control, es el controlador, el cual constituye
el cerebro de un sistema de control automático, este realiza las correcciones necesarias
para lograr que la variable de un proceso o variable controlada, alcance el valor de
referencia.
Actualmente, la mayorı́a de los controladores que se utilizan en la industria siguen

una estructura PID, cuyo significado proviene del tipo de acción que realizan sobre
una señal de error o variable manipilada: P de proporcional, I de integral y D de deri-
vativo. El controlador PID es sumamente robusto para la mayorı́a de las aplicaciones
prácticas, sobretodo en aquellas cuya dinámica no es rápida y de alto orden.
El control PID utiliza como mecanismo motor el error de seguimiento, esto es:
e(t) = r(t) − y(t) (7.1)
Donde e(t) representa la señal de error, r(t) es la señal de referencia y y(t) la señal de
salida del sistema.
La accióm de control se construye según la siguiente estructura básica:

Z t
de(t)
u(t) = Kp e(t) + Ki e(τ )dτ + Kd (7.2)
0 dt
Donde Kp , Ki y Kd representan las ganancias proporcional, integral y derivativa, res-

pectivamente. En la figura 1 se muestra el diagrama de bloques de un controlador PID.
37
38
Figura 7.1: Diagrama de bloques de un control PID.
En el controlador solamente se tienen que garantizar las ganancias Kp , Ki y Kd , de tal

manera que se garantice estabilidad y desempeño, tanto en estado transitorio, como en
estado estable.
El controlador PID en forma de función de transferencia puede escribirse como si-

gue:
Ki
C(s) = Kp + + Kd s (7.3)
s
En algunas ocasiones no es necesario utilizar todo el control PID, puede bastar con
alguna de las siguientes combinaciones:
C(s) = Kp (7.4)
Ki
C(s) = Kp + (7.5)
s
C(s) = Kp + Kd s (7.6)
Según se requiera en la aplicación.
Cada acción dentro del controlador tiene distintas funciones que se describen a conti-
nuación:
39
Acción proporcional. Provee una contribución que depende del valor instantáneo
del error, con este controlador se puede controlar cualquier planta estable, pero
su desempeño es bastante limitado, aunado a que el error en estado estable, no
se garantiza cero.
Acción integral. Provee una salida proporcional al error acumulado, por lo que
es de lenta reacción, además garantiza error en estado estable cero ante entrada
escalón unitario de referencia y perturbaciones.
Acción derivativa. Provee contribución proporcional a la tasa de cambio del error

de referencia, por lo que es de rápida acción, aunque desaparence en presencia
de errores constantes aunado a que puede amplificar el ruido presente en la planta.
Existen diferentes métodos para poder sintonizar las ganancias de un controlador, te-
niendose métodos experimentales y analı́ticos.
7.1.1. Método de Ziegler and Nichols.

Es un método presentado por Ziegler and Nichols, el cual está basado en el registro de
la respuesta del sistema en lazo abierto ante una entrada escalón. Esta respuesta está
caracterizada por los parámetros: τ , td y R. Se asume que la planta puede representarse
por la siguiente estructura:
Hetd s
G(s) = (7.7)
τs + 1
Para obtener los parámetros, se determina el punto donde la pendiente de la curva
de respuesta es máxima y se dibuja una lı́nea tangente en este punto. En el punto
donde intersecta con la lı́nea donde se estabiliza el sistema, se traza una recta vertical,
obteniendo los parámetros td y τ como se muestra en la figura 2.
Figura 7.2: Metodologı́a Ziegler-Nichols.
En donde H es la altura de la respuesta del sistema y td es el retardo. La respuesta es

40
obtenida al evaluar el escalón unitario y se toman los siquientes parámetros:
L = td (7.8)
H
R= (7.9)
τ
Las ganancias del controlador se obtienen para un amortiguamiento de ζ = 0.21, de

acuerdo a la tabla 7.1.
Tabla 7.1: Ganancias Ziegler Nichols para ζ = 0.21

Controlador Ganancias
1
P Kp = RL
0.9 0.27
PI Kp = RL , Ki = RL 2
1.2 0.6 0.6
PID Kp = RL , Ki = RL 2 , Kd = R
7.2. Objetivo
Diseñar un controlador por método gráfico de Ziegler-Nichols para sintonizar paráme-
tros de control en sistemas de primer orden.
7.3. Material
Lápiz y papel
2. Equipo utilizado
7.4. Procedimiento
1. Obtener la planta mediante el método de Ziegler-Nichols para la respuesta a es-
calón unitario de la Figura 7.3.
2. Realizar cálculo para obtener controladores P, PI y PID.

Realizar simulación del modelo sin controlar.
41
Figura 7.3: Respuesta de un sistema de primer orden.
3. Realizar simulación del modelo con cada uno de los controladores y mostrarlos
en la misma gráfica.
4. Obtener la planta mediante el método de Ziegler-Nichols para la respuesta a es-

calón unitario de la figura 7.4. La lı́nea tangente se debe de trazar en el punto de
inflexión.
Figura 7.4: Respuesta de segundo sistema.
5. Realizar cálculo para obtener controladores P, PI y PID.

42

Reportar: las gráficas con las lı́neas que se emplearon, todos los cálculos y las
ganancias obtenidas.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
7.6. Conclusiones
Práctica 8

por método experimental parte II.
8.1. Introducción
Cuando se descubrieron los sistems de control PID, la forma de ajuste de los paráme-
tros se basaba en la experiencia del operador de la planta, por lo que eran procesos
lentos en los cuales se requeria mucho tiempo para llevarlos a cabo e implicaban varios
ensayos de prueba y error para verificar los datos. Los ajustes permitı́an mejorar la
calidad en la producción entre un 5 % y un 25 %.
En 1942 Ziegler y Nichols desarrollaron diversos métodos para sintonizar parámetros

de un controlador, los cuales se obtuvieron a través de numerosas pruebas, sin tener
conocimiento especı́fico del sistema que estaban tratando.
Desarrollaron dos metodologı́as:
Ziegler-Nichols en lazo abierto.
Ziegler-Nichols en lazo cerrado.
Ambos métodos coinciden en el objetivo, el cual es conseguir un valor sobre el impulso

para reducirlo al 25 %.
El método de lazo abierto, fue descrito en la práctica anterior. El método de lazo

cerrado, mejor conocido como método de la última ganancia, se muestra en la figura
8.1 en diagrama de bloques, el cual consiste en obtener una ganancia máxima de forma
proporcional que se le puede asignar a un sistema en lazo cerrado sin que el sistema se
haga inestable. El sistema en retroalimentación comenzará a oscilar, y la ganancia que
provoca el régimen oscilatorio se le denominará Ku y el perı́odo resultante como Pu ,
que es el perı́odo de oscilación, como se muestra en la figura 8.2.
43
44
Figura 8.1: Diagrama de bloques método de última ganancia.
Figura 8.2: Perı́odo de oscilación método de última ganancia.
Una vez que se han obtenido los valores de la última ganancia y el último periodo,
Ziegler y Nichols proponen las fórmulas que se muestran en la tabla 8.1.
Tabla 8.1: Ganancias para el método de última ganancia.

Controlador Ganancias
P Kp = 12 Ku
PI Kp = 0.45Ku , Ki = 0.54 KPu
u
Ku
PID Kp = 0.6Ku , Ki = 1.2 Pu , Kd = 0.075Ku Pu
8.2. Objetivo
Diseñar un controlador por el método última ganancia para sintonizar parámetros de
control.
45
8.3. Material
Lápiz y papel
2. Equipo utilizado
8.4. Procedimiento
1. Obtener la planta de la respuesta del sistema de la figura 8.3 a una entrada es-
calón unitario.
Figura 8.3: Respuesta de un sistema a escalón unitario.
2. Simular la planta obtenida en el punto anterior en Matlab en retroalimentación

con control proporcional y encontrar la ganancia con la que el sistema oscila de
manera estable.
3. Obtener el periodo de oscilación Pu .
4. Sintonizar un controlador P, PI y PID.
5. Simular el sistema con los 3 controladores y graficar en una misma figura los
resultados.
46
6. Simular la siguiente función de transferencia en Simulink:
100
G(s) = (8.1)
s3 + 13s2 + 30s + 100
Obtener la ganancia para la cual el sistema oscila de manera estable.
7. Obtener el periodo de oscilación Pu .
8. Sintonizar un controlador P, PI y PID.
9. Simular el sistema con los 3 controladores y graficar en una misma figura los
resultados.

Reportar las gráficas con las lı́neas que se hicieron, todos los cálculos y las ga-
nancias obtenidas.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
8.6. Conclusiones
Práctica 9

forma analı́tica parte I.
9.1. Introducción
Cuando se cuenta con el modelo matemático de la planta a controlar, se puede em-
plear la sintonización analı́tica, y se puede definir el comportamiento deseado para la
respuesta transitoria y de estado estacionario en lazo cerrado.
Si solamente se desea estabilidad en lazo cerrado, se puede utilizar el criterio de Routh-

Hurwitz para encontrar el rango de variación permitido para las ganancias del contro-
lador, ya sea P, PI o PID.
En el criterio de estabilidad de Routh, la estabilidad de un sistema lineal invarian-

te en el tiempo se tendrá en sentido de una salida acotada a una entrada acotada, lo
cual puede apreciarse mediante la respuesta al impulso. Una función de transferencia
logrará la condición de estabilidad si y solamente si tiene todos sus polos con parte real
negativa.
El criterio de estabilidad de Routh-Hurwitz permite evaluar si una función de transfe-

rencia G(s) tiene todos sus polos estables, es decir, la parte real es negativa.
Se asume que G(s) es una función racional de la forma:
b0 sm + b1 sm−1 + ... + bm
G(s) = (9.1)
sn + a1 sn−1 + ... + an
La estabilidad de G(s) se evalúa en función de las raı́ces del denominador de la función

de transferencia. Enseguida se construye el arreglo de Routh y se evalúa el número de
cambios de signo en la primera columna del arreglo, lo cual indicará la cantidad de
raı́ces que están en el semiplano derecho (inestabilidad).
47
48
El arreglo de Routh para el polinomio, es como a continuación se muestra:
den(s) = sn + a1 sn−1 + a2 sn−2 + ... + an (9.2)
sn 1 a2 a4 ...
sn−1 a1 a3 a5 ...
sn−2 δ1 δ2 δ3 ...
sn−3 β1 β2 β3 ...
. . . . ...
. . . . ...
. . . . ...
s1 * * 0 0
s0 * 0 0 0
49
En donde:

−1 1 a2 a1 a2 − a3
δ1 = = (9.3)
a1 a1 a3 a1

−1 1 a4 a1 a4 − a5
δ2 = = (9.4)
a1 a1 a5 a1

−1 1 a6 a1 a6 − a7
δ3 = = (9.5)
a1 a1 a7 a1

−1 a1 a3 δ1 a3 − δ2 a1
β1 = = (9.6)
δ1 δ1 δ2 δ1

−1 a1 a5 δ1 a5 − δ3 a1
β2 = = (9.7)
δ1 δ1 δ3 δ1

−1 a1 a7 δ1 a7 − δ4 a1
β3 = = (9.8)
δ1 δ1 δ4 δ1
..
.
etc.
9.1.1. Ejemplo.
Considerar la siguiente planta:
1
G(s) = (9.9)
(s − 1)(s − 4)
Obtenga el rango de valores para un controlador proporcional Kp de tal forma que

mantenga estable el sistema retroalimentado.
Primeramente se realiza la retroalimentación unitaria del sistema, multiplicado por

el control proporcional de la siguiente forma:

1 1
L(s) = Kp ÷ 1 + Kp (9.10)
(s − 1)(s + 4) (s − 1)(s + 4)
Kp
= 2 (9.11)
s + 3s − 4 + Kp
Se analizarán las raı́ces del polinomio denominador:
s2 + 3s − 4 + Kp = 0 (9.12)
Por el criterio de Routh-Hurwitz

50
s2 1 Kp −4
s1 3 0
s0 Kp −4 0
Para evitar un cambio de signo y que el sistema se vuelva inestable se debe de cumplir
que:
Kp − 4 ≥ 0 (9.13)
Por lo que Kp debe ser mayor o igual que 4. Como nota, cabe destacar que el valor de
4 es un valor lı́mite, por lo que la ganancia debe ser mayor a este valor.
9.2. Objetivo
Diseñar un controlador criterio de Routh-Hurwitz para sintonizar parámetros de contol
de manera analı́tica.
9.3. Material
Lápiz y papel
2. Equipo utilizado
9.4. Procedimiento
1. Obtener el sistema retroalimentado con un controlador proporcional para la si-
guiente planta:
s+1
G(s) = (9.14)
s(s − 2)(s + 4)
2. Realizar el criterio de Routh-Hurwitz para encontrar el rango de valores de la
ganancia proporcional.
3. Simular el sistema sin control en Simulink.
4. Simular el sistema con control en un rango de valores, empezando por el valor

lı́mite (mı́nimo 4 valores más).
51
5. Obtener en una misma gráfica la respuesta de los controladores y determinar cuál
fue la mejor ganancia.
6. Considerar la siguiente función de transferencia:
1
G(s) = (9.15)
(s + 3)(s + 1)
7. Obtener el sistema retroalimentado con un controlador PI aplicado a la planta.
8. Obtener el rango de las ganancias proporcional e integral mediante el criterio de

estabilidad de Routh-Hurwitz.
9. Simular el sistema sin controlar en Simulink.
10. Simular el sistema con controlador para varios rangos de valores, empezando por
los valores lı́mite (establecer al menos 4 valores más por cada ganancia).
11. Reportar los resultados de los controladores en una misma gráfica y determinar
cuáles fueron las mejores ganancias.

Incluir todo el procedimiento para el cálculo de las ganancias.
Reportar todas las gráficas obtenidas de los sistemas con y sin control.
Establecer un rango de valores para las ganancias y simularlas, a fin de obtener

el mejor control.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
52
9.6. Conclusiones
Práctica 10

forma analı́tica parte II.
10.1. Introducción
Otro método de sintonización analı́tica, que se utiliza cuando se conoce el modelo
matemático de la planta, es la ubicación de los polos dominantes en el lazo cerrado.
Durante el proceso de diseño es posible poder cancelar un polo estable con dinámica
lenta de la planta con un cero del controlador.
Se toman además consideraciones en la respuesta estacionaria para poder definir las

restricciones en las ganancias del controlador. Si se requiere un error de cero en estado
estable, se necesita que el controlador contenga acción integral, es decir PI ó PID.
10.1.1. Ejemplo.
Considerar el siguiente modelo matemático con función de transferencia:
250
G(s) = (10.1)
(s + 5)(s + 12)
Se diseñará un controlador PI que tiene la siguiente forma:

Ki
C(s) = Kp + (10.2)
s
Cual busca que garantice:
Error en estado estable cero ante escalones de referencia.

Error en estado estable menor o igual a uno para una rampa de referencia eer = 1.
53
54
Sobretiro menor o igual al 5 %
Tiempo de asentamiento menor o igual a 1 segundo
Primeramente se analizan los polos del sistema en lazo abierto, teniendo que la función
de transferencia propuesta tiene dos polos estables, uno en s = −5 (polo lento) y otro
en s = −12, y ningún cero.
El controlador PI puede reescribirse de la siguiente forma:
(s + Ki /Kp )
C(s) = Kp (10.3)
s
El cual tiene un polo en s = 0 y un cero en s = −Ki /Kp . Por lo que se puede proponer
cancelar un polo lento de la planta con el cero del controlador, por lo tanto se define
el cero como:
Ki
=5 (10.4)
Kp
Por lo que la planta en lazo abierto queda como:
(s + 5) 250
L(s) = C(s)G(s) = Kp (10.5)
s (s + 5)(s + 12)
250Kp
= (10.6)
s(s + 12)
Debido a que el controlador contiene parte integral, se cumple que el error en estado
estable sea igual a cero. Para la condición de la rampa se utiliza el teorema del valor
final, con lo que se tiene que:
lı́m sL(s) ≥ 1/eer (10.7)

s=0
Por lo que se tiene:
250Kp 250Kp 1
lı́m s = ≥ (10.8)
s=0 s(s + 12) 12 1
Por lo tanto:
12
Kp ≥ = 0.048 (10.9)
250
Siguiendo con los requerimentos del ejercicio, se tiene una condición de sobretiro me-
nor o igual al 5 %, por lo que se tiene una restricción para el amortiguamiento en lazo
55
cerrado como a continuación se muestra:
s
(ln(0.05))2
ζ≥ = 0.69 (10.10)
π 2 + (ln(0.05))2
Enseguida, la función de transferencia en retroalimentación queda:
L(s) 250Kp
H(s) = = 2 (10.11)
1 + L(s) s + 12s + 250Kp
Al comparar los elementos con el polinomio caracterı́stico de segundo orden, se tiene

que:
ωn2 = 250Kp (10.12)

2ζωn = 12 (10.13)
Entonces:
12
ζ= p ≥ 0.69 (10.14)
2 250Kp
De donde:
Kp = 0.472 (10.15)
Ası́ Kp puede tomar cualquier valor entre 0.048 ≤ Kp ≤ 0.472, por lo que se fijará el
valor en Kp = 0.2. Para la última restricción, se requiere que el tiempo de asentamiento
sea menor a 1 segundo, por lo que se tiene:
3.91 3.91
ts = = p = 0.802 (10.16)
ζωn (0.69)( (250)(0.2)
El tiempo obtenido es menor a 1 segundo, pero en caso contrario de no serlo, se habrı́a

de ajustar de nueva cuenta la ganancia.
Finalmente, se ajustó el valor de Kp = 0.2 y sustituyendo este valor en la ecuación

(10.4), se tiene como resultado que:
Ki = 1 (10.17)
Por lo que el controlador queda de la forma:

1
C(s) = 0.2 + (10.18)
s
56
10.2. Objetivo
Diseñar un controlador de modificación de polos de forma analı́tica para sintonización
de parámetros de control.
10.3. Material
Lápiz y papel
2. Equipo utilizado
10.4. Procedimiento
1. Diseñar un controlador PI para la siguiente planta:
325
G(s) = (10.19)
s2+ 7s + 10
2. Se deberá de considerar los siguientes parámetros de diseño:
Error en estado estable al escalón igual a 0

Error en estado estable a rampa menor o igual a 0.5
Sobretiro menor o igual a 7.5 %
El tiempo de asentamiento menor a 2 segundos
3. Simular el sistema sin controlador en Simulink.
4. Simular el sistema con control y graficarlo en la misma figura que el sistema sin
controlador.

Incluir todo el procedimiento para el cálculo de las ganancias.
Reportar todas las gráficas obtenidas de los sistemas con y sin control.
tes secciones:
1. Portada
57
3. Desarrollo
4. Resultados
6. Referencias
10.6. Conclusiones
58
Práctica 11
Creación de planta analógica de

segundo orden.
Un amplificador operacional, es un amplificador diferencial de muy alta ganancia que
tiene como caracterı́stica: alta impedancia de entrada y baja impedancia de salida. Den-
tro de sus aplicaciones podemos encontrar la resolución de operaciones matemáticas
como suma, resta, integración y derivación. Actualmente los amplificadores operacio-
nales son circuitos integrados confiables y baratos con voltajes en CD relativamente
bajos. Un amplificador operacional básico se muestra en la figura 11.1.
Figura 11.1: Circuito amplificador operacional.
Gracias a las distintas configuraciones y operaciones que se pueden implementar en un

amplificador operacional, es posible resolver ecuaciones diferenciales ordinarias, usando
principalmente: amplificador inversor, amplificador sumador y amplificador integrador.
59
60
Los sistemas de segundo orden son ampliamente estudiados en teorı́a de control, y

se caracterizan porque en su función de transferencia contienen dos polos. Estos sis-
temas tienen distintos tipos de respuestas de acuerdo a las raı́ces de su polinomio
caracterı́stico:
s2 + 2ζωn s + ωn2 (11.1)

En donde al igualar a cero este polinomio, es posible obtener diversos sistemas:
Sistemas subamortiguados
Sistemas sobreamortiguados
Sistemas crı́ticamente amortiguados
Sistemas oscilatorios
En el caso de la respuesta subamortiguada, se puede caracterizar la respuesta y con
base a ella seleccionar los parámetros de control para mejorarla y cada término se de-
fine como:
Tiempo de crecimiento: es el tiempo que le lleva a la respuesta al escalón ir del

0-100 % de su respuesta estacionaria.
Tiempo pico: tiempo en llegar la respuesta al escalón máximo.
Sobretiro: valor pico de la respuesta al escalón.
Tiempo de asentamiento: tiempo en que la respuesta al escalón alcanza un error
menor o igual al 2 % de su valor estacionario.
11.2. Objetivo
Construir una planta analógica de segundo orden mediante amplificadores operaciona-
les para poder comprobar los valores de los elementos fı́sicos y variar la constante de
amortiguamiento.
11.3. Material
Lápiz y papel
Amplificadores operacionales
Resistencias
Capacitores
Protoboard
Cable para protoboard
2. Equipo utilizado
Fuentes de poder
61
Generador de funciones
Osciloscopio
Equipo de cómputo con software de simulación de circuitos (Multisim, Ps-
pice, Proteus, etc.)
11.4. Procedimiento
1. Con base a la configuración que se muestra en la figura 11.2 de una planta analógi-
ca de segundo orden, se obtiene la siguiente función de transferencia:
Y (s) k1 k2 kc1 kc2

G(s) = = 2 (11.2)
X(s) s + k3 k4 kc1 s + k5 kc1 kc2
En donde:
R2
k1 = (11.3)
R1
R4
k2 = (11.4)
R3
R8
k3 = (11.5)
R7
R4
k4 = (11.6)
R9
R4
k5 = (11.7)
R10
1
kc1 = (11.8)
R5 C1
1
kc2 = (11.9)
R6 C2
2. Tomando como base la fórmula general de las funciones de transferencia de se-

gundo orden:
K
G(s) = (11.10)
s2 + 2ζωn s + ωn2
Identifique los siguientes parámetros:
K =?
ωn =?
ζ =?
62
Figura 11.2: Planta analógica de segundo orden.
3. Se requiere modificar solamente el coeficiente de amortiguamiento ζ sin alterar

ninguna de las otras caracterı́sticas de la planta, para lograrlo, defina valores
convenientes de los componentes reales dejando solo uno como variable, a fin de
variar ζ y obtener:
Sistema subamortiguado
Sistema sobreamortiguado
Sistema crı́ticamente amortiguado
4. Realice la simulación de los 3 casos en Simulink considerando un tren de pulsos
a la entrada de 5V. Defina un periodo adecuado utilizando un ciclo de trabajo
de 50 %.
5. Realice la simulación de los 3 casos en un simulador de circuitos (Multisim, Pspice,
Proteus, etc.), considerando un tren de pulsos de 5V a la entrada con un ciclo de
trabajo de 50 %. Defina el periodo adecuado.
6. Construya el sistema completo en protoboard con amplificadores operacionales,
resistencias capacitores y obtenga las 3 diferentes configuraciones. Ajuste el gene-
rador de funciones a un tren de pulsos de entrada de 5V con un ciclo de trabajo
de 50 %.
7. Compare los resultados obtenidos con los de simulación.

Incluir todo el procedimiento para el cálculo de los parámetros del circuito.
63
Incluir los diagramas de simulación de Simulink y del circuito.
Incluir todas las gráficas obtenidas tanto en simulación, como experimental.
Realizar la comparativa entre los resultados obtenidos.
En los resultados incluya siempre la señal de referencia.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
11.6. Conclusiones
64
Práctica 12
Implementación de un PID
analógico.
Muchas de las aplicaciones de procesos de control se pueden resolver mediante contro-
ladores PID. En la actualidad la mayorı́a de los controladores que existen son digitales,
aunque una implementación analógica puede ser igual de efectiva y más económica.
Gracias a las diferentes configuraciones de los amplificadores operacionales, es posible

obtener de manera experimental cada una de las etapas de un controlador: proporcio-
nal, integral y derivativa.
12.1.1. Etapa proporcional.

Es la etapa más simple y consiste de un amplificador con una ganancia Kp ajustable.
La función de esta acción es aumentar la velocidad de respuesta y reducir el error en
estado estacionario. El diagrama del circuito se muestra en la figura 12.1.
Donde la ganancia viene dada por:
V sal R2
=− (12.1)
V en R1
12.1.2. Etapa integral.

Esta etapa añade la capacidad de procesamiento temporar del controlador, guardando
la historia de la magnitud del error y contribuir a reducirlo a cero en el estado estacio-
nario. En la figura 12.2 se muestra la etapa integral con amplificador operacional.
En donde la ganancia viene dada por:
65
66
Figura 12.1: Etapa proporcional.
Figura 12.2: Etapa integral.
V sal 1 1
=− (12.2)
V en R Cs
12.1.3. Etapa derivadora.

A pesar de que la etapa integral reduce el error en estado estacionario, vuelve lento la
respuesta del sistema. Al añadir una etapa derivadora permite al controlador mejorar
la etapa de amortiguamiento, permitiendo aumentar la etapa proporcional y por ende
la velocidad de respuesta. Esta etapa se muestra en la figura 12.3.
En donde la ganancia viene dada por:

67
Figura 12.3: Etapa derivadora.
V sal
= −RCs (12.3)
V en
12.2. Objetivo
Implementar un controlador PID utilizando amplificadores operacionales y sintonizarlo
para el sistema de segundo orden descrito en la práctica anterior.
12.3. Material
Lápiz y papel
Amplificadores operacionales
Resistencias
Capacitores
Protoboard
Cable para protoboard
2. Equipo utilizado
Fuentes de poder
Generador de funciones
Osciloscopio
Equipo de cómputo con software de simulación de circuitos (Multisim, Ps-
pice, Proteus, etc.)
68
12.4. Procedimiento
1. Considerando el proceso de segundo orden construido en la práctica anterior co-
mo planta en modo subamortiguado, sintonizar un PID o PI analógico que mejore
alguno de los parámetros de la siguiente lista:
Error en estado estable

Sobretiro
2. Realizar la simulación del controlador sintonizado y la planta considerando una
retroalimentación negativa unitaria en Simulink. Considerar un tren de pulsos a
la entrada, definiendo un periodo adecuado y un ciclo de trabajo del 50 %.
3. Realizar la simulación del controlador sintonizado y la planta en un software
de simulación de circuitos y considerar la misma entrada descrita en el punto
anterior, con una amplitud de 5V.
4. Construir el sistema completo utilizando amplificadores operacionales.
5. Obtener respuestas reales del sistema. Ajustar el generador de funciones a un
tren de pulsos a la entrada con la amplitud y frecuencia simulados.

Incluir todo el procedimiento para el cálculo de las ganancias del controlador.
Incluir los diagramas de simulación de Simulink y del circuito.
Incluir todas las gráficas obtenidas tanto en simulación, como de forma experi-
mental.
Realizar la comparativa entre los resultados obtenidos.
En los resultados incluya siempre la señal de referencia.
tes secciones:
1. Portada
3. Desarrollo
4. Resultados
6. Referencias
69
12.6. Conclusiones
70
Bibliografı́a
[1] K. Ogata. Ingenierı́a de control moderna. Pearson Educación, Madrid, 2010.

ISBN:978-84-8322-660-5.
[2] B. Kuo. Sistemas de control automático. 7ma edición. Pearson Educación, 1996.
ISBN: 9688807230.
[3] N. Matthew, O. Sadiku. Fundamentos de los circuitos eléctricos. Tercera edición.

McGraw Hill, 2006.
[4] L.M. Faulkenberry. Introducción a los amplificadores operacionales con aplicacio-

nes a CI lineales. Limusa, 2005.
71
72
Índice de figuras
2.1. Bloque donde se especifica la función de transferencia. . . . . . . . . . . 7

2.2. Configuración de parámetros de Simulink. . . . . . . . . . . . . . . . . 8
2.3. Diagrama de bloques de lazo de control retroalimentado. . . . . . . . . . 11
2.4. Implementación de lazo de control en Simulink. . . . . . . . . . . . . . 12
2.5. Diagramas de Simulink para comparar sistemas de segundo orden. . . . 12
2.6. Sistema en Simulink con tiempo muerto. . . . . . . . . . . . . . . . . . 13
3.1. Figura de ejemplo 1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

3.2. Sistema de dos masas, dos resosrtes y un amortiguador. . . . . . . . . . 19
3.3. Sistema rotacional masa resorte amortiguador. . . . . . . . . . . . . . . 19
4.1. Circuito RLC en serie. . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

4.2. Circuito RLC en paralelo. . . . . . . . . . . . . . . . . . . . . . . . . . 24
6.1. Obtención de parámetros de un modelo de primer orden. . . . . . . . . 32

6.2. Segundo método para la obtención de parámetros de un modelo de primer
orden. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
7.1. Diagrama de bloques de un control PID. . . . . . . . . . . . . . . . . . 38

7.2. Metodologı́a Ziegler-Nichols. . . . . . . . . . . . . . . . . . . . . . . . . 39
7.3. Respuesta de un sistema de primer orden. . . . . . . . . . . . . . . . . 41
7.4. Respuesta de segundo sistema. . . . . . . . . . . . . . . . . . . . . . . . 41
8.1. Diagrama de bloques método de última ganancia. . . . . . . . . . . . . . 44

8.2. Perı́odo de oscilación método de última ganancia. . . . . . . . . . . . . 44
8.3. Respuesta de un sistema a escalón unitario. . . . . . . . . . . . . . . . 45
11.1. Circuito amplificador operacional. . . . . . . . . . . . . . . . . . . . . . 59

11.2. Planta analógica de segundo orden. . . . . . . . . . . . . . . . . . . . . 62
12.1. Etapa proporcional. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66

12.2. Etapa integral. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
12.3. Etapa derivadora. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
73
74
Índice de cuadros
1.1. Ejemplos de comandos en Matlab . . . . . . . . . . . . . . . . . . . . . 5
2.1. Ejemplos de bloques en Simulink . . . . . . . . . . . . . . . . . . . . . 9

2.2. Parámetros de simulación en Simulink . . . . . . . . . . . . . . . . . . 10
2.3. Instrucciones de Matlab . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.1. Ecuaciones de elementos de sistemas mecánicos traslacionales. . . . . . 16

3.2. Ecuaciones de elementos de sistemas mecánicos rotacionales . . . . . . 17
4.1. Ecuaciones de sistemas eléctricos . . . . . . . . . . . . . . . . . . . . . 21
7.1. Ganancias Ziegler Nichols para ζ = 0.21 . . . . . . . . . . . . . . . . . 40
8.1. Ganancias para el método de última ganancia. . . . . . . . . . . . . . . 44
75
76

Sistemas de Control I

Cargado por

Copyright:

Formatos disponibles

Sistemas de Control I

Cargado por

Información del documento

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

Sistemas de Control I

Cargado por

Copyright:

Formatos disponibles

Prácticas de Sistemas

Roberto Martínez Montejano

Isaac Campos Cantón

Raúl Castillo Meraz

Tecnología y ciencias aplicadas

© Roberto Martínez Montejano, Isaac Campos Cantón, Raúl Castillo Meraz.

Edición a cargo Roberto Martínez Montejano, Isaac Campos Cantón.

ISBN: 978-607-535-122-3 (e-Libro)

San Luis Potosí, México.

2. Uso básico de Simulink 7

3. Modelado matemático de sistemas fı́sicos:parte 1 15

4. Modelado matemático de sistemas fı́sicos: parte II 21

6. Identificación de procesos de primer orden 31

7. Sintonización de controladores PID por método experimental parte

8. Sintonización de controladores PID por método experimental parte

9. Sintonización de controladores PID forma analı́tica parte I. 47

10.Sintonización de controladores PID forma analı́tica parte II. 53

11.Creación de planta analógica de segundo orden. 59

12.Implementación de un PID analógico. 65

Los sistemas de control son de mucha utilidad en el campo de la investigación y en

Estas prácticas de control abarca lo que se ha denominado teorı́a de control clásico, se

El modelado matemático de un sistema dinámico, es el conjunto de ecuaciones que

Un sistema que se denomina lineal e invariante en el tiempo, es aquel en el que se

En un sistema lineal invariante en el tiempo, una función de transferencia se define

En donde u(t) representa la entrada y y(t) la salida. La función de transferencia que

Y (s) bm sm + bm−1 sm−1 + .. + b1 s + b0

A partir de este concepto, es posible representar la dinámica de un sistema mediante

Tabla 1.1: Ejemplos de comandos en Matlab

1.4. Reporte de resultados

Uso básico de Simulink

Las funciones de transferencia de los bloques se escribe dentro de estos, interconec-

Figura 2.1: Bloque donde se especifica la función de transferencia.

Un diagrama de bloques contiene información relacionada con el comportamiento

Se pueden modificar los parámetros de la simulación en Simulink, para ajustarlo a

Figura 2.2: Configuración de parámetros de Simulink.

Tabla 2.1: Ejemplos de bloques en Simulink

Genera una señal de tipo escalón, en Sources

Define una función de transferencia Continuous

Simula el tiempo muerto del proceso. Continuous

Multiplexa las señales de entrada en Signal Routing

Genera una gráfica de la simulación con Sinks

Permite transferir los resultados de la Sinks

Suma señales. El número y signo de las Math Opera-

Tabla 2.2: Parámetros de simulación en Simulink

Tabla 2.3: Instrucciones de Matlab

Figura 2.3: Diagrama de bloques de lazo de control retroalimentado.

La figura 2.4 muestra la implementación en Simulink del sistema de control mos-

Figura 2.4: Implementación de lazo de control en Simulink.

2. Simule un cambio de referencia R(t) de 0 a 1 iniciando en t=1, además de la

Figura 2.5: Diagramas de Simulink para comparar sistemas de segundo orden.

4. Utilice el diagrama de Simulink descrito en el punto anterior para comparar la

6. Grafique los resultados de la simulación en Matlab, utilizando el siguiente código:

7. El diagrama de simulink de la figura 2.6, muestra la manera de incluir un tiem-

Figura 2.6: Sistema en Simulink con tiempo muerto.

8. Construya las siguientes funciones de transferencia y defina los parámetros. Mues-

2.5. Reporte de resultados

Modelado matemático de sistemas

3.1.1. Sistemas mecánicos traslacionales.

En la tabla 3.1, se muestran las ecuaciones de dichos componentes:

Tabla 3.1: Ecuaciones de elementos de sistemas mecánicos traslacionales.