Segundo Orden

Control de Procesos. Julio A. Romero. Dpto. Tecnología

Identificación experimental de sistemas dinámicos por

métodos gráficos usando Matlab.

Introducción

Los métodos gráficos para la identificación de sistemas tienen la ventaja de que, dada su

sencillez, para su aplicación sólo se necesita “lápiz y papel”. Hoy en día, sin embargo, son

habituales los ordenadores con posibilidades gráficas potentes, y además la existencia de

programas especializados para el manejo de datos. En este contexto es importante conocer

como aprovechar las potencialidades de estas herramientas computacionales para la

aplicación de estos métodos de identificación.

Esta guía tiene el objetivo de brindar los procedimientos básicos a seguir para la

aplicación de los métodos gráficos de identificación experimental utilizando el programa

Matlab. Para ello se han desarrollado los siguientes ejemplos:

1. Obtención de un modelo de primer orden.

2. Identificación de un modelo de segundo orden con polos iguales mediante el

método de Strecj.

3. Identificación de un modelo de segundo orden con polos diferentes mediante el

método de gráficas logarítmicas.

Ejemplo 1

Se quiere obtener el modelo de un sistema cuya respuesta ante un escalón unitario está en

el archivo ensayo.dat.

Cargar los datos al Workspace de Matlab:

>> load ensayo.dat

mediante la instrucción “who” podemos comprobar que los datos han sido cargados desde

el archivo:

>> who

Your variables are:

ensayo

La variable ensayo es una matriz de 3 columnas: la primera columna contiene los valores

de tiempo t y la segunda los valores de la entrada u(t)y la tercera columna los valores de

la salida del sistema y(t).

Obtener una gráfica de los datos experimentales:

>> plot(ensayo(:,1),ensayo(:,3))

Con esta instrucción Matlab genera una figura (Figure No. 1) que contiene la gráfica de la

respuesta del sistema: y(t) vs t

Para facilitar la búsqueda de valores en la gráfica añadir divisiones mediante la

instrucción:

>> grid

1

Control de Procesos. Julio A. Romero. Dpto. Tecnología

Sobre esta gráfica se aplican los métodos de identificación gráficos.

Identificación de un modelo de 1er orden

Supongamos que queremos aproximar el comportamiento por un sistema de 1 orden.

G s = k

τs1

Hemos de obtener los parámetros: ganancia estática y constante de tiempo.

Ganancia estática k:

k=

Δy ∞

Δu

de la gráfica Δy ∞=1.5 y conociendo que la entrada ha sido un escalón unitario

Δu=1 . Por tanto k=1.5 .

Constante de tiempo t

Calcular el 63% del valor final de estado estable de la salida 0.63 Δy ∞ :

>> y_tau=0.63*1.5

y_tau =

0.9450

Para este valor se obtiene de la gráfica por inspección visual el valor de t. Para obtener un

valor más preciso se puede utilizar la instrucción ginput, con la cual es posible obtener las

coordenadas (x,y) de un punto seleccionado con el ratón en un gráfico. Una vez

introducida la instrucción en la línea de comando, se selecciona en el gráfico el punto del

cual deseamos conocer las coordenadas:

hacemos clic con el ratón y aparecen las coordenadas del punto seleccionado

2

Hacer coincidir

aproximadamente

con el valor de y_tau

Control de Procesos. Julio A. Romero. Dpto. Tecnología

>>tau=ginput(1)

tau =

2.5461 0.9287

de donde tenemos que τ≈2.55 .

Verificación del modelo

Una vez obtenidos los dos parámetros del modelo podemos comprobar la similitud de la

respuesta del modelo con los datos experimentales.

Primero crearemos el modelo con los parámetros calculados con la instrucción tf, que

permite definir una función de transferencia a partir de dos vectores que contienen los

coeficientes de los polinomios del numerado y denominador en potencias decrecientes de

s:

>> G0=tf(1.5,[2.55 1])

Transfer function:

1.5

----------

2.55 s + 1

Luego

...