Descargar

Introducción a la identificación de sistemas (página 2)

Enviado por Pablo Turmero


Partes: 1, 2
edu.red

Consideraciones prácticas sobre elproceso de Identificación. Elección de las señales a elegir. Elección del tipo de entrada. Elección del periodo de muestreo. Elección del número de muestras a tomar.

DE LA OBTENCIÓN DE LOS DATOS:

edu.red

Consideraciones prácticas sobre elproceso de Identificación. Eliminación de perturbaciones de alta frecuencia. Eliminación de datos erróneos. Tratamiento de valores en continua (remover medias).

DEL PRE TRATAMIENTO DE LOS DATOS:

edu.red

Consideraciones prácticas sobre elproceso de Identificación. Validación en base a la aplicación del modelo (simulación, predicción, diseño de un controlador) . Comprobación de parámetros físicos. Coherencia con el comportamiento de entrada-salida. Reducción del modelo. Simulación. Análisis de residuos.

DE LA VALIDACIÓN DEL MODELO:

edu.red

Problemática existente.

Un paso importante en el proceso de identificación consiste en preparar y ejecutar un experimento que permita obtener información sobre la dinámica total del sistema. Uno de los problemas al que nos enfrentamos al intentar aplicar la identificación de sistemas a nuestro proceso es la limitante de no poder visitar la planta regularmente para la toma de datos, aplicar señales de entrada y registrar los datos de las salidas para así realizar todas las pruebas de identificación .

edu.red

Problemática existente.

Al sacar nuestro modelo matemático del proceso obtendremos un equivalente de lo que sería trabajar con la planta real, con lo que esperamos tener resultados similares y que se aproximen mucho a la realidad.

edu.red

Ventajas del modelado en la identificación. Predecir el comportamiento del sistema y poder hacer cambios al mismo sin que tenga que experimentarse físicamente sobre él.

edu.red

Modelamiento matemático de procesos. Representa una alternativa al diseñador, cuando en la identificación de sistemas es imposible obtener los datos de entrada/salida porque no se dispone físicamente de la planta. El modelo base es una ecuación matemática o un conjunto de ellas que describen el comportamiento del sistema.

edu.red

Aplicación práctica: Modelo matemático base para el sistema de control de temperatura en uno de los cuartos térmicos de una industria cervecera en donde se realiza la refrigeración de la cebada.

edu.red

Diagrama del proceso.

edu.red

Modelado matemático de la planta. La condición de equilibrio de los sistemas térmicos establece que el calor administrado a un sistema es igual al calor almacenado por el sistema más el calor liberado por el sistema y más las pérdidas del mismo.

edu.red

Consideraciones generales de la dinámica del sistema. Transferencia de calor hacia las paredes. O carga de Fuga, es una medición del calor que fluye a través de las paredes del espacio refrigerado del exterior hacia el interior.

edu.red

Consideraciones generales de la dinámica del sistema. Carga del producto o calor cedido por el producto. Cuando el producto entra al espacio de almacenamiento a temperatura mayor que la que se tiene dentro del espacio refrigerado, el producto cederá calor al dicho espacio hasta la temperatura que se tiene en el espacio.

QPRODUCTO = m Ce (? T)

edu.red

Consideraciones generales de la dinámica del sistema. Calor de respiración. Calor de transpiración, es el calor eliminado por los alimentos durante su almacenamiento y refrigeración.

QRESPIRACIÓN = m * Calor de respiración * 24 Hr

edu.red

Consideraciones generales de la dinámica del sistema. Cambios de aire. Cargas varias (personas, alumbrado, equipos eléctricos, etc.).

edu.red

Consideraciones generales de la dinámica del sistema.

BALANCES DE ENERGÍA. Balance de energía en el saladin. Balance de energía en el radiador.

PÉRDIDAS DEL SISTEMA. Transferencia de calor hacia las paredes. Calor de respiración del producto.

edu.red

Consideraciones generales de la dinámica del sistema.

DESPRECIABLES. Cambios de aire. Cargas varias.

PERTURBACIONES. Calor producido por el volteador mecánico que separa los granos. Cambios en la temperatura externa del cuarto.

