Laboratorio 02

Control Clásico – Matlab y Simulink
Laboratorio 02
Programas de Matlab y Simulink

1. Representación en Matlab de los resultados de una
simulación
En la práctica previa se ha utilizado el elemento Scope para representar los valores
que toma una señal cualquiera en un esquema Simulink. Sería interesante disponer,
no sólo de la representación de esas señales sino también de sus valores numéricos,
para poder trabajar posteriormente sobre ellos (por ejemplo, para poder calcular el
valor máximo y mínimo de una señal y los instantes de tiempo en que se producen,
etc).
El elemento Scope permite, al mismo tiempo que representa los resultados de una
simulación, guardar esos resultados en una variable de Matlab, de modo que se
puedan consultar sus valores desde la ventana de comandos.
Como ejemplo volveremos a trabajar sobre el esquema utilizado en la práctica anterior

para simular el comportamiento de un conjunto masa – resorte – amortiguador. El
aspecto del esquema es el que se muestra a continuación:
𝒙
𝒌
F
𝒃 m
Se realizará la simulación de este esquema con los siguientes parámetros:

 Tiempo de simulación: 10 segundos.
 Amplitud del escalón de entrada: 5 unidades.
Docente: Jorge Ortiz Porras 1

El resultado debería ser igual al que muestra la siguiente figura:
Nota: en caso de no obtener exactamente estos resultados es posible que la

configuración del bloque escalón no sea la correcta. En la figura siguiente se muestran
los parámetros del bloque escalón y su significado sobre la representación gráfica del
mismo:

En general, el escalón se producirá en el instante cero (no estará retrasado) con lo que
el parámetro step time deberá ser siempre cero; el valor inicial del escalón (initial
value) también será cero; y el valor final (parámetro final value) deberá indicar la
amplitud del escalón, que en este caso particular es cinco.
Ahora haremos click sobre el botón Configuration Properties del elemento Scope,
para abrir la ventana de configuración:
En la ventana de Configuration Properties deberemos seleccionar logging y en ella

marcaremos la opción ‘Log data to workspace’ o, lo que es lo mismo, ‘registrar
datos en el espacio de trabajo de Matlab’.

Será necesario indicar dos parámetros:

 En primer lugar, tendremos que asignar un nombre a la variable en la que
deseamos guardar los datos. Por defecto, esta variable es ‘ScopeData’, pero
nosotros cambiaremos ese nombre por ‘posicion’ (sin tilde para evitar
problemas) que es el dato que representa ese osciloscopio en nuestro sistema
(posición del extremo del conjunto masa – resorte – amortiguador).
 En segundo lugar, será necesario especificar como guardar formato para los
datos, Save format. Simulink ofrece tres posibilidades: ‘Structure with time’,
‘Structure’ y ‘Array’. Nosotros elegiremos este último formato, que quiere decir
que los datos se guardarán en forma de vector.
Una vez hecho esto, lanzaremos de nuevo la simulación. El resultado sobre la ventana
de Simulink será el mismo, pero en la ventana de comandos de Matlab podremos
comprobar cómo se han creado dos nuevas variables. Para ello bastará teclear el
comando who:
>> who
Your variables are:
posicion tout
Se han creado dos variables: posicion y tout. La única variable que nos interesa es la
primera de ellas.
Si comprobamos el tamaño de la variable con el comando size, veremos que se trata
de un vector de 2 columnas y 58 filas (este último dato puede variar según los
ordenadores):
>> size(out.posicion)
ans =
58 2
La primera de las columnas contiene instantes de tiempo y la segunda contiene los
valores que toma la variable que se representa en el osciloscopio en cada instante. Si

representamos las primeras 10 filas de la variable podremos comprobar esto (nota: se

usa la instrucción ‘format long’ para mostrar más cifras decimales):
>> format long
>> out.posicion (1:10,:)
ans =
0 0
0.000000000000000 0.000000000000000
0.000040190182904 0.000000004038073
0.000241141097425 0.000000145360881
0.001245895670029 0.000003879023912
0.006269668533048 0.000098063591952
0.031388532848144 0.002435533506525
0.132449768341823 0.041378878201275
0.311148914122735 0.202039697904474
0.511148914122735 0.447797118045490
La primera de las columnas, como se ha dicho, corresponde a los instantes de tiempo;
y podemos comprobar cómo no están separados uniformemente. La segunda columna
corresponde a los valores que toma la señal, y podemos comprobar que los datos son
correctos haciendo zoom sobre el gráfico que aparece en el elemento Scope (por
ejemplo, en el instante 0.511 segundos el valor de la señal es 0.448
aproximadamente):

