# -*- 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')
示例#2
0
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')
示例#4
0
# -*- 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')
示例#5
0
# -*- 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")