edu.red

Modelo matemático base. Balance de energía en el radiador.

Tomando las variables de desviación:

edu.red

Modelo matemático base. La ecuación en términos de variables de desviación es la siguiente:

edu.red

Modelo matemático base.

De donde:

faire= 10.41 m3/s, flujo de aire a la salida del ventilador que ingresa al cuarto térmico. ?agua= 1000Kg/m3 , densidad del agua. ?aire= 1.22Kg/m3 , densidad del aire.     4186 J/Kg°C , calor específico del agua. 1004.64 J/Kg°C calor específico del aire.

 

edu.red

Modelo matemático base. Balance de energía en el saladin. Flujo de energía entrada – Flujo de energía salida – Pérdidas = Acumulación de energía. Tomando variables de desviación:

edu.red

Modelo matemático base. Aplicando la transformada de Laplace

edu.red

Modelo matemático base. De donde:

faire= 10.41 m3/s , flujo de aire a la salida del ventilador que ingresa al cuarto térmico. ?aire= 1.22Kg/m3 , densidad del aire.   Ce-aire= 1004.64 J/Kg°C, calor específico del aire. TV= Cambio en la temperatura del aire frío a la salida del radiador. TS= Cambio en la temperatura del sensor dentro del saladin.

edu.red

Modelo matemático base.

Qpared = Pérdida de calor hacia las paredes. Qtecho = Pérdida de calor hacia el techo. Qpiso = Pérdida de calor hacia el piso. Qrespir = Pérdida de calor producida por la cebada al respirar. m = Masa total de la cebada dentro del saladin.

edu.red

Modelo matematico base. Balance de energía total.

La función de transferencia general del sistema es entonces:

edu.red

Modelo matemático base. Entonces la ecuación que describe la dinámica del proceso considerando las perturbaciones debidas a los cambios de temperatura externa y la que entra al ventilador es la siguiente:  

edu.red

Modelo matemático del proceso simulado en Simulink de Matlab. Figura : Modelo de la planta real simulado en Simulink.

edu.red

Modelo matemático del proceso simulado en Simulink de Matlab. Figura : Modelo en Simulink de la función de transferencia del proceso.

edu.red

Validación de la respuesta en base al modelo matemático obtenido comparada con la del proceso real. Figura: Evolución de la temperatura en función del tiempo en el proceso real.

edu.red

Validación de la respuesta en base al modelo matemático obtenido comparada con la del proceso real. Figura: Evolución de la temperatura en función del tiempo en el proceso simulado. Cambio de un grado de la temperatura ocurre en 20 horas

edu.red

Identificación de sistemas. La identificación se simplifica utilizando una herramienta en Matlab llamada System Identification. Esta tiene una serie de funciones programadas que sirven para realizar las siguientes operaciones: Carga archivo de datos, es decir los datos registrados en un osciloscopio. Organiza y filtra los datos, removiendo los valores medios y filtrando el ruido en los mismos. Adicionalmente separa datos para identificación y datos para validación. Realiza proceso de validación, comparando gráficamente los resultados. Entrega análisis de correlación de los residuos con la entrada.

edu.red

Diseño de la señal de entrada haciendo uso de Matlab.

edu.red

La PRBS es una señal de dos niveles consistente en una sucesión de impulsos rectangulares modulados en su ancho, de forma que esta se aproxime a un ruido blanco con un contenido rico de frecuencias.

edu.red

Análisis de auto-correlación es idéntico al de ruido blanco, lo que nos garantiza un barrido alrededor de las frecuencias de interés.

edu.red

La entrada PRBS escogida considera las frecuencias en las que trabaja nuestro proceso.

edu.red

Recolección de datos a la entrada y salida del proceso a identificar. La señal de entrada PRBS diseñada es importada al modelo del proceso en Simulink. Luego se compila para ver el comportamiento de las salidas.

edu.red

Datos de entrada y salida del proceso.

edu.red

Identificación del sistema haciendo uso de la herramienta ident de matlab.

edu.red

Datos de entrada y salida del proceso que fueron importados a ident.

