Вопросы по теме 'odeint'
ODEINT не может вызвать функцию
поэтому я пытаюсь решить систему из трех ODE и разработал следующий код для их решения с помощью ODEint. Но когда я бегу, ODEint не может вызвать мою функцию для системы уравнений.
from scipy.integrate import odeint
#initial parameters
X0 =...
249 просмотров
schedule
22.10.2021
C ++ - подходящее решение для системы значительно различается
Я использую odeint для решения системы из 4 связанных уравнений, которые имитируют вибрацию транспортного средства во время вождения. Я надеялся, что мои результаты будут похожи на те, что я получаю в MATLAB, но, к сожалению, этого не происходит. Я...
82 просмотров
schedule
02.11.2021
ODEINT с несколькими параметрами (зависит от времени)
Я пытаюсь решить одно ОДУ первого порядка с помощью ODEINT. Ниже приведен код. Я ожидаю получить 3 значения y для 3 временных точек. Проблема, с которой я борюсь, - это возможность передать n-е значение mt и nt для вычисления dydt . Я думаю,...
4828 просмотров
schedule
02.10.2021
Ошибка измерения в коде для решения связанного ODE первого порядка с odeint
Я написал этот код для системы ODE первого порядка, но, очевидно, python понимает y как целое число, когда я имел в виду список. Может кто увидит ошибку? Я просто не знаю, где это.
import matplotlib.pyplot as plt
from scipy.integrate import...
114 просмотров
schedule
21.02.2022
Использование встроенных функций Python для связанных ODE
ЭТА ЧАСТЬ - ПРОСТО ФОН, ЕСЛИ ВАМ НУЖНО
Я разрабатываю численный решатель для модели Курамото второго порядка. Функции, которые я использую для нахождения производных тэты и омеги, приведены ниже.
# n-dimensional change in omega
def...
132 просмотров
schedule
03.03.2022
odeint для переменного во времени возбуждения
Я хочу решить связанное дифференциальное уравнение.
Здесь A — возбуждение, а y(t) — реакция, где y=[y1,y2] . Я написал следующий код, который работает с постоянным значением одиночного возбуждения. Но это не удается, когда я ставлю...
20 просмотров
schedule
06.04.2022
Что мне нужно сделать, чтобы функция интеграции odeint компилировалась в другом классе?
Я пытался в течение нескольких дней попытаться скомпилировать функцию ODEINT Integrated_adaptive() внутри класса. Если я беру функцию вне класса и имею общие имена функций, то программа компилируется и работает так, как требуется. (out — это просто...
561 просмотров
schedule
22.05.2022
odeint сбрасывает объект во время итерации
Выполнение следующего кода:
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
class CSystem
{
private:
int counter=0;
public:
void operator() ( const double...
48 просмотров
schedule
28.05.2022
Параметрический график решения 2х2 разн. система на Python, Mathematica
Я реализовал решение следующей системы уравнений
dy/dt = -t*y(t) - x(t)
dx/dt = 2*x(t) - y(t)^3
y(0) = x(0) = 1.
0 <= t <= 20
сначала в Mathematica, а затем в Python.
Мой код в Mathematica:
s = NDSolve[
{x'[t] == -t*y[t] - x[t],...
483 просмотров
schedule
19.06.2022
Способы ускорить работу SciPy?
Я вызываю функцию, которая использует odeint при каждом проходе через цикл for (к сожалению, я не могу ничего сломать из этого цикла). Но дела идут намного медленнее, чем я надеялся. Вот код:
def get_STM(t_i, t_f, X_ref_i, dxdt, Amat):...
565 просмотров
schedule
21.06.2022
Как передать только один аргумент в odeint?
Я пытаюсь использовать scipy odeint для решения некоторых обыкновенных дифференциальных уравнений. Единственная проблема в том, что я хочу определить только один аргумент, и кажется, что для создания кортежа вам нужно как минимум два значения....
754 просмотров
schedule
19.06.2022
Почему нельзя вывести аргумент шаблона?
Я использую odeint из библиотеки boost и как новичок в C++ не понимаю, в чем проблема. Учитывая следующий исходный код в файле с именем ODESystem.hpp
#include "blaze/Math.h"
class ODESystem {
public:
using Complex = std::complex<...
183 просмотров
schedule
17.07.2022
Решение простого ODE с использованием scipy odeint дает прямую линию в 0
Я пытаюсь решить простую ОДУ: dN/dt = N*(rho(t)-бета)/лямбда
Rho — это функция времени, и я сгенерировал ее с помощью linspace. Код работает для других уравнений, но каким-то образом дает ровную прямую линию в точке 0. (Вы можете видеть это на...
171 просмотров
schedule
27.07.2022
OpenMP с ODEINT в функции ODE
Я пытаюсь внутренне распараллелить функцию ODE, которая интегрирована с ODEINT.
Я сделал следующий небольшой пример
#include <iostream>
#include <chrono>
#include <Eigen/Dense>
#include <omp.h>
#include...
62 просмотров
schedule
23.07.2022
Передать массив в Python odeint
Я новичок в Python, поэтому извините меня, если на следующий вопрос есть ответ «дух».
Итак, я пытаюсь решить ODE с помощью odeint и хочу передать массив. Но TypeError: не может умножать последовательность на не-int типа 'float', продолжает...
2757 просмотров
schedule
20.09.2022
python: прерывание odeint при выполнении условия
Я использую функцию odeint из пакета scipy.integrate:
r0 = np.array([1,2,3,4])
t=np.linspace(0,1,20)
def drdt(r,t):
return r # or whatever else
r = odeint(drdt,r0,t)
r0 — это массив numpy, который содержит начальные позиции определенного...
972 просмотров
schedule
21.12.2022
Возникли проблемы с ODEINT в python
Я относительно новичок в Python и пытаюсь использовать его для решения нелинейного дифференциального уравнения второго порядка, в частности уравнения Пуассона-Больцмана в электролите.
phi''(r) + (2/r)*phi'(r) = (k^2)*sinh(phi(r))
По сути, он...
1312 просмотров
schedule
16.04.2023
Использование dopri5 с библиотекой odeint boost
система уравнений
Привет. Я хочу развить эти уравнения во времени от нуля до 10 ^ 16 и начальных условий x (0) = 10 ^ 8 и y (0) = 0,5. Из-за зависимости уравнений от x в знаменателе я думаю, что использование odeint с runge_kutta_dopri5...
264 просмотров
schedule
13.06.2023
Оптимизация написания уравнений для odeint
Я написал следующий код для простой химической сети, которую должен решить odeint :
def chemnet(y,t):
assoc=0.1,0.001,1
oxy=0.001,0.1,0.001
f=zeros(6,float)
f[0]=...
47 просмотров
schedule
13.05.2023
Решение системы ОДУ в обратном направлении во времени в С++
К сожалению, я заметил, что одеинт BOOST не может решить систему ОДУ назад во времени, т.е. когда я изменяю условия так, чтобы
typedef std::vector< double > state_type;
void ode_function(const state_type &x, state_type &dxdt, const...
301 просмотров
schedule
10.03.2023