Desde la ventana de comandos de Matlab, también será posible representar mediante

la instrucción plot la señal guardada en la variable posición. Dado que la primera
columna contiene valores de tiempo y la segunda columna contiene valores de
posiciones, la forma correcta de la instrucción plot debería ser:
>> plot(out.posicion(:,1), out.posicion(:,2))
Supongamos que nos interesa calcular el valor máximo de la señal. Dado que los
valores de la señal están contenidos en la segunda columna de la variable posicion,
bastará con utilizar el comando max de Matlab (se recomienda consultar la ayuda de
Matlab acerca de este comando):
>> [maximo, indice] = max(out.posicion(:,2))

maximo =
0.788347181331777
indice =
12
Los valores pueden ser ligeramente distintos en diferentes ordenadores debido a las
precisiones de los cálculos. En cualquier caso, el comando max recorre la segunda
columna de la variable posicion y nos devuelve su valor máximo y cuál es la fila en la
que se produce ese valor (lo que hemos denominado índice). Si queremos comprobar
a qué instante de tiempo corresponde ese máximo, bastará con comprobar el valor
que toma la primera columna de la variable posicion para ese mismo índice (nota: se
utiliza la instrucción format short para mostrar menos cifras decimales).
>> format short
>> out.posicion(indice, :)
ans =
0.9111 0.7883

En resumen, hemos obtenido como resultado que la posición de la masa móvil (x) en
el ejemplo del sistema muelle-amortiguador alcanza un valor máximo de 0.7883
metros en el instante 0.9111 segundos.
2. Introducción de señales arbitrarias en simulink

En la práctica anterior se comprobó como Simulink dispone de muchas señales de
entrada, válidas para la mayor parte de las aplicaciones: se dispone de las señales
escalón, rampa, señal senoidal, tren de pulsos, etc.
No obstante, en algunas ocasiones se desea utilizar señales de forma arbitraria, no
disponibles directamente en Simulink. Por ejemplo, podríamos desear utilizar una
señal de entrada (fuerza aplicada al conjunto masa – resorte – amortiguador en el
sistema sobre el que se está trabajando en esta práctica) como la que se muestra a
continuación:
𝒙
2𝟎𝟎
𝟏𝟎𝟎
𝟏𝟎 𝟑𝟎 𝟑𝟓 𝒈
La forma más sencilla de introducir una señal como la propuesta es mediante el
bloque From Workspace de la categoría Sources. A continuación, se muestra el
aspecto que debe tener el esquema Simulink de la práctica una vez sustituido el
bloque escalón por el bloque From workspace:

Una vez modificado el esquema Simulink, haremos clic sobre el elemento ‘From
Workspace’ para acceder a sus parámetros de configuración. Aparecerá una ventana
como la que se muestra a continuación:
El único parámetro que nos interesa por el momento es el nombre de la variable de

Matlab de la que Simulink tomará los datos para utilizarlos como entrada. Por
defecto, esta variable es simin.
Nosotros cambiaremos este nombre por otro más explicativo: dado que la entrada, en
el sistema considerado, es la fuerza aplicada sobre el conjunto masa – Resorte –
amortiguador, llamaremos a esta variable fuerza.
El formato que debe tener esta variable es similar al de las variables guardadas en
Matlab a través del bloque Scope: la primera columna debe contener los valores de
los instantes de tiempo y la segunda los valores de la señal.
Un dato importante que hay que fijar es el intervalo entre cada dos instantes de
tiempo. Cuanto más pequeño sea este intervalo, mejor será la definición de la señal.
En este ejemplo utilizaremos un intervalo de 0.1 segundos.
A continuación, pasaremos a definir la señal (véase el gráfico anterior de x(t)). Lo más
cómodo será definirla por tramos, utilizando un tramo para cada trazo recto de la
señal.
Además, y para simplificar la creación de la señal, se usarán variables distintas para el
tiempo y para la señal. Los tiempos se guardarán en la variable tmp y los valores de la
señal (o datos) en la variable dat.
Primer tramo: entre 0 y 10 segundos. La variable tiempo (vector tmp1) deberá tomar
valores entre 0 y 10 a intervalos de 0.1; y la señal (vector dat1) deberá responder a la
ecuación de la recta x=10t.