edu.red

Datos de entrada y salida tratados.

edu.red

Datos para identificación y validación.

edu.red

METODOS DE IDENTIFICACION Dependiendo del tipo de modelo obtenido:

Métodos no paramétricos: modelos no paramétricos. análisis de la respuesta transitoria, análisis de la respuesta en frecuencia, análisis de la correlación, análisis espectral.

Métodos paramétricos: modelos paramétricos. Estructura del modelo. Criterio de ajuste de parámetros. Estimación de los parámetros.

edu.red

Identificación no paramétrica vía análisis de correlación Este método es muy apropiado para obtener una idea rápida de la relación entre distintas señales del sistema, retardos, constantes de tiempo y ganancias estáticas del mismo.

edu.red

Con el análisis de correlación obtenemos los coeficientes del modelo respuesta al impulso finito (FIR), que relaciona la entrada con la salida.

Consideraciones antes del análisis de correlación: La series de tiempo de la entrada y la salida deben ser estacionarias o por lo menos cuasi-estacionarias. La entrada y el ruido no deben estar correlacionados.

Análisis de Correlación.

edu.red

La solución para encontrar los parámetros es:

Donde:

edu.red

Si la señal de entrada es una señal de ruido blanco con varianza entonces la matriz se simplifica:

Y el cálculo de los parámetros queda de la siguiente forma:

edu.red

Este resultado motiva a usar señales de entrada pre-blanqueadas en el análisis de correlación.

Los requisitos para el análisis de correlación se lo ha cumplido ya que nuestros datos de entrada son estacionarios y la señal de ruido se ha escogido independiente de la entrada. Por lo que ahora haremos el análisis de correlación de nuestro proceso.

edu.red

Del análisis de correlación obtenemos el modelo FIR:

1*y(t)+0*y(t-T)+…+ 0*y(t-nT)=9.1539*u(t) -13.5640*u(t-T)+… -4.7292*u(t-45*T) + e(t)

Modelo FIR con 45 parámetros.

edu.red

Observamos que es un modelo FIR con 45 parámetros. Al ser una identificación no paramétrica los parámetros no representan parámetros de nuestro modelo real y tampoco se puede obtener caracteristicas dinámicas y estaticas de nuestro sistema a simple vista. Por lo que a continuación obtendremos la respuesta transiente para determinar y tener una idea de las caracteristicas de nuestro sistema.

1*y(t)+0*y(t-T)+…+ 0*y(t-nT)=9.1539*u(t) -13.5640*u(t-T)+… -4.7292*u(t-45*T) + e(t)

edu.red

Con estas caracterices podemos modelar un sistema de primer orden. Pero como habíamos mencionado antes los parámetros estimados no necesariamente representan los parámetros del proceso real.

En esta respuesta podemos ver que la relación Señal-Ruido es alta que es lo contrario si solo introducimos una señal step al proceso para determinar su respuesta transiente. Análisis de Correlación. Respuesta Transiente.

edu.red

En la siguiente respuesta de frecuencia vemos que existe un alto contenido de varianza a altas frecuencias y bias bajo. Por lo que a altas frecuencias tendremos problemas. Para corregir esto usaremos una ventana de truncamiento.

Análisis de Correlación. Respuesta en Frecuencia.

edu.red

Una ventana de truncamiento común es la de TUKEY:

Donde M es el parámetro de truncamiento de la ventana. Mientras se incrementa M decrece el bias en la estimación y la varianza asociada con la estimación se incrementa. Por lo que escoger una ventana de truncamiento y un parametro de truncamiento que me permita un equilibrio entre “fidelidad” (bias) y “estabilidad” (varianza).

edu.red

En los gráficos observamos que se ha llegado a un compromiso entre la fidelidad y la estabilidad del sistema. Una vez hecho esto se puede determinar una relación para la variable de entrada y salida. Ventana TUKEY y parámetro truncamiento

edu.red

