Modelos Matemáticos de Sistemas de Control

download Modelos Matemáticos de Sistemas de Control

of 71

Transcript of Modelos Matemáticos de Sistemas de Control

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    1/71

    1 - 1

    FUNDAMENTOS Y MODELOS

    MATEMTICOS DE SISTEMAS DE CONTROL

    INTRODUCCIN

    En las ltimas dcadas, los sistemas de control han desempeado un rol vital en el

    desarrollo y avance tecnolgico de la ciencia y la ingeniera, y por ende de la sociedad moderna.

    En la ltima dcada se han convertido en componentes esenciales en el control de vehculos

    espaciales, sistemas robticos, y ms indispensables an, en el procesamiento de productos

    alimenticios, combustibles, industria petroqumica, generacin y distribucin de energa elctrica,

    distribucin y tratamiento de aguas residuales y servidas, electrnica de automviles,

    electrodomsticos, etc. Los sistemas de control estn involucrados de manera implcita en todos

    los aspectos de la vida diaria, siendo su objetivo fundamental el de mantener los ms altos

    estndares de calidad de los productos (composicin, pureza, color, etc.), manteniendo los niveles

    de produccin a mnimo costo y proporcionando adems las condiciones de trabajo adecuadas

    para satisfacer la seguridad industrial y ambiental, con la menor intervencin del ser humano.

    En este captulo se tratarn los aspectos relacionados con el desarrollo del modelo

    matemtico, las estructuras tpicas, los componentes fsicos y dems elementos indispensables

    para el anlisis de un sistema de control en tiempo continuo. Utilizando como herramienta la

    transformada de Laplace, se aplicar el concepto de funcin de transferencia para modelar el

    comportamiento dinmico del sistema, estableciendo su relacin con la respuesta impulso. La

    representacin grfica del sistema utilizando diagrama de bloques, facilitar el modelado de cada

    uno de los componentes del sistema de control. A partir del diagrama de bloques se desarrollar

    el grfico de flujo de seales, que permitir el uso de la frmula de ganancia de Mason, como

    herramienta para evaluar la funcin de transferencia de un sistema con varios lazos de control.

    Finalmente, utilizando el concepto de variables fsicas fundamentales para establecer

    analoga entre sistemas fsicos, se desarrollar el modelo matemtico de un conjunto de sistemas

    tpicos, asociados con sistemas elctricos, sistemas mecnicos de traslacin, sistemas mecnicos

    de rotacin, sistemas electromecnicos, sistemas trmicos y sistemas hidrulicos.

    1.1 CONCEPTOS BSICOS DE SISTEMAS DE CONTROL

    En esta seccin se presentar una visin global y un conjunto de aspectos fundamentales

    relacionados con el propsito, definicin, componentes, estrategias, seales caractersticas

    y campos de aplicacin de los sistemas de control.

    1

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    2/71

    1-2 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    As mismo se formular el problema de control, con el propsito de identificar los

    elementos que deben considerarse en el diseo del controlador, as como otros aspectos

    adicionales que intervienen en la solucin de este problema.

    Definiciones bsicas

    Aunque existen diferentes definiciones clsicas relacionadas con el objetivo y propsito de

    un sistema de control, [Dorf05], [Ogata03a], [Kuo95], [Franklin91], [Phillips00], la

    siguiente definicin incluye dos condiciones que caracterizan a un sistema de control:

    DEFINICION 1.1 Propsito del sistema de control

    Conjunto de componentes interconectados, de modo que puedan ser

    comandados, o regulados por s mismos o por otro sistema en forma

    automtica, para lograr una condicin deseadade una variable fsica.

    Condiciones mnimas de un sistema de control:

    Segn esta definicin, existen dos condiciones mnimasque debe satisfacer un sistema de

    control: la primera se refiere a la capacidad de regulacin de sus componentes

    interconectados, para responder a las especificaciones de uso de acuerdo con la variable

    fsica a controlar o variable controlada. La segunda establece que la regulacin debe ser

    automtica, lo cual implica que no es necesaria la intervencin del ser humano.

    Podramos imaginarnos las consecuencias de tener a una persona ajustando manualmente la

    vlvula de vapor de un sistema de control de temperatura que utiliza un intercambiador decalor. En primer lugar, el alto nivel de riesgo por descuido del operador podra elevar la

    temperatura a valores peligrosos. En segundo lugar, la calidad en la regulacin del sistema

    sera muy pobre, por la imposibilidad de garantizar que el operador est pendiente de las

    variaciones en la temperatura del vapor, todas las horas del da y todos los das del ao.

    Diagrama funcional y modelo del sistema:

    Permite identificar la relacin causa efecto asociada con la seal fsica a regular ovariable controlada, y la condicin esperada para esta variable o valor deseado(setpoint).

    Esta identificacin permite formular la relacin entrada salidadel sistema tal como semuestra en la figura 1.1, la cual establece a su vez el propsito del sistema de control.

    SISTEMADE CONTROL

    Valor deseado

    Causa (entrada)

    Variable controlada

    Efecto (salida)

    Figura 1.1Diagrama funcionalasociado con elpropsito delsistema de control.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    3/71

    1.1 CONCEPTOS BSICOS DE SISTEMAS DE CONTROL 1-3

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    El planteamiento anterior es consistente con la siguiente definicin [Carlson98]:

    DEFINICION 1.2 Sistema dinmico

    Proceso en el cual existe una relacin causa efecto y es posibleformular en forma algebraica o grfica una relacin entrada salida,para evaluar su comportamiento en el tiempo.

    La formulacin de la relacin entrada salida conduce al desarrollo del modelo delsistema, el cual puede ser analtico o grfico y es el elemento bsico para el anlisis y

    diseo del sistema de control.

    Anlisis y diseo:

    El problema de anlisis se refiere a la evaluacin de la respuesta dinmica del sistema

    (salida) para una entrada conocida, referida generalmente como seal de prueba,

    asumiendo que se conoce el modelo del sistema. Por otro lado, el problema de diseo se

    refiere a determinar el modelo del sistema para satisfacer condiciones especficas de una

    relacin entrada salida, referida normalmente como requerimientos de diseo.

    Sistemas SISO y sistemas MIMO:

    Aunque el diagrama de la figura 1.1 muestra un sistema de 1-entrada y 1-salida o SISO

    (Single-Input-Single-Output), la mayor parte de las aplicaciones prcticas corresponden a

    sistemas multivariables o MIMO (Multiple-Input-Multiple-Output), cuyo diagrama

    funcional se muestra en la figura 1.2. La flecha doble se usa para indicar que existen varias

    entradas (causas) y salidas (varios efectos) en el proceso de regulacin del controlador.

    Un ejemplo tpico de un sistema de control multivariable es el caso del sistema de

    regulacin de velocidad de un automvil, mostrado en la figura 1.3, donde la velocidad

    final depende del par mecnico o torque mT ejercido por la inercia total del vehculo y elngulo de posicin del acelerador

    acel fijado por el conductor.

    Figura 1.2Diagrama funcionalde un sistema desistema de controlmultivariable (MIMO).

    SISTEMA DECONTROL

    MULTIVARIABLE

    Entradas Salidas

    Valoresdeseados

    Variablescontroladas

    CONTROL DEVELOCIDADVEHICULO

    mT

    acel

    Velocidad

    Figura 1.3Diagrama funcional delsistema de regulacinde velocidad de unvehculo.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    4/71

    1-4 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Campos de aplicacin de los sistemas de contro l

    Aunque existe una gran variedad de sistemas de control que se pueden encontrar en

    aplicaciones domsticas, comerciales e industriales, es frecuente clasificarlos en dos

    grandes categoras [Johnson02], segn la naturaleza de variable fsica a regular:

    - control de procesos

    - servomecanismos

    El control de procesostrata de forzar a que una variable fsica mantenga un valor constante

    en el tiempo, respecto de un valor deseadoo setpoint. De este modo el control de procesos

    se aplica generalmente en operaciones automticas de control de nivel, temperatura, flujo,

    presin, posicin, relacionadas con procesos domsticos, comerciales e industriales.

    Los servomecanismosobedecen a otro tipo de sistema de control, donde el objetivo es el

    seguimiento o rastreo de una seal fsica, forzando a que se mantenga cercana a valor

    especfico o target. El trmino servomecanismo se debe a la clase de componentes que

    utiliza para lograr el propsito del sistema de control. Ejemplos tpicos son: el control de

    posicin de una antena de un radar, el control de la direccin de un vehculo y el uso de

    robots en aplicaciones industriales y biomedicina, para lograr movimientos precisos en el

    espacio como una funcin del tiempo.

    Otras clasificaciones hacen referencia a aplicaciones ms especficas como el control

    secuencial, utilizado en sistemas electrodomsticos y en procesos de manufactura de

    productos que utilizan mquinas herramientas computarizadas. El control analgicodondela funcin de regulacin es realizada por dispositivos analgicos electrnicos, neumticos o

    hidrulicos y el control digitalque utiliza un microprocesador como unidad de control.

    1.2 ESTRUCTURAS TPICAS DE LOS SISTEMAS DE CONTROL

    En esta seccin se har una descripcin de las estrategias utilizadas para satisfacer el

    requisito de regulacin automticadel sistema de control, sus estructuras, componentes y

    caractersticas. Se har nfasis en elprincipio de control por realimentacin(feedback) por

    su caracterstica de regulacin del error del sistema de control, adems de otras efectos ensu comportamiento dinmico como: estabilidad, capacidad de rechazo a las perturbaciones

    y baja sensibilidadpor cambio en sus parmetros.

    Estrategias de control

    La condicin de regulacin automtica del sistema se puede lograr utilizando diversas

    estrategias de control. Estas estrategias se desarrollan a travs de esquemas de controlque

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    5/71

    1.2 ESTRUCTURAS TPICAS DE LOS SISTEMAS DE CONTROL 1-5

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    ofrecen caractersticas particulares. Los esquemas clsicos utilizados con mayor frecuencia

    en aplicaciones prcticas, son:

    - sistema de control de lazo abierto

    - sistema de control de lazo cerrado o realimentado (feedback)- sistema de control de accin precalculada

    - sistema de control en cascada

    Sistema de control de lazo abierto:

    La figura 1.4 muestra los elementos bsicos de un sistema de control de lazo abierto, donde

    se ha introducido el trmino de compensador, ampliamente utilizado en la teora clsica de

    sistemas de control para hacer referencia al controlador.

    Segn la definicin 1.1, el sistema debe ser capaz de regular la seal de salida o variable

    controlada ( )y t dentro de lmites aceptables del valor deseadoo setpoint, establecido por

    la seal de entrada ( )r t . La seal de control ( )m t es determinada por la accin del

    controlador o compensador y se encarga de ajustar el proceso o planta para garantizar quela variable controlada se mantenga cerca del valor deseado.

    La seal ( )m t acta sobre un componente del proceso referido como el elemento final de

    control (EFC), para regular su funcionamiento, tal como se muestra en la figura 1.4.

    Ejemplos tpicos del EFC son: vlvulas, fuentes de potencia, reguladores, servomotores,

    etc. En esquemas posteriores para efecto de simplificacin del esquema se omitir el EFC.

    La seal ( )p t en la figura 1.4 se utiliza para simular la presencia de perturbacionesen el

    proceso, entendida como una seal o seales que pueden modificar la variable controlada

    o salida del sistema ( )y t . Para cada sistema en particular es posible identificar este tipo de

    seales. Por ejemplo, la temperatura ambiente, la temperatura de entrada del fluido de

    control y la masa de fluido cuya temperatura se desea regular, son ejemplos tpicos de

    perturbacin en un sistema de control de temperatura de un proceso.

    Si el sistema de regulacin de velocidad de un vehculo mostrado en la figura 1.3 es de lazo

    abierto, la perturbacin puede ocurrir cuando se presenta una pendiente en la trayectoria del

    Figura 1.4Componentes de unsistema de control delazo abierto.

    COMPENSADOR OCONTROLADOR

    ( )r t ( )m t PROCESOO PLANTA

    ( )y t

    ( )p t

    EFC

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    6/71

    1-6 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    vehculo, la cual reduce de inmediato la velocidad . A menos que el conductor ajustemanualmenteel ngulo de posicin del acelerador

    acel (seal de referencia o setpoint), la

    variable controlada seguira la trayectoria mostrada en la figura 1.5. Luego, el sistema en

    lazo abierto no es regulado, dado que no es capaz de ajustar automticamente la entrada

    ( )r t para responder a las perturbaciones del proceso ( )p t .

    En este orden de ideas, las caractersticas del sistema de control de lazo abierto, son:

    - es simple y econmico.

    - no responde a cambios en la variable controlada por efecto de perturbaciones.

    - es un sistema no regulado.

    Sistema de control de lazo cerrado:

    Para lograr un sistema regulado, es necesario modificar el esquema de la figura 1.4,

    insertando un lazo de realimentacina travs del cual se pueda informar continuamente al

    controlador o compensador de los cambios que ocurren en la variable controlada por efectode perturbaciones en el proceso. Esta estrategia se logra con el esquema mostrado en la

    figura 1.6, conocido como sistema de control de lazo cerradoo control realimentado.

    En el esquema de la figura 1.6, los cambios en la variable controlada ( )y t por efecto de

    perturbaciones en el proceso ( )p t son transmitidos al controlador a travs de la seal de

    realimentacin ( )b t . El controlador, segn la magnitud de estos cambios ajusta

    aplicacin dem

    T

    t

    Figura 1.5Respuesta de unsistema de control delazo abierto ante unaperturbacin.

    Figura 1.6Sistema de control delazo cerrado o controlrealimentado.

    COMPENSADORCONTROLADOR

    ( )r t ( )m t PROCESOO PLANTA

    ( )y t

    ( )p t

    TRANSMISORMEDIDOR

    ( )b t

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    7/71

    1.2 ESTRUCTURAS TPICAS DE LOS SISTEMAS DE CONTROL 1-7

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    automticamente la seal de control ( )m t para lograr que la variable controlada ( )y t se

    mantenga cerca del valor deseado ( )r t ,

    Retomando el ejemplo del sistema de regulacin de velocidad del automvil de la figura

    1.3, podemos asumir que se instala un sensor de velocidad y un regulador que ajustaautomticamente el ngulo de posicin

    acel del acelerador, cada vez que ocurra un cambio

    en la velocidad del motor. La figura 1.7 muestra la respuesta dinmica de este sistema ante

    una perturbacin, originada por un aumento en torque mecnicom

    T .

    Las caractersticas de un sistema de control de lazo cerrado , se pueden resumir en:

    - transmite continuamente al controlador la informacin sobre el estado actual de la

    variable controlada ( )y t .

    - determina la accin de control ( )m t en funcin de los cambios de la variable

    controlada, respecto del valor deseado ( )r t .

    - es un sistema regulado, porque responde a los cambios en la variable controlada, porefecto de perturbaciones en el proceso.

    - ms complejo y costoso que un sistema de control de lazo abierto.

    - puede ser lento en la respuesta.

    El esquema de la figura 1.6 es el fundamento de los sistemas realimentados de controlo

    feedback y ser el modelo clsico a utilizar en los captulos posteriores relacionados con

    el anlisis y diseo del sistema de control. Sus caractersticas y propiedades sern

    analizadas con detalle ms adelante.

    Sistema de control por accin precalculada:

    Una forma de mejorar la velocidad de reaccin del sistema de control en lazo cerrado, es

    utilizar el esquema de control por accin precalculada, mostrado en la figura 1.8. En esta

    estrategia, el controlador recibe continuamente informacin relacionada con el estado

    actual de las variables de perturbacin del proceso ( )p t . De este modo, el sistema

    determina la accin de control ( )m t necesaria para ajustar el funcionamiento del proceso,

    aplicacin de

    mT

    t

    Figura 1.7

    Respuesta dinmica deun sistema de controlde lazo cerrado.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    8/71

    1-8 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    antes de que ocurran los cambios en la variable controlada ( )y t , mejorando la velocidad de

    respuestadel sistema de control. Esta caracterstica es fundamental en sistemas de control

    de temperatura, donde las constantes de tiempo del proceso son elevadas.

    Sin embargo, el clculo de la accin de control ( )m t es complejo y generalmente se recurre

    a registros histricos del comportamiento de la variable controlada ( )y t respecto de lasvariables de perturbacin ( )p t .

    Las caractersticas del control por accin precalculada, pueden resumirse en:

    - evala continuamente las variables de perturbacin del proceso, para determinar la

    accin de control ( )m t .

    - es ms rpido en la respuesta, que el sistema de control en lazo cerrado.

    - es ms complejo y ms costoso de implementar.

    - no utiliza realimentacin de la variable controlada.

    Sistema de control en cascada:

    El esquema tpico se muestra en la figura 1.9, donde se muestran 2 lazos de control que

    utilizan el principio de realimentacin. El propsito de este esquema es regular una de las

    variables de perturbacin que pueden tener mayor efecto en las alteraciones de la variable

    controlada del proceso. A partir de esta variable se establece el lazo secundario de control,

    logrando as minimizar su efecto sobre la variable controlada del sistema ( )y t .

    Un caso prctico de control en cascada es el sistema mostrado en la figura 1.18, que utiliza

    un intercambiador de calor para regular la temperatura de un fluido. El fluido de controles

    vapor y una de las variables de perturbacin que tienen mayor incidencia en la variable

    controlada (temperatura), son las posibles variaciones en el flujo de vapor aguas arriba de la

    vlvula de control.

    Una forma de minimizar el efecto de estas variaciones de flujo, es instalar un lazo

    secundario de control de flujo, que regule la cantidad de vapor que llega al intercambiador,

    COMPENSADORCONTROLADOR

    ( )r t ( )m t PROCESOO PLANTA

    ( )y t

    ( )p t

    Figura 1.8Sistema de control poraccin precalculada.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    9/71

    1.2 ESTRUCTURAS TPICAS DE LOS SISTEMAS DE CONTROL 1-9

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    tal como se muestra en la figura 1.9. Otra perturbacin, como la temperatura ambiente

    puede minimizarse usando revestimiento trmico apropiado en el intercambiador de calor.

    De acuerdo con la figura 1.9, el sistema de control en cascada utiliza un lazo de control

    primario, cuya accin de control ( )pm t se convierte en el setpoint del controlador

    secundario, el cual responder a su vez con una accin de control ( )s

    m t para regular la

    variable controlada o variable del proceso ( )y t . En aplicaciones prcticas el controlador

    primario se denomina como control maestro y el controlador secundario como control

    esclavoy el conjunto global como sistema de control maestro-esclavo.

    Caractersticas del sistema control realimentado (feedback)

    La razn de utilizar el principio de realimentacincomo estrategia de control se present

    en el sistema de regulacin de velocidad de un automvil, mostrado en la figura 1.3. En

    esta estrategia de control podemos identificar 3 operaciones bsicas que debe realizar el

    sistema de control:

    1. Detectar el valor actual de la variable controlada ( )y t , a travs del sistema de

    medicin y transmitirla al controlador.

    2. Comparar la seal de realimentacin ( )b t

    con el valor deseado ( )r t

    de la variablecontrolada. El resultado de esta comparacin establece el error del sistema, como

    ( ) ( ) ( )e t r t b t (1.1)

    3. A partir del error del sistema ( )e testablecer la accin de control ( )m t necesaria para

    corregir la desviacin de la variable controlada ( )y t .

    Figura 1.9Sistema de control encascada.

    CONTROLADORPRIMARIO

    ( )r t ( )pm t PROCESOO PLANTA

    ( )y t

    ( )p t

    MEDIDOR

    PRIMARIO

    CONTROLADORSECUNDARIO

    ( )s

    m t

    MEDIDORSECUNDARIO

    [Flujo]

    [Temperatura]

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    10/71

    1-10 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    La figura 1.10 muestra los componentes y las seales que intervienen en el desarrollo de las

    3 operaciones que realiza el sistema de control realimentado. La lnea punteada permite

    identificar al controlador o compensador mostrado en la figura 1.6. Esta similitud entre el

    diagrama de las figuras 1.10 y 1.6, hace que el sistema de control con realimentacin, se

    reconozca tambin como sistema de control de lazo cerrado.

    La estrategia para el anlisis del sistema de control mostrado en la figura 1.10 consiste en

    modelar la relacin entrada salida de cada componente por un bloque funcional, usandolas seales que se describen a continuacin:

    ( )y t : variable controlada o variable del proceso, la cual establece el propsito del

    sistema de control.

    ( )r t : valor deseado de la variable controlada ( )y t , seal de referencia, o setpoint.( )b t: valor medido de la variable controlada ( )y t o seal de realimentacin.

    ( )e t: seal de error, como una medida de la desviacin que sufre la variable controlada

    ( )y t , respecto del valor deseado ( )r t .

    ( )m t : seal de control, calculada a partir de la seal de error ( )e t, de acuerdo con el

    modo de accindel controlador.

    ( )p t: perturbacin del proceso, razn de uso del principio de realimentacin.

    El sumador mostrado en la figura 1.10 para evaluar la seal de error, se reconoce como el

    detector de error.

    El problema de control

    Es posible formular el problema de control, en trminos de la siguiente definicin:

    DEFINICION 1.3 El problema de control

    Controlar con un mnimo de precisin un proceso o planta, utilizando el

    principio de realimentacin, a travs del esquema de lazo cerrado.

    Figura 1.10Componentes delsistema de controlrealimentado o de lazocerrado.

    + MODO DECONTROL

    ( )r t ( )m t PROCESOO PLANTA

    ( )y t

    ( )p t

    MEDICIN

    ( )e t

    ( )b t

    CONTROLADOR

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    11/71

    1.2 ESTRUCTURAS TPICAS DE LOS SISTEMAS DE CONTROL 1-11

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    De este modo, el propsito de la realimentacin es el de minimizar el error, y su magnitud

    es una medida de la precisinlograda por el sistema de control. La figura 1.11 muestra la

    seal de error en dos sistemas de control que puede utilizarse para comparar su exactitud.

    La figura 1.12 muestra las 4 etapas [Bishop97] que se utilizan en la Ingeniera de Control

    para la solucin del problema.

    La descripcin detallada de cada etapa se presenta a continuacin:

    t

    ( )e t baja exactitudalta exactitud

    Figura 1.11Exactitud de dossistemas de control apartir de la seal deerror.

    Figura 1.12Fases en la solucindel problema decontrol.

    Propsito del sistema de control:- variables a ser reguladas- variables de perturbacin- requerimientos del sistema

    Modelo del sistema de control:- esquema de control- modelo del sensor y actuador- modelo del proceso o planta

    Diseo del sistema de control:- especificaciones de diseo- ajuste de parmetros- modelo del controlador

    Se cumplen especificaciones?SNo Documentacin

    del proyecto

    Verificacin de resultados:- comprobar especificaciones- anlisis de sensibilidad

    - rechazo a perturbaciones

    1

    2

    3

    4

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    12/71

    1-12 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Fase 1: Establecimiento del propsito del sistema de control

    En esta etapa es necesario identificar las seales a ser reguladas y las seales de

    perturbacin, las cuales establecen el propsito del sistema de control. As mismo es

    necesario formular los requerimientos del sistema en trminos de valores caractersticos

    esperados para la respuesta transitoria y permanente.

    Fase2: Desarrollo del modelo del sistema de control

    Esta fase es la que presenta mayor complejidad en la solucin del problema y se inicia

    estableciendo la estrategia de control y el esquema a ser utilizado: lazo abierto, lazo

    cerrado, cascada, accin precalculada, etc. donde la experiencia prctica del diseador es

    fundamental para lograr una estrategia de control sencilla pero efectiva, segn los

    requerimientos del sistema. Un segundo elemento a considerar en esta fase, es la seleccin

    de sensores para medicin de la seal de campo y de actuadores para modificar el proceso.

    A continuacin es necesario desarrollar el modelo del procesoo planta, del actuador y del

    sensor, aplicando criterios prcticos para lograr una abstraccin del modelo fsico,

    mediante el uso de elementos conceptuales de fsica, qumica, mecnica, etc. para lograr un

    modelo matemtico simplificado, pero que a su vez sea una adecuada representacin de los

    componentes fsicos del proceso o planta.

    Fase 3: Diseo del sistema de control

    Esta fase se inicia formulando las especificaciones de diseoa partir de lo requerimientos

    del sistema presentados en la fase 1. De acuerdo con el esquema de control seleccionado en

    la fase 2, es posible establecer el modelo matemtico del controlador o compensador a

    utilizar y a partir de este calcular el ajuste de sus parmetros, aplicando mtodos clsicos o

    modernos de diseo.

    Fase 4: Verificacin de resultados y documentacin del proyecto

    Una vez diseado el controlador, es necesario verificar el resultadoobtenido, evaluando la

    respuesta dinmica del sistema a la luz de las especificaciones de diseo. En esta fase

    generalmente se recurre al uso de herramientas de simulacin. Adems de verificar el

    cumplimiento de las especificaciones de diseo, se deben evaluar otros aspectos

    relacionados con la sensibilidad por cambio en parmetrosdel sistema y el rechazo a las

    perturbaciones.

    Si el resultado del diseo no es satisfactorio es necesario retornar a la fase 2, para la

    revisin del esquema seleccionado y de las simplificaciones hechas en el desarrollo del

    modelo de los componentes del proceso, hasta lograr un resultado que se ajuste a los

    requerimientos del sistema, formulados en la fase 1.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    13/71

    1.2 ESTRUCTURAS TPICAS DE LOS SISTEMAS DE CONTROL 1-13

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Una vez logrado un resultado satisfactorio, la etapa final consiste en la documentacin del

    proyecto, usando tcnicas y herramientas de la ingeniera de detallepara la descripcin de

    las especificaciones de cada componente del sistema de control, el desarrollo de planos y

    diagramas de control usando simbologa ISA (Instrument Society of America) [ISA92].

    1.3 EJEMPLOS TIPICOS DE SISTEMAS DE CONTROL

    A continuacin se presentar la descripcin de de casos prcticos de sistemas de control,

    algunos de los cuales sern utilizados en los ejemplos de anlisis y diseo en captulos

    posteriores, con el objeto de identificar las 3 acciones operaciones bsicas, la accin de

    control y posible perturbaciones.

    Sistema de control de nivel de un tanque

    La figura 1.13 muestra los componentes de un sistema de control de nivel de lazo cerrado

    para regular el nivel de un tanque que contiene un fluido. Se desea mantener constante el

    nivel (variable controlada), usando como elemento final de control (EFC) la vlvula de

    entrada (VE). La vlvula de salida (VS) se asume que est abierta en una posicin fija.

    Si se presenta un aumento en el caudal de salida ( )sq tel nivel disminuye, generando una

    seal de error ( )e trespecto del valor deseadoo setpoint(SP). La seal de error es utilizada

    por el controlador para generar la seal de control ( )m t que se encarga de abrir la vlvula

    de controlVE, para aumentar el caudal de entrada ( )q ty recuperar as el nivel ( )h t del

    tanque. En este ejemplo la accin de controlconsiste en abrir o cerrar la vlvula VEsi baja

    o sube el nivel del tanque. Una posible perturbacin en este sistema es el cambio en elcaudal del fluido de control, aguas arriba de la vlvula de control: VE.

    Sistema de contro l de temperatura de una cmara de cultivo

    La figura 1.14 muestra el caso tpico de un sistema de control de temperatura de lazo

    cerrado[Phillips00], cuyo propsito es regular la temperatura de una cmara utilizada para

    el cultivo orgnico de plantas. La variable controladaes la temperatura de la cmara y para

    Figura 1.13Sistema de controlde nivel.

    VE

    VS

    SP

    ( )sq t

    ( )q t

    ( )h t

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    14/71

    1-14 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    medirla se usa una resistencia RTD (Resistive Thermal Device), cuyas variaciones ( )ohms

    se convierten en una seal de voltaje ( )mVa travs de un puente de Wheatstone.

    Como la seal de salida del puente suele ser de pocos milivoltios, se usa un amplificador de

    voltaje para acondicionarla a un nivel de 1-5V. La salida del amplificador se compara con

    el valor deseado o setpoint (SP) para generar la seal de error ( )e t, que es utilizada por el

    controlador o compensador para generar la seal de control ( )m t , necesaria para modificar

    la salida de la fuente de potencia que alimenta la resistencia de calefaccin. De este modo,

    la accin de control se traduce en aumentar o disminuir la potencia suministrada a la

    resistencia de calefaccin, cada vez que disminuya o aumente la temperatura interior de la

    cmara de cultivo. Unaposible perturbacinen este sistema de control es la apertura de la

    puerta de entrada, tal como se analizar en el ejemplo 1.4.

    Sistema de control de posicin de una antena

    La figura 1.15 los componentes de este sistema de control, cuyo propsito es regular el

    ngulo de posicin de una antena. El controlador genera una seal de voltaje que es

    utilizada para posicionar el eje de un servomotor, el cual a travs de una caja de engranajes

    regula la posicin de la antena. El sensor reporta una seal de voltaje como una medida del

    ngulo de posicin de la antena, que es utilizada para compararla con el SP para generar la

    seal de error ( )e tentregada a controlador. La accin de controlse reduce a establecer la

    posicin del eje del servomotor y unaposible perturbacines la presin del viento sobre la

    superficie de la antena.

    Figura 1.14Sistema de controlde temperatura.

    Controlador ocompensador

    Servomotor

    Sensor deposicin

    Caja deengranajes

    SP

    +

    voltios grados

    ( )e t ( )m t

    Figura 1.15Sistema de controlde posicin.

    Controlador ocompensador Fuente depotencia

    Cmara decultivo

    Puente deWheatstone

    K

    Resistencia

    RTD

    SP

    Amplificador

    +

    - mV Puerta deacceso

    ( )e t

    ( )m t

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    15/71

    1.3 EJEMPLOS TIPICOS DE SISTEMAS DE CONTROL 1-15

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Sistema de control digital

    Los ejemplos presentados hasta este momento utilizan y procesan seales continuas y en

    este sentido se reconocen como sistemas analgicos de control. Sin embargo, el desarrollo

    de los microprocesadores ha tenido una fuerte incidencia en los sistemas de control,

    permitiendo que la funcin de control pueda efectuarse en forma digital.

    La figura 1.16 muestra los componentes de un sistema de control digital de lazo cerrado

    donde la funcin del controlador analgico se ha sustituido por un microcontrolador ( C).La accin de control [ ]m kse establece a travs de un algoritmo de controlalmacenado en

    el C. En este esquema es necesario incluir un convertidor anlogo digital (A/D) y unconvertidor digital analgico (D/A) como dispositivos de interfase. Detalles relacionados

    con el anlisis y diseo de este modo de control sern tratados en los captulos 5 y 6.

    Sistema de control multivariable

    Los ejemplos anteriores se han referido a sistemas de una entrada y una salida o sistemas

    SISO. Sin embargo, en aplicaciones prcticas, los sistemas de control pueden incluir varios

    lazos de control, que se identifican a travs de la variable controlada de cada uno.

    La figura 1.17 muestra el caso tpico del control de lazo cerrado de un turbogenerador,

    formado por tres componentes: una calderapara la produccin de vapor, una turbinapara

    convertir la energa trmica en energa cintica y un generador de corriente alterna o

    alternador, para transformar la energa cintica en energa elctrica. En este sistema es

    necesario regular simultneamente cuatro variables:

    - composicin de la mezcla en la caldera.

    - temperatura de salida del valor.

    - presin de salida del vapor.

    - frecuencia del voltaje de salida del generador.

    Figura 1.16Sistema de controldigital.

    +( ) [ ]m kC Proceso

    y(t)(t)

    Sensor oTransmisor

    A/D D/A

    Controlador

    ( )m t

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    16/71

    1-16 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    De acuerdo con la estrategia de lazo cerrado, es necesario medir cada una de estas variables

    y enviar su estado actual al controlador, que en este caso en un microprocesador, el cual se

    encarga de generar la respectiva seal de control en funcin de cada seal de error y delmodo de control establecido para regular cada una de las variables controladas.

    Para regular cada una de las variables anteriores es necesario establecer una variable

    manipulada ( )m t para cada lazo de control, que sea capaz mantener su respectiva variable

    controlada cerca del setpoint. Segn la figura 1.17 se identifican cuatro lazos de control:

    - flujo de entrada de aire a la caldera.

    - flujo de entrada de combustible a la caldera.

    - flujo de entrada de agua a la caldera.

    - velocidad nde la turbina en revoluciones por minuto (rpm)

    Para facilitar la identificacin de los lazos de control, en la figura 1.17 las seales

    correspondientes a la accin de control se han dibujado en lneas punteadas. De este modo

    se identifican cuatro lazos de control y se reconoce como un sistema de control

    multivariableo sistema MIMO.

    Figura 1.17Sistema de controlmultivariable para laregulacin de unturbogenerador.

    Combustible

    Aire

    Agua

    Caldera Turbina Generador

    Reguladorde velocidad

    Medidor demezcla

    Medidor detemperatura

    Medidor defrecuencia

    Medidor depresin

    C

    Setpoint de cada lazo de control

    n

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    17/71

    1.3 EJEMPLOS TIPICOS DE SISTEMAS DE CONTROL 1-17

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Simbologa estndar en sistemas de control de procesos

    En los ejemplos anteriores se han utilizando diagramas funcionales, para identificar los

    componentes de cada sistema de control. Sin embargo, en aplicaciones prcticas, en

    particular en los sistemas de control de procesos, se utiliza una simbologa estndar

    [ISA92] desarrollada y aprobada en julio de 1992 entre el American National Standard

    Institute (ANSI) y la Instrument System and Automation Society (ISA), conocida como

    ANSI-ISA S5.1-1984 (R1992).

    Esta simbologa se utiliza en los planos de la ingeniera de detalle del proyecto y la figura

    1.18 muestra el caso de un sistema de control en cascada similar al esquema de la figura

    1.9, para controlar la temperatura de salida2( )T tde un fluido, usando como actuador un

    intercambiador de calor.

    El sistema utiliza un lazo secundario de control para regular el flujo de vapor con el objeto

    de contrarrestar sus posibles variaciones, aguas arriba de la vlvula de control. En el

    apndice F se presenta un resumen de los smbolos y letras utilizadas en este diagrama.

    Intercambiador decalor

    TT25

    TRC

    25

    FY25

    FRC

    25

    FT25

    T

    Vapor

    1( )T t 2( )T t

    SP

    Figura 1.18Sistema de controlen cascada usandosimbologa ISA.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    18/71

    1.4 FUNCION DETRANSFERENCIA Y RESPUESTA IMPULSO 1 - 18

    1 - 18

    1.4 FUNCIN DE TRANSFERENCIA Y RESPUESTA IMPULSO

    En las secciones anteriores se presentaron los elementos necesarios para comprender el

    funcionamiento de un sistema de control. En esta seccin se har referencia a los tres

    modelos clsicosusados en el anlisis y diseo del sistema de control en tiempo continuo:ecuacin diferencial (ED), la respuesta impulso ( )h t y la funcin de transferencia (FT).

    Para el desarrollo del modelo de FT se utilizar como herramienta la transformada de

    Laplace(TL), cuyos fundamentos bsicos se presentan en el apndice B.

    Respuesta dinmica como solucin de una ecuacin d iferencial

    Los 3 modelos clsicos utilizados en el anlisis de un sistema de control se derivan de la

    teora de sistemas lineales y se reconocen como [ReySoto08]:

    - ecuacin diferencial- respuesta impulso

    - funcin de transferencia

    La ecuacin diferencial ordinaria de coeficientes constantes (ED) permite formular la

    relacin entrada-salida representada simblicamente como ( ) ( )x t y t , del sistema linealinvariante en el tiempo (LIT) mostrado en la figura 1.19.

    Asumiendo que el sistema es de orden-2, la ED normalizada se formula como

    2

    1 0 1 02

    d y dy dx a a y b b x

    dt dt dt (1.2)

    En el sistema anterior 2n y 1m . Como m n se puede demostrar [ReySoto08] que es

    causal o realizable. Si el sistema LIT es modelado por una ED, su solucin permite obtener

    la respuesta dinmica del sistema ( )y t a una entrada arbitraria ( )x t. Aunque existen

    mtodos clsicos para esto [ReySoto08], es preferible utilizar mtodos de transformacin

    basados en la transformada de Laplace, tal como se muestra en el siguiente ejemplo.

    EJEMPLO 1.1: Obtener la respuesta dinmica de un sistema LIT modelado por la siguiente

    ED, ante una entrada escaln unitario.

    SISTEMA

    CONTINUO

    ( )x t ( )y tFigura 1.19Relacin entrada

    salida de un sistemaLIT continuo.

    Respuesta

    dinmica como

    solucin de la ED

    de un sistema LIT

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    19/71

    1.4 FUNCIN DE TRANSFERENCIA Y RESPUESTA IMPULSO 1-19

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    2

    2

    (0) 12 6 4 10 ( )

    '(0) 2

    yd y dy y x t

    ydt dt

    Solucin: Comenzamos normalizando la ED

    2

    2 3 2 ( ) 5 ( )d y dy y t x t

    dt dt

    Llevando al dominio-spara ( ) ( )x t u t y aplicando la tabla B.1, obtenemos

    2 ( ) (0) '(0) 3[ ( ) (0)] 2 ( ) ( )s Y s sy y sY s y Y s X s

    Resolviendo para ( )Y s

    2

    2 2

    5 1 5( ) ( )

    ( 3 2) 3 2 ( 1)( 2)

    s s sY s X s

    s s s s s s s

    respuesta naturalrespuesta forzada

    La expresin anterior permite identificar las 2 componentes de la solucin: la

    respuesta natural ( )NY s o respuesta a entrada cero y la respuesta forzada

    ( )FY so respuesta a estado cero.

    Para la respuesta natural ( )N

    y t, usando la tabla B.1, obtenemos

    21 1( ) ( ) , 0

    ( 1)( 2) 2

    t

    N N

    sY s y t e t

    s s s

    Para la respuesta forzada, descomponemos ( )FY sen fracciones parciales

    5 2.5 5 2.5( )

    ( 1)( 2) 1 2FY s

    s s s s s s

    que puede verificarse usando matemtica simblica de MATLAB

    :

    syms s, Ys=5/ s/ ( s 2+3*s+2) ; FPYFs=di f f ( i nt ( Ys) )

    FPYFs = 5/ 2/ s- 5/ ( s+1) +5/ 2/ ( s+2)

    Utilizando la tabla B.1, obtenemos finalmente,

    2( ) 2.5 5 2.5 , 0t tF

    y t e e t

    Luego la respuesta completa es

    2( ) ( ) ( ) 2.5 5 1.5 , 0t t

    N Fy t y t y t e e t

    La solucin completa y las dos componentes pueden obtenerse usando la

    funcin dsolve() del Toolbox de Matemtica Simblica (TBMS) de

    MATLAB:

    yF=dsol ve( ' D2y+3*Dy+2*y=5, y( 0) =0, Dy( 0)=0' )

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    20/71

    1-20 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    yF = 5/ 2+5/ 2*exp( - 2*t ) - 5*exp( - t )

    yN=dsol ve( ' D2y+3*Dy+2*y=0, y( 0)=- 1, Dy( 0)=2' )

    yN = - exp( - 2*t )

    y=dsol ve( ' D2y+3*Dy+2*y=5, y( 0)=- 1, Dy( 0)=2' )

    y = 5/ 2+3/ 2*exp( - 2*t ) - 5*exp( - t )

    Respuesta impulso

    El segundo modelo de un sistema LIT es la respuesta impulso: ( )h t, que se obtiene

    [ReySoto08] asumiendo que la entrada del sistema de la figura 1.19 es ( ) ( )x t t . Comoen este caso la entrada ( ) 0x t para 0t , no es posible incluir condiciones iniciales en lasolucin de la ED y se reconoce como un sistema en reposo. (Respuesta forzada).

    EJEMPLO 1.2: Obtener la respuesta impulso del sistema LIT de ejemplo 1.1, cuya ED es

    2

    2 3 2 5 ( )

    d y dy y x t

    dt dt

    (0) 1

    '(0) 2

    y

    y

    Solucin: Como se asume que el sistema est en reposo, llevamos la ED al dominio-s,

    asumiendo (0) '(0) 0y y , para ( ) ( )x t t y ( ) ( )y t h t :

    2 ( ) 3 ( ) 2 ( ) 5s H s sH s H s

    Resolviendo para ( )H s, obtenemos ( )h tcomo

    2

    2

    5 5 5( ) ( ) 5 5 , 0

    3 2 1 2

    t tH s h t e e t

    s s s s

    Usando convolucin lineales posible [ReySoto08] obtener la respuesta del sistema a una

    entrada arbitraria como

    ( ) ( ) ( ) ( ) ( )y t h x t d h t x t

    (1.3)

    Desplazando la funcin ( )h tse obtiene una expresin alterna para la convolucin

    ( ) ( ) ( ) ( ) ( )y t x h t d x t h t

    (1.4)

    Por lo tanto ( ) ( ) ( ) ( )h t x t x t h t que se reconoce como lapropiedad conmutativa.

    Respuesta

    impulso a partir

    de la ED de un

    sistema LIT

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    21/71

    1.4 FUNCIN DE TRANSFERENCIA Y RESPUESTA IMPULSO 1-21

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Si el sistema y la seal de entrada son causales solo existen para 0t y las expresiones(1.3) y (1.4) se reducen a

    0 0

    ( ) ( ) ( ) ( ) ( )

    t t

    y t h x t d x h t d (1.5)

    Como la respuesta impulso implica que el sistema est en reposo, la convolucin (1.3),

    (1.4) o (1.5) conducen a la respuesta forzada del sistema, suficiente para el anlisis y

    diseo del sistema de control.

    De este modo, un sistema de control LIT puede ser modelado a travs de la respuesta

    impulso ( )h t, tal como se muestra en la figura 1.20, para evaluar su respuesta dinmica

    ante una entrada arbitraria ( )x t.

    Existen tres aspectos importantes en el uso de este modelo:

    1. Como el sistema est en reposo, la solucin ( )y t corresponde a la respuesta forzada.

    2. El uso de la convolucin lineal para obtener la respuesta del sistema puede presentar

    dificultades algebraicas.

    3. Para facilitar la evaluacin de la respuesta se recurre a mtodos de transformacin.

    Funcin de transferencia

    Aunque es posible utilizar la ED o la respuesta impulso para evaluar la respuesta dinmica

    de un sistema de control, es ms prctico recurrir al tercer modelo conocido como la

    funcin de transferencia(FT), representando el sistema de la figura 1.20 en el dominio de

    Laplace, asumiendo condiciones iniciales cero (reposo).

    Si el sistema se modela en funcin de su respuesta impulso ( )h t, aplicando la propiedad de

    convolucin de la transformada de Laplace (tabla B.2), obtenemos:

    ( ) ( ) ( ) ( ) ( ) ( )y t h t x t Y s H s X s (1.6)

    Por lo tanto, el problema se reduce a resolver la ecuacin algebraica (1.6), donde ( )H ses

    la transformada de Laplace de la respuesta impulso ( )h t, es decir

    1( ) ( )} ( ) { ( )}H s h t h t H s L { L (1.7)

    ( )h t( )x t ( )y t

    (Reposo)

    ( )X s ( )Y s

    Figura 1.20Respuesta impulsocomo modelo de unsistema continuo enreposo.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    22/71

    1-22 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Despejando ( )H sen la ecuacin (1.6) es posible definir la FT en los siguientes trminos.

    DEFINICION 1.4 Funcin de transferencia

    La funcin de transferencia (FT) es la funcin caracterstica de un

    sistema LIT en reposo en el dominio-sy se obtiene como la relacin entre

    la TL de la salida ( )y t y TL de la entrada ( )x t:

    ( )( )

    ( )reposo

    Y sH s

    X s (1.8)

    El siguiente ejemplo muestra el clculo de la FT y de la respuesta impulso, a partir del

    modelo bsico de ED de un sistema LIT, y su aplicacin para obtener la respuesta dinmica

    del sistema ante una entrada escaln.

    EJEMPLO 1.3: Obtener la respuesta escaln del sistema LIT del ejemplo 1.2, cuya ED es

    2

    2 3 2 5 ( )

    d y dy y x t

    dt dt

    Solucin: Para obtener la FT llevamos la ED al dominio-s, asumiendo el sistema en reposo

    2 ( ) 3 ( ) 2 ( ) 5 ( )s Y s sY s Y s X s

    Utilizando la definicin (1.8) de FT

    2

    ( ) 5( )

    ( ) 3 2

    Y sH s

    X s s s

    Para obtener la respuesta escaln aplicamos (1.6), asumiendo ( ) 1/X s s

    2

    5 1 5( )

    3 2 ( 1)( 2)Y s

    s s s s s s

    Usando fracciones parciales obtenemos la TIL

    22.5 5 2.5( ) ( ) 2.5 5 2.5 , 01 2

    t tY s y t e e t

    s s s

    FT y respuesta

    escaln de

    sistema LIT.

    Figura 1.21Respuesta escalnde un sistemacontinuo en reposo.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    23/71

    1.4 FUNCIN DE TRANSFERENCIA Y RESPUESTA IMPULSO 1-23

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    La grfica de ( )y t se muestra en la figura 1.21 y corresponde a la respuesta

    forzadade la solucin que se obtuvo en el ejemplo 1.1.

    Estabilidad a partir de la respuesta impulso y la FT

    El resultado del ejemplo 1.3 muestra que la respuesta escaln es suficiente para evaluar la

    estabilidad de un sistema LIT. Lo anterior se basa en el concepto de estabilidad acotada

    que se define a continuacin:

    DEFINICION 1.5 Estabilidad acotada o estabilidad BIBO

    Un sistema LIT se considera que posee estabilidad acotada, si para una

    entrada ( )x tacotada, su salida ( )y t tambin es cotada. Esta condicin se

    define como estabilidad BIBO (Bounded-Input Bounded-Output) la cual

    establece que para todo valor de t, debe existir un valor finito M, tal que

    ( ) ( )x t M y t M (1.9)

    Un sistema LIT que satisface (1.9) se dice que es absolutamente estable. Como

    consecuencia de lo anterior, si el sistema es modelado por su respuesta impulso ( )h t,

    aplicando convolucin (1.3) para ( ) ( )x t t se puede demostrar [ReySoto08] que elsistema tiene estabilidad absoluta, si ( )h tes absolutamente integrable, es decir

    0

    ( )h t dt M

    (1.10)

    De acuerdo con (1.7), la respuesta impulso ( )h t puede evaluarse como la transformada

    inversa de la FT ( )H s. Para esto, asumiendo que ( )H ses causaly tiene npolos simples,

    reales o complejos conjugados, aplicando el mtodo de fracciones parciales obtenemos:

    1 21 21 2

    1 2

    ( ) ( ) np t p t p t nn

    n

    k k kH s h t k e k e k e

    s p s p s p

    (1.11)

    Luego, cada polokp a jb genera en (1.11) un trmino de la forma ( )a jb t at jbt ke ke e .

    Por lo tanto, para lograr la condicin de estabilidad absoluta (1.10) todos los trminos en(1.11) deben tener parte real negativa, es decir

    0kp R e (1.12)

    Lo anterior se consigue si los polos de ( )H sse ubican en el semiplano izquierdo (SPI) del

    plano-s. La expresin algebraica que permite evaluar los polos del sistema se reconoce

    como la ecuacin caractersticay se obtiene a partir del denominador de la FT.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    24/71

    1-24 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    EJEMPLO 1.4: Evaluar la estabilidad del sistema LIT del ejemplo 1.3, donde

    2

    5( )

    3 2H s

    s s

    Solucin: Para evaluar los polos de ( )H sdebemos resolver la ecuacin caracterstica:

    23 2 0s s

    El resultado es 1 21, 2p p , que se ubican en el semiplano izquierdo (SPI)del plano-s y por lo tanto el sistema es absolutamente estable. Tambin es

    posible evaluar la estabilidad a partir de la respuesta impulso ( )h t, que

    corresponde a la TIL de ( )H s. Utilizando el resultado del ejemplo 1.2:

    2( ) 5 5 , 0t th t e e t

    Considerando que ( )h ten (1.10) es real

    2 20 0

    0 0

    ( ) 5 5 5 5 2.5t t t t

    h t dt e e dt e e

    El resultado anterior se puede verificar utilizando MATLAB:

    pk=r oot s( [ 1 3 2] )

    pk = - 2 - 1

    syms t , ht =5*exp( - t ) - 5*exp( - 2*t ) ; Sh=i nt ( abs( ht ) , 0, i nf )

    Sh = 5/ 2

    Propiedades de la func in de transferencia

    El anlisis anterior nos permite enumerar las siguientes propiedades de la funcin de

    transferencia (FT) cuando se utiliza para modelar un sistema LIT:

    P1.La FT es la transformada de Laplace de la respuesta impulso de un sistema LIT.

    P2.Es posible obtener la ED del sistema LIT a partir de su FT, aplicando la propiedad

    de derivacin de la transformada de Laplace, asumiendo condiciones iniciales cero.

    P3.Los componentes de un sistema de control pueden ser modelados por una FT

    obtenida a partir de su ED, asumiendo condiciones iniciales cero.

    Aplicando la propiedad P3 al esquema de la figura 1.10, obtenemos el modelo del sistema

    de control de lazo cerradomostrado en la figura 1.22.

    Estabilidad a

    partir de la FT y

    la respuesta

    impulso.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    25/71

    1.4 FUNCIN DE TRANSFERENCIA Y RESPUESTA IMPULSO 1-25

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Los componentes de este modelo son:

    ( )c

    G s: FT del controlador o compensador, segn el modo de control utilizado.

    ( )pG s: FT del proceso o planta.

    ( )H s: FT del sistema de medicin y transmisin.

    La figura 1.22 constituye la base de la teora de diagramas de bloque, cuyos fundamentos

    sern presentados en la seccin 1.5. Las seales que intervienen en este modelo son:

    ( )Y s: variable controlada. Establece el propsito del sistema de control.

    ( )B s: seal de realimentacin. Es una medida de la variable controlada.

    ( )R s: seal de referencia, setpoint o valor deseado de la variable controlada.

    ( ):E s seal de error. Es una medida de las desviaciones de la variable controlada.

    ( )M s: accin de control. Calculada en funcin de la seal de error.

    En algunas aplicaciones prcticas ( ) 1H s y se refiere como un sistema de control conrealimentacin unitariasimplemente de lazo cerrado unitario.

    EJEMPLO 1.5: Obtener la ED del proceso del sistema de control de lazo cerrado de la

    figura 1.22, asumiendo que su FT es:

    2

    3 2 ( )( )

    2 5 ( )p

    s Y sG s

    s s M s

    Solucin: Utilizando descomposicin directa de la FT, obtenemos

    2( 2 5) ( ) (3 2) ( )s s Y s s M s

    Como el concepto de FT implica que el sistema est en reposo, aplicamos la

    propiedad de la derivada, llevamos la expresin anterior al dominio-t,

    asumiendo condiciones iniciales cero:

    2

    2 2 5 ( ) 3 2 ( )

    d y dy dm y t m t

    dt dt dt

    que es la ED del proceso del sistema de control de lazo cerrado.

    Figura 1.22Componentes de unsistema de control

    de lazo cerrado.

    +( )c

    G s( )R s ( )M s

    ( )pG s( )Y s

    ( )H s

    ( )E s

    ( )B s

    ED asociada

    con una FT.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    26/71

    1.5 FUNCION DE MATLABEN SISTEMAS DE CONTROL 1 - 26

    1 - 26

    1.5 FUNCIONES DE MATLABEN SISTEMAS DE CONTROL

    Utilizando una estrategia similar a la del ejemplo 1.5, podra obtenerse la ED a partir de la

    FT de cada componente del esquema de control de lazo cerrado de la figura 1.22. Sin

    embargo, es conveniente revisar antes los comandos y funciones que ofrece MATLAB

    parala representacin y manipulacin algebraica de una funcin racional de variable compleja

    ( )G s, las cuales se incluyen como parte del Toolbox de Control [MWorks92].

    Creacin de objetos LIT

    La versin 4.2 del Toolbox de Control (TBC) de MATLAB incorpora el uso de objetos

    para la creacin de modelos de sistemas LIT. Se trata de una estructura basada en arreglos

    de celdas, que permite encapsular en una sola variable diferentes tipos de datos de un

    modelo continuo o discreto. Existen tres tipos de objeto LIT [Hanselman97]:

    - tipo TF: modelo de FT representada por la relacin de dos polinomios.

    - tipo ZPK: modelo de FT en forma de factores de polos, ceros y constante de ganancia.

    - tipo SS: modelo de estado (ME).

    Estos tres tipos de objetos LIT pueden crearse usando las siguientes funciones:

    Gs=t f ( nGs, dGs)

    donde nGsy dGsson arreglos del polinomio del numerador y denominador de ( )G s.

    Gs=zpk(z, p, k)

    donde zy pson arreglos de los ceros y polos; el escalar kes la constante de ganancia de ( )G s. meC=ss( A, B, C, D)

    donde A, B, C, Dson arreglos correspondientes a las matrices del modelo de estado.

    Recuperacin de datos en objetos LIT

    Una vez creado el objeto LIT en cualquiera de las tres formas anteriores, es posible obtener

    los datos asociados con cada una, utilizando las siguientes funciones:

    [ num, den] =t f dat a( sys, ' v' )

    Devuelve el polinomio del numerador y denominador de la forma TF del objeto sys.

    [ z, p, k] =zpkdat a( sys, ' v' )

    Devuelve los ceros, polos y constante de ganancia de la forma ZPK del objeto sys.

    [ a, b, c, d] =ssdat a( sys)

    Devuelve las matrices del modelo de estado de la forma SS del objeto sys.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    27/71

    1.5 FUNCIONES DE MATLABEN SISTEMAS DE CONTROL 1-27

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    El siguiente ejemplo muestra el uso de estas funciones para el desarrollo del modelo de un

    sistema continuo LIT tipo SISO. En el captulo 5 se presentarn las variantes de estas

    funciones para la creacin de objetos discretos LIT. El tutorial de MATLABdel apndice D

    incluye el uso de estas y otras funciones de objetos tipo MIMO.

    EJEMPLO 1.6: Obtener los objetos LIT de un sistema SISO cuya FT viene dada por:

    3 2

    2 3( )

    2 5

    sG s

    s s s

    Solucin: Para crear ( )G scomo objeto en la forma TF,

    num=[ 2 3] ; den=[ 1 2 5 0] ; G1s=t f ( num, den)

    Tr ansf er f unct i on:

    2 s + 3

    - - - - - - - - - - - - - - - - -s 3 + 2 s 2 + 5 s

    Para crear ( )G sen la forma ZPK, calculamos los polos y ceros de ( )G s, como

    las races del polinomio del numerador y denominador

    z=r oot s( num) ; p=r oot s( den) ; k=2; G2s=zpk(z, p, k)

    Zer o/ pol e/ gai n:

    2 ( s+1. 5)- - - - - - - - - - - - - - - -s ( s 2 + 2s + 5)

    El resultado anterior muestra que las funciones tf()y zpk()crean modelos de

    funcin de transferencia, que son equivalentes. La diferencia est en que la

    forma TF representa a ( )G scomo la relacin de dos polinomios, mientras que

    la forma ZPK lo hace en forma factorizada asociada con sus ceros, polos y

    constante de ganancia. Cuando los polos o ceros son complejos conjugados la

    forma ZPK incluye un polinomio de orden-2. La creacin de modelos de estado

    como objetos LIT en la forma SS ser tratada en el captulo 7.

    Una vez creado un objeto de FT es posible cambiarlo a cualquiera de las dos

    formas de FT, aplicando la misma funcin utilizada para la creacin. Por

    ejemplo, a partir del modelo TF obtenemos el modelo ZPK como

    G2sm=zpk( G1s)

    Zer o/ pol e/ gai n:2 ( s+1. 5)

    - - - - - - - - - - - - - - - -s ( s 2 + 2s + 5)

    Objetos LIT en

    las forma TF,

    ZPK y SS.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    28/71

    1-28 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    De modo similar, a partir de la forma ZPK, obtenemos

    G1sm=t f ( G2s)

    Tr ansf er f unct i on:

    2 s + 3

    - - - - - - - - - - - - - - - - -

    s 3 + 2 s 2 + 5 s

    Como era de esperar, se obtiene el mismo resultado. Para la recuperacin de

    datos asociados con las formas TF y ZPK

    [ nGs, dGs] =t f dat a(G1s, ' v' )

    nGs = 0 0 2 3dGs = 1 2 5 0

    Para los datos de la forma ZPK

    [ z, p, k]=zpkdat a( G2s, ' v' )

    z = p = k =

    - 1. 5000 0 2- 1. 0000 + 2. 000i- 1. 0000 2. 000i

    Es posible obtener datos de la forma TF a partir de la forma ZPK

    [ num, den] =t f dat a( G2s, ' v' )

    num = 0 0 2 3den = 1 2 5 0

    o datos de la forma ZPK a partir de la forma TF

    [ zm, pm, km] =zpkdat a( G1s, ' v' )

    zm = pm = k =- 1. 5000 0 2

    - 1. 0000 + 2. 0000i- 1. 0000 - 2. 0000i

    Utilizando la funcin pzmap() del TBC, es posible capturar directamente los

    polos y ceros de una FT creada como objeto TF o ZPK:

    [ p1, z1] =pzmap( G1s) %a part i r de f or ma TF

    p1 = z1 =

    0 - 1. 5000- 1. 0000 + 2. 0000i- 1. 0000 - 2. 0000i

    [ p2, z2] =pzmap( G2s) %a part i r de f or ma ZPK

    p1 = z1 =0 - 1. 5000

    - 1. 0000 + 2. 0000i- 1. 0000 - 2. 0000i

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    29/71

    1.5 FUNCIONES DE MATLABEN SISTEMAS DE CONTROL 1-29

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Comentarios:

    1. Las funciones tf() y zpk() devuelven modelos equivalentes de funcin de

    transferencia en las formas TF y ZPK. Estas mismas funciones pueden utilizarse para

    cambiar de una forma a otra.

    2. Para recuperar datos de las formas TF y ZPK es necesario usar la cadena 'v'paraindicar que se requieren los valores guardados como arreglo de celdas.

    3. Aunque los modelos de FT creados en estas formas son equivalentes, el modelo ZPK

    ofrece mejor precisin numrica.

    4. En operaciones con objetos LIT existe un orden de precedencia: SS ZP TF(ver ejemplo 1.7). Lo anterior implica que si se combinan dos modelos en forma TF y

    ZPK, el resultado se dar en la forma ZPK.

    5. Existen otras funciones del TBC, asociadas con objetos LIT que sern incorporadas a

    medida que sean requeridas.

    Ms detalles sobre el uso de las funciones tf(), zpk() y ss()pueden conseguirse en el

    Tutorial de MATLABque se presenta en el apndice D.

    1.6 DIAGRAMAS DE BLOQUE Y GRAFICO DE FLUJO DE SEALES

    En la figura 1.22 se utiliz el concepto de funcin de transferencia (FT) para representar

    cada componente del sistema de control de lazo cerrado, en forma de bloques funcionales

    entrada salida. En esta seccin se utilizarn los conceptos de diagrama de bloques(DB)

    y grfico de flujo de seales (GFS) para desarrollar el modelo grfico del sistema decontrol. Aplicando lafrmula de ganancia de Mason(FGM) al GFS, ser posible evaluar la

    funcin de transferencia equivalentede lazo cerrado de un sistema de control.

    Elementos del diagrama de bloques

    La figura 1.23 muestra los elementos utilizados en la construccin de un diagrama de

    bloques(DB), donde la FT se utiliza para representar la gananciade cada bloque.

    Usando estos 3 elementos bsicos se construy el diagrama del sistema de control de lazo

    cerrado de la figura 1.22.

    ( )Y s( )G s

    ( )X s

    ( )X s

    ( )Y s

    ( ) ( )X s Y s

    +

    ( )X s ( )X s

    ( )X sBloque de ganancia

    Sumador Derivacin o toma

    Figura 1.23Elementos de undiagrama debloques.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    30/71

    1-30 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Aunque existe un conjunto de reglas algebraicas [Dorf05] que pueden aplicarse para

    simplificar un diagrama de bloques (DB), como se demostrar posteriormente es ms

    prctico hacerlo usando la frmula de ganancia de Mason (FGM). Sin embargo, tres de

    estas reglas son de utilidad en la manipulacin de DB aplicados a sistemas de control y se

    describen a continuacin:

    R1. Bloques en serie o en cascada:

    La figura 1.24 muestra dos bloques conectados en serie o en cascada, los cuales

    pueden reducirse a un bloque equivalente, como el producto de sus ganancias.

    Para obtener la FT equivalente, formulamos la relacin entrada salida del segundobloque como: 2( ) ( ) ( )Y s G s M s . De modo similar, en el primer bloque obtenemos:

    1( ) ( ) ( )M s G s X s . Sustituyendo ( )M sen la expresin anterior:

    2 1( ) ( ) ( ) ( ) ( ) ( )eY s G s G s X s G s X s

    donde,

    2 1( ) ( ) ( )eG s G s G s (1.13)

    es la ganancia oFT equivalentemostrada en la figura 1.23. Si se trata de un sistemaSISO, los valores de 1( )G sy 2( )G sson escalares y el orden del producto no afecta el

    resultado de (1.13).

    R2. Bloques en paralelo:

    La figura 1.25 muestra dos bloques conectados en paralelo, que pueden reducirse a un

    equivalente, como la suma de sus ganancias.

    Para desarrollar la expresin equivalente, calculamos la salida de cada bloque como

    1 1( ) ( ) ( )Y s G s X s y 2 2( ) ( ) ( )Y s G s X s . Sumando estas dos seales, la salida ( )Y ses

    ( )Y s

    2 1( ) ( )G s G s

    ( )X sFigura 1.24Reduccin debloques en serie oen cascada.

    1( )G s

    ( )Y s( )X s

    2( )G s( )M s

    Figura 1.25Reduccin debloques enparalelo.

    1( )G s

    ( )Y s( )X s

    2( )G s

    +

    ( )Y s1 2( ) ( )G s G s ( )X s

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    31/71

    1.6 DIAGRAMAS DE BLOQUE Y GRFICO DE FLUJO DE SEALES 1-31

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    1 2 1 2( ) ( ) ( ) ( ) ( ) ( ) ( )Y s Y s Y s G s X s G s X s

    Finalmente, factorizando ( )X s

    1 2( ) ( ) ( )eG s G s G s (1.14)

    que es la ganancia oFT equivalentemostrada en la figura 1.25.

    R3. Bloques en realimentacin:

    La figura 1.26 muestra la conexin de dos bloques en realimentacin, conocida como

    forma cannica de lazo cerrado, por su similitud con el modelo del sistema de control

    en lazo cerrado de la figura 1.22.

    Para conseguir la expresin equivalente, calculamos la salida como

    ( ) ( ) ( )Y s G s E s

    Por otro lado, la seal ( )E spuede expresarse como

    ( ) ( ) ( ) ( ) ( ) ( )E s R s B s R s H s Y s

    Sustituyendo en la expresin anterior

    ( ) ( ) ( ) ( ) ( )Y s G s R s H s Y s

    Reagrupando trminos y simplificando

    ( )( ) ( )

    1 ( ) ( )

    G sY s R s

    G s H s

    (1.15)

    Por lo tanto

    ( )( )

    1 ( ) ( )e

    G sG s

    G s H s

    (1.16)

    que es la gananciao FT equivalentemostrada en la figura 1.26.

    Como el signo del sumador de la figura 1.26 es negativo, se refiere como realimentacin

    negativa. Si el signo de ( )B s en la figura 1.26 es positivo, el denominador de (1.16) se

    convierte en 1 ( ) ( )G s H s y se refiere como realimentacin positiva. Sin embargo, como sedemuestra ms adelante, en aplicaciones prcticas de control la realimentacin positiva

    genera inestabilidad en el sistema y por lo tanto no es de inters prctico.

    Figura 1.26FT equivalente dela forma cannicade lazo cerrado.

    ( )Y s( )

    1 ( ) ( )

    G s

    G s H s

    ( )R s

    ( )G s( )Y s( )R s

    ( )H s

    + ( )E s

    ( )B s

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    32/71

    1-32 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Funcin de transferencia de lazo abierto y de lazo cerrado

    Asociadas con la forma cannica de la figura 1.26, se definen dos expresiones que son

    fundamentales en el anlisis y diseo de un sistema de control realimentado (feedback). La

    primera se reconoce como lafuncin de transferencia de lazo cerrado(FTLC): ( )T s, que

    permite establecer la relacin ( ) ( )R s Y s y puede obtenerse aplicando (1.16). Para elcaso de realimentacin negativa, obtenemos

    ( ) ( )( )

    ( ) 1 ( ) ( )

    Y s G s T s

    R s G s H s

    (1.17)

    La segunda es la funcin de transferencia de lazo abierto (FTLA): ( )F s, utilizada para

    formular la relacin ( ) ( )E s B s , como

    ( )( ) ( ) ( )

    ( )

    B sF s G s H s

    E s

    (1.18)

    Tomando el denominador de (1.17) es posible identificar la ecuacin caracterstica (EC)

    del sistema de control modelado por laforma cannica, como

    1 ( ) ( ) 0G s H s (1.19)

    De acuerdo con (1.12) la estabilidad absoluta del sistema de control se logra si las races de

    esta ecuacin estn en el semi-plano izquierdo (SPI) del plano-s.

    EJEMPLO 1.7: Aplicando las reglas bsicas del lgebra de bloques, reducir el siguiente

    diagrama de bloques a laforma cannicade lazo cerrado y a partir de este

    obtener la FTLA y FTLC. Evaluar adems su estabilidad.

    1

    1( )

    2

    sG s

    s

    2( )

    1

    sG s

    s

    1( ) 2H s

    2

    1( )H s

    s

    Solucin: Para el lazo de realimentacin2 2( ) ( )G s H saplicamos (1.16)

    23

    2 2

    ( )( )

    1 ( ) ( )

    G sG s

    G s H s

    Simplificacin

    del DB para

    calcular la

    FTLA, FTLC y

    EC en la forma

    cannica.

    2( )G s

    2( )H s

    1( )H s

    + +1( )G s

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    33/71

    1.6 DIAGRAMAS DE BLOQUE Y GRFICO DE FLUJO DE SEALES 1-33

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Utilizando los datos del enunciado

    2

    31( )

    1 ( 1) 21

    1

    s

    s ssG ss s s s s

    s s

    Este resultado puede verificarse usando objetos LIT de MATLAB

    G1s=t f ( [ 1 1] , [ 1 2] ) ; G2s=zpk(0, - 1, 1) ; H1s=2; H2s=t f ( 1, [ 1 0] ) ; G3s=mi nreal ( G2s/ ( 1+G2s*H2s) )

    Zer o/ pol e/ gai n:s

    - - - - -( s+2)

    Para los bloques en cascada 1( )G sy 3( )G saplicamos (1.13)

    4 1 3 2

    1 ( 1)( ) ( ) ( )

    2 2 ( 2)

    s s s s G s G s G s

    s s s

    Utilizando objetos LTI de MATLAB

    G4s=mi nreal ( G1s*G3s)

    Zer o/ pol e/ gai n:s ( s+1)- - - - - - -( s+2) 2

    Luego, el diagrama anterior se reduce a la siguiente forma estndar

    2

    ( 1)( )

    ( 2)

    s sG s

    s

    ( ) 2H s

    Aplicando (1.17) obtenemos la FTLC como

    forma TF forma ZPK

    22

    2 2

    2

    1 1 13 3 3

    4 43 3

    ( 1)

    ( 1)( ) ( 2)( )( 1)1 ( ) ( ) 2 2

    1 2( 2)

    s s

    s s s s G s sT ss sG s H s s s s s

    s

    La forma ZPK muestra el denominador como un polinomio de orden-2, lo cual

    permite identificar que los polos de ( )T s son complejos conjugados. Usando

    objetos LTI de MATLAB:

    2

    +2

    ( 1)

    ( 2)

    s s

    s

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    34/71

    1-34 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Gs=G4s, Hs=H1s, Ts=f eedback( Gs, Hs)

    Zer o/ pol e/ gai n:0. 33333 s ( s+1)

    - - - - - - - - - - - - - - - - - -( s 2 + 2s + 1. 333)

    Aunque la expresin de ( )T spuede calcularse usando operaciones con objetos

    LIT, tal como se hizo al evaluar a3( )G s, el TBC incluye la funcin

    feedback()que facilita el clculo directo de la FTLC, a partir de ( )G sy ( )H s

    correspondientes a laforma cannicade lazo cerrado.

    Aplicando (1.18), la FTLA es

    2 2

    ( 1) 2 ( 1)( ) ( ) ( ) 2

    ( 2) ( 2)

    s s s s F s G s H s

    s s

    La ecuacin caracterstica(EC) necesaria para evaluar la estabilidad es

    22

    2 2

    2 ( 1) ( 2) 2 ( 1)1 0 3 6 4 0

    ( 2) ( 2)

    s s s s s s s

    s s

    que corresponde al denominador de ( )T s. Calculando sus races, los polos del

    sistema en lazo cerrado son: 1,2 1 0.5774p j . Como se encuentran en el SPI,el sistema es absolutamente estable.

    Utilizando MATLAB:

    pk=pol e(Ts)pk = - 1. 0000 + 0. 5774i

    - 1. 0000 0. 5774i

    Comentarios:

    1. En el ejemplo anterior las funciones1( )G sy

    2( )H sse crearon como objetos LIT en

    la forma TF, mientras que2( )G sse cre en la forma ZPK.

    2. Sin embargo, al combinar estos objetos para obtener 3( )G saparece en la forma ZPK,

    como consecuencia del orden de precedenciaen las operaciones con objetos.3. La funcin minreal() del TBC se utiliza para simplificar en un resultado, los

    factores comunes del numerador y denominador.

    4. En aplicaciones prcticas es conveniente usar la funcin feedback()para evaluar la

    funcin de transferencia de lazo cerrado correspondiente a laforma cannica.

    5. La funcin pole()del TBC permite calcular los polos de una FT. En el ejemplo 1.7

    se utiliz para calcular las races de la ecuacin caracterstica, a partir de la FTLC.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    35/71

    1.6 DIAGRAMAS DE BLOQUE Y GRFICO DE FLUJO DE SEALES 1-35

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Fundamentos del grfico de flujo de seales

    Una forma prctica para evaluar la FTLC de un sistema que incluye varios lazos de control,

    sin necesidad de reducirlo a la forma cannica, consiste en utilizar el grfico de flujo de

    seales (GFS), introducido por S.J. Mason en 1953 para la representacin causaefectode un sistema LIT modelado a travs de ecuaciones algebraicas. El GFS es un modelo

    grficode la funcin de transferencia de un sistema LIT y puede considerarse como una

    versin simplificadadel diagrama de bloques (DB) presentado en el prrafo anterior.

    El GFS est asociado a un DB y utiliza los 3 elementos mostrados en la figura 1.27, donde

    las seales se representan mediante nodos interconectados por ramas orientadas. A cada

    rama se le asigna la ganancia G correspondiente a la FT del bloque. Esta ganancia se

    reconoce como la transmitanciade la rama, para indicar que la seal ( )X ses transmitida

    desde un nodo hacia otro nodo, para crear la seal ( )Y s.

    Un nodo se considera como un sumideroo depsito en el cual se guarda el valor de una

    seal. El propsito del GFS es formular una ecuacin en el dominio-spara expresar la seal

    en un nodo cualquiera, como relacin causaefectoo entradasalida: ( ) ( )X s Y s De este modo, se logra un conjunto de ecuaciones estndarde la forma

    ( ) ( ) ( )j jk kY s G s X s (1.20)

    donde ( )jkG ses la gananciao transmitanciaentre un nodo de entrada ( )kX sy un nodo de

    salida ( )jY s. En este sentido, en un GFS se identifican 3 tipos de nodos:

    Nodo de entrada: solo transmite seales; solo pueden existir ramas que salen.

    Nodo de salida: solo recibe seales; a este nodo solo pueden llegar ramas.

    Nodo mixto: recibe y transmite seales; pueden existir ramas entrando y saliendo.

    Al recorrer un GFS entre un par de nodos en el sentido de sus ramas, es posible identificar

    los siguientes elementos:

    Trayectoria directa: trayectoria abierta desde un nodo de entrada hasta un nodo de

    salida, sin pasar por un nodo ms de una vez. El producto de las

    ( )Y s( )X s G( ) ( )X s Y s

    1

    ( )X s

    ( )Y s

    ( )X s

    ( )X s

    ( )X s

    Figura 1.27Elementos de ungrfico de flujo deseales.

    Ganancia

    Sumador Derivacin o toma

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    36/71

    1-36 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    ganancias de las ramas recorridas establece la ganancia de la

    trayectoria directa.

    Lazo: trayectoria cerrada que inicia y termina en el mismo nodo, sin pasar por un nodo

    ms de una vez. El producto de las ganancias de las ramas recorridas establece la

    ganancia del lazo.

    Trayectorias conexas: trayectorias directas o lazos que se tocan o que tienen un

    elemento en comn (nodo o rama).

    La figura 1.28 muestra un ejemplo tpico de un GFS utilizado para formular la relacin

    ( ) ( )R s Y s , donde existe un nodo de entrada: ( )R s, un nodo de salida: ( )Y sy 4 nodosmixtos:

    1( )E s,

    2( )E s,

    3( )E s,

    4( )E s.

    As mismo se identifican los siguientes elementos:

    2 trayectorias directas para la relacin ( ) ( )R s Y s : 1 1 2 3 4T G G G G , 2 4 5T G G 2 lazos:

    1 1 1L G H , 2 2L H

    En cuanto a la relacin entre lazos y trayectorias, se observa que:

    Los lazos 1L y 2L no son conexos.La trayectoria 1Tes conexacon los 2 lazos.

    La trayectoria 2Tes conexasolo con 1L .

    Estas caractersticas son indispensables al aplicar la frmula de ganancia de Mason.

    Frmula de ganancia de Mason

    Esta frmula o algoritmo fue desarrollada por S. J. Mason en el ao 1956 como parte de un

    procedimiento grfico para determinar la funcin de transferencia (FT) entre un nodo de

    entrada y un nodo de salida, de un grfico de flujo de seales (GFS) o de su diagrama debloques (DS) equivalente. La base de este algoritmo es la regla de Crammer, utilizada para

    resolver un sistema simultneo de ecuaciones algebraicas.

    Asumiendo un sistema continuo LIT, la frmula de ganancia de Mason (FGM) permite

    calcular la ganancia entre un nodo de entraday un nodo de salidacomo:

    1G( )R s

    1H

    2G1E 2E 3G

    2H

    3E 4E ( )Y s

    5G

    4G

    Figura 1.28Ejemplo tpico deun grfico de flujode seales.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    37/71

    1.6 DIAGRAMAS DE BLOQUE Y GRFICO DE FLUJO DE SEALES 1-37

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    1

    ( ) ( )( )

    ( )( ) ( )

    NTD

    k k

    k

    T s sY s

    T sX s s

    (1.21)

    donde:

    ( )T s ganancia o FT equivalente para la relacin ( ) ( )X s Y s .( )s determinante del sistema, asociado con los lazos del GFS.( )

    kT s ganancia de trayectorias directaspara la relacin ( ) ( )X s Y s .

    ( )ks cofactorde cada trayectoria directa, evaluado a partir de ( )s .

    Para facilitar la aplicacin de la FGM, se recomienda seguir un orden especfico al evaluar

    los trminos que conforman la ecuacin (1.21), usando el siguiente procedimiento:

    1. Obtener el determinante del sistema: ( )s , como:

    2 3( ) 1 i i is L L L (1.22)donde

    iL representa la ganancia de los lazos simples(cada lazo); 2iL la ganancia de los

    lazos no conexosen grupos de2;3iL la ganancia de los lazos no conexos en grupos de

    3, y as sucesivamente. Como el determinante ( )s solo depende de los lazos del GFS,se convierte en un valornico del grafo, independiente del par de nodos entrada-salida

    considerados al aplicar (1.21).

    2. Calcular la ganancia de las trayectorias directas: ( )k

    T s.

    Las trayectorias directas dependen nicamente de la relacin ( ) ( )X s Y s para la cual

    se quiere calcular la FT, aplicando (1.21). Se obtienen recorriendo el grafo desde el nodo

    de entrada ( )X shasta el nodo de salida ( )Y s.

    3.Evaluar los cofactores asociados con cada trayectoria directa: ( )ks .

    El cofactor de una trayectoria directa se obtiene anulando en el determinante del sistema

    ( )s , la ganancia de los lazos conexos con dicha trayectoria.

    El siguiente ejemplo muestra la aplicacin de la FGM para determinar la ganancia o FT

    equivalente del modelo grfico de un sistema continuo LIT, utilizando el procedimiento

    anterior.

    EJEMPLO 1.8: Aplicar la FGM para obtener la ganancia o FT equivalente del sistema

    continuo LIT cuyo modelo grfico se muestra en la figura 1.28.

    Solucin: Aplicando el procedimiento anterior, obtenemos

    1.Determinante del sistema:

    FGM para

    calcular

    ganancia

    equivalente.

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    38/71

    1-38 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    Existen 2 lazos cuyas ganancias son 1 1 1L G H y 2 2L H . Como no sonconexos, conforman un grupo de 2 cuya ganancia es

    1 2 1 1 2L L G H H . Por lo

    tanto, aplicando (1.22), obtenemos

    1 2 1 2 1 1 2 1 1 2( ) 1 ( ) 1s L L L L G H H G H H 2. Trayectorias directas:

    Para la relacin ( ) ( )R s Y s existen 2 trayectorias directas con ganancias:

    1 1 2 3 4T G G G G

    2 4 5T G G

    3. Cofactores de trayectorias directas:

    La trayectoria1

    T es conexa con los 2 lazos; para obtener su cofactor

    anulamos 1L y 2L en ( )s , obteniendo 1( ) 1s . Por otro lado, como latrayectoria

    2T es conexa solo con

    1L debemos anular este lazo en ( )s ,

    obteniendo 2 2 2( ) 1 1s L H .

    4. Ganancia equivalente:aplicando la FGM (1.21)

    1 1 2 2 1 2 3 4 4 5 2

    1 1 2 1 1 2

    (1 )( )

    1

    T T G G G G G G H T s

    G H H G H H

    Estabilidad a partir de la FGM

    Adems de facilitar el clculo de la FT de un sistema de control con mltiples lazos decontrol, la FGM permite evaluar su estabilidad, si se considera que los polos de ( )T sen la

    expresin (1.21), se pueden obtener a partir de las races del denominador, como

    ( ) 0s (1.23)

    De este modo la expresin (1.23) es la ecuacin caractersticade un sistema de control que

    utiliza un esquema arbitrario. La ecuacin (1.19) es un caso particular que solo puede ser

    aplicado a un sistema cuyo esquema corresponde a laforma cannicade la figura 1.25.

    Del diagrama de bloques al grfico de flujo de seales

    Como se mencion anteriormente la FGM puede aplicarse a un diagrama de bloques (DB).

    Sin embargo, es ms fcil identificar lazos y trayectorias en el GFS. Para lograr un GFS

    con el nmero mnimo de nodos, es conveniente tomar en cuenta las siguientes sugerencias:

    - identificar y dar nombre en el DB, a todas las seales de entrada a cada bloque

    (podra trabajarse con las seales de salida. Ver ejemplo 1.16.).

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    39/71

    1.6 DIAGRAMAS DE BLOQUE Y GRFICO DE FLUJO DE SEALES 1-39

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    - seleccionar el nmero mnimo de nodos igual al nmero de seales anteriores, ms

    las seales de entrada y salida del sistema.

    - desarrollar el GFS a partir de las relaciones algebraicas causa efecto del DB.

    EJEMPLO 1.9: Aplicar la FGM para obtener FTLC del sistema del ejemplo 1.7.

    Solucin: Antes de utilizar la FGM desarrollamos el GFS, a partir del DB. Utilizando la

    sugerencia anterior, el nmero mnimo de nodos:

    - seales de entrada y salida del sistema: ( )R s, ( )Y s

    - seales de entrada a cada bloque: , ,1 2( ) ( ) ( )E s E s Y s

    No. mnimo de nodos 2 2 4

    Usando estos 4nodos desarrollamos el GFS considerando las relaciones causa

    efecto del DB.

    Aplicando el procedimiento propuesto para la relacin ( ) ( )R s Y s :

    1.Determinante del sistema.

    Dos lazos simples:1 2 2L G H y

    2 1 2 1L G G H

    No existen grupos no conexos.

    Aplicando (1.22), obtenemos

    1 2 2 2 1 2 1( ) 1 ( ) 1s L L G H G G H

    2. Trayectorias directas.Una trayectoria directa para la relacin ( ) ( )R s Y s :

    1 1 2T G G

    GFS a partir

    del DB y FGM

    para calcular

    FTLC. 2( )G s

    2( )H s

    1( )H s

    + +1( )G s

    ( )R s 1( )E s 2( )E s ( )Y s

    2H

    1G( )R s 2G1( )E s 2( )E s ( )Y s

    1H

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    40/71

    1-40 Captulo 1 FUNDAMENTOS Y MODELOS MATEMTICOS DE SISTEMAS DE CONTROL

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    3. Cofactores de trayectorias directas.

    1Tes conexa con

    1L y

    2L , luego

    1 1

    4. Ganancia o FT equivalente.

    1 1 1 2

    2 2 1 2 1

    ( )1

    T G GT s

    G H G G H

    Sustituyendo los datos del ejemplo 1.7

    2

    2

    1

    2 1 2( )1 1 1 2 3 6 4

    1 2 11 2 1 1 2

    s s s

    s ss s sT ss s s s s s

    s s s s s s

    Utilizando MATLAB

    G1s=t f ( [ 1 1] , [ 1 2] ) ; G2s=zpk(0, - 1, 1) ; H1s=2; H2s=t f ( 1, [ 1 0] ) ; Ts=mi nreal ( G1s*G2s/ ( 1+G2s*H2s+G1s*G2s*H1s) )

    Zer o/ pol e/ gai n:0. 33333 s ( s+1)

    - - - - - - - - - - - - - - - - - -( s 2 + 2s + 1. 333)

    que corresponde a la forma ZPK del resultado que se obtuvo usando la FGM.

    EJEMPLO 1.10: Obtener el GFS del sistema continuo LIT cuyo DB se muestra a

    continuacin y aplicar la FGM para obtener la ganancia o FT equivalente

    de lazo cerrado.

    Solucin: Para determinar el nmero mnimo de nodos:

    - seales de entrada y salida del sistema: ( )R s, ( )Y s

    - seales de entrada a cada bloque:1( )E s,

    2( )E s,

    3( )E s,

    4( )E s

    GFS a partir de

    un DB y FGM

    para calcular

    FT equivalente

    de lazo cerrado.

    2E

    2H

    +

    1G 2G 3G 4G( )Y s

    4E3E1E( )R s +

    1H

    3H

    ++

    5G

    1Y

  • 7/26/2019 Modelos Matemticos de Sistemas de Control

    41/71

    1.6 DIAGRAMAS DE BLOQUE Y GRFICO DE FLUJO DE SEALES 1-41

    ANALISIS Y DISEO DE SISTEMAS DE CONTROL CON APLICACIONES EN MATLAB Carlos Alberto Rey Soto - 2009

    No. mnimo de nodos 2 4 6

    Utilizando estos 6 nodos desarrollamos el GFS considerando las relaciones

    causa efecto del DB.

    Se observa que no es necesario incluir el nodo1( )Y sdel DB en el GFS, dado

    que es la seal de salida del bloque 3G, se evalu como 1 3 3Y G E .

    Aplicando el procedimiento propuesto para la relacin ( ) ( )R s Y s :

    1.Determinante del sistema.

    3 lazos simples:1 1 1

    L G H ,2 1 2 3 2L G G G H y

    3 3 3L G H

    1 grupo de 2 no conexos: 1 3 1 3 1 3L L G G H H

    Aplicando (1.22), obtenemos

    1 2 3 1 3 1 1 1 2 3 2 3 3 1 3 1 3( ) 1 ( ) 1s L L L L L G H G G G H G H G G H H

    2. Trayectorias directas.

    2 trayectorias directas para la relacin ( ) ( )R s Y s :

    y1 1 2 3 4 2 1 4 5

    T G G G G T G G G

    3. Cofactores de trayectorias directas.

    1Tes conexa con

    1L ,

    2L y

    3L , luego

    1 1

    2Tes conexa con 1L y 2L , luego 2 3 31G H

    4. Ganancia o FT equivalente.

    1 2 3 4 1 4 5 3 3

    1 1 1 2 3 2 3 3 1 3 1 3

    (1 )( )

    1

    G G G G G G G G HT s

    G H