Проблема с начальными значениями odeint

@Koala1996 · Регистрация: 04.12.2019

Студворк — интернет-сервис помощи студентам

Система ДУ подразумевает описание рост клеток от 0 и до 150 дней. И должна соответствовать графику (привожу пример во вложении). Проблема заключается в том, что начальные значение x1 т.е ([s0]) должен прибавлять +10. У меня такое чувство, что odeint не принимает +10 как константу. Помогите решить проблему.

Начальные значения x1, x2

Python
1
2
3
4
5
a1 = 3; a2= 4; g11 = 0.5; g12 = 1; g21 = -0.5; g22 = 0; b1 = 0.2; b2 = 0.02
gamma = g12*g21 - (1-g11)*(1-g22)
x1 = pow((b1/a1),((1-g22)/gamma))*pow((b2/a2),(g21/gamma))
x2 = pow((b1/a1),(g21/gamma))*pow((b2/a2),((1-g11)/gamma))
print (gamma, x1, x2)

Диф.уравнение

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import odeint 
 
def f(s,t):
    a1 = 3; a2= 4; g11 = 0.5; g12 = 1; g21 = -0.5; g22 = 0; b1 = 0.2; b2 = 0.02
    dx1dt = a1*s[0]**g11*s[1]**g21-b1*s[0]
    dx2dt = a2*s[0]**g12*s[1]**g22-b2*s[1]
    return [dx1dt,dx2dt]
s0 = [1.0606601717798212 + 10, 3.651483716701107]
t = np.linspace(0, 150, 100)                          
s = odeint(f, s0, t)
plt.plot(t,s[:, 0],'r-',linewidth=2.0,label='s1')
plt.plot(t,s[:, 1],'g-',linewidth=2.0,label='s2')
plt.xlabel("t")
plt.ylabel("s")
plt.legend()
plt.grid()
plt.show()

@u235 · 24.01.2022, 22:13

Сообщение от Koala1996

У меня такое чувство, что odeint не принимает +10 как константу

Все принимает, проверил.

Python
1
t = np.linspace(0, 150, 100)

лучше заменить на

Python
1
t = np.linspace(0, 150, 1000)

не будет изломов на графике.

@sabbra_caddabra · 30.01.2022, 23:02

Возможно вам будет полезно. Сами разработчики scipy рекомендуют использовать эту функцию вместо odeint https://docs.scipy.org/doc/sci... e_ivp.html

Новые блоги и статьи Все статьи Все блоги /
Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	Киев стоит - украинская песня zorxor 28.01.2026 wfWdiRqdTxc О Господи, Вечный, Ты . . . Я помоги, Бесконечный. . . Я прошу Ты. . . Я погибаю, спаси. . . Я прошу Тебя Вечный. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .	Загрузка PNG с альфа-каналом на SDL3 для Android: с помощью SDL3_image 8Observer8 27.01.2026 Содержание блога SDL3_image - это библиотека для загрузки и работы с изображениями. Эта пошаговая инструкция покажет, как загрузить и вывести на экран смартфона картинку с альфа-каналом, то есть с. . .
влияние грибов на сукцессию anaschu 26.01.2026 Бифуркационные изменения массы гриба происходят тогда, когда мы уменьшаем массу компоста в 10 раз, а скорость прироста биомассы уменьшаем в три раза. Скорость прироста биомассы может уменьшаться за. . .	Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .	Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .

@sabbra_caddabra 1 / 1 / 0 Регистрация: 14.10.2021 Сообщений: 4
	30.01.2022, 23:02
	Возможно вам будет полезно. Сами разработчики scipy рекомендуют использовать эту функцию вместо odeint https://docs.scipy.org/doc/sci... e_ivp.html 1