# -*- coding: utf8 -*- from metodos import eq_diferenciales from metodos import transformada_fourier import math import graficas #definición de las funciones del sistema de ecuaciones diferenciales def f(x, y, t): return y def g(x, y, t): return 0.516 * math.cos(1.2 * t) - x - 0.52 * y - 1.09 * (x * x) h = 0.01 # configuración del paso de integración num_periodos = 8 #periodos de la señal visualizados x_0 = 0.0 # condición inicial x(t_0) = 0 y_0 = 0.0 # condición inicial dx/dt(t_0) = 0 t_0 = 0.0 # condición inicial t_0 = 0 #cálculo y representación gráfica de la señal y la espectro de potencias graficas.mostrar_graficas(h, num_periodos, f, g, x_0, y_0, t_0, 'oscilador no lineal forzado y amortiguado - S1', freq=1.2, fichero='memoria/caso_anarmonico_s1.pdf')
import math import graficas # definición de las funciones del sistema de ecuaciones diferenciales def f(x, y, t): return y def g(x, y, t): return -x h = 0.01 # configuración del paso de integración num_periodos = 5 # periodos de la señal visualizados x_0 = 0.0 # condición inicial x(t_0) = 0 y_0 = 0.0 # condición inicial dx/dt(t_0) = 0 t_0 = 0.0 # condición inicial t_0 = 0 # cálculo y representación gráfica de la señal y la espectro de potencias graficas.mostrar_graficas( h, num_periodos, f, g, x_0, y_0, t_0, "oscilador lineal no forzado ni amortiguado", fichero="memoria/caso_lineal.pdf", )
# -*- coding: utf8 -*- from metodos import eq_diferenciales from metodos import transformada_fourier import math import graficas def f(x, y, t): return y def g(x, y, t): return -x - 0.52 * y h = 0.01 num_periodos = 5 x_0 = 1.0 y_0 = 0.0 t_0 = 0.0 graficas.mostrar_graficas(h, num_periodos, f, g, x_0, y_0, t_0, 'oscilador lineal no forzado pero amortiguado')
# -*- coding: utf8 -*- from metodos import eq_diferenciales from metodos import transformada_fourier import math import graficas #definición de las funciones del sistema de ecuaciones diferenciales def f(x, y, t): return y def g(x, y, t): return math.cos(t) - x h = 0.01 # configuración del paso de integración num_periodos = 5 #periodos de la señal visualizados x_0 = 0.0 # condición inicial x(t_0) = 0 y_0 = 0.0 # condición inicial dx/dt(t_0) = 0 t_0 = 0.0 # condición inicial t_0 = 0 #cálculo y representación gráfica de la señal y la espectro de potencias graficas.mostrar_graficas(h, num_periodos, f, g, x_0, y_0, t_0, 'oscilador lineal forzado no amortiguado', fichero='memoria/caso_forzado.pdf')
# -*- coding: utf8 -*- from metodos import eq_diferenciales from metodos import transformada_fourier import math import graficas def f(x, y, t): return y def g(x, y, t): return -x h = 0.01 num_periodos = 5 x_0 = 1.0 y_0 = 0.0 t_0 = 0.0 graficas.mostrar_graficas(h, num_periodos, f, g, x_0, y_0, t_0, "oscilador lineal no forzado ni amortiguado")