Вопросы по теме 'odeint'

Повторяющаяся дельта-функция в odeint
У меня есть следующая программа odeint: #include <boost/array.hpp> #include <boost/numeric/odeint.hpp> using namespace std; typedef boost::array< double , 1 > state_type; void eqsystem(const state_type &x, state_type...
1585 просмотров
schedule 03.04.2023

Оптимизация констант в дифференциальных уравнениях в Python
Итак, как мне подойти к написанию кода для оптимизации констант a и b в дифференциальном уравнении, например dy/dt = a*y^2 + b, с помощью curve_fit? Я бы использовал odeint для решения ODE, а затем curve_fit для оптимизации a и b. Если бы вы могли,...
3987 просмотров
schedule 26.04.2022

Есть ли аналог lsode в Boost OdeInt?
Существует буст 1.53 и OdeInt в нем. Он имеет диапазон числовые решатели оды в нем. Тем не менее мне интересно, есть ли аналог lsode в это - для случаев, когда у вас есть гарантированный аналитический раствор?
440 просмотров
schedule 07.07.2023

Повышение класса odeint с производной и якобианом
Я намерен использовать библиотеку Boost odeint в подпрограмме MCMC для оценки параметров в модели ODE. Поскольку эти ОДУ могут быть жесткими, мне нужно иметь возможность передать якобиан в решатель с производной. Я хотел бы создать класс, который...
1209 просмотров
schedule 15.06.2022

DDE с использованием boost odeint
Можно ли решить дифференциальные уравнения с запаздыванием, используя библиотеку C++ Boost-odeint? Для примера ниже уравнения: x'(t) = r*x(t)*(1 - x(t-tau)), где tau — постоянная величина временной задержки.
645 просмотров
schedule 11.04.2023

Нелинейные изменяющиеся во времени системные входы с boost :: odeint
Я работаю над бесконечномерными алгоритмами оптимизации с использованием методов оптимального управления (генерация траектории и оптимизация). Системы, в которых я хотел бы провести эту оптимизацию траектории, являются нелинейными и имеют вид $ \ dot...
182 просмотров
schedule 06.09.2022

С++ odeint выводит не число
Я работаю с C++ boost/odeint для интеграции сетей химических реакций. Все работает нормально, за исключением случаев, когда я устраняю некоторые реакции. Затем иногда в выводе есть -nan. вот небольшая часть этого вывода. первый столбец - время, а...
483 просмотров
schedule 22.07.2023

решатель оды с фиксированным размером шага
Я хотел бы знать, фиксирован ли размер шага в odeint или нет. В степпер Основная концепция шагового двигателя. Базовый степпер, следующий этой концепции степпера, может выполнить один шаг решения x(t) ОДУ, чтобы получить x(t+dt),...
352 просмотров
schedule 03.06.2023

Eigen Vector3d как вектор состояния с symplectic_rkn_sb3a_mclachlan odeint
Я реализую симуляцию n тел, определяя отдельные «частицы» с переменными (частица к частице, не зависящими от времени) свойствами, которые влияют на динамику, а затем определяя «Систему» ​​этих частиц как вектор, который определяет различные операции...
155 просмотров
schedule 03.02.2023

Как заставить SciPy.integrate.odeint останавливаться, когда путь закрыт?
edit: Прошло пять лет, SciPy.integrate.odeint еще не научился останавливаться? Приведенный ниже сценарий объединяет линии магнитного поля вокруг замкнутых путей и останавливается, когда оно возвращается к исходному значению в пределах...
1696 просмотров
schedule 24.04.2022

python — scipy.integrate.odeint возвращает неверные результаты
Я пытался интегрировать прямоугольную волну, используя python 3.5 и функцию scipy.integrate.odeint , но результаты не имеют никакого смысла и сильно различаются в зависимости от выбранного массива временных точек. Прямоугольная волна имеет период...
659 просмотров
schedule 03.11.2022

Использование Eigen::MatrixXd в качестве состояния в степпере ошибок odeint
Я новичок здесь, поэтому я надеюсь задать свой вопрос должным образом! :) Я пытаюсь реализовать решатель для диссипативного уравнения Линдблада с помощью odeint на С++, используя библиотеку Eigen для обработки матриц. С постоянным размером...
299 просмотров
schedule 20.07.2022

Установка odeint и boost
Я действительно неопытен в программировании. У меня проблемы с установкой odeint (и boost) для решения системы дифференциальных уравнений. Я просмотрел много вопросов об установке boost, перепробовал много из них, но пока не смог использовать...
1096 просмотров
schedule 10.06.2022

Печать только окончательного вывода odeint
Извините, это может показаться довольно глупым вопросом, но мне нужно спросить, можно ли напечатать только конечное выходное значение при решении связанных дифференциальных уравнений в odeint? На самом деле я пытаюсь решить два связанных...
402 просмотров
schedule 07.10.2023

Проблема с получением odeINT для принятия уравнений, сгенерированных циклом
Я хочу использовать odeINT Python для интеграции нескольких наборов уравнений, которые генерируются циклом. Все уравнения связаны между собой, поэтому их необходимо интегрировать одновременно с помощью одного вызова odeINT. Проблема в том, что...
96 просмотров

boost odeint слишком много аргументов для работы
Я использую boost/odeint для решения оды первого порядка. Определение системы в raman_system.h выглядит следующим образом: class raman_signal_system { double alpha_s; double dB_to_neper; double g_pump_b; double P_p; public:...
333 просмотров
schedule 02.11.2022

Как сделать два ползунка в matplotlib
Я хотел бы сделать два ползунка в matplotlib, чтобы вручную изменять значения N и P в моей модели хищник-жертва: import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint def lotka(x,t,params): N, P = x alpha,...
874 просмотров
schedule 18.02.2024

Python: переполнение встречается в двойных скалярах
У меня есть система из 5 дифференциальных уравнений, которые я решаю с помощью odeint, и я столкнулся с ошибкой деления на ноль в одном из моих дифференциальных уравнений. Чтобы обойти это, я попробовал два метода: (а) один создает функцию и...
778 просмотров
schedule 31.05.2023

Как решить модель пространства состояний с помощью Odeint?
Я пытаюсь реализовать численное моделирование модели пространства состояний с помощью Eigen и Odeint. Моя проблема в том, что мне нужно ссылаться на данные управления U (предопределено перед интеграцией) , чтобы правильно решить часть Ax + Bu...
603 просмотров
schedule 28.04.2022

python odeint с кусочной функцией
У меня возникли проблемы с использованием odeint для решения поштучных од, когда ограничения по времени очень велики. Я включил здесь минималистский пример: import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint t0...
398 просмотров
schedule 30.09.2022