<?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>0122-3461</journal-id>
<journal-title><![CDATA[Ingeniería y Desarrollo]]></journal-title>
<abbrev-journal-title><![CDATA[Ing. Desarro.]]></abbrev-journal-title>
<issn>0122-3461</issn>
<publisher>
<publisher-name><![CDATA[Fundación Universidad del Norte]]></publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id>S0122-34612014000100008</article-id>
<title-group>
<article-title xml:lang="es"><![CDATA[Filtro de Kalman dual de estadística robusta para la identificación de un generador eòlico en presencia de datos espurios]]></article-title>
<article-title xml:lang="en"><![CDATA[Robust statistic dual Kalman filter for wind generator identification in presence of outliers]]></article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Zuluaga Ríos]]></surname>
<given-names><![CDATA[Carlos David]]></given-names>
</name>
<xref ref-type="aff" rid="A01"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Giraldo Suárez]]></surname>
<given-names><![CDATA[Eduardo]]></given-names>
</name>
<xref ref-type="aff" rid="A02"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname><![CDATA[Álvarez López]]></surname>
<given-names><![CDATA[Mauricio Alexander]]></given-names>
</name>
<xref ref-type="aff" rid="A03"/>
</contrib>
</contrib-group>
<aff id="A01">
<institution><![CDATA[,Universidad Tecnológica de Pereira  ]]></institution>
<addr-line><![CDATA[Pereira Risaralda]]></addr-line>
<country>Colombia</country>
</aff>
<aff id="A02">
<institution><![CDATA[,Universidad Tecnológica de Pereira  ]]></institution>
<addr-line><![CDATA[Pereira Risaralda]]></addr-line>
<country>Colombia</country>
</aff>
<aff id="A03">
<institution><![CDATA[,Universidad Tecnológica de Pereira  ]]></institution>
<addr-line><![CDATA[Pereira Risaralda]]></addr-line>
<country>Colombia</country>
</aff>
<pub-date pub-type="pub">
<day>00</day>
<month>01</month>
<year>2014</year>
</pub-date>
<pub-date pub-type="epub">
<day>00</day>
<month>01</month>
<year>2014</year>
</pub-date>
<volume>32</volume>
<numero>1</numero>
<fpage>115</fpage>
<lpage>137</lpage>
<copyright-statement/>
<copyright-year/>
<self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_arttext&amp;pid=S0122-34612014000100008&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_abstract&amp;pid=S0122-34612014000100008&amp;lng=en&amp;nrm=iso"></self-uri><self-uri xlink:href="http://www.scielo.org.co/scielo.php?script=sci_pdf&amp;pid=S0122-34612014000100008&amp;lng=en&amp;nrm=iso"></self-uri><abstract abstract-type="short" xml:lang="es"><p><![CDATA[Las turbinas eólicas son una de las tecnologías más avanzadas dentro del conjunto de las energías limpias para obtener energía eléctrica. Debido al comportamiento estocástico del viento, se debe realizar el control de estas turbinas para maximizar la potencia de salida. Con el fin de solucionar el problema de control óptimo se debe conocer la dinámica de la máquina, por ello, los esquemas de identificación de sistemas adquieren importancia en este contexto. Este estudio presenta la identificación de un generador eólico que emplea un filtro de Kalman dual de estadística robusta. Aunque el filtro de Kalman dual (DKF) ha sido utilizado en la identificación de sistemas, su rendimiento está asociado a la ausencia de datos espurios. Para realizar la identificación de este sistema se utilizó una tarjeta de adquisición NI USB-6009. Pese a sus buenas características de muestreo, esta tarjeta incorpora datos espurios que afectan la identificación del sistema. El estudio propone solucionar este problema combinando el DKF y el filtro de Kalman de estadística robusta. La metodología propuesta o fue implementada en Matlab, mostrando la identificación del generador eólico en tiempo real ante variaciones en las entradas del sistema, con un comportamiento robusto en presencia de datos espurios, comparado con otra metodología estándar de identificación.]]></p></abstract>
<abstract abstract-type="short" xml:lang="en"><p><![CDATA[Among the technologies for clean energy generation, wind turbines are one of the most advanced. Due to stochastic behavior of the wind, a proper control action must be exerted over these turbines in order to get a better use of the output power. A propel control action usually implies the knowledge of the system dynamics, and that is why the identification stage is also crucial. In this paper, the authors propose a robust statistics dual Kalman filter for identification of a wind generator. Although dual Kalman filter (DKF) has been used in system identification, its performance heavily depends on the absence of outliers in the measured data. However, outliers might be easily produced during the data acquisition stage. In this study, the authors show how by combining the DKF and the robust statistic Kalman filter, the problem of outliers in systems identification can be avoided. The method was implemented in Matlab R2009a®. Results for the identification of a wind generator working online are shown and compared to the performance of dual Kalman filter.]]></p></abstract>
<kwd-group>
<kwd lng="es"><![CDATA[datos espurios]]></kwd>
<kwd lng="es"><![CDATA[filtro de Kalman]]></kwd>
<kwd lng="es"><![CDATA[generador eólico]]></kwd>
<kwd lng="es"><![CDATA[identificación de sistemas]]></kwd>
<kwd lng="en"><![CDATA[outliers]]></kwd>
<kwd lng="en"><![CDATA[Kalman filter]]></kwd>
<kwd lng="en"><![CDATA[wind generator]]></kwd>
<kwd lng="en"><![CDATA[system identification]]></kwd>
</kwd-group>
</article-meta>
</front><body><![CDATA[   <font face="verdana" size="2">     <p><b>ART&Iacute;CULO DE INVESTIGACI&Oacute;N / </b>RESEARCH ARTICLE</b></p> <font size="4">    <p align="center"><b>Filtro de Kalman dual de estad&iacute;stica robusta para la identificaci&oacute;n de un generador eòlico en presencia de datos espurios</b></p></font> <font size="3">    <p align="center"><b>Robust statistic dual Kalman filter for wind generator identification in presence of outliers</b></p></font>      <p>Carlos David Zuluaga R&iacute;os*    <br>  Eduardo Giraldo Su&aacute;rez**    <br> Mauricio Alexander &Aacute;lvarez L&oacute;pez***    <br> <i>Universidad Tecnol&oacute;gica de Pereira</i></p>     <p>* Universidad Tecnol&oacute;gica de Pereira, Programa de Ingenier&iacute;a El&eacute;ctrica, joven investigador, mag&iacute;ster en Ingenier&iacute;a El&eacute;ctrica. Pereira, Risaralda, Colombia. <a href="mailto:cardazu@utp.edu.co"><i>cardazu@utp.edu.co</i></a></p>     <p><i>** </i>Universidad Tecnol&oacute;gica de Pereira, Programa de Ingenier&iacute;a El&eacute;ctrica, profesor asociado, Ph.D. en Ingenier&iacute;a.Pereira, Risaralda, Colombia.  <a href="mailto:egiraldos@utp.edu.co"><i>egiraldos@utp.edu.co</i></a></p>     ]]></body>
<body><![CDATA[<p>*** Universidad Tecnol&oacute;gica de Pereira, Programa de Ingenier&iacute;a El&eacute;ctrica, profesor asociado, Ph.D. en Ciencias de la Computaci&oacute;n. Pereira, Risaralda, Colombia. <i><a href="mailto:malvarez@utp.edu.co">malvarez@utp.edu.co</a></i></p>     <p><b>Correspondencia: </b>Carlos David Zuluaga R&iacute;os, Universidad Tecnol&oacute;gica de Pereira, vereda La Julita, Edificio de Ingenier&iacute;a El&eacute;ctrica, Oficina E-230, Pereira, Risaralda, Colombia. Tel: 57 (6) 3137 122 ext. 109.</p>      <p><b>Fecha de recepci&oacute;n:</b> 5 de agosto de 2013    <br> <b>Fecha de aceptaci&oacute;n:</b> 4 de diciembre de 2013</p>  <hr>      <p><b>Resumen</b></p>      <p>Las turbinas e&oacute;licas son una de las tecnolog&iacute;as m&aacute;s avanzadas dentro del conjunto de las energ&iacute;as limpias para obtener energ&iacute;a el&eacute;ctrica. Debido al comportamiento estoc&aacute;stico del viento, se debe realizar el control de estas turbinas para maximizar la potencia de salida. Con el fin de solucionar el problema de control &oacute;ptimo se debe conocer la din&aacute;mica de la m&aacute;quina, por ello, los esquemas de identificaci&oacute;n de sistemas adquieren importancia en este contexto. Este estudio presenta la identificaci&oacute;n de un generador e&oacute;lico que emplea un filtro de Kalman dual de estad&iacute;stica robusta. Aunque el filtro de Kalman dual (DKF) ha sido utilizado en la identificaci&oacute;n de sistemas, su rendimiento est&aacute; asociado a la ausencia de datos espurios. Para realizar la identificaci&oacute;n de este sistema se utiliz&oacute; una tarjeta de adquisici&oacute;n NI USB-6009. Pese a sus buenas caracter&iacute;sticas de muestreo, esta tarjeta incorpora datos espurios que afectan la identificaci&oacute;n del sistema. El estudio propone solucionar este problema combinando el DKF y el filtro de Kalman de estad&iacute;stica robusta. La metodolog&iacute;a propuesta o   fue implementada en Matlab, mostrando la identificaci&oacute;n del generador e&oacute;lico en tiempo real ante variaciones en las entradas del sistema, con un comportamiento robusto en presencia de datos espurios, comparado con otra metodolog&iacute;a est&aacute;ndar de identificaci&oacute;n.</p>      <p><b>Palabras clave: </b>datos espurios, filtro de Kalman, generador e&oacute;lico, identificaci&oacute;n de sistemas. </p>  <hr>      <p><b>Abstract</b></p>      <p>Among the technologies for clean energy generation, wind turbines are one of the most advanced. Due to stochastic behavior of the wind, a proper control action must be exerted over these turbines in order to get a better use of the output power. A propel control action usually implies the knowledge of the system dynamics, and that is why the identification stage is also crucial.</p>     <p>In this paper, the authors propose a robust statistics dual Kalman filter for identification of a wind generator. Although dual Kalman filter (DKF) has been used in system identification, its performance heavily depends on the absence of outliers in the measured data. However, outliers might be easily produced during the data acquisition stage. In this study, the authors show how by combining the DKF and the robust statistic Kalman filter, the problem of outliers in systems identification can be avoided. The method was implemented in Matlab R2009a&reg;. Results for the identification of a wind generator working online are shown and compared to the performance of dual Kalman filter.</p>     ]]></body>
<body><![CDATA[<p><b>Keywords: </b>outliers, Kalman filter, wind generator, system identification.</p> <hr>      <p><b>1. INTRODUCCI&Oacute;N</b></p>     <p>En la actualidad, la energ&iacute;a el&eacute;ctrica se ha convertido en un componente importante en el desarrollo de cualquier pa&iacute;s. La electricidad con base en combustibles ha pasado de ser una ayuda energ&eacute;tica a resolver problemas sociales y ambientales, por lo que se sugiere la penetraci&oacute;n de las energ&iacute;as limpias, como la energ&iacute;a e&oacute;lica. Para aprovechar la energ&iacute;a por medio del viento se han desarrollado diferentes tipos de turbinas, que incluyen los generadores de inducci&oacute;n y los generadores s&iacute;ncronos. Una de las complicaciones dentro de estas tecnolog&iacute;as consiste en la amplia variabilidad en la velocidad del viento &#91;1&#93;, y es por ello que se deben tener estrategias de control que permitan un mayor aprovechamiento de la energ&iacute;a proveniente del viento. Con el fin de implementar un esquema de control &oacute;ptimo es necesario conocer la din&aacute;mica del sistema, por lo cual se hace necesario recurrir a t&eacute;cnicas de identificaci&oacute;n de sistemas.</p>     <p>Para realizar la identificaci&oacute;n y control de este tipo de m&aacute;quinas se deben obtener medidas de velocidad angular, voltajes y corrientes en terminales de la red. Sin embargo, dado que estas medidas se adquieren por medio de sistemas sensoriales, es posible que contengan informaci&oacute;n err&oacute;nea o datos an&oacute;malos, por lo que se tendr&iacute;a una operaci&oacute;n posterior defectuosa. Se debe entonces contar con una identificaci&oacute;n que sea capaz de mitigar los efectos asociados a los datos an&oacute;malos o espurios.</p>     <p>Por otro lado, una t&eacute;cnica eficiente para realizar la estimaci&oacute;n de estados en sistemas din&aacute;micos lineales es el filtro de Kalman &#91;2&#93;. Sus ra&iacute;ces se encuentran en modelos de espacio de estados perturbados por ruidos gaussianos, proporcionando soluciones recursivas al problema de estimaci&oacute;n &oacute;ptima de estas variables; uno de sus inconvenientes, sin embargo, est&aacute; relacionado con su sensibilidad a datos espurios entregados por los sistemas de medida, haciendo que su rendimiento se vea comprometido.</p>     <p>En la literatura se han presentado diferentes versiones del filtro de Kalman con el fin de volverlo robusto a datos espurios, incluyendo: el filtro de Kalman umbralizado, propuesto en &#91;3&#93;, el cual utiliza un an&aacute;lisis comparativo para determinar si hay o no datos espurios. No obstante, el m&eacute;todo entrega valores estimados incorrectos cuando la salida estimada se aleja mucho del valor observado. Otra versi&oacute;n es el filtro de Kalman con ponderacio-nes Bayesianas, propuesto en &#91;4&#93;, donde se usa un enfoque similar al de m&iacute;nimos cuadrados ponderados, pero estas ponderaciones son modeladas probabil&iacute;sticamente; sin embargo este filtro entrega resultados inadecuados cuando se tienen se&ntilde;ales sin tendencia. Otra versi&oacute;n es el filtro de Kalman robusto a datos espurios, propuesto en &#91;5&#93;, en el que se supone que el ruido de medida es una variable oculta que se modela probabil&iacute;sticamente.</p>      <p>Adem&aacute;s de realizar la estimaci&oacute;n de estados, estas metodolog&iacute;as presentan estructuras recursivas que permiten llevar a cabo la estimaci&oacute;n robusta de los par&aacute;metros de sistemas din&aacute;micos. No obstante, estos m&eacute;todos asumen par&aacute;metros constantes con el tiempo, introduciendo problemas de estimaci&oacute;n cuando las se&ntilde;ales presentan cambios fuertes de nivel; igualmente, realizan la estimaci&oacute;n de par&aacute;metros asumiendo no tener una relaci&oacute;n entre las entradas controladas y las salidas observadas. Adicionalmente, los m&eacute;todos anteriores no han sido utilizados para la identificaci&oacute;n de sistemas multivariables.</p>     <p>Este estudio propone la combinaci&oacute;n de tres filtros secuenciales, uno de ellos corresponde al filtro de Kalman de estad&iacute;stica robusta (rsKF) propuesto en &#91;6&#93;, los otros dos filtros componen el filtro de Kalman dual (DKF) &#91;2&#93;. El DKF son dos filtros de Kalman que operan secuencialmente en la estimaci&oacute;n de los estados y las matrices de un modelo de espacio de estado de un sistema din&aacute;mico lineal variante en el tiempo. Debido a que los par&aacute;metros de las m&aacute;quinas el&eacute;ctricas cambian con su operaci&oacute;n, el DKF utiliza cada dato individual tan pronto como est&eacute; disponible para actualizar tanto la estimaci&oacute;n del modelo como los estados del sistema. Esta caracter&iacute;stica hace que estos algoritmos sean atractivos para aplicaciones <img src="img/revistas/inde/v32n1/v32n1a08f02.jpg">&#91;2&#93;. Por otra parte, el rsKF es una variaci&oacute;n del filtro de Kalman, que desarrolla una ponderaci&oacute;n para cada dato de entrada basada en un enfoque de m&aacute;xima verosimilitud argumentada en la teor&iacute;a de estad&iacute;stica robusta &#91;7&#93;.</p>     <p>Empleando las t&eacute;cnicas anteriormente mencionadas se obtuvieron resultados experimentales sobre un ejemplo de prueba que se realiz&oacute; a manera de simulaci&oacute;n: los datos sint&eacute;ticos fueron generados a partir de un sistema conocido con el fin de comprobar la eficiencia del rsKF; por otra parte, se aplic&oacute; la combinaci&oacute;n propuesta (DKF y rsKF) para la identificaci&oacute;n en tiempo real de un generador e&oacute;lico, en donde se utiliz&oacute; una tarjeta de adquisici&oacute;n NI-USB6009, fabricada por <img src="img/revistas/inde/v32n1/v32n1a08f03.jpg">para obtener lecturas de la se&ntilde;al de velocidad y voltaje en terminales, en donde la identificaci&oacute;n de sistemas realizada por la metodolog&iacute;a propuesta no se ve afectada por los datos espurios entregados por el sistema de adquisici&oacute;n.</p>     <p>El documento se encuentra organizado de la siguiente manera. La secci&oacute;n 2 muestra el modelado del generador s&iacute;ncrono mediante una representaci&oacute;n en modelo de espacio de estados. La secci&oacute;n 3 presenta la teor&iacute;a b&aacute;sica de los filtros de Kalman utilizados. La secci&oacute;n 4 muestra la metodolog&iacute;a empleada para realizar la identificaci&oacute;n del generador e&oacute;lico. La secci&oacute;n 5 expone los resultados y la discusi&oacute;n de la utilizaci&oacute;n de los filtros mencionados para datos sint&eacute;ticos y la identificaci&oacute;n de un generador e&oacute;lico. Finalmente, se presentan las conclusiones en la secci&oacute;n 6.</p>     ]]></body>
<body><![CDATA[<p><b>2. GENERADOR S&Iacute;NCRONO</b></p>     <p>Dentro de los elementos m&aacute;s importantes de un sistema de potencia est&aacute; el generador s&iacute;ncrono, ya que con esta m&aacute;quina se puede controlar la estabilidad del sistema de potencia, conociendo su modelo. Una de las tareas importantes al realizar el control de estas m&aacute;quinas es tener una estimaci&oacute;n apropiada de los par&aacute;metros de dicha m&aacute;quina. Debido a que sus par&aacute;metros cambian en el tiempo, se hace necesario tener una etapa de identificaci&oacute;n de sistemas. El generador e&oacute;lico puede ser representado por un generador s&iacute;ncrono descrito como &#91;8&#93;,</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f04.jpg"></p>      <p>donde <img src="img/revistas/inde/v32n1/v32n1a08f05.jpg"> es el voltaje en eje directo del generador, <img src="img/revistas/inde/v32n1/v32n1a08f06.jpg"> es el voltaje en eje en cuadratura, <img src="img/revistas/inde/v32n1/v32n1a08f07.jpg"> es la corriente en eje directo, <img src="img/revistas/inde/v32n1/v32n1a08f08.jpg"> es la corriente en eje en cuadratura, <img src="img/revistas/inde/v32n1/v32n1a08f09.jpg">es la velocidad angular de la m&aacute;quina, <img src="img/revistas/inde/v32n1/v32n1a08f10.jpg"> es el torque mec&aacute;nico, <img src="img/revistas/inde/v32n1/v32n1a08f11.jpg"> es la posici&oacute;n angular, y <img src="img/revistas/inde/v32n1/v32n1a08f12.jpg"> es una constante rot&oacute;rica que depende de la constante del rotor y la inductancia mutua de la m&aacute;quina.</p>     <p><b>3. FILTRO DE KALMAN</b></p>     <p>El filtro de Kalman es un m&eacute;todo que nace de la representaci&oacute;n de sistemas din&aacute;micos lineales en modelo de espacio de estados, donde interact&uacute;an diferentes variables como: salidas, estados, entradas controladas y las entradas no controladas.</p>     <p>La representaci&oacute;n de un sistema din&aacute;mico en espacio de estados est&aacute; dada por,</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f13.jpg"></p>      <p>donde <img src="img/revistas/inde/v32n1/v32n1a08f14.jpg"> es la matriz de transici&oacute;n de estados y en ella est&aacute; consignada toda la din&aacute;mica del sistema, <img src="img/revistas/inde/v32n1/v32n1a08f15.jpg"> es la matriz de relaci&oacute;n de entradas y estados, <img src="img/revistas/inde/v32n1/v32n1a08f16.jpg"> es la matriz de salida, <img src="img/revistas/inde/v32n1/v32n1a08f17.jpg">es la matriz que relaciona las entradas con las salidas. Se asume que <img src="img/revistas/inde/v32n1/v32n1a08f18.jpg"> y <img src="img/revistas/inde/v32n1/v32n1a08f19.jpg"> son ruidos gaussianos de media cero que no est&aacute;n correlacionados, es decir, <img src="img/revistas/inde/v32n1/v32n1a08f20.jpg"> donde <img src="img/revistas/inde/v32n1/v32n1a08f21.jpg"> y <img src="img/revistas/inde/v32n1/v32n1a08f22.jpg">. <img src="img/revistas/inde/v32n1/v32n1a08f23.jpg"> son matrices de covarianza diagonales para el ruido de estado y observaci&oacute;n, respectivamente. De aqu&iacute; en adelante se denominar&aacute;n los par&aacute;metros del modelo como <img src="img/revistas/inde/v32n1/v32n1a08f24.jpg"></p>     <p>Dentro de los sistemas din&aacute;micos hay variables a las cuales se puede tener acceso, como las variables de salida <img src="img/revistas/inde/v32n1/v32n1a08f25.jpg"> y la se&ntilde;al de control <img src="img/revistas/inde/v32n1/v32n1a08f26.jpg"> ; los ruidos siempre van a existir en los sistemas pr&aacute;cticos reales como variables a las que no se tiene acceso. Otras variables de estos sistemas son los estados <img src="img/revistas/inde/v32n1/v32n1a08f27.jpg"> y par&aacute;metros <img src="img/revistas/inde/v32n1/v32n1a08f28.jpg"> que no pueden ser observados. Pese a que los estados no se observan, son las variables m&aacute;s importantes del sistema din&aacute;mico, debido a que pueden ser utilizadas para representar la evoluci&oacute;n din&aacute;mica del sistema. Uno de los m&eacute;todos con mayor &eacute;xito para estimar los estados en sistemas din&aacute;micos lineales es el filtro de Kalman. El objetivo de esta t&eacute;cnica, desde un punto de vista probabil&iacute;stico, consiste en obtener el mejor estado estimado dado los par&aacute;metros  <img src="img/revistas/inde/v32n1/v32n1a08f11.jpg"> (asumiendo que  <img src="img/revistas/inde/v32n1/v32n1a08f29.jpg">) y las observaciones  <img src="img/revistas/inde/v32n1/v32n1a08f30.jpg">&#91;9&#93;. El mejor estimador  <img src="img/revistas/inde/v32n1/v32n1a08f31.jpg"> para en el instante  <img src="img/revistas/inde/v32n1/v32n1a08f32.jpg"> se puede obtener maximizando la probabilidad conjunta de los estados  <img src="img/revistas/inde/v32n1/v32n1a08f33.jpg"> y las observaciones  <img src="img/revistas/inde/v32n1/v32n1a08f34.jpg"> en el instante  <img src="img/revistas/inde/v32n1/v32n1a08f32.jpg">, dado los par&aacute;metros  <img src="img/revistas/inde/v32n1/v32n1a08f11.jpg"> y las observaciones hasta el instante  <img src="img/revistas/inde/v32n1/v32n1a08f35.jpg">. Del mismo modo, el mejor estimador para  <img src="img/revistas/inde/v32n1/v32n1a08f33.jpg"> se puede obtener minimizando el logaritmo negativo de la probabilidad conjunta de  <img src="img/revistas/inde/v32n1/v32n1a08f33.jpg"> y  <img src="img/revistas/inde/v32n1/v32n1a08f36.jpg">, lo cual conduce al siguiente problema de minimizaci&oacute;n,</p>      ]]></body>
<body><![CDATA[<p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f37.jpg"></p>      <p>donde <img src="img/revistas/inde/v32n1/v32n1a08f38.jpg"> es el n&uacute;mero de estados; <img src="img/revistas/inde/v32n1/v32n1a08f78.jpg">es el n&uacute;mero de salidas; <img src="img/revistas/inde/v32n1/v32n1a08f39.jpg">  son la media y la matriz de covarianza de <img src="img/revistas/inde/v32n1/v32n1a08f31.jpg">  dadas las observaciones <img src="img/revistas/inde/v32n1/v32n1a08f30.jpg"> ; los <img src="img/revistas/inde/v32n1/v32n1a08f40.jpg">  son los elementos del vector <img src="img/revistas/inde/v32n1/v32n1a08f41.jpg"> ; son elementos del vector <img src="img/revistas/inde/v32n1/v32n1a08f42.jpg"> es la <img src="img/revistas/inde/v32n1/v32n1a08f43.jpg">  componente del vector <img src="img/revistas/inde/v32n1/v32n1a08f44.jpg">  y <img src="img/revistas/inde/v32n1/v32n1a08f45.jpg"> es la i-esima fila de la matriz <img src="img/revistas/inde/v32n1/v32n1a08f46.jpg"> .</p>     <p>Resolviendo el problema de minimizaci&oacute;n se obtiene &#91;9&#93;. </p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f47.jpg"></p>      <p> donde, </p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f48.jpg"></p>      <p>La matriz <img src="img/revistas/inde/v32n1/v32n1a08f49.jpg"> en la ecuaci&oacute;n 8 se conoce como la ganancia de Kalman &#91;2&#93;.</p>     <p>El filtro de Kalman es un algoritmo eficiente que utiliza las observaciones ruidosas disponibles para estimar los estados de un sistema din&aacute;mico lineal. El filtro cuenta con dos etapas esenciales: una de predicci&oacute;n y una de actualizaci&oacute;n. En la etapa de predicci&oacute;n se calcula una estimaci&oacute;n a priori de los estados <img src="img/revistas/inde/v32n1/v32n1a08f50.jpg"> y de la matriz de covarianza <img src="img/revistas/inde/v32n1/v32n1a08f51.jpg"> Luego, con estas dos cantidades se realiza la etapa de actualizaci&oacute;n, calculando inicialmente la ganancia de Kalman <img src="img/revistas/inde/v32n1/v32n1a08f52.jpg"> y seguidamente actualizando los estados por medio de la ecuaci&oacute;n 7. La matriz de covarianza de <img src="img/revistas/inde/v32n1/v32n1a08f53.jpg"> tambi&eacute;n debe actualizarse (para m&aacute;s detalles ver &#91;2&#93;).</p>     <p><b>Filtro de Kalman dual </b>(DKF)</p>     <p>En el procedimiento de estimaci&oacute;n de estados del filtro de Kalman, se asumi&oacute; que los par&aacute;metros <img src="img/revistas/inde/v32n1/v32n1a08f54.jpg"> eran conocidos para todos los valores de <img src="img/revistas/inde/v32n1/v32n1a08f55.jpg">Debido a que los par&aacute;metros de las m&aacute;quinas el&eacute;ctricas son de dif&iacute;cil acceso y cambian con su operaci&oacute;n, es necesario contar con un m&eacute;todo que sea capaz de estimar estos par&aacute;metros cambiantes. Entre los m&eacute;todos que han sido empleados para la estimaci&oacute;n de los par&aacute;metros de m&aacute;quinas el&eacute;ctricas se encuentran el m&eacute;todo de los m&iacute;nimos cuadrados &#91;10&#93;, el m&eacute;todo de los m&iacute;nimos cuadrados recursivos &#91;11&#93;, enjambres de part&iacute;culas &#91;12&#93;, filtro de Kalman &#91;13&#93;, redes neuronales artificiales &#91;14&#93;, observadores en modos deslizantes &#91;15&#93;, entre otros. Algunos de los inconvenientes que presentan los anteriores m&eacute;todos es que demandan altos costos computa-cionales o no se pueden usar en aplicaciones en l&iacute;nea. Algunos presentan errores cuadr&aacute;ticos medios considerables en la estimaci&oacute;n de par&aacute;metros, otros &uacute;nicamente han sido probados en ambientes controlados, y algunos m&aacute;s no estiman todos los par&aacute;metros de la m&aacute;quina. Por otra parte, el filtro de Kalman dual (DKF) es una herramienta &uacute;til y apropiada que, por su enfoque recursivo de estimaci&oacute;n simult&aacute;nea de estados y par&aacute;metros de sistemas din&aacute;micos, sobresale cuando se compara contra muchas otras metodolog&iacute;as importantes que hacen parte de la identificaci&oacute;n de sistemas. Por consiguiente, a continuaci&oacute;n se introduce el DKF para realizar la estimaci&oacute;n dual simult&aacute;nea de los estados del sistema din&aacute;mico y las matrices <img src="img/revistas/inde/v32n1/v32n1a08f54.jpg"></p>       ]]></body>
<body><![CDATA[<p>Con el fin de modelar la ecuaci&oacute;n 5 por medio de una estructura matricial, se definir&aacute;n los siguientes t&eacute;rminos,</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f57.jpg"></p>      <p>Por lo tanto, el sistema mostrado en la ecuaci&oacute;n 5 puede ser descrito como</p>     <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f58.jpg"></p>       <p>Para llevar a cabo la estimaci&oacute;n de los par&aacute;metros <img src="img/revistas/inde/v32n1/v32n1a08f59.jpg"> se asume que los estados son conocidos. Se plantea para ello un m&eacute;todo de estimaci&oacute;n dual que trabaja usando el modelo para estimar la se&ntilde;al y luego se emplea la se&ntilde;al para estimar el modelo. Este proceso se puede realizar de forma se-cuencial tan pronto como est&eacute; disponible cada muestra individual, lo cual hace que estos algoritmos sean atractivos para aplicaciones en l&iacute;nea &#91;2&#93;. Por consiguiente, el filtro de Kalman puede ser usado para estimar par&aacute;metros y estados en una nueva representaci&oacute;n de espacio de estados dada por </p>       <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f60.jpg"></p>     <p>donde el modelo descrito por las Ecs. (11) y (12) corresponde a la din&aacute;mica de los par&aacute;metros y los estados, dentro de un mismo enfoque; <img src="img/revistas/inde/v32n1/v32n1a08f61.jpg"> es el ruido presente en la evoluci&oacute;n de los par&aacute;metros.</p>     <p>De manera similar se puede plantear una funci&oacute;n de costo para el modelo de las ecuaciones 11 y 12; y se puede demostrar que la minimizaci&oacute;n de esta funci&oacute;n de costo conduce a un problema de minimizaci&oacute;n similar al planteado en la ecuaci&oacute;n 6, para m&aacute;s detalle ver &#91;2&#93;.</p>     <p><b>Filtro de Kalman de estad&iacute;stica robusta</b></p>     <p>Con el objetivo de obtener un filtro de Kalman robusto, el filtro de Kalman de estad&iacute;stica robusta (rsKF) introduce una matriz de ponderaci&oacute;n <img src="img/revistas/inde/v32n1/v32n1a08f63.jpg">, para cada observaci&oacute;n, tal que la matriz de covarianza de las observaciones <img src="img/revistas/inde/v32n1/v32n1a08f64.jpg"> sea ponderada mediante <img src="img/revistas/inde/v32n1/v32n1a08f63.jpg">. El c&aacute;lculo de esta matriz de ponderaci&oacute;n se basa en la teor&iacute;a de estimaci&oacute;n por m&aacute;xima verosimilitud para estad&iacute;stica robusta &#91;16&#93;. Un estimador robusto para <img src="img/revistas/inde/v32n1/v32n1a08f65.jpg"> se obtiene minimizando &#91;6&#93;,</p>      ]]></body>
<body><![CDATA[<p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f66.jpg"></p>       <p>donde los t&eacute;rminos <img src="img/revistas/inde/v32n1/v32n1a08f67.jpg"> y fueron definidos en la ecuaci&oacute;n 6; <img src="img/revistas/inde/v32n1/v32n1a08f68.jpg"> es una funci&oacute;n de p&eacute;rdida dada como &#91;16&#93;,</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f69.jpg"></p>       <p>La funci&oacute;n de p&eacute;rdida en la ecuaci&oacute;n 6, es igual a la funci&oacute;n identidad. La constante <img src="img/revistas/inde/v32n1/v32n1a08f70.jpg">es seleccionada seg&uacute;n el grado de penalizaci&oacute;n en la funci&oacute;n de p&eacute;rdida <img src="img/revistas/inde/v32n1/v32n1a08f68.jpg">.</p>     <p>Se puede demostrar &#91;6&#93; que la ganancia de Kalman presentada en la ecuaci&oacute;n 8 puede ser calculada como</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f71.jpg"></p>       <p>donde <img src="img/revistas/inde/v32n1/v32n1a08f72.jpg">. Las ponderaciones <img src="img/revistas/inde/v32n1/v32n1a08f73.jpg"> en el instante <img src="img/revistas/inde/v32n1/v32n1a08f55.jpg"> se calculan como</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f74.jpg"></p>       <p>con <img src="img/revistas/inde/v32n1/v32n1a08f75.jpg"> igual a la funci&oacute;n de p&eacute;rdida de Huber &#91;16&#93;,</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f76.jpg"></p>  evistas/inde/v32n1/v32n1a08f60.jpg"></p>       ]]></body>
<body><![CDATA[<p>donde <img src="img/revistas/inde/v32n1/v32n1a08f77.jpg">  es la funci&oacute;n signo. Para los experimentos realizados en este estudio el valor de <img src="img/revistas/inde/v32n1/v32n1a08f79.jpg"> fue seleccionado como <img src="img/revistas/inde/v32n1/v32n1a08f80.jpg"> para una contaminaci&oacute;n del 5% de los datos &#91;6&#93;.</p>     <p>La aplicaci&oacute;n del rsKF se hace de manera similar que el filtro de Kalman est&aacute;ndar (sKF), lo &uacute;nico en que se diferencia es en el c&aacute;lculo de la ganancia de Kalman.</p>     <p><b>4. METODOLOG&Iacute;A</b></p>     <p>El desarrollo de este estudio cont&oacute; con la utilizaci&oacute;n de dos tipos de datos: sint&eacute;ticos y datos a partir de un generador e&oacute;lico real implementado por medio de un generador s&iacute;ncrono. A continuaci&oacute;n se describir&aacute;n de forma detallada los datos, c&oacute;mo se realiz&oacute; la identificaci&oacute;n mediante el enfoque propuesto y, finalmente, se explican las m&eacute;tricas de validaci&oacute;n para observar la influencia de los datos espurios en las observaciones.</p>     <p><b>Base de datos sint&eacute;ticos a partir de un modelo conocido</b></p>     <p>A partir de un modelo conocido (funci&oacute;n coseno) se crearon 1.000 observaciones. A la se&ntilde;al de salida (q = 1) se le agreg&oacute; ruido blanco gaussiano con varianza <img src="img/revistas/inde/v32n1/v32n1a08f81.jpg"> . Los datos espurios se introdujeron en una cantidad igual al 5% de todos los datos. Los datos espurios se ubicaron aleatoriamente como parte de las primeras 500 observaciones. Para el sistema conocido se asumi&oacute; un &uacute;nico estado (n = 1). Adem&aacute;s, se asumi&oacute; <img src="img/revistas/inde/v32n1/v32n1a08f82.jpg">&#91;4&#93; . Igualmente se asumi&oacute; que las matrices B =D = 0.</p>     <p><b>Base de datos de un generador e&oacute;lico</b></p>     <p>En la <A href="#f1">figura 1</A> se muestra el montaje que se utiliz&oacute; para llevar a cabo la identificaci&oacute;n del generador e&oacute;lico. Se emple&oacute; un variador de velocidad (VV) (ver ap&eacute;ndice C) que alimenta un motor de inducci&oacute;n para simular el elemento encargado de impulsar el generador e&oacute;lico a una velocidad espec&iacute;fica. En este contexto, el motor de inducci&oacute;n es llamado primo motor. Igualmente, se utiliz&oacute; un m&oacute;dulo de estado s&oacute;lido (SSM) (ver ap&eacute;ndice D) para alimentar el campo del generador. Se asumi&oacute; que el grupo completo es de orden ocho, donde las entradas corresponden al voltaje en el VV y al voltaje en el SSM, donde se hicieron variaciones en sus valores para tener condiciones de velocidad y voltaje variable. Las salidas corresponden a la velocidad del grupo y al voltaje en terminales del generador. Para adquirir los datos a partir del generador e&oacute;lico se utiliz&oacute; una tarjeta NI-USB6009 de  <img src="img/revistas/inde/v32n1/v32n1a08f83.jpg">(ver ap&eacute;ndice E) con un tiempo de muestreo ts = 0,1s. A cada salida se le agreg&oacute; un ruido blanco gaussiano, con varianza de <img src="img/revistas/inde/v32n1/v32n1a08f84.jpg"> para la velocidad; y para la tensi&oacute;n en terminales una varianza de <img src="img/revistas/inde/v32n1/v32n1a08f85.jpg"></p>      <p align="center"><a name="f1"></a><img src="img/revistas/inde/v32n1/v32n1a08f86.jpg"></p>      <p><b>Procedimiento general para la identificaci&oacute;n robusta del generador e&oacute;lico</b></p>     ]]></body>
<body><![CDATA[<p>La secci&oacute;n 2 expuso de forma b&aacute;sica la estructura matem&aacute;tica de las t&eacute;cnicas que se emplearon en este estudio. En esta subsecci&oacute;n se hablar&aacute; de c&oacute;mo se combinan el DKF y el rsKF. El modelo que combina el DKF y el rsKF se denominar&aacute; DKF-rsKF. Para explicar c&oacute;mo se emplearon estos dos m&eacute;todos en la identificaci&oacute;n de sistemas se presenta un diagrama de bloques en la <a href="#f2">figura 2</a>.</p>      <p align="center"><a name="f2"></a><img src="img/revistas/inde/v32n1/v32n1a08f87.jpg"></p>      <p>Como se ha mencionado, las dos variables que se pueden medir en un sistema din&aacute;mico son la se&ntilde;al de entrada o de control u y la se&ntilde;al de salida y. La se&ntilde;al de salida se pasa por el rsKF para obtener una se&ntilde;al que no presente alteraciones frente a datos espurios; esta se&ntilde;al filtrada se convierte en las observaciones disponibles que, junto con la se&ntilde;al de control del sistema, se usa como entrada al DKF que estima las matrices <img src="img/revistas/inde/v32n1/v32n1a08f88.jpg">, y los estados <img src="img/revistas/inde/v32n1/v32n1a08f89.jpg"> . Luego de obtener los par&aacute;metros y estados estimados se calcula la se&ntilde;al de salida estimada.</p>     <p><b>Validaci&oacute;n</b></p>     <p>Con la finalidad de ilustrar la influencia de los datos espurios presentes sobre la predicci&oacute;n realizada por los m&eacute;todos (sKF y rsKF) se calcularon dos estad&iacute;sticas mu&eacute;strales: la skewness muestral <img src="img/revistas/inde/v32n1/v32n1a08f90.jpg"> y la kurtosis muestral <img src="img/revistas/inde/v32n1/v32n1a08f91.jpg">definidas como &#91;17&#93;</p>      <p align="center"><img src="img/revistas/inde/v32n1/v32n1a08f92.jpg"></p>      <p>donde <img src="img/revistas/inde/v32n1/v32n1a08f93.jpg"> es la media muestral. El valor de -3 que aparece en la ecuaci&oacute;n 18 hace que su valor sea id&eacute;nticamente cero para el caso en el que los datos sigan una funci&oacute;n de distribuci&oacute;n gaussiana. Otra m&eacute;trica usada para observar el rendimiento de los m&eacute;todos es el error residual <img src="img/revistas/inde/v32n1/v32n1a08f94.jpg"> entre la predicci&oacute;n <img src="img/revistas/inde/v32n1/v32n1a08f95.jpg"> y su valor real sin datos espurios <img src="img/revistas/inde/v32n1/v32n1a08f96.jpg">(datos limpios) por medio de <img src="img/revistas/inde/v32n1/v32n1a08f97.jpg"> . La kurtosis, la skewness y el error residual solo se calcularon para los datos sint&eacute;ticos.</p>     <p>La influencia de los datos espurios en las observaciones puede ser demostrada a partir de la skewness y kurtosis, puesto que ambas cantidades son sensibles a datos anormales. Ambos momentos estad&iacute;sticos han constituido la base para proponer m&eacute;todos que detectan datos an&oacute;malos &#91;17&#93;. Lo quese pretende es comparar el rendimiento del sKF y del rsKF usando estas m&eacute;tricas.</p>     <p><b>5. RESULTADOS Y DISCUSI&Oacute;N</b></p>     <p>Primero se presenta un experimento con datos sint&eacute;ticos, en el que se muestra el desempe&ntilde;o de un filtro de Kalman de estad&iacute;stica robusta comparado con un filtro de Kalman est&aacute;ndar.</p>     ]]></body>
<body><![CDATA[<p>En un segundo experimento se emple&oacute; el sistema DKF-rsKF, que se mostr&oacute; en la <a href="#f2">figura 2</a>, con el prop&oacute;sito de realizar la identificaci&oacute;n de un generador e&oacute;lico.</p>     <p><b>Datos sint&eacute;ticos a partir de un sistema conocido</b></p>     <p>En la <a href="#f3">figura 3</a> se muestran los datos sint&eacute;ticos, a partir de una funci&oacute;n coseno ruidosa, que se emplearon en el experimento.</p>      <p align="center"><a name="f3"></a><img src="img/revistas/inde/v32n1/v32n1a08f98.jpg"></p>      <p>La <a href="#f4">figura 4</a> presenta la respuesta del sKF y rsKF ante el problema de tratamiento de datos espurios en las observaciones.</p>      <p align="center"><a name="f4"></a><img src="img/revistas/inde/v32n1/v32n1a08f99.jpg"></p>      <p>La <a href="#f4">figura 4</a> muestra c&oacute;mo el sKF es sensible a estos datos. Por otra parte, el rsKF pondera fuertemente la varianza de las observaciones cuando se tiene un error de estimaci&oacute;n grande y, de esta manera, las predicciones del rsKF no se ven perturbadas por la presencia de datos espurios. Tambi&eacute;n en la <a href="#f4">figura 4</a> se puede observar que en el caso de no tener datos espurios, es decir en las &uacute;ltimas 500 muestras, ambos filtros entregan estimaciones adecuadas sin importar el comportamiento de las observaciones. El rsKF obtuvo un error residual de 0,0087 &plusmn; 0,0023, mientras el sKF un error de 0,1728 &plusmn; 0,1945.</p>     <p>Como se mencion&oacute; anteriormente, la kurtosis y la skewness son susceptibles a datos espurios. Por ejemplo, un valor alto de kurtosis o un valor de skewness diferente de cero implican la existencia de datos producidos con distribuciones de colas pesadas en las observaciones. Usualmente, estas distribuciones de colas pesadas son las que mejor describen la incertidum-bre en fuentes generadoras de datos espurios &#91;18&#93;. Para este experimento se obtuvo una kurtosis para los datos limpios de -0,7765&plusmn;0,5423, para el rsKF de -0,8586 &plusmn; 0,5318 y finalmente para el sKF de 5,5647 &plusmn; 12,0013. Se observa entonces que el sKF se ve afectado por los datos espurios presentesen las observaciones, mientras el rsKF ajusta sus predicciones a los valores similares de la se&ntilde;al sin datos an&oacute;malos.</p>     <p>Similarmente, con respecto a la skewness, para los datos limpios se tiene un valor de -0,007 &plusmn; 0,46; para el rsKF de -0,0048 &plusmn; 0,5722, y para el sKF se obtiene una skewness de 0,0199 &plusmn; 1,5275.</p>     <p><b>Resultados en tiempo real para un generador e&oacute;lico</b></p>     ]]></body>
<body><![CDATA[<p>Los resultados anteriores se basaron en un modelo con datos sint&eacute;ticos, obteniendo propiedades relevantes para el caso de tratamiento de datos espurios. Posteriormente, con el prop&oacute;sito de observar c&oacute;mo es el comportamiento del algoritmo de identificaci&oacute;n propuesto en un ambiente real, se realiz&oacute; la identificaci&oacute;n de un generador e&oacute;lico. Su descripci&oacute;n se hizo en la secci&oacute;n 4.</p>     <p>Como primera prueba se analiz&oacute; la respuesta del DKF ante cambios aleatorios de las se&ntilde;ales de entrada. La <a href="#f5">figura 5</a> muestra el comportamiento del m&eacute;todo de identificaci&oacute;n DKF para el generador e&oacute;lico, en donde se observa que el algoritmo se ve afectado por el ruido presente en la se&ntilde;al y los datos espurios entregados por el m&oacute;dulo de adquisici&oacute;n, como se nota en las muestras 410 y 756 para la velocidad (<a href="#f5">figura 5</a> superior) y para las muestras 437, 465 y 862 del voltaje en terminales de la m&aacute;quina s&iacute;ncrona (<a href="#f5">figura 5</a> inferior). En las primeras muestras de esta figura se observa que las se&ntilde;ales tienen peque&ntilde;as diferencias; a medida que se presentan m&aacute;s muestras, la se&ntilde;al estimada por el DKF se hace muy similar a las observaciones entregadas por la tarjeta de adquisici&oacute;n; y al final se presenta un traslape entre la se&ntilde;al medida y la estimada por el m&eacute;todo. Para mostrar lo anterior se realiz&oacute; un acercamiento de la <a href="#f5">figura 5</a> entre las muestras 750 y 760 para la velocidad del generador y entre las muestras 850 y 880 para el voltaje en terminales, con el fin de observar dicha sobreposici&oacute;n (ver <a href="#f6">figura 6</a>). De la <a href="#f6">figura 6</a> es claro que el DKF sigue los datos espurios presentes en las observaciones, lo cual es una propiedad negativa del proceso de identificaci&oacute;n.</p>      <p align="center"><a name="f5"></a><img src="img/revistas/inde/v32n1/v32n1a08f100.jpg"></p>     <p align="center"><a name="f6"></a><img src="img/revistas/inde/v32n1/v32n1a08f101.jpg"></p>     <p>Despu&eacute;s de haber mostrado el efecto que tienen los datos espurios en las observaciones usando el DKF para la identificaci&oacute;n, se procedi&oacute; a estudiar la identificaci&oacute;n del generador empleando el DKF-rsKF: la se&ntilde;al adquirida por el sistema sensorial es analizada inicialmente por el rsKF para obtener se&ntilde;ales no perturbadas por los datos espurios, y a continuaci&oacute;n se aplica el DKF para obtener los estados, salidas y matrices estimadas, como se muestra en la <a href="#f2">figura 2</a>. Igualmente que en el caso anterior, se efectuaron cambios aleatorios en la se&ntilde;al de control y se le agreg&oacute; ruido blanco gaussiano a cada salida. En la <a href="#f7">figura 7</a> se muestra el comportamiento del voltaje en terminales (<a href="#f7">figura 7</a> inferior) y la velocidad de la m&aacute;quina s&iacute;ncrona (<a href="#f7">figura 7 </a>superior). En estas gr&aacute;ficas se puede observar que los datos espurios y el ruido presente en la se&ntilde;al no afectan la estimaci&oacute;n de las se&ntilde;ales de salida. Adem&aacute;s, se ve que para las muestras 300, 694 y 989 en la velocidad del generador, el algoritmo aten&uacute;a el efecto de los datos espurios y entrega una estimaci&oacute;n apropiada del valor de la se&ntilde;al en esas muestras. De igual forma, para el voltaje en terminales del generador, en las muestras 612, 1.055 y 1.100, arrojadas por el equipo de adquisici&oacute;n, el algoritmo de identificaci&oacute;n multivariable estima de manera adecuada estos valores, brindando robustez en la estimaci&oacute;n de la se&ntilde;al de salida del sistema din&aacute;mico. Para observar con claridad la identificaci&oacute;n hecha con el DKF-rsKF se realiz&oacute; un acercamiento de la <a href="#f7">figura 7</a> entre las muestras 280 y 370 para la velocidad del generador y entre las muestras 570 y 660 para el voltaje en terminales, como se muestra en la <a href="#f8">figura 8</a>.</p>     <p align="center"><a name="f7"></a><img src="img/revistas/inde/v32n1/v32n1a08f102.jpg"></p>     <p align="center"><a name="f8"></a><img src="img/revistas/inde/v32n1/v32n1a08f103.jpg"></p>     <p>Para obtener una idea de cu&aacute;nto se tarda el m&eacute;todo propuesto en realizar el procedimiento de identificaci&oacute;n, se analiz&oacute; el tiempo de c&aacute;lculo promedio para cada iteraci&oacute;n; para ello se obtuvo el valor promedio de 1.500 iteraciones y su respectiva desviaci&oacute;n est&aacute;ndar. El DKF-rsKF present&oacute; un tiempo de c&aacute;lculo de 0,0120 &plusmn; 0,0039 segundos, mientras el DKF, un tiempo de c&aacute;lculo de 0,0109&plusmn;7,67e-04 segundos. El DKF-rsKF tarda 1,1 milisegundos m&aacute;s que el DKF para identificar el generador e&oacute;lico en presencia de datos espurios, lo cual no se aleja mucho de los tiempos de ejecuci&oacute;n aplicando metodolog&iacute;as est&aacute;ndar.</p>      <p><b>CONCLUSIONES</b></p>     <p>En este estudio se present&oacute; una metodolog&iacute;a para la identificaci&oacute;n de un generador e&oacute;lico implementado a trav&eacute;s de un generador s&iacute;ncrono. El sistema anterior corresponde a un sistema multivariable acoplado, el m&eacute;todo propuesto permiti&oacute; realizar la estimaci&oacute;n del modelo sin verse afectado por los datos espurios, lo que s&iacute; pasa con el filtro de Kalman dual est&aacute;ndar. Tambi&eacute;n se pudo demostrar que la identificaci&oacute;n del generador e&oacute;lico que exhibe comportamientos no lineales se realiza satisfactoriamente empleando una combinaci&oacute;n de m&eacute;todos lineales que incluyen el DKF y el rsKF.</p>     ]]></body>
<body><![CDATA[<p>En particular, se mostr&oacute; que el filtro de Kalman est&aacute;ndar se ve afectado por datos espurios, obteniendo m&eacute;tricas de validaci&oacute;n err&oacute;neas, como el caso de la skewness, la kurtosis y el error residual.</p>     <p>Igualmente, se demostr&oacute; que el rsKF es un m&eacute;todo capaz de atenuar el efecto de los datos espurios en las observaciones, alcanzando medidas para la kurtosis y skewness similares en comparaci&oacute;n con las se&ntilde;ales sin datos espurios. Tambi&eacute;n logr&oacute; un error residual de 0,87%, el cual es menor que el error obtenido por el sKF (17,82%).</p>      <p><b>Agradecimientos</b></p>     <p>Este trabajo se desarroll&oacute; en el marco del proyecto de investigaci&oacute;n: "Identificaci&oacute;n de par&aacute;metros de sistemas multivariables basado en el filtro de Kalman y filtro de part&iacute;culas dando robustez a outliers", en el programa de j&oacute;venes investigadores "Virginia Guti&eacute;rrez de Pineda" y fue financiado por medio de Colciencias y la Universidad Tecnol&oacute;gica de Pereira.</p>     <p><b>Ap&eacute;ndice A, motor de inducci&oacute;n</b></p>     <p>Modelo IM-100 de <img src="img/revistas/inde/v32n1/v32n1a08f104.jpg"> motor jaula de ardilla con estator bobinado, trif&aacute;sico, cuatro polos, 240V, 1.5A, 60Hz y 1725 RPM.</p>     <p><b>Ap&eacute;ndice B, generador s&iacute;ncrono</b></p>     <p>Modelo IM-100-3 de <img src="img/revistas/inde/v32n1/v32n1a08f104.jpg"> generador trif&aacute;sico, cuatro polos, 240V, 1.5A, 60Hz, 0.33HP y 1800 RPM.</p>     <p><b>Ap&eacute;ndice C, variador de velocidad-Varispeed-606PC3</b></p> </font>    <p><font size="2" face="verdana">Equipo que controla la velocidad de m&aacute;quinas, especialmente motores; el elaborado por <img src="img/revistas/inde/v32n1/v32n1a08f106.jpg"> proporciona un modo remoto de operaci&oacute;n. Sirve para suministrar un voltaje de salida con una frecuencia proporcional a un voltaje DC de referencia.</font></p> <font face="verdana" size="2">    ]]></body>
<body><![CDATA[<p><b>Ap&eacute;ndice D, m&oacute;dulo de estado s&oacute;lido </b>H-IEC-DC2</p>     <p>Es un sistema que realiza el control por medio de rectificadores controlados de silicio (SCR). Este m&oacute;dulo, fabricado por <img src="img/revistas/inde/v32n1/v32n1a08f104.jpg"> es un equipo para uso acad&eacute;mico. El m&oacute;dulo tambi&eacute;n permite entregar un voltaje variable entre 0-125V DC proporcional a un voltaje de referencia de 0-10V.</p>     <p><b>Ap&eacute;ndice E, m&oacute;dulo de adquisici&oacute;n</b></p>     <p>La tarjeta NI-USB6009 de <img src="img/revistas/inde/v32n1/v32n1a08f107.jpg"> brinda funcionalidad de adquisici&oacute;n de datos b&aacute;sica para aplicaciones como registro de datos simple y medidas port&aacute;tiles. Algunas de sus caracter&iacute;sticas incluyen: 8 entradas anal&oacute;gicas, 2 salidas anal&oacute;gicas, rango de voltaje m&aacute;ximo de entrada &#91;-10 10&#93;V y rango de voltaje m&aacute;ximo de salida &#91;0 5&#93;V.</p> <hr>     <p><b>REFERENCIAS</b></p>     <!-- ref --><p>&#91;I&#93; S. Thomsen, K. Rothenhagen, and F. Fuchs, "Online parameter identification methods for doubly fed induction generators", in Power Electronics Specialists Conference, 2008. PESC 2008. IEEE, junio, pp. 2735-2741.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000128&pid=S0122-3461201400010000800001&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;2&#93;    S. Haykin, Kalman Filtering and Neural Networks. John Wiley &amp; Sons, 2001.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000130&pid=S0122-3461201400010000800002&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;3&#93; I. C. Schick and S. K. Mitter, "Robust recursive estimation in the presence of heavy-tailed observation noise", pp. 1045-1080, 1994.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000132&pid=S0122-3461201400010000800003&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;4&#93; J.-A. Ting, E. Theodorou, and S. Schaal, "A Kalman filter for robust outlier detection", Intelligent Robots and Systems, 2007. IROS 2007. IEEE/RSJ International Conference on, octubre-noviembre 2007, pp. 1514-1519.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000134&pid=S0122-3461201400010000800004&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;5&#93; G. Agamennoni, J. Nieto, and E. Nebot, "An outlier-robust Kalman filter", in Robotics and Automation (ICRA), 2011 IEEE International Conference on, mayo 2011, pp. 1551-1558.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000136&pid=S0122-3461201400010000800005&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;6&#93; T. Cipra and R. Romera, "Kalman filter with outliers and missing observations", TEST: An Official Journal of the Spanish Society of Statistics and Operations Research, vol. 6, n.&deg; 2, pp. 379-395, 1997.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000138&pid=S0122-3461201400010000800006&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;7&#93;    P. J. Huber, "Robust estimation of a location parameter", The Annals of Mathematical Statistics, vol. 35, n.&deg; 1, pp. 73-101, 1964.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000140&pid=S0122-3461201400010000800007&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;8&#93; J. Chiasson, Modeling and High-Performance Control of Electric Machines. IEEE Press/Wiley-Interscience, 2005.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000142&pid=S0122-3461201400010000800008&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;9&#93; A. Nelson, "Nonlinear estimation and modeling of noisy time-series by dual Kalman filtering methods", Ph.D. dissertation, septiembre 2000.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000144&pid=S0122-3461201400010000800009&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;10&#93; A. Mansouri and F. Krim, "Genetic algorithms based parameters identification of induction machine armax model", in Power Engineering and Optimization Conference (PEOCO), 2011 5th International, junio 2011, pp. 182-189.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000146&pid=S0122-3461201400010000800010&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;II&#93; Y. He, Y. Wang, Y. Feng, and Z. Wang, "Parameter identification of an induction machine at standstill using the vector constructing method", Power Electronics, IEEE Transactions on, vol. 27, n.&deg; 2, pp. 905-915, febrero 2012. DOI: 10.1109/TPEL.2010.2089699.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000148&pid=S0122-3461201400010000800011&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;12&#93; P. Kou, J. Zhou, C. Wang, H. Xiao, H. Zhang, and C. Li, "Parameters identification of nonlinear state space model of synchronous generator", Engineering Applications of Artificial Intelligence, vol. 24, n.&deg; 7, pp. 1227-1237, 2011. DOI: 10.1016/j.engappai.2011.05.012.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000150&pid=S0122-3461201400010000800012&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;13&#93; M. Huang, W. Li, and W. Yan, "Estimating parameters of synchronous generators using square-root unscented Kalman filter", Electric Power System Research, vol. 80, n.&deg; 9, pp. 1137-1144, 2010.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000152&pid=S0122-3461201400010000800013&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;14&#93; H. Sediki, A. Bechouche, D. O. Abdeslam, and S. Haddad, "Adaline approach for induction motor mechanical parameters identification", Mathematics and Computers in Simulation, n.&deg; 0, 2012. DOI: 10.1016/j.matcom.2012.05.003.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000154&pid=S0122-3461201400010000800014&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;15&#93; A. Meghnous, M. Pham, and X. Lin-Shi, "Dynamic identification of a synchronous machine using an extended sliding mode observer", Mathematics and Computers in Simulation, n.&deg; 0, 2012. DOI: 10.1016/j.matcom.2012.05.020.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000156&pid=S0122-3461201400010000800015&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;16&#93;  P. Huber, Robust Statistics, ed. New York: Wiley Series in Probability and Statistics, 2009.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000158&pid=S0122-3461201400010000800016&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;17&#93; R. Pearson, "Exploring process data", <img src="img/revistas/inde/v32n1/v32n1a08f121.jpg"> vol. 11, n.&deg; 2, pp. 179-194, 2001. DOI: 10.1016/S0959-1524(00)00046-9.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000160&pid=S0122-3461201400010000800017&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>     <!-- ref --><p>&#91;18&#93; C. M. Bishop, Pattern Recognition and Machine Learning. Secaucus, NJ, USA: Springer Verlag New York, Inc., 2006.    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[&#160;<a href="javascript:void(0);" onclick="javascript: window.open('/scielo.php?script=sci_nlinks&ref=000162&pid=S0122-3461201400010000800018&lng=','','width=640,height=500,resizable=yes,scrollbars=1,menubar=yes,');">Links</a>&#160;]<!-- end-ref --></p>  </font>      ]]></body><back>
<ref-list>
<ref id="B1">
<label>I</label><nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Thomsen]]></surname>
<given-names><![CDATA[S]]></given-names>
</name>
<name>
<surname><![CDATA[Rothenhagen]]></surname>
<given-names><![CDATA[K]]></given-names>
</name>
<name>
<surname><![CDATA[Fuchs]]></surname>
<given-names><![CDATA[F]]></given-names>
</name>
</person-group>
<source><![CDATA[parameter identification methods for doubly fed induction generators"]]></source>
<year></year>
<conf-name><![CDATA[ Power Electronics Specialists Conference]]></conf-name>
<conf-date>2008</conf-date>
<conf-loc> </conf-loc>
</nlm-citation>
</ref>
<ref id="B2">
<label>2</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Haykin]]></surname>
<given-names><![CDATA[S]]></given-names>
</name>
</person-group>
<source><![CDATA[Kalman Filtering and Neural Networks]]></source>
<year>2001</year>
<publisher-name><![CDATA[John Wiley & Sons]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B3">
<label>3</label><nlm-citation citation-type="">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Schick]]></surname>
<given-names><![CDATA[I. C]]></given-names>
</name>
<name>
<surname><![CDATA[Mitter]]></surname>
<given-names><![CDATA[S. K]]></given-names>
</name>
</person-group>
<source><![CDATA["Robust recursive estimation in the presence of heavy-tailed observation noise"]]></source>
<year>1994</year>
<page-range>1045-1080</page-range></nlm-citation>
</ref>
<ref id="B4">
<label>4</label><nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Ting]]></surname>
<given-names><![CDATA[J.-A]]></given-names>
</name>
<name>
<surname><![CDATA[Theodorou]]></surname>
<given-names><![CDATA[E]]></given-names>
</name>
<name>
<surname><![CDATA[Schaal]]></surname>
<given-names><![CDATA[S]]></given-names>
</name>
</person-group>
<source><![CDATA["A Kalman filter for robust outlier detection"]]></source>
<year></year>
<conf-name><![CDATA[ Intelligent Robots and Systems, 2007. IROS 2007. IEEE/RSJ International Conference on]]></conf-name>
<conf-date>octubre-noviembre 2007</conf-date>
<conf-loc> </conf-loc>
<page-range>1514-1519</page-range></nlm-citation>
</ref>
<ref id="B5">
<label>5</label><nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Agamennoni]]></surname>
<given-names><![CDATA[G]]></given-names>
</name>
<name>
<surname><![CDATA[Nieto]]></surname>
<given-names><![CDATA[J]]></given-names>
</name>
<name>
<surname><![CDATA[Nebot]]></surname>
<given-names><![CDATA[E]]></given-names>
</name>
</person-group>
<source><![CDATA["An outlier-robust Kalman filter"]]></source>
<year></year>
<conf-name><![CDATA[ 2011 IEEE International Conference on]]></conf-name>
<conf-date>mayo 2011</conf-date>
<conf-loc> </conf-loc>
<page-range>1551-1558</page-range></nlm-citation>
</ref>
<ref id="B6">
<label>6</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Cipra]]></surname>
<given-names><![CDATA[T]]></given-names>
</name>
<name>
<surname><![CDATA[Romera]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Kalman filter with outliers and missing observations"]]></article-title>
<source><![CDATA[TEST: An Official Journal of the Spanish Society of Statistics and Operations Research]]></source>
<year>1997</year>
<volume>6</volume>
<numero>2</numero>
<issue>2</issue>
<page-range>379-395</page-range></nlm-citation>
</ref>
<ref id="B7">
<label>7</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Huber]]></surname>
<given-names><![CDATA[P. J]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Robust estimation of a location parameter"]]></article-title>
<source><![CDATA[The Annals of Mathematical Statistics]]></source>
<year>1964</year>
<volume>35</volume>
<numero>1</numero>
<issue>1</issue>
<page-range>73-101</page-range></nlm-citation>
</ref>
<ref id="B8">
<label>8</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Chiasson]]></surname>
<given-names><![CDATA[J]]></given-names>
</name>
</person-group>
<source><![CDATA[Modeling and High-Performance Control of Electric Machines]]></source>
<year>2005</year>
<publisher-name><![CDATA[IEEE Press/Wiley-Interscience]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B9">
<label>9</label><nlm-citation citation-type="">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Nelson]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
</person-group>
<source><![CDATA["Nonlinear estimation and modeling of noisy time-series by dual Kalman filtering methods"]]></source>
<year></year>
</nlm-citation>
</ref>
<ref id="B10">
<label>10</label><nlm-citation citation-type="confpro">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Mansouri]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
<name>
<surname><![CDATA[Krim]]></surname>
<given-names><![CDATA[F]]></given-names>
</name>
</person-group>
<source><![CDATA["Genetic algorithms based parameters identification of induction machine armax model"]]></source>
<year></year>
<conf-name><![CDATA[ 2011 5th International]]></conf-name>
<conf-date>junio 2011</conf-date>
<conf-loc> </conf-loc>
<page-range>182-189</page-range></nlm-citation>
</ref>
<ref id="B11">
<label>II</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[He]]></surname>
<given-names><![CDATA[Y]]></given-names>
</name>
<name>
<surname><![CDATA[Wang]]></surname>
<given-names><![CDATA[Y]]></given-names>
</name>
<name>
<surname><![CDATA[Feng]]></surname>
<given-names><![CDATA[Y]]></given-names>
</name>
<name>
<surname><![CDATA[Wang]]></surname>
<given-names><![CDATA[Z]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Parameter identification of an induction machine at standstill using the vector constructing method"]]></article-title>
<source><![CDATA[Power Electronics, IEEE Transactions on]]></source>
<year>febr</year>
<month>er</month>
<day>o </day>
<volume>27</volume>
<numero>2</numero>
<issue>2</issue>
<page-range>905-915</page-range></nlm-citation>
</ref>
<ref id="B12">
<label>12</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Kou]]></surname>
<given-names><![CDATA[P]]></given-names>
</name>
<name>
<surname><![CDATA[Zhou]]></surname>
<given-names><![CDATA[J]]></given-names>
</name>
<name>
<surname><![CDATA[Wang]]></surname>
<given-names><![CDATA[C]]></given-names>
</name>
<name>
<surname><![CDATA[Xiao]]></surname>
<given-names><![CDATA[H]]></given-names>
</name>
<name>
<surname><![CDATA[Zhang]]></surname>
<given-names><![CDATA[H]]></given-names>
</name>
<name>
<surname><![CDATA[Li]]></surname>
<given-names><![CDATA[C]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Parameters identification of nonlinear state space model of synchronous generator"]]></article-title>
<source><![CDATA[Engineering Applications of Artificial Intelligence]]></source>
<year>2011</year>
<volume>24</volume>
<numero>7</numero>
<issue>7</issue>
<page-range>1227-1237</page-range></nlm-citation>
</ref>
<ref id="B13">
<label>13</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Huang]]></surname>
<given-names><![CDATA[M]]></given-names>
</name>
<name>
<surname><![CDATA[Li]]></surname>
<given-names><![CDATA[W]]></given-names>
</name>
<name>
<surname><![CDATA[Yan]]></surname>
<given-names><![CDATA[W]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Estimating parameters of synchronous generators using square-root unscented Kalman filter"]]></article-title>
<source><![CDATA[Electric Power System Research]]></source>
<year>2010</year>
<volume>80</volume>
<numero>9</numero>
<issue>9</issue>
<page-range>1137-1144</page-range></nlm-citation>
</ref>
<ref id="B14">
<label>14</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Sediki]]></surname>
<given-names><![CDATA[H]]></given-names>
</name>
<name>
<surname><![CDATA[Bechouche]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
<name>
<surname><![CDATA[Abdeslam]]></surname>
<given-names><![CDATA[D. O]]></given-names>
</name>
<name>
<surname><![CDATA[Haddad]]></surname>
<given-names><![CDATA[S]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Adaline approach for induction motor mechanical parameters identification"]]></article-title>
<source><![CDATA[Mathematics and Computers in Simulation]]></source>
<year>2012</year>
<numero>0</numero>
<issue>0</issue>
</nlm-citation>
</ref>
<ref id="B15">
<label>15</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Meghnous]]></surname>
<given-names><![CDATA[A]]></given-names>
</name>
<name>
<surname><![CDATA[Pham]]></surname>
<given-names><![CDATA[M]]></given-names>
</name>
<name>
<surname><![CDATA[Lin-Shi]]></surname>
<given-names><![CDATA[X]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Dynamic identification of a synchronous machine using an extended sliding mode observer"]]></article-title>
<source><![CDATA[Mathematics and Computers in Simulation]]></source>
<year>2012</year>
<numero>0</numero>
<issue>0</issue>
</nlm-citation>
</ref>
<ref id="B16">
<label>16</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Huber]]></surname>
<given-names><![CDATA[P]]></given-names>
</name>
</person-group>
<source><![CDATA[Robust Statistics]]></source>
<year>2009</year>
<publisher-loc><![CDATA[New York ]]></publisher-loc>
<publisher-name><![CDATA[Wiley Series in Probability and Statistics]]></publisher-name>
</nlm-citation>
</ref>
<ref id="B17">
<label>17</label><nlm-citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Pearson]]></surname>
<given-names><![CDATA[R]]></given-names>
</name>
</person-group>
<article-title xml:lang="en"><![CDATA["Exploring process data"]]></article-title>
<source><![CDATA[<img border=0 width=161 height=17 src="img/revistas/inde/v32n1/v32n1a08f121.jpg">]]></source>
<year>2001</year>
<volume>11</volume>
<numero>2</numero>
<issue>2</issue>
<page-range>179-194</page-range></nlm-citation>
</ref>
<ref id="B18">
<label>18</label><nlm-citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname><![CDATA[Bishop]]></surname>
<given-names><![CDATA[C. M]]></given-names>
</name>
</person-group>
<source><![CDATA[Pattern Recognition and Machine Learning]]></source>
<year>2006</year>
<publisher-loc><![CDATA[Secaucus^eNJNew York^eUSA NJUSA]]></publisher-loc>
<publisher-name><![CDATA[Springer Verlag]]></publisher-name>
</nlm-citation>
</ref>
</ref-list>
</back>
</article>
