Я пытаюсь нарисовать модель SIR на python.Тем не менее, я не смог найти хорошее решение для рисования моей модели, которая приведена ниже:
dS / dt = -aSI… восприимчивый
dI / dt = aSI-bI… Зараженный
dR / dt = bI… выздоровел
Следует отметить, что третье уравнение не влияет на первые два.В результате его можно игнорировать.Используя только первые два уравнения, разработайте скрипт Python, который будет отображать фазовое пространство модели.Вы должны получить что-то вроде этого в качестве результата.
Изображение модели
Я уже выглядел так:
from pylab import *
g,L = 1,1
xvalues, yvalues = meshgrid(arange(-3, 3.1, 0.1), arange(-3, 3.1, 0.1))
vx = xvalues
vy = yvalues - xvalues
streamplot(xvalues, yvalues, vx, vy)
show()
, но результат:
Изображение результата