Карта Пуанкаре в питоне - осциллятор - PullRequest
Новая
       17

Карта Пуанкаре в питоне - осциллятор

0 голосов
/ 04 ноября 2019

Я хочу построить карту Пуанкаре для осциллятора. Я построил кривую t vs x, где x - решение данного уравнения осциллятора. Чтобы построить график Пуанкаре, мне нужно сохранить те значения х, которые появляются в интервале 2 * пи. Кто-нибудь может мне с этим помочь? Я использовал Odeint Scipy, чтобы решить это. 

from scipy.integrate import odeint
import matplotlib.pyplot as plt

def mod(r,t):
    w0 = -1
    b = 1
    a = 0.5
    w = 1
    f = 0.33

    x = r[0]
    dx = r[1]
    rdot = [[],[]]
    rdot[0] = dx
    rdot[1] = f*np.sin(w*t) - a*dx - x*w0 - b*x**3
    return rdot

r0 = [0,-1]
t = np.arange(0,20,2*np.pi)


sol = odeint(mod,r0,t)
...