Вопросы по теме 'ode'
Возможная ошибка во взаимодействии odeint ‹-› interp1d?
Я относительно новичок в python и scipy, будучи конвертером из MATLAB. Я быстро тестировал функцию odeint в scipy.integrate и наткнулся на эту потенциальную ошибку. Рассмотрим следующий фрагмент:
from scipy import *
from scipy.integrate import...
611 просмотров
schedule
29.06.2023
Решение Maple ODE, графика не завершена. Почему?
Итак, мы пытаемся решить простую ОДУ. Мы используем текст book (предупреждение - русский) . Это код, представленный в качестве примера в книге (стр. 18):
> with(DEtools); with(plots);
> eq1 := diff(x(t), t) = y(t); eq2 := diff(y(t), t) =...
173 просмотров
schedule
15.12.2022
Использование boost::numeric::odeint внутри класса
Для моделирования я использую boost::numeric::odeint, но у меня возникла проблема. Я использую функцию интеграции внутри метода одного из моих классов, и у меня возникает ошибка «нет соответствующей функции для вызова интеграции». Чтобы было...
1672 просмотров
schedule
24.05.2023
Может ли параметр tspan = [to: very_small_step: tf] повлиять на размер шага решателя ode45?
Я знаю, что решатель ode45 имеет адаптивный размер шага, управляемый самой программой Matlab. Приведенное ниже описание предоставлено веб-сайтом Matlab :
Указание tspan с более чем двумя элементами не влияет на внутренние временные шаги,...
466 просмотров
schedule
03.09.2023
цикл for для значений параметров, зависящих от времени, в решателе ode работает только для некоторых значений t
Я использую простой цикл if для изменения значений параметров в моем сценарии ode. Вот пример сценария, который я написал, который демонстрирует ту же проблему. Итак, первая версия, которая работает:
function aah = al(t,x)
if (t>10000...
958 просмотров
schedule
29.05.2024
Могу ли я запустить ode45 назад во времени?
Я использую ode45 для решения дифференциальных уравнений, но tspan всегда должен быть таким, чтобы система работала вперед во времени. Могу ли я заставить ode45 запустить систему в обратном направлении до отрицательного t?
Конкретная проблема, с...
6480 просмотров
schedule
03.10.2023
Можно ли использовать анонимную функцию в качестве функции события при решении ОДУ в Matlab
Можно ли использовать анонимную функцию в качестве функции события в Matlab. То, что я хочу сделать, это в основном
opt = odeset('Events', @(t,y) (deal(y(end)-t^2,1,0)));
[T,Y] = ode45(@odefun,[tstart tend],y0,opt);
Однако это возвращает...
527 просмотров
schedule
26.09.2022
Задачи Mathematica с некоторыми дифференциальными уравнениями и DSolve
У меня проблемы с началом. Я учусь на программе финансового инжиниринга и пытаюсь использовать книгу, написанную в 2003 году, чтобы помочь мне смоделировать дифференциальные уравнения в частных производных, модель Блэка Шоулза и т. д.
Но во...
633 просмотров
schedule
16.12.2023
std::bad_alloc и решатель GSL ODE
Я пытаюсь решить (большую) систему ОДУ с помощью решателей GSL. Когда я использую метод драйвера, я получаю сообщение об ошибке could not allocate space for gsl_interp_accel , когда я определяю управление, ошибку и степпер вручную, я получаю...
297 просмотров
schedule
08.07.2023
Matlab: как использовать массив в функции dsolve?
У меня есть система ОДУ из двух уравнений, но я хочу минимизировать ее, используя только одно уравнение с результатом другого.
1)
t=linspace(0,2,3);
syms x(t) y(t);
inits='x(0)=2,y(0)=0';
[x,y]=dsolve('Dx=y','Dy=(y*2)-x', inits)
x =...
879 просмотров
schedule
27.10.2022
Matlab: нахождение коэффициентов системы ОДУ
У меня есть все данные и система ОДУ из трех уравнений с 9 неизвестными коэффициентами (a1, a2,..., a9).
dS/dt = a1*S+a2*D+a3*F
dD/dt = a4*S+a5*D+a6*F
dF/dt = a7*S+a8*D+a9*F
t = [1 2 3 4 5]
S = [17710 18445 20298 22369 24221]
D = [1357.33 1431.92...
2519 просмотров
schedule
02.04.2023
Matlab: возможно ли численно решить систему од со смесью начальных и конечных условий?
Я пытаюсь использовать ode45 для решения системы ОДУ:
[X,Y]= ode45(@sys,[0, T],y0);
куда,
function dy = sys(t,y)
dy(1) = f_1(y)
dy(2) = f_2(y)
dy(3) = f_3(y)
end
Проблема в том, что функция ode45...
1827 просмотров
schedule
12.04.2023
Обновите переменные ODE в Matlab
Я пытаюсь воспроизвести результаты статьи ( Immune Response ). Короче говоря, он смотрит на концентрацию патогенов, антител и плазматических клеток, а также на влияние иммунного ответа на орган. В зависимости от состояния органа будет меняться...
141 просмотров
schedule
29.12.2022
Случайное событие Matlab/Octave ode45
У меня есть некоторые проблемы с пониманием того, как реализовать события в октаве/матлабе при решении дифференциальных уравнений.
Рассмотрим, например, этот простой код для решения дифференциального уравнения y' = -y:
function dy =...
993 просмотров
schedule
25.03.2022
Абсолютная погрешность методов ODE45 и Рунге-Кутта по сравнению с аналитическим решением
Буду признателен, если кто-то поможет со следующей проблемой. У меня есть следующий ODE:
dr/dt = 4*exp(0.8*t) - 0.5*r ,r(0)=2, t[0,1] (1)
Я решил (1) двумя разными способами. С помощью метода Рунге-Кутты (4-й порядок) и с помощью...
6659 просмотров
schedule
20.12.2022
(Как) scipy.integrate.odeint ускоряет оценку функции?
Обычно чистый Python работает примерно в 50 раз медленнее, чем собственный код (C, Fortran), если он состоят из плотного цикла с простой арифметикой. Когда вы используете scipy.odeint , как описано для Например, в этом руководстве вы просто...
888 просмотров
schedule
07.11.2023
решение двух несвязанных ODE в цикле с использованием python и scipy.integrate.ode
У меня возникла проблема с решением двух очень простых несвязанных ODE с использованием scipy.integrate.ode. Например, следующий простой код:
from scipy.integrate import ode
def f(t, y, r_r=1.68,mu_ext=0. ,tau_m=0.020, tau_s=0.005, gs=...
674 просмотров
schedule
15.06.2022
Эмуляция MATLAB ode15s в Python
Я работаю над переводом модели из MATLAB в Python. Суть модели заключается в ode15 от MATLAB. В исполнении MATLAB ode15s имеет стандартные опции:
options = odeset()
[t P] = ode15s(@MODELfun, tspan, y0, options, params)
Для справки, y0...
3655 просмотров
schedule
11.03.2022
Вывод ступенчатой функции Хевисайда в качестве параметра для другой функции Matlab
Я пытаюсь реализовать математическую модель, описанную в этой ссылке:
Чад Лю, Чуан-Юань Ли, Фан Юань. Математическое моделирование восходящего пути Феникса. Plos One, Computational Biology.v.10, 2014.
Я хотел бы, чтобы выходная...
867 просмотров
schedule
13.05.2022
Остановить интеграцию по истечении заданного промежутка времени в Matlab
Я хочу прекратить решение дифференциального уравнения в Matlab, если оно занимает больше заданного времени. Я пробовал следующее, но это не работает...
options = odeset('AbsTol',1e-8,'RelTol',1e-5);
RUNTIME=5;
timerID = tic;
while...
1431 просмотров
schedule
02.03.2023