SOLUCIÓN GRÁFICA DE ECUACIONES DIFERENCIALES EN SIMULINK
1. OBJETIVO.
i. Simular la solución de ecuaciones diferenciales de segundo orden con el entorno
Simulink de Matlab.
2. INTRODUCCIÓN.
Simulink es una extensión gráfica de Matlab que permite construir y simular modelos de
sistemas físicos y sistemas de control mediante diagramas de bloques. El comportamiento
de dichos sistemas se define mediante funciones de transferencia, operaciones
matemáticas, elementos de Matlab y señales predefinidas de todo tipo. Simulink dispone
de una serie de utilidades que facilitan la visualización, análisis y guardado de los
resultados de simulación.
Graficaremos la solución de las siguientes ecuaciones diferenciales.
1 1 1
2𝑥̈ − 𝑥̇ + 𝑥 = sin 2𝑡 (1)
2 4 4
Con las condiciones iniciales.
𝑥(0) = 0
𝑥̇ (0) = 1
𝑥̈ + 𝜔2 𝑥 + 𝛽𝑥 3 = 𝐹 sin 8𝑡 (2)
Donde 𝜔 y 𝛽 son constantes; con las condiciones iniciales
𝑥̇ (0) = 2
𝑥(0) = 0
3. MATERIALES.
Programa Matlab R2015a
Laptop HP bs019la
Especificaciones técnicas:
Procesador: Intel Core i7-7500U CPU @ 2.70GHz 2.90GHz
Memoria RAM: 8.00 GB
Tipo de sistema: Sistema operativo Windows 10 de 64 bits, procesador x64
Tarjeta gráfica: AMD Radeon 530 de 4 GB
4. PROCEDIMIENTO.
1º Iniciamos Simulink desde Matlab, haciendo click en Simulink Library.
2º Hacemos click en New Model y construimos los diagramas de bloques de las
ecuaciones 1 y 2 arrastrando los bloques desde la librería hasta la ventana de
trabajo (ver Ilustraciones 1 y 2).
Ilustración 1: Diagrama de bloques de la ecuación diferencial 1
Ilustración 2: Diagrama de bloques de la ecuación diferencial 2
3º Colocamos las respectivas condiciones iniciales de cada ecuación en los
Integrators haciendo doble click en ellos y colocando su valor en Initial
Condition, así también colocamos el valor de la amplitud y la frecuencia en la
función senoidal. En la Ilustración 2 se colocaron los valores de 𝜔 y 𝛽
(constantes) iguales a la unidad.
4º Posteriormente se analizaron las gráficas 𝑥 𝑣𝑠 𝑡, 𝑥̇ 𝑣𝑠 𝑡 y 𝑥̇ 𝑣𝑠 𝑥 para distintos
valores de las constantes.
5. DISCUSIÓN.
I. Ecuación diferencial 1.(𝑡 = 60𝑠)
En la Ilustración 3, 4 y 5 observamos que la partícula realiza un movimiento de
derecha a izquierda o en sentido horario, parte con una rapidez igual a la unidad
desde el origen dadas por las condiciones iniciales hasta alcanzar una rapidez de
1500 m/s en la posición x = 4800 m aproximadamente. La partícula pasa cuatro
veces por el origen. A partir de 𝑡 = 35 𝑠 la partícula comienza a ganar valores
muy altos en la rapidez y la posición.
Ilustración 3: Gráfica 𝑥̇ 𝑣𝑠 𝑥 para la ecuación 1
Ilustración 4: Gráfica 𝑥 𝑣𝑠 𝑡 para la ecuación diferencial 1
Ilustración 5: Gráfica 𝑥̇ 𝑣𝑠 𝑡 para la ecuación diferencial 1
II. Ecuación diferencial 2. (𝑡 = 20𝑠)
Ilustración 6: Gráfica 𝑥 𝑣𝑠 𝑡 para la ecuación 2
Ilustración 7: Gráfica 𝑥̇ 𝑣𝑠 𝑡 para la ecuación 2
Ilustración 8: Gráfica 𝑥̇ 𝑣𝑠 𝑥 para la ecuación 2
Hemos considerado los valores de las constantes 𝜔, 𝛽 y 𝐹 iguales a la unidad
para las Ilustraciones 6,7 y 8. Al disminuir el tiempo para el que se capturaron
las gráficas anteriores se puede observar que la partícula se mueve también hacia
la derecha y luego hacia la izquierda o en sentido horario, completando su
primera vuelta a los 3.9 s aproximadamente; la solución 𝑥(𝑡) es oscilatoria
alcanzando un desplazamiento máximo de 1.44 m aproximadamente, pero en sus
puntos máximos y mínimos no es suave, de igual forma la gráfica 𝑥̇ 𝑣𝑠 𝑡, se
alcanza una rapidez máxima de 2.08 m/s aproximadamente en los puntos
máximos y una rapidez mínima cercana a 2.11 m/s en los puntos mínimos.
Ahora variamos los valores de las constantes, pondremos 𝜔 = 0.1 y los otros
valores los dejamos iguales; la partícula alcanza un desplazamiento máximo
cercano a 1.716 m y una rapidez máxima de 2.25 m/s arriba y de 2.21 m/s abajo
y se tiene que la gráfica se deforma más en los máximos y mínimos; las curvas
en el espacio de fase (gráfica 𝑥̇ 𝑣𝑠 x) tienden a separarse.
Ilustración 9: Gráfica 𝑥̇ 𝑣𝑠 𝑥 de la ecuación 2 para 𝜔 = 0.1
Si disminuimos el valor de 𝐹, por ejemplo igual a 0.1 para 𝜔 = 0.1, se tiene más
suavidad en las tres gráficas, especialmente en la gráfica 𝑥̇ 𝑣𝑠 t, alcanzando una
rapidez máxima y mínima de 2 m/s.
Al hacer 𝛽 = 0.1 y los otros valores iguales a 0.1, se produce un mayor
desplazamiento, máximo y mínimo, igual a 3m; esto hace que la gráfica 𝑥̇ 𝑣𝑠 x
se haga más ancha.
Al incrementar todos los valores de tal forma que sean mayores a 1, se va
perdiendo la simetría en todas las gráficas, se adquiere un comportamiento
caótico que se puede “regular” si hacemos el valor de F menor a uno.
Ilustración 10: Gráfica 𝑥̇ 𝑣𝑠 𝑡 para 𝜔 = 4, 𝛽 = 2, 𝐹 = 5 de la ecuación 2
6. CONCLUSIONES.
1) Para valores de las constantes 𝜔, 𝛽 y 𝐹 menores o iguales a la unidad se observa
un movimiento periódico.
2) Un valor de 𝜔 mayor que uno provoca un menor desplazamiento en la partícula.
3) Valores de 𝜔 y 𝛽 menores que uno provocan un mayor desplazamiento en la
partícula.
4) Para valores de las constantes 𝜔, 𝛽 y 𝐹 mayores o iguales a la unidad se pierde
la simetría en las gráficas.
5) Se puede retomar la simetría si hacemos a 𝐹 tan pequeño como queramos.
7. REFERENCIAS.
i. [Link]
132199 visitado el 29 de Abril de 2019
ii. [Link] visitado el 29 de Abril de 2019