Esto lo conseguimos con las siguientes sentencias de Matlab:
>> tmp1 = [0:0.1:10]

>> dat1 = 10*tmp1
Segundo tramo: entre los instantes 10.1 segundos y 30 segundos. la variable tiempo
(vector tmp2) deberá tomar valores entre 10.1 y 30 a intervalos de 0.1; y la señal
(vector dat2) deberá responder a la ecuación de la recta x=100. Esto lo conseguimos
con las siguientes sentencias de Matlab:
>> tmp2 = [10.1:0.1:30]

>> dat2(1:200)=100
La segunda instrucción crea un vector de 200 elementos, tantos elementos como tiene
el vector tmp2, todos ellos con valor 100.
Tercer tramo: entre 30.1 y 35 segundos. La variable tiempo (vector tmp3) deberá
tomar valores entre 30.1 y 35 a intervalos de 0.1; y la señal ( vector dat3) deberá
responder a la ecuación de la recta x=20t-500. Esto lo conseguimos con las siguientes
sentencias de Matlab:
>> tmp3 = [30.1:0.1:35]

>> dat3 = 20*tmp3-500
El último paso es crear las variables tmp y dat como concatenación de las variables
tmp1, tmp2, tmp3 y dat1, dat2, dat3 disponibles. Dado que tanto tmp como dat
deben ser vectores columna, se deberán combinar los vectores y trasponer el
resultado, tal y como indican las siguientes sentencias de Matlab (atención a la comilla
final que indica la operación de trasposición):
>> tmp = [tmp1, tmp2, tmp3]

>> dat = [dat1, dat2, dat3]
Para comprobar que la señal se ha generado correctamente, lo más sencillo es utilizar

el comando plot de Matlab: la instrucción plot(tmp,dat) debería producir el siguiente
resultado:

Una vez comprobado que las variables tmp y dat contienen valores correctos,
generaremos a partir de ellas la variable fuerza necesaria para el bloque From
Workspace. La primera columna deberá contener los tiempos y la segunda los datos,
con lo que las instrucciones Matlab a utilizar serán:
>> fuerza(:,1)=tmp
>> fuerza(:,2)=dat
Con la variable fuerza creada, es posible lanzar la simulación. Dado que la señal está
definida durante 35 segundos, el tiempo de simulación se fijará exactamente en 35
segundos (desde la opción Simulation parameters del menú Simulation). El
resultado sobre el osciloscopio se muestra en la figura siguiente, y se debe apreciar
cómo la salida del sistema (posición) reproduce aproximadamente los valores de la
entrada (fuerza aplicada) con pequeñas oscilaciones.


Ejercicio: COMPORTAMIENTO DINÁMICO DE UN SISTEMA
Trabajaremos sobre un montaje ya estudiado en otro laboratorio: el servomecanismo

de posición accionado mediante un motor eléctrico. Recordemos que el objetivo es
actuar sobre el motor o aplicar al motor una tensión v(t) de forma que la pieza móvil se
sitúe en el lugar x(t) que se le indique mediante la referencia x0(t):
El diagrama de bloques correspondiente a este montaje es el siguiente:
Donde los valores de las constantes del sistema son las siguientes:
 R = 1.25 (resistencia de los devanados del motor).
 I = 0.8 (momento de inercia del conjunto).
 B = 0.5 (coeficiente de rozamiento viscoso).
 KP = 3 (constante de par del motor).
 KV = 0.01 (constante de velocidad del motor).
La constante KC representa el controlador empleado: la acción sobre el motor v(t) será

proporcional al error o diferencia entre la referencia x0(t) y la posición real x(t). La
constante de proporcionalidad es KC, y ajustar el controlador equivale en este caso a
elegir el valor de KC que hace que el comportamiento del sistema sea el deseado.
Se pide:

 Representar el esquema anterior en Simulink y obtener la respuesta ante