IDENTIFICACIÓN PARAMÉTRICA ELECCION DE LA ESTRUCTURA: Si no se tiene conocimiento previo del proceso a identificar la única forma de saber que estructura vamos a elegir es observando cual se adapta y describe mejor los datos de entrada y salida. Nosotros hemos probado algunas estructuras y hemos elegido a tres candidatas.

Modelo OE (output error) Modelo ARX. Modelo FIR.

edu.red

Estructura ARX:

Para determinar los ordenes de A y B en la estructura ARX primero sobre-parametrizamos el orden de los dos polinomios con esto se puede aproximar arbitrariamente bien cualquier sistema lineal ya que el bias es muy bajo. Pero nosotros debemos buscar el mínimo orden de los polinomios tal que se llegue a un equilibrio entre la “fidelidad” y la “varianza”.

edu.red

El modelo ARX con sus respectivos ordenes para los polinomios que hemos elegido es:

ARX(112) del cual: El orden del polinomio A es 1. El orden del polinomio B es 1. El retardo es 2, lo cual indica un retardo de una vez el tiempo de muestreo.

edu.red

MODELO ARX:

edu.red

Estructura OE:

Una de las ventajas de esta estructura es que independientemente se parametriza el modelo de la planta y del ruido. Y tampoco requiere una estructura para el modelo del ruido.

edu.red

El modelo OE con sus respectivos ordenes para los polinomios que hemos elegido es:

OE(112) del cual: El orden del polinomio B es 1. El orden del polinomio F es 1. El retardo es 2, lo cual indica un retardo de una vez el tiempo de muestreo.

edu.red

MODELO OE:

edu.red

Estructura FIR:

El orden apropiado para el polinomio B viene dado por el tiempo de muestreo y el tiempo de estabilización del proceso, el resultado usualmente es alto. Cuando la entrada esta constantemente excitando al proceso, el ruido es estacionario y la entrada y el ruido no están correlacionados los coeficientes de respuesta al impulso finito estimado no tendrán bias pero tendrán alta varianza. Se puede llegar a un equilibrio entre bias y varianza realizando operaciones de suavizado.

edu.red

MODELO FIR:

edu.red

CRITERIO DE AJUSTE DE PARAMETROS.

Para la estimación de los parámetros de los polinomios de los modelos antes mencionados se utilizó el método de mínimos cuadrados. El cual, como se muestra a continuación tiene como objetivo estimar los parámetros del modelo tal que la suma cuadrática del error entre los datos del modelo estimado y el real tienda a cero.

edu.red

VALIDACIÓN DE LOS MODELOS:

Análisis de polos y ceros. Análisis de residuos. Model OUTPUT.

edu.red

Análisis de polos y ceros. Polos y ceros modelo OE. Polos y ceros modelo ARX. Polos y ceros modelo FIR. Si hubiera superposición de polos y ceros entonces se podría reducir el orden de los modelos.

edu.red

Análisis de residuos. Se conocen como residuos de un sistema a los errores de predicción obtenidos según la expresión:

siendo q el vector de parámetros del modelo, y(t) la respuesta real del sistema e ye(t) la respuesta estimada por el modelo para la misma entrada.

Idealmente, estos residuos deben ser independientes de la entrada. Si no sucede así, significa que hay componentes en e(t) que proceden de la entrada u(t), lo cual a su vez significa que el modelo no es capaz de describir completamente la dinámica del sistema.

edu.red

Análisis de residuos y correlación modelo OE Si hay componentes en e(t) que proceden de la entrada u(t), significa que el modelo no es capaz de describir completamente la dinámica del sistema.

edu.red

Análisis de residuos y correlación del modelo ARX

edu.red

Análisis de residuos y correlación modelo FIR

edu.red

Model OUTPUT Con el model OutPut vemos las respuestas de los tres sistemas a la entrada PRBS que hemos diseñado, las tres se comparan con la señal de salida “real“ y se escoge la que mejor se adapta y describa la dinámica del proceso.

Para saber cuál es el modelo que mejor describe y se adapta a los datos de entrada y salida del proceso real, cogeremos el que tiene el BEST FIT más alto. En nuestro caso el modelo con mayor best fit es el OE, con un best fit de 92.99. Esto significa que la variación en la entrada explica 93.0 % de la variación en la salida.

