Я пытаюсь вычислить градиент функционала решения стохастического дифференциального уравнения (SDE) с учетом конкретной реализации шума. Я могу успешно вычислить эти градиенты, если оставлю шум неопределенным, как показано в DiffEqFlux.jl:...

using DifferentialEquations f(u,p,t) = 1.01*u u0 = 1/2 tspan = (0.0,1.0) prob = ODEProblem(f,u0,tspan) sol = solve(prob, Tsit5(), reltol=1e-8, abstol=1e-8) using Plots plot(sol,linewidth=5,title="Solution to the linear ODE with a thick...

В Julia я хочу использовать пакет DifferentialEquations.jl для решения \ ddot {u} + f (u, \ dot {u}, p) = g (t) где g (t) задается как вектор значений в эквидистантные моменты времени t. Эта ситуация отличается от ситуации в...

Я использую DifferentialEquations.jl для решения системы ODE, как показано ниже. Результат на самом деле не имеет значения, поскольку p содержит только тестовые параметры для целей создания MWE, но ключевым моментом является то, что я вижу много...

Недавно я начал знакомиться с пакетом DiffEqPhysics и DifferentialEquations от Джулии. Мне интересно, есть ли способ (например, функции обратного вызова) завершить решатель с определенным условием сходимости вывода. Например, если результат,...

Подобно этому вопросу , я пытаюсь решить это ODE с зависящим от времени входным параметром . Он состоит из серии дискретные обратные вызовы . В определенные моменты изменяется параметр (не состояние!). Время и значения хранятся в nx2 Array ....

Мне приходится решать систему дифференциальных уравнений много раз, перебирая параметр. Для этого я запускаю цикл по списку параметров и сохраняю решение (оцененное в массиве значений времени) для каждого параметра. Итак, у меня есть 2D-массив, в...