Articles of Ode

Graficando ODEs, Isoclines usando Python

Estoy buscando un paquete de Python que me permita trazar algo similar al applet de Java que se ve a continuación: http://math.mit.edu/mathlets/mathlets/isoclines/ ¿Alguien sabe algún paquete de trazado ODE para esto? Puedo codificar algo desde cero usando Numpy, Matplotlib, pero primero quería preguntar. Gracias,

codigo python para multiples odas

Quiero escribir un cuarto orden de Adams Bashforth para resolver el sistema. Y lo siguiente es lo que tengo: El sistema se encuentra en el siguiente enlace: Sistema que tenemos. def AdamsBashforth4( f, x0, t ): “”” Fourth-order Adams-Bashforth method:: u[n+1] = u[n] + dt/24.*(55.*f(u[n], t[n]) – 59*f(u[n-1], t[n-1]) + 37*f(u[n-2], t[n-2]) – 9*f(u[n-3], t[n-3])) […]

resolviendo dos ecuaciones diferenciales de dimensión en python con scipy

Soy un novato en Python. Tengo un sistema diferencial simple, que consta de dos variables y dos ecuaciones diferenciales y condiciones iniciales x0=1, y0=2 : dx/dt=6*y dy/dt=(2t-3x)/4y Ahora estoy tratando de resolver estas dos ecuaciones diferenciales y elijo odeint . Aquí está mi código: import matplotlib.pyplot as pl import numpy as np from scipy.integrate import […]

sistemas ODE complejos en scipy

Estoy teniendo problemas para analizar la ecuación del bloch óptico, que es un sistema ODE de primer orden con valores complejos. He encontrado que Scipy puede resolver este tipo de sistema, pero su página web ofrece muy poca información y casi no puedo entenderla. Tengo 8 EDO de primer orden acopladas, y debo generar una […]

Usar scipy.integrate.complex_ode en lugar de scipy.integrate.ode

Estoy tratando de usar el método complex_ode en lugar del método ode en scipy.integrate. La página de ayuda para complex_ode no proporciona un ejemplo, así que puedo haber hecho algo mal. Este código funciona correctamente con scipy.integrate.ode: from scipy.integrate import ode y0, t0 = [1.0j, 2.0], 0 def f(t, y, arg1): return [1j*arg1*y[0] + y[1], […]

¿Funciona scipy.integrate.ode.set_solout?

La interfaz scipy.integrate.ode a las rutinas de integración proporciona un método para detener la integración si se viola una restricción en cualquier paso, set_solout . Sin embargo, no puedo hacer que este método funcione, incluso en los ejemplos más simples. Aquí hay un bash: import numpy as np from scipy.integrate import ode def f(t, y): […]

Odd SciPy ODE Integration error

Estoy implementando un modelo muy simple de recuperación de la infección infectada con una población estable para un proyecto del lado ocioso, normalmente una tarea bastante trivial. Pero me estoy encontrando con errores de solucionador usando PysCeS o SciPy, los cuales usan lsoda como su solucionador subyacente. Esto solo sucede con los valores particulares de […]

Pase args para solve_ivp (nueva API de SciPy ODE)

Para resolver EDOs simples usando SciPy, solía usar la función odeint, con la forma: scipy.integrate.odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, hmin=0.0, ixpr=0, mxstep=0, mxhnil=0, mxordn=12, mxords=5, printmessg=0)[source] donde una función simple para integrarse podría incluir argumentos adicionales de la forma: def dy_dt(t, y, arg1, arg2): # processing […]

Optimizar constantes en ecuaciones diferenciales en Python.

Bueno, ¿cómo me acercaría a escribir un código para optimizar las constantes a y b en una ecuación diferencial, como dy / dt = a * y ^ 2 + b, usando curve_fit? Estaría usando odeint para resolver la EDO y luego curve_fit para optimizar a y b. Si pudiera darnos su opinión sobre esta […]

Usando tamaños de pasos adaptables con scipy.integrate.ode

La documentación (breve) de scipy.integrate.ode dice que dos métodos ( dopri5 y dop853 ) tienen control de tamaño grande y salida densa. Mirando los ejemplos y el código en sí, solo puedo ver una forma muy sencilla de obtener resultados de un integrador. Es decir, parece que simplemente hace que el integrador avance un cierto […]