edu.red

Model OUTPUT de los tres modelos

edu.red

Diseño del controlador. En los sistemas en lazo cerrado la variable controlada retroalimenta la entrada, comparándose continuamente el error existente entre la referencia y esta variable controlada. Nuestro objetivo es diseñar un controlador discreto por lo que para eso diseñaremos primero un controlador continuo y luego lo discretizaremos.

edu.red

Proceso sin Controlador.

edu.red

Respuesta transitoria en lazo cerrado para el sistema sin controlador.

edu.red

Diseño del controlador continuo PI Real. De donde tenemos que Kp=19. La Función de transferencia del controlador PI es la siguiente: Ahora vamos a calcular la posición del cero y el polo del regulador para que se cumpla la restricción del error en régimen permanente que nos dicen (5%).

La nueva constante de error en posición, al introducir el regulador PI es:

Ahora bien, la constante de error de posición que necesitamos es:

edu.red

Luego:

z=8.184*10-05

De esta forma:

p=(1/19)*z=0.000004307

Por tanto, el regulador PI que vamos a utilizar es el siguiente:

edu.red

Proceso con Controlador Continuo.

edu.red

Respuesta transitoria en lazo cerrado para el sistema continuo con controlador.

edu.red

Diseño del Controlador Discreto. Los parámetros KP, KI y KD del controlador PID discreto están relacionados con los parámetros Kp, Ki y Kd del controlador PID continuo, de la siguiente forma:

edu.red

Con estas relaciones determinamos los parámetros del controlador discreto, también utilizamos un muestreador a T=1000 segundos.

Los parámetros para el controlador discreto son: KP=0.9625 KI=0.075 KD=0

edu.red

Al discretizar todos los componentes de proceso nos queda de la siguiente forma:

Proceso dincretizado.

edu.red

Respuesta transitoria en lazo cerrado para el sistema discreto con controlador.

edu.red

Vemos que el efecto del muestreo o discretización, es desestabilizante, reflejándose este en un sobrepaso mayor y un tiempo de establecimiento mas largo. Es importante decir que el aumento en la taza de muestreo disminuye este efecto y la reducción de esta, aumenta la inestabilidad del sistema.

Sin embargo debemos tomar en cuenta que la tasa de muestreo no es una variable de libre elección, y depende de factores tales como la velocidad de procesamiento, la cantidad de lazos de control utilizados o la velocidad del proceso a controlar. Una regla práctica indica que se debe tomar una frecuencia de muestreo con un valor de 10 veces el ancho de banda del sistema de control en lazo cerrado.

edu.red

Conclusiones. Realizar una identificación en sistemas reales implica grandes costos debido a los paros de producción que requiere la experimentación, por tanto, para fines académicos, es de gran ayuda trabajar con un modelo matemático base que represente la dinámica del proceso. Se demostró que el proceso de identificación de sistemas nos ofrece una alternativa práctica y eficaz para determinar un modelo satisfactorio en control y su aplicación no sólo se limita al control automático.

edu.red

Conclusiones. El proceso de identificación va desde el diseño del experimento, la adquisición y tratamiento de los datos, la elección de la estructura del modelo y la selección de los parámetros, hasta concluir con la fase de validación del modelo. Para conseguir una coherencia de estimación, todos los pasos citados deben tratarse con el mismo orden y rigor; y descuidos en consideraciones sobre ellos, provocarán inexactitudes y comportamientos erróneos del modelo final.

edu.red

Recomendaciones. Validar varios modelos antes de escoger el modelo final, considerando la mayor cantidad de análisis posibles (análisis de residuos, análisis de las entradas-salidas, salidas-perturbaciones, análisis de estabilidad en base al criterio de polos y ceros, etc.).

edu.red

Recomendaciones. Revisar datos de entrada-salida para asegurar que estos proporcionen información de la dinámica del sistema. Revisar si la estructura escogida proporciona una buena descripción del modelo. Verificar si el ajuste de parámetros es el adecuado.

Partes: 1, 2
 Página anterior Volver al principio del trabajoPágina siguiente