escalón de 20 unidades para los 3 valores siguientes de la constante KC:
KC =0.1, KC =0.3 y KC =0.5. El sistema debe simularse durante 30 segundos,
y se debe ajustar el parámetro ‘max step size’ a 0.01 segundos para que la
representación de la respuesta sea suficientemente detallada.
 Anotar para cada valor de KC la sobre oscilación que se produce en la señal de

salida x(t) (sobreoscilación = exceso porcentual sobre el valor final) y el tiempo
de pico (tiempo de pico = instante en el que se alcanza el máximo valor). Debe
comprobarse cómo la sobreoscilación aumenta al aumentar KC y el tiempo de
pico disminuye al aumentar KC.
 Sobre Simulink, probar con distintos valores de KC hasta encontrar uno que
cumpla las siguientes especificaciones:
 Sobreoscilación < 28%.

 Tiempo de pico < 4 segundos.
 Una vez elegido el valor adecuado para KC, calcular con la ayuda de Matlab la
situación de los polos del sistema resultante. Se deberán emplear las
siguientes instrucciones de Matlab: tf para crear las funciones de transferencia;
series y feedback para reducir el diagrama de bloques; y finalmente pole para
determinar la posición de los polos en el sistema resultante.
Solución
R = 1.25;
I = 0.8;
B = 0.5;
Kp = 3;
Kv = 0.01;
G = tf(1,[I,B])
sys1 = 1/R*Kp*G
sys2 = feedback(sys1,Kv)
sys3 = 0.25*sys2*tf(1,[1 0])
sys4 = feedback(sys3,1)
step(sys4,30)
t = 0:0.1:30;
n = length(t)
u = 20./ones(1,n)
lsim(sys4,u,t)
pole(sys4)

Para KC =0.1

Laboratorio 02

Cargado por

Copyright:

Formatos disponibles

Laboratorio 02

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

Laboratorio 02

Cargado por

Copyright:

Formatos disponibles

Control Clásico – Matlab y Simulink

Programas de Matlab y Simulink

Como ejemplo volveremos a trabajar sobre el esquema utilizado en la práctica anterior

Se realizará la simulación de este esquema con los siguientes parámetros:

Docente: Jorge Ortiz Porras 1

El resultado debería ser igual al que muestra la siguiente figura:

Nota: en caso de no obtener exactamente estos resultados es posible que la

Docente: Jorge Ortiz Porras 2

En la ventana de Configuration Properties deberemos seleccionar logging y en ella

Docente: Jorge Ortiz Porras 3

Será necesario indicar dos parámetros:

Docente: Jorge Ortiz Porras 4

representamos las primeras 10 filas de la variable podremos comprobar esto (nota: se

Docente: Jorge Ortiz Porras 5

Desde la ventana de comandos de Matlab, también será posible representar mediante

>> plot(out.posicion(:,1), out.posicion(:,2))

>> [maximo, indice] = max(out.posicion(:,2))

Docente: Jorge Ortiz Porras 6

2. Introducción de señales arbitrarias en simulink

Docente: Jorge Ortiz Porras 7

El único parámetro que nos interesa por el momento es el nombre de la variable de

Docente: Jorge Ortiz Porras 8

Esto lo conseguimos con las siguientes sentencias de Matlab:

>> tmp1 = [0:0.1:10]

>> tmp2 = [10.1:0.1:30]

>> tmp3 = [30.1:0.1:35]

>> tmp = [tmp1, tmp2, tmp3]

Para comprobar que la señal se ha generado correctamente, lo más sencillo es utilizar

Docente: Jorge Ortiz Porras 9

Docente: Jorge Ortiz Porras 10

Docente: Jorge Ortiz Porras 11

Ejercicio: COMPORTAMIENTO DINÁMICO DE UN SISTEMA

Trabajaremos sobre un montaje ya estudiado en otro laboratorio: el servomecanismo

El diagrama de bloques correspondiente a este montaje es el siguiente:

La constante KC representa el controlador empleado: la acción sobre el motor v(t) será

Docente: Jorge Ortiz Porras 12

 Representar el esquema anterior en Simulink y obtener la respuesta ante

 Anotar para cada valor de KC la sobre oscilación que se produce en la señal de

 Sobreoscilación < 28%.

Docente: Jorge Ortiz Porras 13

Docente: Jorge Ortiz Porras 14

También podría gustarte