Informe_Pendulo_Invertido__Copy__(2)
Informe_Pendulo_Invertido__Copy__(2)
Proyecto no guiado
I Semestre 2026
1. Resumen
En este proyecto se realizó el análisis, identificación y diseño de control para una planta de
péndulo invertido sobre carro. Este sistema representa un caso clásico de control automático debido
a que es naturalmente inestable alrededor de la posición vertical superior y, además, es un sistema
subactuado, ya que se cuenta con una única entrada de control aplicada al motor del carro, mientras
se desea regular simultáneamente la posición del carro y el ángulo del péndulo. Para obtener el
modelo del sistema se utilizó un enfoque experimental basado en identificación de sistemas. En lugar
de desarrollar un modelo físico completo a partir de masas, inercias, longitudes y coeficientes de
fricción, se obtuvieron funciones de transferencia para la posición del carro y para el ángulo del
péndulo a partir de datos medidos en la planta. Posteriormente, estas funciones de transferencia
fueron transformadas a una representación en variables de estado y se construyó un modelo SIMO,
con una entrada y dos salidas. A partir del modelo obtenido se diseñó una estrategia de control
por realimentación de estados, considerando tanto ubicación de polos como control LQR. El diseño
final se basó en un regulador LQR con sistema aumentado, de manera que se pudiera penalizar
fuertemente el ángulo del péndulo, controlar el desplazamiento del carro y mantener un compromiso
razonable entre rapidez de respuesta y esfuerzo de control. Durante el desarrollo se identificaron
limitaciones importantes de la planta real, tales como backlash en la transmisión, fricción no lineal,
zona muerta del motor, ruido en los encoders y posibles corrimientos del cero angular. Estos factores
explican por qué la respuesta experimental puede diferir de la respuesta simulada.
2. Introducción
El péndulo invertido sobre carro es uno de los sistemas más utilizados para estudiar técnicas
de control moderno, debido a que combina inestabilidad natural, dinámica acoplada, restricciones
físicas y sensibilidad a perturbaciones. A diferencia de un péndulo convencional, cuya posición esta-
ble se encuentra hacia abajo, el péndulo invertido debe mantenerse alrededor de la posición vertical
superior. Esta condición genera una dinámica inestable en lazo abierto, por lo que cualquier pe-
queña desviación angular tiende a crecer si no existe una acción de control adecuada. La planta de
laboratorio consiste en un carro que se desplaza horizontalmente sobre un riel, impulsado por un
motor de corriente directa mediante una transmisión mecánica. Sobre el carro se encuentra montado
el péndulo, cuyo movimiento angular se mide mediante un encoder. También se mide la posición
horizontal del carro mediante otro encoder. La entrada de control corresponde a la tensión aplicada
al motor, mientras que las variables de interés son la posición del carro y el ángulo del péndulo.
Una característica importante del sistema es que se trata de una planta subactuada. Esto signifi-
ca que el sistema posee más variables relevantes que actuadores disponibles. En este caso, solo se
puede actuar directamente sobre el movimiento del carro, pero el ángulo del péndulo no tiene un
actuador independiente. Por tanto, la estabilización angular debe lograrse indirectamente moviendo
el carro de forma adecuada. Esta condición hace que el diseño del controlador sea más exigente, ya
que el movimiento necesario para estabilizar el péndulo puede entrar en conflicto con el objetivo de
mantener el carro cerca del centro del riel. Durante el trabajo se observaron varias no idealidades
de la planta real. Entre ellas se encuentran la fricción, el backlash de la transmisión, la zona muerta
del motor, el ruido en las mediciones y posibles offsets del encoder del ángulo. Estas no idealidades
dificultan que un modelo lineal reproduzca exactamente el comportamiento físico. En particular, se
observó que el sistema real puede amortiguarse más rápido que el modelo, probablemente por efectos
1
de fricción no lineal que no son capturados por una función de transferencia lineal e invariante en
el tiempo. Por esta razón, el modelo utilizado en este proyecto se obtuvo mediante identificación
experimental. Esta decisión permite representar de manera más directa la dinámica dominante ob-
servada en la planta, aunque no incluye de forma explícita todos los parámetros físicos del sistema.
Posteriormente, el modelo identificado fue utilizado para diseñar controladores en variables de es-
tado, verificando su comportamiento mediante simulación antes de su implementación en la planta.
El informe se organiza de la siguiente manera. Primero se presentan los objetivos del proyecto y las
características principales de la planta. Luego se describe el proceso de identificación del sistema,
el modelo obtenido y la construcción de la representación en variables de estado. Posteriormente se
presenta el diseño del controlador, la estrategia LQR seleccionada y la implementación experimental.
Finalmente, se discuten las preguntas del proyecto, las limitaciones del modelo, las conclusiones y
recomendaciones.
El sistema de péndulo invertido se basa en dos funciones de transferencia que describen la po-
sición y el ángulo. Se pueden describir matemáticamente, y encontrar cada parámetro como masa,
inercia, fricción viscosa, o realizar una aproximación no paramétrica a partir de experimentos donde
se registre la respuesta del sistema. En nuestro caso escogimos la segunda, ya que nos aportaba un
modelo más cercano a la realidad al ser difícil medir con precisión las variables del modelo matemá-
tico. En el caso de la posición esta se puede describir con un integrador, debido a que la posición
se acumula a medida que avanza el carro, y no se pierde, simulando el efecto de integrador junto
al sistema de primer orden. El ángulo se puede definir como un sistema de segundo orden, lo cual
describe correctamente la oscilación de un péndulo y su amortiguamiento. Este cuenta con un cero,
que describe las aceleraciones del carrito o cambios bruscos, asemejando un derivador.
El sistema de pendulo invertido es similar al de una grua o pendulo no invertido, la diferencia
radia en los polos del ángulo, al ser invertido el ángulo debe tener un polo en el semiplano positivo,
el cual representa la inestabilidad del mismo. A continuación las ecuacines teóricas del péndulo
invertido
El modelo SISO de un sistema en forma canónica controlable es el siguiente
kp kp
P osicin(s) = = 2 (1)
s(s + ax ) s + sax + 0
ṗ 0 1 p 0
ẋ = = + u (2)
p̈ 0 −ax ṗ kp
p
(3)
p= 1 0
ṗ
Se aumenta las dimensiones de la matriz para luego unificar el modelo SISO del ángulo y posición
en un SIMO.
ṗ 0 0 10 p 0
θ̇ 0 0 00 θ 0
= + u (4)
p̈ 0 0 −ax
0 ṗ kp
θ̈ 0 0 00 θ̇ 0
p
θ
(5)
p= 1 0 0 0 ṗ
θ̇
2
En el caso del ángulo, el signo del último término del sistema subamortiguado es el que representa
el cambio de sistema pendulo o grúa a péndulo invertido.
ka (s + za ) ka s + ka · za
ngulo(s) = = 2 (6)
s2
+ 2ζωn s + ωn2 s + 2ζωn s + (−ωn2 )
θ̇ 0 1 θ 0
ẋ = = + u (7)
θ̈ −(−ωn2 ) −2ζωn θ̇ 1
θ
(8)
θ = ka · za ka
θ̇
En forma canónica controlable
θ̇ 0 1 θ ka
ẋ = = + u (9)
θ̈ −(−ωn2 ) −2ζωn θ̇ ǩ ∗ za
θ
(10)
θ= 1 0
θ̇
Se aumenta las dimensiones del sistema.
ṗ 0 0 0 0 p 0
θ + ka u
θ̇ 0 0 0 1
(11)
=
p̈ 0 0 0 0 ṗ 0
θ̈ 0 −(−ωn2 ) 0 −2ζωn θ̇ ǩ ∗ za
p
θ
(12)
θ= 0 1 0 0
ṗ
θ̇
Se unen ambos sistemas SISO expandidos en uno SIMO.
ṗ 0 0 1 0 p 0
θ + ka u
θ̇ 0 0 0 1
(13)
=
p̈ 0 0 ax 0 ṗ kp
2
0 ωn 0 −2ζωn
θ̈ θ̇ ǩ ∗ za
p
p θ
1 0 0 0
(14)
=
θ 0 1 0 0 ṗ
θ̇
Esta es la representación del sistema en ecuaciones de estado, siendo la posición y el ángulo las
salidas.
===============================================================
3
3. Objetivo
3.1. Objetivo general
Identificar experimentalmente el modelo dinámico de una planta de péndulo invertido sobre carro
y diseñar un controlador por realimentación de estados que permita estabilizar el péndulo alrededor
de la posición vertical superior, manteniendo el carro dentro de su rango físico de operación.
2. Obtener modelos experimentales para la dinámica de posición del carro y del ángulo del pén-
dulo.
4
4.1. Caracterización de los componentes principales
Debido a que el proyecto corresponde a un experimento no guiado, fue necesario caracterizar los
componentes principales de la planta antes de realizar la identificación y el diseño del controlador.
Esta caracterización permite comprender la función de cada elemento dentro del sistema, así como
las limitaciones físicas que pueden afectar el desempeño experimental.
5
Por tanto, aunque el motor tiene una tensión nominal de 12 VDC, experimentalmente se utilizó
un límite de operación cercano a ±16 V para mejorar la autoridad de control del sistema. Esta
condición permitió obtener una respuesta más rápida del carro, necesaria para estabilizar el péndulo
invertido, pero también implica operar por encima de la condición nominal, por lo que debe evitarse
una saturación prolongada.
La caja reductora permite disminuir la velocidad de giro del motor y aumentar el torque dispo-
nible en el eje de salida. Esto es útil para mover el carro con suficiente fuerza; sin embargo, también
introduce efectos no ideales como fricción interna, holgura mecánica y backlash.
E5-900-250-NE-S-D-3-B. (21)
Además, se observaron códigos individuales como M212922 y M1605R6. Estos códigos se inter-
pretan como identificadores de fabricación, lote, serie o revisión, y no como el modelo principal del
encoder.
Uno de los encoders se utiliza para medir la posición horizontal del carro sobre el riel. Esta
señal permite conocer el desplazamiento del carro y verificar que se mantenga dentro de la carrera
disponible. El segundo encoder se utiliza para medir el ángulo del péndulo respecto a la vertical.
Esta segunda medición es fundamental, ya que el controlador utiliza el ángulo para estabilizar el
péndulo invertido.
6
El valor 900 en el código del encoder indica una resolución nominal de 900 ciclos por revolución.
Al utilizar decodificación en cuadratura, considerando los cuatro flancos de las señales A y B, la
resolución efectiva puede expresarse como:
360◦
∆θ = = 0,1◦ /cuenta. (23)
3600
En radianes:
2π
∆θ = = 1,745 × 10−3 rad/cuenta. (24)
3600
Aunque esta resolución es adecuada para el experimento, la calidad de la medición angular puede
verse afectada por vibraciones, backlash, ruido eléctrico y posibles desplazamientos del cero del
encoder. Estos factores deben considerarse al interpretar la respuesta experimental del controlador.
Figura 1: Planta de péndulo invertido sobre carro. Se puede observar el riel horizontal, el carro
móvil, el péndulo articulado, la transmisión de banda, los encoders en ambos extremos y la tarjeta
de control montada sobre el carro.
7
Figura 2: Fuente de alimentación de corriente directa GW Instek GPE-3323 configurada a 16 V y
2,99 A para alimentar la planta. También se observa la tarjeta de comunicación que implementa el
protocolo UART mediante Bluetooth, encargada de la transmisión de datos entre la computadora y
la planta en tiempo real.
8
Figura 3: Circuito driver de puente H utilizado para controlar el motor de corriente directa. El
módulo cuenta con dos puentes H independientes, aunque para este proyecto solo se utilizó uno
de ellos. Este circuito recibe la señal PWM generada por PASCAL y la convierte en la tensión de
potencia aplicada al motor.
9
4.3. Descripción del experimento
Para la caracterización de la planta se realizaron tres pruebas principales. Estas pruebas permitie-
ron obtener información sobre el modelo dinámico, la zona muerta del actuador y el comportamiento
inestable del sistema como péndulo invertido.
10
y mantener el carro exactamente fijo al mismo tiempo no siempre es posible ante perturbaciones,
porque el movimiento del carro es precisamente el mecanismo que permite estabilizar el ángulo. Por
esta razón, el diseño del controlador debe dar prioridad al ángulo, pero sin descuidar la posición del
carro.
11
4.7. Offset del encoder y ruido angular
Otra consideración importante es el offset del encoder del ángulo. En algunas pruebas, el cero
del encoder no necesariamente coincidía con la posición vertical real del péndulo. Esto implica que el
controlador puede interpretar una pequeña inclinación como si fuera cero, o puede intentar corregir
un error que en realidad corresponde a un desplazamiento de referencia.
Además, el backlash y la fricción pueden producir pequeñas oscilaciones en la medición angular.
Estas oscilaciones afectan especialmente el cálculo de θ̇, ya que derivar una señal ruidosa amplifica
sus componentes de alta frecuencia. Por esta razón, el uso de filtrado o una estimación cuidadosa de
la velocidad angular resulta fundamental.
12
El diagrama está compuesto por dos encoders de cuadratura (QDEC_0 y QDEC_1), dos deri-
vadores (DER_0 y DER_1), una constante (CTE_0), un restador (REST_0), el bloque de retro-
alimentación de estados (RE_0), un saturador (SAT_0) y un generador PWM (PWM_0). La señal
de referencia también se incluye como entrada al bloque RE.
La ganancia de cada encoder se determina a partir de la relación física entre las cuentas del
encoder y la magnitud física medida.
13
5.2. Constante y restador (CTE_0 y REST_0)
14
Figura 7: Configuración de los derivadores DER_0 y DER_1 (ambos con los mismos parámetros).
Estos derivadores no son diferenciadores puros, sino filtros derivativos de segundo orden. Un
diferenciador puro amplifica el ruido de alta frecuencia de forma ilimitada, lo que resulta especial-
mente problemático para la señal del ángulo, que como se mencionó anteriormente, presenta ruido
de cuantificación del encoder y vibraciones mecánicas. El filtro derivativo limita esta amplificación
al introducir una dinámica de segundo orden que atenúa las componentes de alta frecuencia.
La frecuencia natural ωn = 60 rad/s define hasta qué frecuencia el bloque se comporta apro-
ximadamente como un diferenciador puro. Por encima de esta frecuencia, la ganancia se atenúa
progresivamente. Un valor más alto de ωn permite derivar señales más rápidas con menos retraso,
pero también deja pasar más ruido. Un valor más bajo filtra mejor el ruido, pero introduce más
retraso en la estimación de velocidad, lo que puede afectar la estabilidad del lazo cerrado.
El amortiguamiento ζ = 0,85 determina el comportamiento transitorio del filtro. Un valor cercano
a 1 produce una respuesta sin oscilaciones en la transición entre la zona de derivación y la zona
de atenuación. Un amortiguamiento menor genera una resonancia alrededor de ωn , lo que puede
amplificar ciertas frecuencias en lugar de atenuarlas. El valor de 0.85 representa un compromiso
entre suavidad de la respuesta y velocidad de seguimiento.
15
5.4. Retroalimentación de estados (RE_0)
(33)
K = −122,51926701 −1231,51732237 −263,30596057 −131,28575213 .
Los cuatro estados corresponden en orden a la posición del carrito p, el ángulo del péndulo θ, la
velocidad del carrito ṗ y la velocidad angular θ̇.
16
5.5. Experimento de zona muerta
Para caracterizar la zona muerta del motor se diseñó un experimento específico en PASCAL. La
zona muerta corresponde al rango de tensiones de entrada para el cual el motor no genera movimiento
apreciable en el carrito, debido a la fricción estática de la transmisión, la caja reductora y la banda
en V.
El experimento consistió en aplicar una señal de tipo trapecio como entrada al motor, configurada
únicamente con tiempo de rampa y sin tiempo en alto, de modo que la señal sube linealmente desde
cero hasta el valor final a lo largo de 60 segundos. Este tiempo de subida largo garantiza que la
tensión aumente muy lentamente, lo que permite identificar con precisión el valor exacto en que el
motor comienza a mover el carrito antes de que la inercia o la velocidad acumulada distorsionen la
medición.
17
Figura 10: Configuración de la señal trapecio para el lado positivo (valor final +7 V) y negativo
(valor final −7 V).
Se realizaron dos pruebas independientes: una con rampa positiva hasta +7 V para caracterizar
el sentido de avance, y otra con rampa negativa hasta −7 V para caracterizar el sentido de retroceso.
En ambos casos se registró simultáneamente la señal de entrada y la posición del carrito medida por
el encoder QDEC_0.
Figura 11: Experimento de zona muerta lado positivo. La posición del carrito permanece constante
hasta que la entrada alcanza 5,63 V.
18
Figura 12: Experimento de zona muerta lado negativo. La posición del carrito permanece constante
hasta que la entrada alcanza −5,296 V.
19
5.6. Saturador (SAT_0)
El saturador limita la señal de control entre −16 V y +16 V. Este límite se seleccionó para
proporcionar al motor suficiente autoridad de control, es decir, suficiente tensión para que el carrito
pueda alcanzar las velocidades y aceleraciones que el controlador requiere para estabilizar el péndulo.
Aunque la tensión nominal del motor es de 12 VDC, operar con un límite de 16 V aumenta la
velocidad y aceleración disponibles, lo cual es necesario dado que el péndulo invertido requiere
correcciones rápidas ante perturbaciones.
20
5.7. PWM (PWM_0)
El bloque PWM_0 convierte la señal de control en una señal de modulación por ancho de pulso
para el driver del motor. Se configuró en modo estándar con un periodo de 500 µs y una amplitud
de 16. Adicionalmente, se activó la corrección de zona muerta con límite superior de 4 e inferior de
−4.
La zona muerta del motor corresponde al rango de tensiones en el que el motor no genera
movimiento apreciable del carrito, debido a la fricción estática de la transmisión y la caja reductora.
Para caracterizar este valor se realizó un experimento específico: se aplicó una señal de rampa de
crecimiento lento sobre la entrada del motor y se observó el valor a partir del cual el carrito comenzaba
a desplazarse de forma apreciable. Este procedimiento permitió identificar empíricamente el umbral
de la zona muerta, que resultó corresponder a un valor de 4 en las unidades internas de la plataforma.
La corrección de zona muerta en el bloque PWM compensa este efecto sumando automáticamente el
valor umbral a cualquier señal no nula, de modo que el motor recibe siempre una tensión suficiente
para vencer la fricción estática.
21
matemáticos que reproduzcan la relación entre entrada y salida. En este proyecto se identificaron
dos funciones de transferencia principales: una para la posición del carro y otra para el ángulo del
péndulo. Se escogió este enfoque porque obtener un modelo físico completo requeriría conocer con
precisión parámetros como masas, inercias, longitud efectiva del péndulo, fricción, constantes del
motor y características de la transmisión. En la planta real, varios de estos parámetros son difíciles
de medir con exactitud y pueden cambiar según el estado mecánico del sistema. Por tanto, la iden-
tificación experimental permite obtener un modelo práctico que represente la dinámica dominante
observada.
Esta función de transferencia representa la relación entre la entrada aplicada al motor y la posición
horizontal del carro.
22
Figura 15: Validación del modelo de posición. Ajuste del 97.18 % entre el modelo identificado y los
datos experimentales.
Figura 16: Validación del modelo angular. Ajuste del 83.28 % entre el modelo identificado y los datos
experimentales.
El modelo de posición presenta un ajuste muy alto, lo cual indica que la estructura de integrador
con dinámica de primer orden representa adecuadamente el comportamiento del carro ante la entrada
aplicada. El modelo angular presenta un ajuste menor, lo cual es esperable dado que el péndulo exhibe
efectos no lineales como fricción en el pivote, backlash en la transmisión y ruido en el encoder, ninguno
de los cuales está incluido en el modelo lineal.
Para complementar la validación, se realizó una simulación adicional mediante lsim en MATLAB,
aplicando un escalón de 7,5 V tanto al modelo identificado como comparándolo contra los datos
experimentales del ángulo.
Figura 17: Comparación entre la respuesta simulada del modelo angular y los datos experimentales
ante un escalón de 7,5 V.
23
Como se observa en la figura, la frecuencia de oscilación del modelo coincide bien con la frecuen-
cia experimental, lo que indica que los parámetros ωn y ζ identificados capturan adecuadamente la
dinámica oscilatoria dominante del péndulo. Sin embargo, el modelo presenta un amortiguamiento
ligeramente menor que el experimental: la amplitud de las oscilaciones simuladas decae más lenta-
mente que en la planta real. Esta diferencia es atribuible a la fricción no lineal presente en el pivote
del péndulo, el riel y la transmisión, la cual no está incluida explícitamente en el modelo lineal
identificado pero produce una disipación adicional de energía en la planta real.
En conjunto, los resultados de validación indican que el modelo obtenido es suficientemente
representativo de la dinámica dominante del sistema para ser utilizado en el diseño del controlador,
aunque no reproduce con exactitud todos los efectos físicos presentes en la planta.
24
1 0 0 0
C= , (45)
0 1 0 0
0
D= . (46)
0
Este modelo representa la dinámica acoplada entre la posición del carro y el ángulo del péndulo.
Aunque se obtuvo a partir de dos modelos SISO, la representación final permite analizar el sistema
como una planta SIMO, donde una única entrada afecta simultáneamente dos salidas.
25
8.3. Sistema aumentado
Para considerar el error acumulado de posición, se utilizó un sistema aumentado con un estado
integral. Este estado permite penalizar el error de posición del carro y ayuda a mantenerlo cerca de
la referencia deseada. El sistema aumentado se define como:
A 0
As = , (49)
−Cpos 0
B
Bs = . (50)
0
La acción integral se aplicó sobre la posición del carro y no directamente sobre el ángulo. Esto se
debe a que el ángulo es una señal muy sensible al ruido, al backlash y al offset del encoder. Si
se integrara directamente el error angular, el controlador podría acumular errores producidos por
mediciones imperfectas y generar una acción de control innecesaria o incluso perjudicial.
La matriz Q penaliza los estados del sistema, mientras que R penaliza el esfuerzo de control. Si
un estado tiene un peso mayor en Q, el controlador tratará de reducirlo con mayor prioridad. Si R
aumenta, el controlador se vuelve menos agresivo porque se penaliza más el uso de la entrada. Para
este proyecto se utilizaron los siguientes pesos:
qp = 10, (52)
qθ = 980000, (53)
qṗ = 400000, (54)
qθ̇ = 500, (55)
qi = 100. (56)
Por tanto:
Q = diag (10, 980000, 400000, 500, 100) , (57)
R = 1. (58)
La selección de estos pesos refleja la prioridad del diseño. El ángulo tiene un peso muy alto porque
es la variable más crítica para evitar la caída del péndulo. La velocidad de posición también se
penaliza fuertemente, con el objetivo de evitar movimientos excesivamente bruscos del carro. La
velocidad angular se penaliza con un peso menor que el ángulo, porque esta señal puede contener
ruido derivativo. Finalmente, el peso integral se mantiene relativamente bajo para evitar que la
acción integral domine el control.
26
8.5. Ganancia obtenida
Con los pesos anteriores se obtuvo la siguiente ganancia de realimentación de estados:
(59)
KLQR = −122,51926701 1231,51732237 263,30596057 131,28575213 .
Es importante aclarar que el signo de implementación debe verificarse directamente con la convención
utilizada en la planta. En MATLAB, la forma usual de realimentación es u = −Kx. Por esta razón, no
necesariamente se debe multiplicar el vector completo de ganancias por −1 antes de implementarlo.
Si el programa en PASCAL ya aplica el signo negativo de la realimentación, invertir nuevamente el
signo de las ganancias puede producir una acción de control incorrecta. Además, el signo efectivo
depende de la polaridad del motor, del sentido positivo de movimiento del carro y de la convención
del encoder angular. Por tanto, la validación experimental del signo es indispensable.
Ganancia Valor
kp -122.51926701
kθ 1231.51732237
kṗ 263.30596057
kθ̇ 131.28575213
Cabe mencionar que para el experimento físico se tenían todas las ganancias negativas, esto hizo
que el sistema real con sus referencais de signo recibiera la correcta acción correctiva para estabilizar
el sistema.
9. Verificación en simulación
Antes de implementar el controlador en la planta real, se verificó el comportamiento del sistema
mediante simulación. Esta etapa es necesaria porque permite revisar la estabilidad del lazo cerrado
sin riesgo de dañar la planta o provocar la caída repetida del péndulo. La simulación permite observar
27
la evolución de la posición, el ángulo y la acción de control. Además, permite verificar si la ganancia
seleccionada estabiliza el polo inestable del sistema y si la respuesta se mantiene dentro de los rangos
esperados. A continuación el diagrama de simulink.
28
mejarse a la dinámica real, como la fricción viscosa, el ruido, factores no lineales como la gravedad
dentro del sistema aproximado o el backlash. Se puede observar un tiempo de estabilización muy
grande y nula oscilación, lo cual no fue observado en el experimento, así mismo, el sistema se com-
portaba mucho más agresivo en la planta real, para poder estabilizar el péndulo invertido venciendo
el efecto continuo de la gravedad. Finalmente, no se aplicó al final el uso de la parte integral en el
experimento, por razones explicadas posteriormente.
A los 38 s se colocó una cinta métrica sobre el péndulo y adicionalmente se golpeó la planta,
combinando una perturbación de masa con un impulso mecánico.
29
A los 116 s se colocó un peso calibrado sobre el péndulo, correspondiente a una masa de
referencia metrológica utilizada en balanzas de precisión.
A los 193 s se retiró el peso calibrado y de forma inmediata se colocó un martillo, representando
una perturbación de masa significativamente mayor y con distribución de masa diferente.
Estas perturbaciones permiten observar cómo el controlador responde ante cambios en la diná-
mica del péndulo, variaciones de masa y excitaciones impulsivas.
La señal RE corresponde a la salida del bloque de realimentación de estados antes de pasar por
el saturador y el PWM. Esta señal representa la acción de control calculada por el controlador a
partir de la combinación lineal de los estados: posición del carro, ángulo del péndulo, velocidad del
carro y velocidad angular.
Durante los primeros segundos después del inicio, la señal RE se mantiene cercana a cero, lo que
indica que el péndulo se encontraba aproximadamente en la posición de equilibrio y el controlador
no requería acción significativa. A partir de la primera perturbación en el segundo 38, se observan
picos pronunciados en la señal RE, los cuales corresponden a la respuesta del controlador ante la
combinación de masa adicional y golpe impulsivo.
Los picos negativos más profundos que se observan alrededor del segundo 50 y posteriormente
en el intervalo entre los segundos 116 y 193 indican que el controlador requirió acciones de control
intensas para compensar perturbaciones de mayor magnitud. Esto es coherente con el hecho de que
en esos intervalos se aplicaron masas más pesadas o perturbaciones más enérgicas.
30
11.3. Velocidad angular del péndulo
La velocidad angular del péndulo se obtuvo a partir de la derivación numérica de la señal del
encoder angular. Debido a que la medición del ángulo contiene ruido proveniente del encoder, de
vibraciones mecánicas y del backlash de la transmisión, la señal derivada amplifica considerablemente
estas perturbaciones. Como resultado, la velocidad angular presenta un nivel de ruido elevado durante
todo el experimento, incluso en los tramos donde el péndulo se encuentra cerca del equilibrio.
Este comportamiento es esperado y se menciona explícitamente porque el controlador utiliza esta
señal como uno de sus estados de realimentación. Un ruido elevado en θ̇ implica que la ganancia kθ̇
no puede ser demasiado alta, ya que amplificaría ese ruido directamente en la acción de control. Por
ello, en el diseño LQR se asignó un peso relativamente bajo a la velocidad angular en comparación
con el ángulo mismo.
En los intervalos donde se aplicaron perturbaciones, especialmente entre los segundos 116 y 193
con el peso calibrado, y entre 193 y 213 con el martillo, se observan excursiones más amplias en la
velocidad angular. Estas excursiones reflejan las oscilaciones más enérgicas del péndulo ante masas
considerablemente mayores, situación ante la cual el controlador debía actuar con mayor intensidad
para evitar la caída.
31
11.4. Velocidad del carrito
La velocidad del carrito refleja la respuesta dinámica del carro ante las acciones de control
calculadas. Un aspecto fundamental del funcionamiento del péndulo invertido es que el carrito debe
moverse en la dirección del péndulo para corregir su error angular. Esto puede parecer contraintuitivo
a primera vista, pero es el mecanismo físico que permite la estabilización.
Cuando el péndulo se inclina hacia un lado, el controlador ordena al carro que se desplace en
esa misma dirección. Este movimiento genera una aceleración horizontal en la base del péndulo, lo
que produce una fuerza de inercia que tiende a empujar la parte superior del péndulo en sentido
contrario al movimiento del carro. De esta forma, el péndulo es empujado de regreso hacia la posición
vertical. Si el carro se moviera en dirección contraria al péndulo, la aceleración generada aumentaría
la desviación angular en lugar de corregirla.
Este principio explica por qué la posición del carro no permanece fija durante el experimento:
el carro se mueve constantemente para generar las aceleraciones necesarias que mantienen el ángulo
cerca de cero. La velocidad del carrito presenta picos en los mismos instantes que las perturbaciones,
confirmando que el controlador reacciona ante cada excitación moviendo activamente el carro.
32
11.5. Ángulo del péndulo
El ángulo del péndulo es la variable más crítica del sistema, ya que el objetivo principal del
controlador es mantenerla cercana a cero. Como se puede observar en la gráfica, durante la mayor
parte del experimento el ángulo se mantiene en un rango pequeño alrededor de cero, lo que indica
que el controlador logra estabilizar el péndulo ante las perturbaciones aplicadas.
Las desviaciones más notorias ocurren en los instantes inmediatamente posteriores a cada pertur-
bación, especialmente cuando se aplicaron masas o golpes. Después de cada excitación, el controlador
recupera el equilibrio en un tiempo relativamente corto. Esto confirma que el diseño LQR con el peso
elevado asignado al ángulo (qθ = 980000) permite priorizar correctamente la estabilización angular.
Se observa además que la señal del ángulo presenta cierto nivel de ruido, aunque menor que el
de la velocidad angular. Esto es consistente con la naturaleza de la medición: el encoder registra
directamente la posición angular con una resolución de 0,1◦ por cuenta, por lo que el ruido en el
ángulo es el cuantificado por el encoder. Al derivar esta señal para obtener θ̇, el ruido se amplifica
considerablemente, tal como se describió en la subsección anterior.
En los intervalos con masa pesada, particularmente entre los segundos 116 y 193 con el peso
calibrado y entre 193 y 213 con el martillo, el ángulo presenta oscilaciones de mayor amplitud.
Esto es esperable porque la masa adicional cambia la dinámica del péndulo respecto al modelo
identificado, lo que representa un error de modelado que el controlador debe compensar en tiempo
real.
33
11.6. Posición del carrito
Prueba 1: [Link]
Prueba 2: [Link]
34
11.8. Comparación entre simulación y experimento
La respuesta experimental difiere de la simulación en varios aspectos. En primer lugar, la señal
angular experimental presenta más ruido del que el modelo predice, debido a los efectos del en-
coder, el backlash y las vibraciones mecánicas. En segundo lugar, la respuesta real puede mostrar
mayor amortiguamiento que la simulada, lo cual es atribuible a la fricción no lineal presente en el
riel, el pivote del péndulo y la transmisión mecánica. Estos efectos disipativos no están incluidos
explícitamente en el modelo lineal identificado.
En tercer lugar, las perturbaciones de masa aplicadas durante el experimento representan cambios
reales en la dinámica del sistema que no estaban contemplados en el modelo de diseño. El controlador
LQR fue diseñado para el péndulo sin masa adicional, por lo que al colocar objetos sobre él, el
punto de equilibrio y la dinámica angular cambian. Pese a esto, el controlador logra mantener la
estabilización en la mayoría de los casos, lo que indica que el diseño tiene un margen de robustez
suficiente para tolerar estas variaciones.
Finalmente, la zona muerta del motor implica que acciones de control pequeñas no producen
movimiento real en el carro. Esto puede observarse en los tramos tranquilos del experimento, donde
la señal RE es pequeña pero el carro prácticamente no se mueve. Esta no linealidad no está incluida
en el modelo y representa una diferencia estructural entre la simulación y la planta real.
35
backlash y fricción. Estas no idealidades pueden afectar la calidad de los datos. Por esta razón,
aunque el modelo obtenido resulta útil para el diseño del controlador, debe interpretarse como una
aproximación válida alrededor de la región donde se realizó la identificación.
36
implementable. Un modelo de orden mayor puede ajustarse mejor a ciertos detalles de los datos
experimentales, pero también puede capturar ruido, vibraciones mecánicas o efectos no repetibles
de la planta. Además, un modelo de mayor orden aumenta la complejidad del controlador. En
una implementación experimental, esto puede complicar el cálculo de estados, la estimación de
velocidades y el filtrado de señales. También puede generar un controlador más sensible a errores de
modelado. Para este sistema, la posición puede representarse razonablemente con un integrador y una
dinámica de primer orden, mientras que el ángulo se representa con una dinámica de segundo orden
inestable. Esta estructura captura los elementos físicos más importantes: acumulación de posición,
oscilación angular e inestabilidad del péndulo invertido.
37
12.10. ¿Cuál es el compromiso entre rapidez y desempeño?
Una respuesta rápida puede parecer deseable, pero en la planta real no siempre es la mejor
opción. Para lograr una respuesta más rápida se requieren ganancias más altas, lo que aumenta
la acción de control. Esto puede provocar saturación del motor, movimientos bruscos del carro y
mayor sensibilidad al ruido. En el péndulo invertido, una acción muy agresiva puede excitar las
holguras mecánicas de la transmisión. El backlash produce pequeños movimientos no lineales que el
controlador puede interpretar como errores reales, generando oscilaciones. Además, si la velocidad
angular se estima a partir de una señal ruidosa, ganancias altas pueden amplificar ese ruido. Por
tanto, el controlador debe ser suficientemente rápido para evitar que el péndulo caiga, pero no tan
agresivo como para saturar el actuador o amplificar las no idealidades. Este compromiso fue una de
las razones para utilizar LQR y ajustar cuidadosamente los pesos.
38
después de ajustar mecánicamente la planta. Esto permitiría obtener un modelo más consistente y
posiblemente mejorar el desempeño del controlador.
13. Conclusiones
Se obtuvo un modelo experimental para la planta de péndulo invertido mediante identificación
de sistemas, utilizando una función de transferencia para la posición del carro y otra para el
ángulo del péndulo.
La planta corresponde a un sistema subactuado, ya que solo existe una entrada de control
aplicada al motor del carro, mientras se desea regular tanto la posición del carro como el
ángulo del péndulo.
El cálculo de backlash mostró una holgura angular aproximada de 6,61◦ , valor significativo
para el control del péndulo, ya que puede afectar la medición angular y generar una zona
muerta mecánica.
La acción integral se consideró sobre la posición del carro y no directamente sobre el ángulo,
debido a que la señal angular puede estar afectada por ruido, backlash y offset del encoder.
Las diferencias entre simulación y experimento se explican por efectos físicos no incluidos
completamente en el modelo lineal, como fricción no lineal, zona muerta, backlash, saturación
del actuador y ruido en los encoders.
14. Recomendaciones
Revisar el estado mecánico de la planta antes de cada sesión experimental, verificando tornillos,
arandelas, separadores, tensión de la banda y sujeción del pivote del péndulo.
Validar el cero del encoder angular antes de cada prueba para evitar que el controlador trabaje
con un punto de equilibrio incorrecto.
Evitar diseños excesivamente agresivos, ya que pueden producir saturación del motor y au-
mentar la sensibilidad al backlash y al ruido.
39
Repetir la identificación después de corregir problemas mecánicos, con el fin de obtener un
modelo más representativo de la planta real.
Referencias
[1] R. C. Dorf y R. H. Bishop, Sistemas de Control Moderno, 10.ª ed. Prentice Hall, 2005.
[2] K. Ogata, Ingeniería de Control Moderna, 4.ª ed. Pearson Prentice Hall, 2003.
[3] K. Ogata, Problemas de Ingeniería de Control utilizando MATLAB. Prentice Hall, 1999.
40