<?xml version="1.0" encoding="ISO-8859-1"?><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
<front>
<journal-meta>
<journal-id>1815-5901</journal-id>
<journal-title><![CDATA[Ingeniería Energética]]></journal-title>
<abbrev-journal-title><![CDATA[Energética]]></abbrev-journal-title>
<issn>1815-5901</issn>
<publisher>
<publisher-name><![CDATA[Universidad Tecnológica de La Habana José Antonio Echeverría, Cujae]]></publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id>S1815-59012012000300002</article-id>
<title-group>
<article-title xml:lang="es"><![CDATA[Aplicación de algoritmos de control clásico, adaptable y robusto a sistemas dinámicos de parámetros variables]]></article-title>
<article-title xml:lang="en"><![CDATA[Classic, adaptable and robust control algorithm application, to variant parameter dynamic system]]></article-title>
</title-group>
<aff id="A01">
<institution><![CDATA[,Universidad de Camagüey  ]]></institution>
<addr-line><![CDATA[ ]]></addr-line>
<country>Cuba</country>
</aff>
<pub-date pub-type="pub">
<day>00</day>
<month>12</month>
<year>2012</year>
</pub-date>
<pub-date pub-type="epub">
<day>00</day>
<month>12</month>
<year>2012</year>
</pub-date>
<volume>33</volume>
<numero>3</numero>
<fpage>184</fpage>
<lpage>195</lpage>
<copyright-statement/>
<copyright-year/>
<self-uri xlink:href="http://scielo.sld.cu/scielo.php?script=sci_arttext&amp;pid=S1815-59012012000300002&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://scielo.sld.cu/scielo.php?script=sci_abstract&amp;pid=S1815-59012012000300002&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://scielo.sld.cu/scielo.php?script=sci_pdf&amp;pid=S1815-59012012000300002&amp;lng=en&amp;nrm=iso"></self-uri><abstract abstract-type="short" xml:lang="es"><p><![CDATA[Existen múltiples sistemas dinámicos cuyos modelos matemáticos se caracterizan por ser de primer orden y parámetros variables con el tiempo. En estos casos las herramientas clásicas no siempre logran un sistema de control que sea estable, posea un buen desempeño dinámico y rechace adecuadamente las perturbaciones, cuando el modelo de la planta se desvía del nominal, para el cual se realizó el diseño. En este trabajo se evalúa el comportamiento de tres estrategias de control en presencia de variación de parámetros. Estas son: control clásico, control adaptable y control robusto. Se realiza un estudio comparativo de las mismas en cuanto a complejidad del diseño, costo computacional de la implementación y sensibilidad ante variaciones en los parámetros y/o presencia de disturbios. Se llega a conclusiones que permiten disponer de criterios para la elección más adecuada, en dependencia de los requerimientos dinámicos que la aplicación demande, así como de los medios técnicos de que se disponga.]]></p></abstract>
<abstract abstract-type="short" xml:lang="en"><p><![CDATA[Many dynamic systems have first order mathematic models, with time variable parameters. In these cases, the classical tools do not satisfy at all control system stability, good performance and perturbation rejection, when the plant model differs from the nominal one, for which the controller was designed. In this article, three control strategies are evaluated in parameter variations and disturbance presence. The strategies are the followings: classical control, adaptive control and robust control. A comparative study is carried out, taking into account the design complexity, the computational cost and the sensitivity. The obtained conclusions helps to provide the criterion to choose the most adequate control strategy, according to the necessary dynamic, as well as the available technical means.]]></p></abstract>
<kwd-group>
<kwd lng="es"><![CDATA[control adaptable]]></kwd>
<kwd lng="es"><![CDATA[control robusto]]></kwd>
<kwd lng="es"><![CDATA[identificación de sistemas]]></kwd>
<kwd lng="en"><![CDATA[adaptable control]]></kwd>
<kwd lng="en"><![CDATA[robust control]]></kwd>
<kwd lng="en"><![CDATA[system identification]]></kwd>
</kwd-group>
</article-meta>
</front><body><![CDATA[ <p align="right"><font size="2" face="Verdana"><b>TRABAJO TEORICO EXPERIMENTAL</b>    </font></p>     <p>&nbsp; </p>     <P><font size="4" face="Verdana"><b>Aplicaci&oacute;n de algoritmos de control    cl&aacute;sico, adaptable y robusto a sistemas din&aacute;micos de par&aacute;metros    variables</b></font><font size="2" face="Verdana"> </font>     <P>&nbsp;     <P><font size="3" face="Verdana"><b>Classic, adaptable and robust control algorithm    application, to variant parameter dynamic system</b></font><font size="2" face="Verdana">    </font>     <P>&nbsp;     <P>&nbsp;     <P>      <P><font size="2" face="Verdana"><b>MSc. Orlando Regal&oacute;n Anias, Dra. </b></font><b><font size="2" face="Verdana">Milagros    Diez Rodr&iacute;guez, MSc. Vlad&iacute;mir Rodr&iacute;guez Diez, Ing. Ra&uacute;l    B&aacute;ez Prieto</font> </b>      <p><font size="2" face="Verdana"> Departamento Ing. El&eacute;ctrica de la Facultad    Electromec&aacute;nica. Universidad de Camag&uuml;ey, Cuba. </font>      ]]></body>
<body><![CDATA[<p>&nbsp;     <P>      <P>      <P>&nbsp; <hr>     <P>      <P>      <P>      <P><font size="2" face="Verdana"><b>RESUMEN</b></font>     <P>      <P><font size="2" face="Verdana">Existen m&uacute;ltiples sistemas din&aacute;micos    cuyos modelos matem&aacute;ticos se caracterizan por ser de primer orden y par&aacute;metros    variables con el tiempo. En estos casos las herramientas cl&aacute;sicas no    siempre logran un sistema de control que sea estable, posea un buen desempe&ntilde;o    din&aacute;mico y rechace adecuadamente las perturbaciones, cuando el modelo    de la planta se desv&iacute;a del nominal, para el cual se realiz&oacute; el    dise&ntilde;o. En este trabajo se eval&uacute;a    el comportamiento de tres estrategias de control en presencia de variaci&oacute;n    de par&aacute;metros. Estas son: control cl&aacute;sico, control adaptable y    control robusto. Se realiza un estudio comparativo de las mismas en cuanto a    complejidad del dise&ntilde;o, costo computacional de la implementaci&oacute;n    y sensibilidad ante variaciones en los par&aacute;metros y/o presencia de disturbios.    Se llega a conclusiones que permiten disponer de criterios para la elecci&oacute;n    m&aacute;s adecuada, en dependencia de los requerimientos din&aacute;micos que    la aplicaci&oacute;n demande, as&iacute; como de los medios t&eacute;cnicos    de que se disponga. </font>      ]]></body>
<body><![CDATA[<P><font size="2" face="Verdana"><b>Palabras clave:</b> control adaptable, control    robusto, identificaci&oacute;n de sistemas.</font>  <hr>     <P><font size="2" face="Verdana"><b>ABSTRACT</b></font>     <P>      <P><font size="2" face="Verdana">Many dynamic systems have first order mathematic    models, with time variable parameters. In these cases, the classical tools do    not satisfy at all control system stability, good performance and perturbation    rejection, when the plant model differs from the nominal one, for which the    controller was designed. In this article, three control strategies are evaluated    in parameter variations and disturbance presence. The strategies are the followings:    classical control, adaptive control and robust control. A comparative study    is carried out, taking into account the design complexity, the computational    cost and the sensitivity. The obtained conclusions helps to provide the criterion    to choose the most adequate control strategy, according to the necessary dynamic,    as well as the available technical means. </font>      <P><font size="2" face="Verdana"><b>Keywords:</b> adaptable control, robust control,    system identification.</font> <hr>     <P>&nbsp;     <P><font size="2" face="Verdana"> </font>      <P>      <P>      <P><font size="3" face="Verdana"><b>INTRODUCCI&Oacute;N</b></font>     ]]></body>
<body><![CDATA[<P>&nbsp;      <P>      <P><font size="2" face="Verdana">Una gran parte de los sistemas de control son    dise&ntilde;ados a partir de un modelo din&aacute;mico de la planta a controlar,    teniendo en cuenta las caracter&iacute;sticas de esta en condiciones nominales;    sin embargo en la pr&aacute;ctica los sistemas tienen un comportamiento no lineal    lo que hace que dicho modelo sea solo v&aacute;lido alrededor de un punto de    operaci&oacute;n. Adem&aacute;s, al cambiar las condiciones bajo las cuales    se obtuvo el modelo, los par&aacute;metros de las ecuaciones diferenciales que    describen el mismo cambian, lo que hace que el proceso se comporte como un sistema    de par&aacute;metros variables con el tiempo. Tal es el caso de las articulaciones    de robots industriales, los accionamientos asincr&oacute;nicos que emplean la    t&eacute;cnica de control vectorial, as&iacute; como los procesos de fermentaci&oacute;n    y secado que ocurren en la industria biotecnol&oacute;gica. Cuando    el objeto de control no presenta variaci&oacute;n con el tiempo o alinealidades,    las herramientas matem&aacute;ticas que brinda la teor&iacute;a de control cl&aacute;sico    se bastan para lograr un dise&ntilde;o del sistema de control, que garantice    estabilidad en lazo cerrado y buenos &iacute;ndices de desempe&ntilde;o en estado    transitorio [1]. Si por el contrario, la planta a controlar fuese no lineal    o de din&aacute;mica variable, como es el caso de un fermentador discontinuo,    entonces con un controlador cl&aacute;sico (P, PI, PID) se logra un desempe&ntilde;o    din&aacute;mico muy malo, entendiendo por mal desempe&ntilde;o cuando los indicadores    de su comportamiento temporal se desv&iacute;an de las especificaciones de dise&ntilde;o.    A&uacute;n cuando la planta es lineal e invariante con el tiempo, caso que logra    ser gobernada con &eacute;xito por controladores cl&aacute;sicos, se pueden    evaluar otras estrategias de control si se desea mejorar la raz&oacute;n de    rechazo a perturbaciones externas. Un ejemplo de ello es el control de posici&oacute;n    de la antena de un radar que est&eacute; expuesta a frecuentes perturbaciones    provocadas por el viento. En casos en que    la teor&iacute;a de control cl&aacute;sico no resuelve totalmente los problemas    de desempe&ntilde;o din&aacute;mico se han aplicado, entre otras t&eacute;cnicas,    las de control adaptable y control robusto. </font>      <P>&nbsp;      <P>      <P><font size="2" face="Verdana"><b>DESARROLLO</b> </font>      <P><font size="2" face="Verdana">Inicialmente se analizar&aacute; el objeto de    control para seguidamente evaluar la aplicaci&oacute;n de cada una de las estrategias    de control propuestas. </font>     <P>      <P><font size="2" face="Verdana"><b>El objeto de control</b> </font>      <P><font size="2" face="Verdana">Para expresar la idea de aproximaci&oacute;n    de un modelo de orden superior a un modelo m&aacute;s sencillo, se ejemplificar&aacute;    a trav&eacute;s del control de velocidad de un motor de corriente directa con    excitaci&oacute;n independiente, el cual cuenta con un modelo muy sencillo y    bien conocido [2]. La velocidad de un motor de corriente directa con excitaci&oacute;n    independiente puede ser controlada manteniendo constante el flujo y variando    el voltaje aplicado a la armadura. Para variar este voltaje es com&uacute;n    utilizar un convertidor CA-CD controlado o semicontrolado. El modelo matem&aacute;tico    del sistema convertidor-motor-carga est&aacute; determinado por la siguiente    <a href="#e1">ecuaci&oacute;n (1)</a>, [3]: </font>      ]]></body>
<body><![CDATA[<P>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0102312.gif" width="522" height="64"></font>    <a name="e1"></a>      
<P><font size="2" face="Verdana">donde: </font>     <br>   &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <font size="2" face="Verdana">W    : velocidad del motor </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;E<sub>m</sub>    : voltaje de mando del convertidor AC-DC </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;K<sub>c</sub>    : ganancia del convertidor AC-DC </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;T<sub>c</sub>    : constante de tiempo del convertidor AC-DC </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;K<sub>&#216;</sub>    : constante de par del motor </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;K<sub>b</sub>    : constante de fuerza contraelectromotriz del motor </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;L<sub>a</sub>    : inductancia de armadura </font>     ]]></body>
<body><![CDATA[<br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;R<sub>a</sub>    : resistencia de armadura del motor </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;J    : inercia del conjunto motor-carga </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;f    : fricci&oacute;n del conjunto motor-carga </font>      <P><font size="2" face="Verdana">Se puede apreciar que el modelo din&aacute;mico    que describe a esta planta es de tercer orden. En la pr&aacute;ctica sucede    que T<sub>c</sub> y L<sub>a</sub> son muy peque&ntilde;as, lo que permite aproximar    este modelo a uno de primer orden, el que de forma general puede ser expresado    en la <a href="#e2">ecuaci&oacute;n (2)</a>: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0202312.gif" width="328" height="64"></font>    <a name="e2"></a>      
<P><font size="2" face="Verdana">donde: </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;C    : se&ntilde;al de salida (en este caso velocidad del motor) </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;R    : se&ntilde;al de entrada (en este caso voltaje de mando del convertidor AC-DC)    </font>     <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;K    : ganancia del sistema (en este caso del conjunto convertidor-motor-carga) </font>        <br>   <font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;T    : constante de tiempo del sistema (en este caso tambi&eacute;n se refiere a    dicho conjunto) </font>      ]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">De forma similar puede ser tratado el modelo    matem&aacute;tico de la articulaci&oacute;n de un robot. Tambi&eacute;n se conoce    que el modelo matem&aacute;tico de un motor asincr&oacute;nico de corriente    alterna, controlado por campo orientado, se puede considerar seg&uacute;n el    modelo descrito por la <a href="#e2">ecuaci&oacute;n (2)</a>, [4]. Por otro    lado, los procesos de fermentaci&oacute;n y secado para la fabricaci&oacute;n    de productos biotecnol&oacute;gicos, al ser linealizados se comportan como modelos    PORT (primer orden con retardo de transporte) [5]. En muchas ocasiones se tiene    que el retardo de transporte es muy peque&ntilde;o en comparaci&oacute;n con    la constante de tiempo y se puede despreciar, obteniendo un modelo resultante    aproximado de primer orden como el descrito en <a href="#e2">ecuaci&oacute;n    (2)</a>. En los ejemplos mencionados anteriormente    la aproximaci&oacute;n es v&aacute;lida si se considera que los par&aacute;metros    del modelo permanecen constantes y no var&iacute;an con el tiempo; pero esto    sucede ya sea por el cambio en la carga que manipula el robot, por las caracter&iacute;sticas    intr&iacute;nsecas del control por campo orientado o por la naturaleza no lineal    de los procesos biotecnol&oacute;gicos. </font><font size="2" face="Verdana">Por    lo antes expuesto se investigar&aacute; c&oacute;mo se comportan las estrategias    de control cl&aacute;sica, adaptable y robusta, ante una planta de primer orden    con par&aacute;metros variables. </font>      <P><font size="2" face="Verdana">El sistema de control bajo estudio puede ser    representado de la forma expuesta en la <a href="#fig1">figura 1</a>. </font>     <P align="center"><img src="/img/revistas/rie/v33n3/f0102312.gif" width="404" height="131">    <a name="fig1"></a>     
<P>      <P><font size="2" face="Verdana">Para comparar el desempe&ntilde;o del sistema    de control se trabajar&aacute; con un modelo nominal de la planta con ganancia    unitaria y constante de tiempo igual a 0.2 segundos, es decir: </font>     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0302312.gif" width="390" height="56"></font>    <a name="e3"></a>      
<P><font size="2" face="Verdana">Y para realizar un estudio de la sensibilidad    de este sistema de control ante variaciones de los par&aacute;metros, as&iacute;    como ante la presencia de disturbios, se escogen los siguientes modelos con    par&aacute;metros alterados: </font>     <P align="center"><img src="/img/revistas/rie/v33n3/e0402312.gif" width="399" height="69">    <a name="e4"></a>     
<P align="center">&nbsp;     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0502312.gif" width="407" height="60"></font>    <a name="e5"></a>      
]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">La elecci&oacute;n de este rango est&aacute;    acorde a lo planteado en [5]. Aqu&iacute; se considera este intervalo de variaci&oacute;n    de par&aacute;metros como el caso m&aacute;s cr&iacute;tico para procesos fermentativos;    estos a su vez son los que presentan mayores alinealidades de las plantas de    inter&eacute;s mencionadas en la introducci&oacute;n. Teniendo en cuenta que    el dise&ntilde;o del sistema de control se hace para el caso nominal, los casos    de peor desempe&ntilde;o se tendr&aacute;n en los extremos del intervalo de    variaci&oacute;n de par&aacute;metros. </font>     <P>      <P><font size="2" face="Verdana"><b>Control cl&aacute;sico </b></font>      <P><font size="2" face="Verdana">En este tipo de control se emplean algoritmos    con estructura sencilla como redes de atraso-adelanto o reguladores de tipo    proporcional-integral-derivativo (PID) [3]. Por su parte el dise&ntilde;o est&aacute;    basado fundamentalmente en m&eacute;todos frecuenciales, donde son muy utilizados    los diagramas de Bode. Las especificaciones    de dise&ntilde;o dependen en lo fundamental de la aplicaci&oacute;n en s&iacute;;    para este trabajo se ha elegido realizar un dise&ntilde;o de forma tal que el    sistema resultante posea un margen de fase <i>&#947;</i>=60&#176; y una frecuencia    de cruce &#969;<sub>c</sub>=1 rad/s. El controlador que se propone emplear es    de tipo PI, cuya funci&oacute;n transferencial es la siguiente: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0602312.gif" width="408" height="63"><a name="e6"></a></font>      
<P><font size="2" face="Verdana">Este, a decir de la literatura, es muy recomendado    para este tipo de plantas, ya que la acci&oacute;n integral logra una respuesta    al paso, en lazo cerrado, con cero error de seguimiento, aspecto este muy deseable    en la inmensa mayor&iacute;a de los sistemas de control [3, 6]. El empleo de    controladores con acci&oacute;n derivativa no siempre es recomendable, debido    a la alta influencia del ruido que puede acompa&ntilde;ar las mediciones. Los    m&eacute;todos de dise&ntilde;o para calcular los par&aacute;metros del controlador    escogido pueden ser m&uacute;ltiples, por ejemplo m&oacute;dulo &oacute;ptimo,    m&eacute;todos anal&iacute;ticos, etc. [4, 7]. En este trabajo se emplea la    t&eacute;cnica de s&iacute;ntesis en el dominio de la frecuencia, donde la ganancia    y el cero del regulador se modifican iterativamente hasta cumplir las especificaciones    de dise&ntilde;o [3]. Con ella se lleg&oacute; al siguiente resultado: </font>      <P align="center"><img src="/img/revistas/rie/v33n3/e0a02312.gif" width="194" height="61">      
<P>      <P><font size="2" face="Verdana">de aqu&iacute;, teniendo en cuenta la funci&oacute;n    transferencia de un regulador PI dada por <a href="#e6">ecuaci&oacute;n (6)</a>,    se tiene que: </font>      <P align="center"><font size="2" face="Verdana"> <i>T<sub>I</sub></i> = 0.05&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<i>K<sub>c</sub>/T<sub>I</sub></i>    =&nbsp;7.07&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&#8756; <i>K<sub>c</sub></i></font>    = 0.354     ]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">La respuesta al paso del sistema de control a    lazo cerrado, empleando el modelo nominal de la planta, se comporta con cero    error de seguimiento, sobrepaso del 12%, dada por el margen de fase escogido,    y tiempo de establecimiento de 0.93 segundos, el cual queda determinado por    el ancho de banda que fue seleccionado. N&oacute;tese la relativa facilidad    con que el m&eacute;todo de dise&ntilde;o garantiza las especificaciones de    lazo abierto, no as&iacute; las de lazo cerrado las que son determinadas de    forma indirecta, aspecto que caracteriza a la teor&iacute;a de control cl&aacute;sica.    </font>      <P><font size="2" face="Verdana">En la <a href="#fig2">figura 2(a,b)</a>, se muestra    el estudio de sensibilidad de la estrategia de control cl&aacute;sica dise&ntilde;ada,    ante la variaci&oacute;n de par&aacute;metros definida por las plantas <i>G<sub>1</sub></i>(s)    y <i>G<sub>2</sub></i>(s). </font>     <P align="center"><img src="/img/revistas/rie/v33n3/f0202312.gif" width="541" height="187">    <a name="fig2"></a>     
<P>      <P><font size="2" face="Verdana">Obs&eacute;rvese el deterioro que sufre la respuesta    del sistema, la que en un caso cr&iacute;tico pudiese perder la estabilidad.    Se nota como a medida que los par&aacute;metros aumentan el sistema se vuelve    m&aacute;s oscilatorio, debido a que el margen de fase disminuye, y viceversa.    Si se requieren altas prestaciones din&aacute;micas,    el control cl&aacute;sico con reguladores de estructura sencilla y par&aacute;metros    fijos, no resuelve del todo el problema de control de esta planta variable con    el tiempo. Este problema dio pi&eacute; a innumerables investigaciones producto    de las cuales se han aplicado a esta planta m&uacute;ltiples esquemas de control,    a decir de la literatura: control con aprendizaje, control neuronal, control    borroso, etc. [8]. Estas variantes de control inteligente tienen la dificultad    de que para su implementaci&oacute;n, es indispensable disponer de potentes    medios de c&aacute;lculo, lo que encarece el sistema de control. Otras    variantes m&aacute;s sencillas desde el punto de vista de su implementaci&oacute;n    tambi&eacute;n han sido estudiadas, como son los casos de control adaptable    y control robusto, las cuales se analizan a continuaci&oacute;n. </font>      <P>      <P><font size="2" face="Verdana"><b>Control adaptable </b></font>      <P><font size="2" face="Verdana">Un regulador adaptable posee la cualidad de modificar    su operaci&oacute;n, en respuesta a cambios en la din&aacute;mica del proceso    o la presencia de disturbios. Existen varias t&eacute;cnicas que posibilitan    esta caracter&iacute;stica entre las que se encuentran: la sinton&iacute;a propia,    el inventario de ganancia y la basada en modelo de referencia. Particularmente    esta &uacute;ltima, es decir, los sistemas de control adaptables con modelo    de referencia (S.C.A.M.R.), son de amplio uso para estabilizar las caracter&iacute;sticas    din&aacute;micas de un sistema de control realimentado ante la variaci&oacute;n    de par&aacute;metros [3, 9]. El modelo    de referencia es el encargado de especificar la respuesta deseada que ha de    tener el sistema ajustable. De acuerdo al comportamiento de las entradas y los    estados de este se conforma el &iacute;ndice de desempe&ntilde;o, quedando para    el mecanismo de adaptaci&oacute;n la tarea de minimizar un funcional de la diferencia    entre la salida (o los estados del sistema ajustable) y los correspondientes    al modelo de referencia. Para dar soluci&oacute;n a esta problem&aacute;tica    el mecanismo de adaptaci&oacute;n modifica los par&aacute;metros del sistema    ajustable y/o genera una se&ntilde;al de entrada auxiliar. </font>      <P><font size="2" face="Verdana">M&eacute;todos basados en la variaci&oacute;n    de par&aacute;metros del controlador son reportados por [9]. Tienen la dificultad    que requieren resolver un gran n&uacute;mero de ecuaciones diferenciales (o    de diferencia). Esto hace que la adaptaci&oacute;n de par&aacute;metros sea    muy costosa computacionalmente y por tanto, poco factible. Por el contrario    las variantes que emplean s&iacute;ntesis de se&ntilde;al auxiliar de entrada    resultan generalmente en sistemas casi tan sencillos como los mismos reguladores.    A continuaci&oacute;n se realiza un estudio    de la estrategia de control adaptable reportada en [6]. Esta variante posee    la cualidad de que s&oacute;lo es necesaria la medici&oacute;n de una variable    de estado para construir la ley de control realimentada, lo que la hace muy    simple y poco costosa de implementar [2]. </font>      <P><font size="2" face="Verdana">Sea la planta y el modelo de referencia elegido,    tales que est&eacute;n descritos por: </font>     ]]></body>
<body><![CDATA[<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0702312.gif" width="411" height="60"></font>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0802312.gif" width="409" height="51"></font>      
<P><font size="2" face="Verdana">donde una elecci&oacute;n adecuada del modelo    de referencia ser&iacute;a tal que <i>x<sub>p</sub>&isin;R<sup>n</sup></i>,    <i>x<sub>m</sub>&isin;R<sup>n</sup></i>, <i>u<sub>p</sub>&isin;R<sup>p</sup></i>,    <i>u<sub>m</sub>&isin;R<sup>p</sup></i>, <i>y<sub>p</sub>&isin;R<sup>q</sup></i>,    <i>y<sub>m</sub>&isin;R<sup>q</sup></i> </font>      <P><font size="2" face="Verdana">La entrada de control se propone como: </font>     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0902312.gif" width="400" height="40"></font>    <a name="e9"></a>      
<P><font size="2" face="Verdana">donde <i>e<sub>0</sub></i> es el error de salida, el cual    es descrito por: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1002312.gif" width="414" height="46"></font>      
<P><font size="2" face="Verdana">Los vectores de ganancia se proponen como: </font>     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1102312.gif" width="412" height="49"></font>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1202312.gif" width="412" height="41"></font>      
]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">donde D es una matriz de ganancia, <i>K<sub>x0</sub>,    K<sub>u0</sub> y K<sub>e0</sub></i> denotan las ganancias de la parte lineal    del control de seguimiento y <i>&#916;K<sub>e</sub></i>(<i>v,t</i>) es la ganancia    de adaptaci&oacute;n aumentada para reducir el efecto de la variaci&oacute;n    de los par&aacute;metros. Para ganar en facilidad de implementaci&oacute;n se    propone una ley de adaptaci&oacute;n del tipo PI: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1302312.gif" width="412" height="60"></font>    <a name="e13"></a>      
<P><font size="2" face="Verdana">donde M<sub>i</sub> y R<sub>i</sub>, i=1, 2 son    matrices de ponderaci&oacute;n con dimensiones adecuadas. M<sub>1</sub>, R<sub>1</sub>    y R<sub>2</sub> son matrices definidas positivas y M2 es una matriz definida    positiva o no negativa. En la <a href="#fig3">figura 3</a>, se muestra un diagrama    de bloques del sistema hasta aqu&iacute; descrito. </font>      <P align="center"><img src="/img/revistas/rie/v33n3/f0302312.gif" width="463" height="414">    <a name="fig3"></a>     
<P>      <P><font size="2" face="Verdana">Despreciando el vector de ganancia de adaptaci&oacute;n    &#916;K<sub>e</sub>(v,t) en <a href="#e13">ecuaci&oacute;n (13)</a>, se puede    demostrar que es posible encontrar el control de seguimiento de modelo perfecto    eligiendo K<sub>e</sub> tal que (A<sub>p</sub>-B<sub>p</sub>K<sub>e</sub>C<sub>p</sub>)    sea una matriz de Hurwitz [9-10]. Adem&aacute;s se debe hacer: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1402312.gif" width="418" height="43"></font>    <a name="e14"></a>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1502312.gif" width="420" height="43"></font>    <a name="e15"></a>      
<P><font size="2" face="Verdana">donde B<sub>p</sub><sup>+</sup>=(B<sup>T</sup>B<sub>p</sub>)<sup>-1</sup>B<sub>p</sub><sup>T</sup>, la que es conocida como    pseudoinversa de la matriz B<sub>p</sub>. Debido a las variaciones de las condiciones de    operaci&oacute;n, los verdaderos par&aacute;metros del sistema pueden diferir    de los usados en el dise&ntilde;o del sistema adaptable, efecto que ser&aacute;    contrarrestado a trav&eacute;s de la ley de adaptaci&oacute;n. </font>     <P><font size="2" face="Verdana">En este caso, un sistema realimentado equivalente    para el sistema de control adaptable con modelo de referencia puede construirse    de las ecuaciones (7) a la (15) como: </font>     ]]></body>
<body><![CDATA[<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1602312.gif" width="406" height="50"></font>    <a name="e16"></a>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1702312.gif" width="406" height="49"></font>    <a name="e17"></a>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1802312.gif" width="528" height="50"></font>    <a name="e18"></a>      
<P><font size="2" face="Verdana">donde <a href="#e16">ecuaci&oacute;n (16)</a>    y <a href="#e17">ecuaci&oacute;n (17)</a> son bloques lineales de control de    adelanto, y <a href="#e18">ecuaci&oacute;n (18)</a> denota realimentaci&oacute;n    variable con el tiempo y no lineal. La estabilidad de este sistema se estudia    usando el teorema de Popov [6]. Obs&eacute;rvese    que debido a la parte integral que inherentemente existe en el mecanismo de    adaptaci&oacute;n de <a href="#e13">ecuaci&oacute;n (13)</a>, se puede contrarrestar    el efecto de cualquier disturbio que aparezca. </font>      <P><font size="2" face="Verdana">El modelo de referencia que se propone emplear    es tal que se iguale a la funci&oacute;n transferencial de la planta con el    controlador PI calculado, o sea: </font>     <P align="center"><img src="/img/revistas/rie/v33n3/e0b02312.gif" width="280" height="132">      
<P>      <P>      <P><font size="2" face="Verdana">De esta forma en presencia de condiciones nominales    tanto la planta como el modelo responder&aacute;n id&eacute;nticamente ante    una misma excitaci&oacute;n, por lo que el error ser&aacute; cero, no siendo    necesaria la generaci&oacute;n de la se&ntilde;al auxiliar de entrada. Por otro    lado seg&uacute;n se plantea en la <a href="#e14">ecuaci&oacute;n (14)</a>,    K<sub>x</sub> tendr&aacute; valor cero con lo que seg&uacute;n <a href="#e9">ecuaci&oacute;n    (9)</a>, no ser&aacute; necesario el c&aacute;lculo de la parte de la se&ntilde;al    de mando correspondiente a los estados del modelo, esto en definitiva disminuye    la complejidad computacional del sistema propuesto. De la <a href="#e15">ecuaci&oacute;n    (15)</a>, se tiene que K<sub>u</sub>=1, y para que (A<sub>p</sub>-B<sub>p</sub>K<sub>e</sub>C<sub>p</sub>)    sea una matriz de Hurwitz, se hace K<sub>e</sub>=1. </font><font size="2" face="Verdana">Las    matrices de ponderaci&oacute;n se han elegidos como sigue: D=5, M<sub>i</sub>=2,    i=1,2. </font>      <P><font size="2" face="Verdana">En la <a href="#fig4">figura 4(a,b)</a>, se muestran    las respuestas transitorias de los modelos nominal G<sub>p</sub>(s) y de par&aacute;metros    alterados G<sub>1</sub>(s) y G<sub>2</sub>(s).dado por las ecuaciones (<a href="#e3">3</a>,    <a href="#e4">4</a> y <a href="#e5">5</a>). Se evidencia la notable mejor&iacute;a    en la operaci&oacute;n del sistema tanto ante cambios en la entrada de referencia    como en los disturbios. Los resultados se deben comparar con los de la <a href="#fig2">figura    2</a>. El efecto integral en la ganancia    de adaptaci&oacute;n K<sub>x</sub>(v,t) depende del cuadrado del error por lo    que tendr&aacute; caracter&iacute;stica de funci&oacute;n mon&oacute;tona creciente    a medida que transcurre el tiempo. Esto implica que aun cuando se alcance un    valor tal que obligue a que el error tienda a cero, si se encuentra en una aplicaci&oacute;n    real no exenta de ruido en la medici&oacute;n, el valor de la integral continuar&aacute;    creciendo lo que en un momento u otro desestabilizar&aacute; el sistema. De    esta forma es necesario introducir un mecanismo de desconexi&oacute;n de la    adaptaci&oacute;n basado en alg&uacute;n criterio pr&aacute;ctico. Esto implica    un incremento en la cantidad de ecuaciones a solucionar y en las decisiones    a tomar en tiempo real, lo que va en detrimento del costo computacional. </font>      ]]></body>
<body><![CDATA[<P align="center"><img src="/img/revistas/rie/v33n3/f0402312.gif" width="554" height="170">    <a name="fig4"></a>      
<P>      <P>      <P><font size="2" face="Verdana"><b>Control Robusto</b> </font>      <P><font size="2" face="Verdana">Las estrategias conocidas como &quot;control    moderno&quot; no resultaron ser del todo satisfactorias, sin embargo dieron    pi&eacute; al surgimiento de la teor&iacute;a de control robusto. Como su nombre    lo indica, con ella se pretende obtener un dise&ntilde;o tal que el sistema    sea insensible a las incertidumbres en la modelaci&oacute;n y a los disturbios.    Esta parte de la teor&iacute;a del control    se encuentra en la actualidad en proceso de desarrollo, aunque ya se han obtenido    resultados concretos como para su aplicaci&oacute;n pr&aacute;ctica. Tal es    el caso de la variante de dise&ntilde;o de un controlador de dos grados de libertad    por un m&eacute;todo sistem&aacute;tico que no utiliza m&aacute;s que una simple    parametrizaci&oacute;n del regulador [11]. Esta variante es evaluada aqu&iacute;    para su comparaci&oacute;n con las estrategias antes presentadas, en cuanto    a desempe&ntilde;o y sensitividad, as&iacute; como en lo que respecta a costo    computacional y complejidad de dise&ntilde;o. Inicialmente    se revisar&aacute;n algunos conceptos. Sea el sistema de control mostrado en    la <a href="#fig5">figura 5(a,b)</a>. </font>      <P><font size="2" face="Verdana">Las distintas variables que aparecen en este    diagrama de bloques se denotan como:    <br>   &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </font><font size="2" face="Verdana">C<sub>A</sub>(s)    y C<sub>B</sub>(s) : controladores (el caso m&aacute;s conocido es cuando C<sub>B</sub>(s)=0).    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;P(s)    : planta.    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;r    : entrada de comando (referencia).    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y    : salida controlada.    ]]></body>
<body><![CDATA[<br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;d    : disturbios en el mando.    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;w    : disturbios en la salida.    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;n    : ruido en la medici&oacute;n. </font>      <P><font size="2" face="Verdana">Se definen: </font>     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e0c02312.gif" width="445" height="54"></font>      
<P><font size="2" face="Verdana">not&aacute;ndose que: </font>     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e1902312.gif" width="307" height="39"></font>    <a name="e19"></a>     
<P><font size="2" face="Verdana">La funci&oacute;n de transferencia S(s) se conoce    como funci&oacute;n de sensitividad y T(s) como funci&oacute;n de sensitividad    complementaria. Seg&uacute;n lo que se evidencia en la <a href="#e19">ecuaci&oacute;n    (19)</a>, si S(s) se minimiza en un rango de frecuencias, entonces T(s) debe    tener m&oacute;dulo cercano a 1 en ese mismo intervalo de frecuencias, y viceversa.    De ah&iacute; que el problema de control robusto sea lograr un regulador tal,    que optimice las funciones de sensitividad y sensitividad complementaria en    determinados rangos de frecuencias, de forma que el sistema sea insensible a    las incertidumbres de la modelaci&oacute;n y a la presencia de disturbios. Debido    a que los disturbios son por lo general de baja frecuencia, la funci&oacute;n    S(s) debe tener un m&oacute;dulo peque&ntilde;o en las bajas frecuencias y como    el ruido presente en la medici&oacute;n posee casi todas sus componentes ubicadas    en las altas frecuencias, la funci&oacute;n de sensitividad complementaria T(s)    debe dise&ntilde;arse de forma tal que sea peque&ntilde;a en las altas frecuencias.    Por el contrario |T(s)| debe ser lo m&aacute;s cercano a 1 en las bajas frecuencias    para lograr buenas condiciones de seguimiento en estado estable (cero error    ante una entrada constante por tramos). Por tanto, S(s) debe comportarse como    un filtro paso alto mientras que T(s) ha de poseer caracter&iacute;sticas de    filtro paso bajo. </font>      <P align="center"><img src="/img/revistas/rie/v33n3/f0502312.gif" width="451" height="377"><a name="fig5"></a>     
<P><font size="2" face="Verdana">Una de las posibles estructuras de un regulador    de dos grados de libertad se muestra en la <a href="#fig5">figura 5(a)</a>,    se pueden encontrar otras combinaciones [9, 12]. Estos controladores obtienen    la se&ntilde;al de mando u de la suma algebraica de dos componentes, en este    caso particular a partir del error de seguimiento y de la salida. La fundamentaci&oacute;n    matem&aacute;tica de ellos ha sido ampliamente estudiada, sin embargo es dif&iacute;cil    encontrar ejemplos de aplicaci&oacute;n satisfactoria a sistemas pr&aacute;cticos    [3]. El m&eacute;todo estudiado aqu&iacute;    logra un dise&ntilde;o robusto con magn&iacute;ficas cualidades y facilidad    de implementaci&oacute;n en comparaci&oacute;n con otros algoritmos de igual    prop&oacute;sito [11]. </font>      ]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">Para comenzar el estudio de este controlador    se introducen las siguientes definiciones: </font>     <P><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;R_(s)    : Anillo de funciones racionales propias y estables [13].    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;R<sub>s</sub>_(s)    : Anillo de funciones racionales estrictamente propias y estables.    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;P(s)    : Modelo real de la planta.    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;P<sub>n</sub>(s)    : Modelo nominal de la planta.    <br>   </font><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;G<sub>ry</sub>(s)    : Modelo deseado de la entrada de referencia r(s) a la salida y(s). </font>      <P><font size="2" face="Verdana">En este caso el sistema es de una entrada y una    salida. Adem&aacute;s, de forma general, se considera que tanto P(s) como P<sub>n</sub>(s)    son estrictamente propias, la entrada de referencia r(s) y los disturbios d(s)    tienen los mismos polos inestables y que el orden relativo de las funciones    racionales en R<sub>s</sub>_(s) es mayor que el de la planta. Las    factorizaciones coprimas de P<sub>n</sub>(s), r(s) y d(s) son definidas como:    P<sub>n</sub>(s)=N<sub>p</sub>/D<sub>p</sub>, r(s)=G<sub>r</sub>/F, y d(s)=G<sub>d</sub>/F,    donde Np&isin;Rs_ (s) y D<sub>p</sub>, G<sub>r</sub>, G<sub>d</sub>, F&isin;R<sub>s</sub>_(s).    </font>      <P><font size="2" face="Verdana">Usando las definiciones de arriba se obtienen    las siguientes expresiones para el controlador de dos grados de libertad [11]:    </font>     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2002312.gif" width="417" height="63"></font>    <a name="e20"></a>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2102312.gif" width="422" height="63"></font>    <a name="e21"></a>      
]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">Las condiciones de dise&ntilde;o est&aacute;n    dadas por: </font>     <P><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;A    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;G<sub>ry</sub>(s)&isin;R<sub>s</sub>_(s), G<sub>ry</sub>(s)/N<sub>p</sub>&isin;R_(s),    y (1-G<sub>ry</sub>(s))/F es una unidad en R_(s). Es decir, que el grado relativo    del modelo sea mayor &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;que    el de la planta. </font>      <P><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;B    &nbsp;&nbsp;&nbsp;&nbsp;Q&isin;R<sub>s</sub>_(s) y Q(s)/N<sub>p</sub>, (1-Q(s))/D<sub>p</sub>&isin;R_(s).    O sea, todos los ceros y polos inestables de la planta deben incluirse en Q(s)    y 1-Q(s) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;respectivamente,    tal que se garantice estabilidad interna del sistema de control y que el grado    relativo de Q(s) tiene que ser mayor o &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;igual    que el grado relativo de P(s) [13]. </font>      <P><font size="2" face="Verdana">La parametrizaci&oacute;n del controlador posee    dos elementos libres: G<sub>ry</sub>(s) y Q(s). Por su parte, G<sub>ry</sub>(s) debe ser escogida    tal que se cumpla la condici&oacute;n de dise&ntilde;o A, mientras que Q(s)    es el m&aacute;s significativo, ya que es quien determina la robustez y el desempe&ntilde;o    del sistema ante disturbios y se elige para que cumpla con B. </font>     <P><font size="2" face="Verdana">En la <a href="#fig5">figura 5(b)</a>, se muestra    la estructura equivalente del sistema de control derivado de ecuaciones (<a href="#e20">20</a>    y <a href="#e21">21</a>), donde: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2202312.gif" width="416" height="65"></font>      
<P><font size="2" face="Verdana">Bajo condiciones nominales y en ausencia de disturbios,    el error e=u'-u es igual a cero. La respuesta a la entrada de referencia coincide,    por tanto, con la especificada por G<sub>ry</sub>(s) independientemente del par&aacute;metro    Q(s). Si la planta difiere de su modelo nominal, o se a&ntilde;aden disturbios    al sistema, entonces e&ne;0 pero es forzado a converger a cero a trav&eacute;s    del lazo de realimentaci&oacute;n formado por Q(s), de aqu&iacute; la importancia    que reviste una adecuada elecci&oacute;n de este par&aacute;metro. </font>     <P><font size="2" face="Verdana">Para el sistema de la figura 5(a), con C<sub>A</sub>(s)    y C<sub>B</sub>(s) dados por las ecuaciones (<a href="#e20">20</a> y <a href="#e21">21</a>),    la funci&oacute;n de sensitividad S(s) y la funci&oacute;n de respuesta a disturbios    G<sub>dy</sub>(s) est&aacute;n dadas por: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2302312.gif" width="413" height="44"></font>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2402312.gif" width="416" height="46"></font>      
]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">Es v&aacute;lido recordar que S(s) es la funci&oacute;n    transferencial entre el disturbio d(s) y la salida y(s) y representa la robustez    contra variaciones en el sistema. Se supone que estas se describen mediante    incertidumbre multiplicativa [13]. Entonces se tiene que la respuesta a la entrada    de comando G<sub>ry</sub>(s) y el desempe&ntilde;o en cuanto a eliminaci&oacute;n de disturbios    G<sub>dy</sub>(s) ser&aacute;n cambiados por: </font>     <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2502312.gif" width="419" height="60"></font>      
<P><font size="2" face="Verdana">Esto implica que la norma |S(j&#969; )| debe    ser suficientemente peque&ntilde;a para alcanzar una elevada robustez ante variaci&oacute;n    de par&aacute;metros. Claro est&aacute;, que seg&uacute;n lo planteado en la    igualdad <a href="#e19">ecuaci&oacute;n (19)</a>, es imposible mantener dicha    norma peque&ntilde;a en todo el rango de frecuencias, por lo que, teniendo en    cuenta que los disturbios son por lo general de baja frecuencia, se dise&ntilde;a    Q(s) tal que S(s) posea caracter&iacute;sticas de filtro paso alto. </font>      <P><font size="2" face="Verdana">La otra caracter&iacute;stica a ser considerada    es la funci&oacute;n de sensitividad complementaria, la cual s obtiene a partir    de <a href="#e19">ecuaci&oacute;n (19)</a> y toma la forma: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2602312.gif" width="429" height="45"></font>      
<P><font size="2" face="Verdana">Esta tendr&iacute;a caracter&iacute;sticas de    filtro paso bajo, lo cual es deseable para rechazar el ruido en la medici&oacute;n.    La pregunta ser&iacute;a: &#191;C&oacute;mo elegir el par&aacute;metro Q(s)?.    Innumerables respuestas pudiese tener la    misma. En [11] se demuestra, aplicando optimizaci&oacute;n de una norma H&infin;,    que si el ancho de banda del par&aacute;metro Q(s) es mucho mayor que el correspondiente    al comportamiento elegido para el sistema en lazo cerrado G<sub>ry</sub>(s),    entonces se cumplen las siguientes expresiones: </font>      <P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2702312.gif" width="425" height="46"></font>    <a name="e27"></a>      
<P align="center"><font size="2" face="Verdana"> <img src="/img/revistas/rie/v33n3/e2802312.gif" width="426" height="49"></font>    <a name="e28"></a>      
<P><font size="2" face="Verdana">Esto implica que las funciones de sensitividad    y sensitividad complementaria pueden ser determinadas de una manera directa,    dise&ntilde;ando el par&aacute;metro Q(s). A    continuaci&oacute;n se aplica esta metodolog&iacute;a al caso bajo estudio,    donde el modelo nominal de la planta es: </font>      <P align="center"><img src="/img/revistas/rie/v33n3/img/e0d02312.gif" width="151" height="58">      
]]></body>
<body><![CDATA[<P>      <P><font size="2" face="Verdana">Se elige G<sub>ry</sub>(s) como la funci&oacute;n    de transferencia del sistema en lazo cerrado formado por la planta y el controlador    PI calculado inicialmente (igual al modelo de referencia del control adaptable),    y se han escogido as&iacute; intencionalmente para facilitar la comparaci&oacute;n    con las estrategias antes evaluadas. De esta forma: </font>     <P align="center"><img src="/img/revistas/rie/v33n3/e0e02312.gif" width="250" height="57">      
<P>      <P><font size="2" face="Verdana">el cual cumple con la condici&oacute;n de dise&ntilde;o    A. El ancho de banda de este sistema es    aproximadamente igual a 6 rad/s. Se escoge Q(s) tal que posea un ancho de banda    diez veces mayor, con el objetivo de garantizar el cumplimiento de las ecuaciones    (<a href="#e27">27</a> y <a href="#e28">28</a>). En la elecci&oacute;n de Q(s)    debe tenerse en cuenta adem&aacute;s la condici&oacute;n de dise&ntilde;o B,    o sea: </font>      <P align="center"><img src="/img/revistas/rie/v33n3/e0f02312.gif" width="238" height="64">      
<P>      <P>      <P><font size="2" face="Verdana">Cumpliendo con lo anterior, se escoge: </font>     <P align="center"><img src="/img/revistas/rie/v33n3/e0g02312.gif" width="241" height="45">      
]]></body>
<body><![CDATA[<P>      <P><font size="2" face="Verdana">la que no es m&aacute;s que un filtro paso bajo    de Butterworth, con ancho de banda de 60 rad/s. Con    las ecuaciones (<a href="#e20">20</a> y <a href="#e21">21</a>), se calculan    C<sub>A</sub>(s) y C<sub>B</sub>(s), resultando: </font>      <P align="center"><img src="/img/revistas/rie/v33n3/e0h02312.gif" width="296" height="105">      
<P>      <P>      <P><font size="2" face="Verdana">La <a href="#fig6">figura 6(a,b)</a>, muestra    los resultados de simulaci&oacute;n. Se brinda la respuesta ante un paso en    la entrada de referencia 6(a) y en los disturbios 6(b), de los modelos nominal    G<sub>p</sub>(s) y de par&aacute;metros alterados G<sub>1</sub>(s) y G<sub>2</sub>(s)    dado por las ecuaciones (<a href="#e3">3</a>, <a href="#e4">4</a> y <a href="#e5">5</a>).    </font>      <P align="center"><img src="/img/revistas/rie/v33n3/f0602312.gif" width="546" height="166">    <a name="fig6"></a>     
<P>      <P><font size="2" face="Verdana">Se evidencia la notable mejor&iacute;a en la    operaci&oacute;n del sistema, si se compara con los resultados del control cl&aacute;sico,    tanto ante cambios en la entrada de referencia como en los disturbios, los que    son casi eliminados por completo. Los resultados se deben comparar con los de    las figuras (<a href="#fig2">2</a> y <a href="#fig4">4</a>). Las escalas se    han igualado intencionalmente a las de estas para facilitar la comparaci&oacute;n.    </font>      <P>      ]]></body>
<body><![CDATA[<P><font size="2" face="Verdana"><b>RESULTADOS Y DISCUSI&Oacute;N</b> </font>      <P><font size="2" face="Verdana">A continuaci&oacute;n se realizar&aacute; una    comparaci&oacute;n de los resultados de simulaci&oacute;n expuestos para cada    estrategia, atendiendo a &iacute;ndices de comportamiento que permitan obtener    criterios sobre su operaci&oacute;n. Se ha elegido en el caso de la respuesta    al paso en la entrada de comando, comparar en cuanto a sobrepaso M<sub>p</sub>    y tiempo de establecimiento t<sub>s</sub> (seg&uacute;n el criterio de 2%),    par&aacute;metros estos que se utilizaron como especificaciones indirectas en    el dise&ntilde;o del control cl&aacute;sico. En el caso de la respuesta a disturbios    se utilizan como criterio la desviaci&oacute;n m&aacute;xima que sufre la se&ntilde;al    de salida, y el tiempo de recuperaci&oacute;n (tiempo necesario para que el    sistema entre en una gama del 2% luego de presentarse el disturbio). Los    resultados se muestran en forma tabular (tablas <a href="#tabla1">1</a>,<a href="#tabla2">2</a>,<a href="#tabla3">3</a>),    present&aacute;ndose los correspondientes a los modelos nominal G<sub>p</sub>(s)    y de par&aacute;metros alterados G<sub>1</sub>(s) y G<sub>2</sub>(s). Tambi&eacute;n    se incluye la cantidad de operaciones de suma algebraica, multiplicaci&oacute;n,    e integraci&oacute;n, que fueron necesarios incluir en los modelos de simulaci&oacute;n.    Estos datos dan una idea de la complejidad de implementaci&oacute;n, ya que    independientemente de las caracter&iacute;sticas de la misma (electr&oacute;nica    anal&oacute;gica o programada), la cantidad de hardware o software a incluir    ser&aacute; proporcional a la cantidad de operaciones matem&aacute;ticas de    la simulaci&oacute;n. Por otro lado se dan criterios cualitativos de la operaci&oacute;n    en cuanto a robustez y desempe&ntilde;o del sistema. </font>      <P>      <P align="center"><img src="/img/revistas/rie/v33n3/t0102312.gif" width="504" height="204">    <a name="tabla1"></a>      
<P>&nbsp;     <P align="center"><img src="/img/revistas/rie/v33n3/t0202312.gif" width="506" height="203">    <a name="tabla2"></a>     
<P>&nbsp;     <P align="center"><img src="/img/revistas/rie/v33n3/t0302312.gif" width="504" height="207">    <a name="tabla3"></a>     
<P>      <P><font size="2" face="Verdana">Del an&aacute;lisis de las tablas anteriores    se observa que: </font>     ]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&nbsp;&#8226; &nbsp;En presencia    de condiciones nominales, las estrategias de control adaptable y robusto mejoran    considerablemente la operaci&oacute;n del sistema &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ante    disturbios.    <br>   &nbsp;&nbsp;&nbsp; </font><font size="2" face="Verdana">&#8226; &nbsp;Cuando    los par&aacute;metros var&iacute;an, las estrategias adaptable y robusta logran    mantener los indicadores de desempe&ntilde;o din&aacute;mico Mp y ts cercano    al &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;valor de condiciones nominales.    <br>   &nbsp;&nbsp;&nbsp;&nbsp;</font><font size="2" face="Verdana">&#8226; &nbsp;Aun    con variaci&oacute;n de par&aacute;metros, es mejor la respuesta ante disturbios    de los sistemas adaptable y robusto, que en el caso nominal del control &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;cl&aacute;sico.    </font>      <P><font size="2" face="Verdana">&nbsp;Como fue reflejado en los comentarios conclusivos    de la evaluaci&oacute;n de cada estrategia: </font>      <P><font size="2" face="Verdana">&nbsp;&nbsp;&nbsp;&#8226; &nbsp;El control adaptable    garantiza buen desempe&ntilde;o, pero se puede perder en presencia de inestabilidad.    <br>   &nbsp;&nbsp; </font><font size="2" face="Verdana">&#8226; &nbsp;El control robusto    garantiza tanto robustez como excelente desempe&ntilde;o, sin riesgos de inestabilidad.    </font>      <P><font size="2" face="Verdana">Cu&aacute;l estrategia elegir, depende de la    relaci&oacute;n de compromiso entre los requerimientos que se planteen para    el sistema, as&iacute; como de las posibilidades de implementaci&oacute;n. O    sea, es la aplicaci&oacute;n quien determina. </font>     <P>&nbsp;      <P>      <P><font size="3" face="Verdana"><b>CONCLUSIONES</b></font><font size="2" face="Verdana">    </font>      ]]></body>
<body><![CDATA[<P><font size="2" face="Verdana">Mediante el an&aacute;lisis de los resultados    aqu&iacute; expuestos, se concluye que la t&eacute;cnica de control cl&aacute;sico    se vuelve incapaz de garantizar un elevado desempe&ntilde;o din&aacute;mico    y robustez, en plantas con modelos de par&aacute;metros variables. Con la aplicaci&oacute;n    del control adaptable se consigue insensibilidad ante variaciones de par&aacute;metros    de la planta, aunque con riesgos de estabilidad, que al tratar de ser evitados,    elevan el costo computacional de la implementaci&oacute;n. Por &uacute;ltimo,    el control robusto permite un dise&ntilde;o de elevada calidad, garantizando    por un lado la robustez a cambios en los par&aacute;metros del sistema, y por    el otro un alto desempe&ntilde;o din&aacute;mico. Estas caracter&iacute;sticas    se consiguen adem&aacute;s con un costo computacional m&aacute;s bajo que con    la aplicaci&oacute;n del control adaptable. </font>     <P>&nbsp;      <P>      <P><font size="3" face="Verdana"><b>REFERENCIAS</b> </font>     <P>&nbsp;      <P><font size="2" face="Verdana">1. Nelles, O., &quot;Nonlinear system identification:    from classical approaches to neural networks and fuzzy models&quot;. Berlin:    Springer-Verlag, 2000,     <br>   ISBN 978-3-540-67369-9. </font>      <!-- ref --><P><font size="2" face="Verdana">2. Sira, H., &quot;Fast Adaptive Trajectory Tracking    Control for a Completely Uncertain DC Motor via Output Feedback&quot;. Computaci&oacute;n    y Sistemas, 2009, vol.12, n.4, p. 397-408, ISSN 1405-5546.     </font>      <!-- ref --><P><font size="2" face="Verdana">3. Ogata, K., &quot;Ingenier&iacute;a de Control    Moderna&quot;. 4ta Edici&oacute;n, Madrid: Pearson Educaci&oacute;n S. A., 2003,    </font> <font size="2" face="Verdana">ISBN 84-205-3678-1.     </font>      <!-- ref --><P><font size="2" face="Verdana">4. Purwanto, E., <i>et al.</i>, &quot;Application    of vector control method for developing ANFIS observer as speed sensor for induction    motor speed control in electric vehicle&quot;. Academic Research International,    2012, vol.21, n.1, p. 22-29, ISSN 2223-9944.     </font>      <!-- ref --><P><font size="2" face="Verdana">5. Wang, H.P., <i>et al.</i>, &quot;Composed    Adaptive Control for a second-order nonlinear model of a biotechnological process&quot;.    En: 19th Mediterranean Conference on Control and Automation. Corfu, Greece,    2011. p. 1140-1143, ISBN 978-1-4577-0124-5.     </font>      <!-- ref --><P><font size="2" face="Verdana">6. Lin, F.; Liaw, C., &quot;Reference Model and    Adaptive Control for Induction Motors Drives&quot;. IEEE Transaction on Automatic    Control, 1993, vol.38, n.10, p. 1594-1600, ISSN 0018-9286.     </font>      <!-- ref --><P><font size="2" face="Verdana">7. L. Wang, &quot;Model predictive control system    design and implementation using MATLAB&quot;. London: Springer-Verlag, 2009,    ISBN 978-1-84882-330-3.     </font>      <!-- ref --><P><font size="2" face="Verdana">8. Angelov, P., &quot;Evolving Intelligent Systems:    Methodology and Applications&quot;. New Jersey: Wiley-IEEE Press, 2010, ISBN    978-0-470-28719-4.     </font>      <!-- ref --><P><font size="2" face="Verdana">9. Tao, G., &quot;Adaptive Control: Design and    Analysis&quot;. New Jersey: John Wiley &amp; Sons Inc., 2003, </font><font size="2" face="Verdana">ISBN    0-471-27452-6.     </font>      <!-- ref --><P><font size="2" face="Verdana">10. Zhou, H., &quot;State feedback sliding mode    control without chattering constructing Hurwitz matrix for AUV movement&quot;. International Journal of Automation and Computing, 2011,    vol.8. n.2, p. 262-268, </font><font size="2" face="Verdana">ISSN 1476-8186.        </font>      <!-- ref --><P><font size="2" face="Verdana">11. Umeno, T., &quot;Robust Servosystem Design    with Two Degrees of Freedom and is Application to Novel Motion Control of Robot    Manipulators&quot;. IEEE Transaction on Automatic Control, 1993, vol.40. n.5,    p. 473-485, </font><font size="2" face="Verdana">ISSN 0278-0046.     </font>      <!-- ref --><P><font size="2" face="Verdana">12. Kishimoto, Y., &quot;Disturbance attenuation    property for random disturbances in two-degree-of freedom optimal servo control    systems&#168;. En: 2011 International Conference on Advanced Mechatronic Systems    (ICAMechS), Zhengzhou, China, 2011. p. 321-324, ISBN 978-1-4577-1698-0.     </font>      <!-- ref --><P><font size="2" face="Verdana">13. Lin, F., &quot;Robust control design: an    optimal approach&quot;. Chichester: John Wiley &amp; Sons Ltd., 2007, </font><font size="2" face="Verdana">ISBN    978-0-470-03191-9.     </font>      <P>&nbsp;     <P>&nbsp;     <P><font size="2" face="Verdana">Recibido: Febrero del 2012    <br>   </font><font size="2" face="Verdana">Aprobado: Julio del 2012 </font>     <P>&nbsp;     <P>&nbsp;     <P><font size="2" face="Verdana">Orlando Regal&oacute;n Anias, Ingeniero en Control    Autom&aacute;tico, M&aacute;ster en Ingenier&iacute;a en Control Autom&aacute;tico,    Profesor Asistente. Departamento Ing. El&eacute;ctrica de la Facultad Electromec&aacute;nica.    Universidad de Camag&uuml;ey, Cuba. e-mail: <a href="mailto:orlando.regalon@reduc.edu.cu">orlando.regalon@reduc.edu.cu</a>    </font>       ]]></body><back>
<ref-list>
<ref id="B1">
<label>1</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Nelles]]></surname>
<given-names><![CDATA[O]]></given-names>
</name>
</person-group>
<source><![CDATA[Nonlinear system identification: from classical approaches to neural networks and fuzzy models]]></source>
<year>2000</year>
<publisher-loc><![CDATA[Berlin ]]></publisher-loc>
<publisher-name><![CDATA[Springer-Verlag]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B2">
<label>2</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Sira]]></surname>
<given-names><![CDATA[H]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[Fast Adaptive Trajectory Tracking Control for a Completely Uncertain DC Motor via Output Feedback]]></article-title>
<source><![CDATA[Computación y Sistemas]]></source>
<year>2009</year>
<volume>12</volume>
<numero>4</numero>
<issue>4</issue>
<page-range>397-408</page-range></nlm-citation>
</ref>
<ref id="B3">
<label>3</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Ogata]]></surname>
<given-names><![CDATA[K]]></given-names>
</name>
</person-group>
<source><![CDATA[Ingeniería de Control Moderna]]></source>
<year>2003</year>
<edition>4ta</edition>
<publisher-loc><![CDATA[Madrid ]]></publisher-loc>
<publisher-name><![CDATA[Pearson Educación S. A.]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B4">
<label>4</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Purwanto]]></surname>
<given-names><![CDATA[E]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[Application of vector control method for developing ANFIS observer as speed sensor for induction motor speed control in electric vehicle]]></article-title>
<source><![CDATA[Academic Research International]]></source>
<year>2012</year>
<volume>21</volume>
<numero>1</numero>
<issue>1</issue>
<page-range>22-29</page-range></nlm-citation>
</ref>
<ref id="B5">
<label>5</label><nlm-citation citation-type="">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Wang]]></surname>
<given-names><![CDATA[H.P]]></given-names>
</name>
</person-group>
<source><![CDATA[Composed Adaptive Control for a second-order nonlinear model of a biotechnological process]]></source>
<year>2011</year>
<page-range>1140-1143</page-range><publisher-loc><![CDATA[Corfu ]]></publisher-loc>
</nlm-citation>
</ref>
<ref id="B6">
<label>6</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Lin]]></surname>
<given-names><![CDATA[F]]></given-names>
</name>
<name>
<surname><![CDATA[Liaw]]></surname>
<given-names><![CDATA[C]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[Reference Model and Adaptive Control for Induction Motors Drives]]></article-title>
<source><![CDATA[IEEE Transaction on Automatic Control]]></source>
<year>1993</year>
<volume>38</volume>
<numero>10</numero>
<issue>10</issue>
<page-range>1594-1600</page-range></nlm-citation>
</ref>
<ref id="B7">
<label>7</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[L]]></surname>
<given-names><![CDATA[Wang]]></given-names>
</name>
</person-group>
<source><![CDATA[Model predictive control system design and implementation using MATLAB]]></source>
<year>2009</year>
<publisher-loc><![CDATA[London ]]></publisher-loc>
<publisher-name><![CDATA[Springer-Verlag]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B8">
<label>8</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Angelov]]></surname>
<given-names><![CDATA[P]]></given-names>
</name>
</person-group>
<source><![CDATA[Evolving Intelligent Systems: Methodology and Applications]]></source>
<year>2010</year>
<publisher-loc><![CDATA[New Jersey ]]></publisher-loc>
<publisher-name><![CDATA[Wiley-IEEE Press]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B9">
<label>9</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Tao]]></surname>
<given-names><![CDATA[G]]></given-names>
</name>
</person-group>
<source><![CDATA[Adaptive Control: Design and Analysis]]></source>
<year>2003</year>
<publisher-loc><![CDATA[New Jersey ]]></publisher-loc>
<publisher-name><![CDATA[John Wiley & Sons Inc.]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B10">
<label>10</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Zhou]]></surname>
<given-names><![CDATA[H]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[State feedback sliding mode control without chattering constructing Hurwitz matrix for AUV movement]]></article-title>
<source><![CDATA[International Journal of Automation and Computing]]></source>
<year>2011</year>
<volume>8</volume>
<numero>2</numero>
<issue>2</issue>
<page-range>262-268</page-range></nlm-citation>
</ref>
<ref id="B11">
<label>11</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Umeno]]></surname>
<given-names><![CDATA[T]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA[Robust Servosystem Design with Two Degrees of Freedom and is Application to Novel Motion Control of Robot Manipulators]]></article-title>
<source><![CDATA[IEEE Transaction on Automatic Control]]></source>
<year>1993</year>
<volume>40</volume>
<numero>5</numero>
<issue>5</issue>
<page-range>473-485</page-range><page-range>0278-0046</page-range></nlm-citation>
</ref>
<ref id="B12">
<label>12</label><nlm-citation citation-type="">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Kishimoto]]></surname>
<given-names><![CDATA[Y]]></given-names>
</name>
</person-group>
<source><![CDATA[Disturbance attenuation property for random disturbances in two-degree-of freedom optimal servo control systems]]></source>
<year>2011</year>
<page-range>321-324</page-range><publisher-loc><![CDATA[Zhengzhou ]]></publisher-loc>
</nlm-citation>
</ref>
<ref id="B13">
<label>13</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Lin]]></surname>
<given-names><![CDATA[F]]></given-names>
</name>
</person-group>
<source><![CDATA[Robust control design: an optimal approach]]></source>
<year>2007</year>
<publisher-loc><![CDATA[Chichester ]]></publisher-loc>
<publisher-name><![CDATA[John Wiley & Sons Ltd.]]></publisher-name>
</nlm-citation>
</ref>
</ref-list>
</back>
</article>
