Вопросы по теме 'odeint'
Объект слишком глубокий для желаемого массива - scipy.integrate.odeint
Я только вчера начал работать с Python и получаю сообщение об ошибке при использовании scipy.integrate.odeint .
Я определил функцию
def SIR(x, t, beta, gamma, mu, M):
который принимает numpy.array объекты x , t и M ; а скалярные...
7305 просмотров
schedule
03.04.2022
Интеграция ODE на GPU с помощью boost и python
Не так давно я писал здесь о модели, которую я пытаюсь построить с помощью pycuda, которая решает около 9000 связанных ODE. Однако моя модель слишком медленная, и член SO предположил, что, вероятно, виновата передача памяти с хоста на GPU.
Прямо...
1502 просмотров
schedule
21.08.2022
Двухточечная граница с одеинтом
Я пытаюсь решить проблему двухточечной границы с помощью odeint. Мое уравнение имеет вид
y'' + a*y' + b*y + c = 0
Это довольно тривиально, когда у меня граничные условия y(x_1) = y_1 , y'(x_2) = y_2 , но когда граничные условия y(x_1) =...
724 просмотров
schedule
13.10.2023
odeint: Как мне регистрировать промежуточные результаты при интеграции?
Я хочу знать, как я могу регистрировать значения, отличные от состояний, во время интеграции с помощью odeint. У меня есть симуляция динамики спутника, которая описывается как дифференциальные уравнения полного углового момента L и количества...
165 просмотров
schedule
25.03.2022
Не могу скомпилировать С++, который использует odeint из boost
Я использую Ubuntu 12.04, и в /usr/include уже есть некоторые файлы повышения. я сделал
sudo apt-get install libboost-all-dev
и это тоже установило много файлов. Я не хочу удалять это ускорение и устанавливать из исходного кода, так как...
1042 просмотров
schedule
30.03.2022
odeint (С++) - даунсэмплинг наблюдений
Извините, если это простой вопрос, но есть ли «лучшая практика» для понижения частоты эволюции переменных состояния в odeint?
Ниже я скопировал хороший пример создания «наблюдателя» для регистрации переменных состояния, представленных в этой статье...
185 просмотров
schedule
18.11.2022
Как построить собственные значения при решении связанных матричных дифференциальных уравнений в PYTHON?
Допустим, у нас есть три комплексные матрицы и система связанных дифференциальных уравнений с этими матрицами.
import numpy, scipy
from numpy import (real,imag,matrix,linspace,array)
from scipy.integrate import odeint
import matplotlib.pyplot as...
3137 просмотров
schedule
24.07.2022
Сравнение runge_kutta4 odeint с ode45 Matlab
Я хочу использовать метод runge_kutta4 в библиотеке odeint C++ . Я решил проблему в Matlab. Мой следующий код в Matlab для решения x'' = -x - g*x' с начальными значениями x1 = 1 , x2 = 0 выглядит следующим образом
main.m
clear all...
3963 просмотров
schedule
12.04.2022
Проблемы с функцией и odeint в Python
В течение нескольких месяцев я начал работать с python, учитывая его огромные преимущества. Но недавно я использовал odeint от scipy для решения системы дифференциальных уравнений. Но в процессе интегрирования реализованная функция работает не так,...
1166 просмотров
schedule
23.12.2022
ODEint: адаптивная интеграция с произвольной точностью
Может ли ODEint использовать адаптивные процедуры интеграции с произвольной арифметикой точности? Например, я хотел бы использовать библиотеки мультиточности Boost с функцией Integrated_adaptive() с управляемым степпером. В документации ODEint...
548 просмотров
schedule
03.04.2023
Расширение подписи типа System при использовании `boost::numeric::odeint`
Система оды $f'(x,t) = f(x,t)$ должна иметь следующую подпись, как указано здесь
void sys( const state_type & /*x*/ , state_type & /*dxdt*/ , const double /*t*/ )
{
// ...
}
Его можно изменить следующим образом
void sys(...
57 просмотров
schedule
31.01.2023
Odeint непостоянен в течение нескольких периодов, моделируя управляемый маятник
Кажется, есть много вопросов, касающихся моделирования маятника и одеинта. Я считаю, что этот вопрос достаточно конкретен, чтобы стоять сам по себе. Он связан с передачей массива времени в odeint.
Я моделирую ведомый демпфированный маятник. Я...
108 просмотров
schedule
13.07.2023
python - разная длина массива по оси интерполяции?
Я пытаюсь использовать функцию интерполяции Python, чтобы получить значение y для данного x, но я получаю сообщение об ошибке «raise ValueError («массивы x и y должны быть равны по длине вдоль оси интерполяции», хотя мои массивы имеют равные размер и...
1186 просмотров
schedule
29.04.2023
Сохранение выходных данных решателя жестких уравнений odeint
Мне было интересно, есть ли способ сохранить вывод следующего код в матрицу повышения или вектор.
size_t num_of_steps = integrate_const( make_dense_output< rosenbrock4< double > >( 1.0e-6 , 1.0e-6 ) ,
make_pair(...
132 просмотров
schedule
28.02.2024
не может получить те же результаты, что и Matlab
Взгляните на следующую передаточную функцию:
С Matlab Simulink:
Результат
В представлении в пространстве состояний система может быть смоделирована следующим образом:
В Matlab мы можем смоделировать систему в...
88 просмотров
schedule
17.01.2023
Использование Eigen::VectorXd (Eigen 3.3.4) в качестве типа состояния в boost::numeric::odeint (Boost 1.65.1)
Во время моей работы мне требовалось использовать Eigen::VectorXcd в качестве типа состояния для решения огромной линейной системы ОДУ. В этом проекте матрица в системе ODE является разреженной. Умножение it на плотный вектор может быть вычислено...
697 просмотров
schedule
02.01.2023
Функция odeint Python, похоже, не работает
Я хочу изучить движение заряженной частицы во время путешествия через магнитное поле, моделируя его с помощью Python. Я попытался использовать функцию odeint из scipy.integrate, и, похоже, она работает не так, как я ожидал. Вот что я ожидал,...
335 просмотров
schedule
01.05.2022
odeint: переменный временной шаг для математического моделирования
В качестве примера я использую зомби-код Sci Py.
Проблема с ODEint заключается в том, что я интегрирую в течение длительного времени (в данном случае от 0 до 500, но для других проектов до тысяч лет), и он выполняет эту интеграцию за один раз. Я...
167 просмотров
schedule
24.03.2023
Как избежать расходящихся решений с помощью odeint? способ съемки
Я пытаюсь решить уравнение в Python. В основном, что я хочу сделать, это решить уравнение:
(1/x^2)*d(Gam*dL/dx)/dx)+(a^2*x^2/Gam-(m^2))*L=0
Это уравнение Клейна-Гордона для массивного скалярного поля в пространстве-времени Шварцшильда....
210 просмотров
schedule
06.06.2024
Система двух многомерных связанных ОДУ
Я пытаюсь решить следующую проблему связанных ОДУ, используя odeint() из scipy. Система выглядит так:
X'_k = среднее (Y_k) + F
Y'_{k,j} = X_k - Y_{k,j}
Это система с 3 переменными X, и для каждой переменной X есть еще 3...
94 просмотров
schedule
25.03.2023