Описание тега desolve

Описание тега Вопросы с тегом
deSolve — это пакет R с функциями для решения обыкновенных дифференциальных уравнений первого порядка (ODE), уравнений в частных производных (PDE), дифференциальных алгебраических уравнений (DAE) и дифференциальных уравнений с запаздыванием. Этот тег следует использовать для вопросов о том, как создавать или отлаживать статистические модели, выполняемые с помощью функций, принадлежащих deSolve.
1 ответ

Добавление оператора if then к начальному значению условия в системе ODE; deSolve

Я пытаюсь добавить оператор if then для определения начального значения одной из моих переменных состояния и использую deSolve. По сути, я хочу ввести 3-е ОДУ (в данном случае 3-й вид в популяцию) после начала моделирования. Вот как выглядит код без…
18 мар '21 в 06:44
0 ответов

Как я могу навязать начальные условия дискретному PDE при использовании deSolve?

Я пытаюсь решить систему дифференциальных уравнений с одним ОДУ и одним УЧП второго порядка. Я дискретизировал PDE, чтобы решить его с помощью deSolve ode.1D, но теперь я не могу понять, как включить мои граничные условия для пространственной произв…
13 апр '21 в 19:00
1 ответ

Как получить результаты с фиксированным шагом?

У меня есть система дифференциальных уравнений, она немного сложна, но самая близкая аналогия, к которой я могу прийти, - это каталитическая химическая реакция, при которой катализатор разлагается со временем, например A + B -> C + B (rate k1) B …
27 апр '21 в 12:17
2 ответа

Ошибка интеграции deSolve ODE, я использую неправильную функцию?

Я пытаюсь решить систему уравнений, связанных с биологическими процессами. Одно уравнение (около 5) предназначено для фармакокинетической (ФК) кривой формы C = Co(exp(k1*t)-exp(k2*t). Необходимо одновременно решить производную этого уравнения вместе…
30 апр '21 в 16:49
0 ответов

Пакет R deSolve: запускайте различные события, используя комбинацию функции и фрейма данных

Я использую функцию ODE в пакете deSolve. Я понимаю, что могу использовать аргумент «события» для запуска события в определенный момент времени с помощью фрейма данных или функции ( https://tpetzoldt.github.io/deSolve-forcing/deSolve-forcing.html#ev…
04 май '21 в 17:39
1 ответ

Различия между Rs deSolve и Pythons odeint

В настоящее время я изучаю систему Лоренца с помощью и и заметил тонкие различия в пакетах. odeint из и оба говорят, что используют их для расчета своих производных. Однако с помощью lsodaкоманда для обоих, кажется, дает совершенно разные результаты…
20 май '21 в 19:03
1 ответ

Установка начальных (состояний) значений для системы ODE в скомпилированной модели (deSolve, Rcpp)

Я борюсь с, вероятно, незначительной проблемой, вызывая скомпилированные ODE для решения с помощью пакета R, и я обращаюсь за советом к более опытным пользователям. Задний план Мне нужно решить пару систем ODE. Я определил ODE в отдельных функциях C…
28 май '21 в 15:14
1 ответ

Реалистичная возрастная структурированная модель с использованием ODE из пакета deSolve

Я пытаюсь смоделировать реалистичную структурированную по возрасту модель, в которой все люди могут перейти в следующую возрастную группу в конце временного шага (а не стареть непрерывно с заданной скоростью), используя ODE из пакета deSolve. Рассма…
20 июл '21 в 17:34
1 ответ

в модели R desolve SIR с петлей и если еще

У меня есть простая модель SIR, я пытаюсь реализовать подход к вакцинации (V), где сначала проверяется, если инфицированные выше порогового значения (100), и если остается достаточно восприимчивых людей (50), он будет вакцинировать определенное коли…
22 авг '21 в 06:22
1 ответ

lsoda создает вымышленные ценности?

В настоящее время я разрабатываю модель для интеграции микробиологии и геохимии, в которой lsoda используется для решения большого количества дифференциальных уравнений. Модель слишком велика, чтобы размещать ее здесь, потому что она состоит из неск…
01 сен '21 в 16:09
1 ответ

Цикл по параметрам для достижения равновесия с deSolve

С петлями борюсь интуитивно. У меня есть простая модель потребительских ресурсов, и я хочу прокрутить значения скорости роста ресурсов. gчтобы получить значения конечного состояния, чтобы затем построить график равновесия как функцию значений параме…
10 сен '21 в 23:03
2 ответа

R [Shiny]: Как сделать реактивные блестящие приложения, отображающие модели динамических систем?

Я хочу создать реактивное приложение Rshiny, которое отображает результаты динамической модели, которая решается пакетом deSolve. Код примера был скопирован из Джима Дагганса System Dynamics Modeling с R. Вот код без R-Shiny, это экономическая модел…
25 окт '21 в 23:06
2 ответа

Как решить систему ОДУ с параметрами, зависящими от времени в R?

Я пытаюсь решить эту систему ODE через deSolve, dX / dt = -X*a + (YX) b + c и dY / dt = -Ya + (XY)*b для времени [0,200], a= 0,30, b= 0,2, но c равно 1 для времени [50,70] и 0 в противном случае. Код, который я использовал, time <- seq(0, 200, by…
04 ноя '21 в 19:57
1 ответ

Дифференциальное уравнение задержки ошибки deSolve (dede)

Я пишу дифференциальное уравнение с задержкой в ​​deSolve (R) и получаю сообщение об ошибке. Я не знаю, как решить. Итак, для некоторой предыстории. У меня есть система с 12 дифференциальными уравнениями, и 3 из них имеют запаздывание. Мне удалось н…
17 ноя '21 в 13:56
1 ответ

deSolve: дифференциальные уравнения с двумя последовательными динамиками

Я моделирую кольцевую трубку с текущей водой и температурным градиентом, используя deSolve::ode(). Кольцо моделируется как вектор, в котором каждый элемент имеет значение температуры и положение. Я моделирую формулу рассеивания тепла: 1) Но я также …
03 дек '21 в 14:49
0 ответов

Смоделируйте температуру вдоль трубы с движущейся водой с помощью deSolve и ReacTran

Я моделирую кольцевую трубку с текущей водой и температурным градиентом, используя deSolve::ode(). Кольцо моделируется как вектор, в котором каждый элемент имеет значение температуры и положение. Я моделирую распространение тепла через: Которая испо…
06 дек '21 в 20:11
1 ответ

deSolve: Не могу понять, как досрочно остановить решатель од с корневыми функциями

Я не понимаю, как остановить решатель при выполнении определенного условия. Я подготовил фиктивную модель SIR, которая должна остановиться, когда I-отсек достигнет определенного значения. Но в моем коде решатель просто продолжает: library(deSolve) l…
16 дек '21 в 12:40
1 ответ

Возможная проблема с индексированием в ODE?

Следуя предложению @tpetzoldt, я открываю это как вопрос после предыдущего обсуждения (значения параметров как функция другого вектора. DeSolve ). Я пытаюсь достичь того, чтобы иметь возможность интегрировать модель на каждом временном шаге по векто…
21 дек '21 в 16:23
3 ответа

Значения параметров как функция другого вектора. deSolve

Я хочу создать модель динамики популяции, в которой каждое значение параметра соответствует температуре этого дня. например Простая модель library(deSolve) set.seed(1) pars <- c(alpha = 1, beta = 0.2, delta = 0.5, gamma = 0.2) lv_model <- func…
20 дек '21 в 18:43
1 ответ

Приток в определенное время с изменением временного шага в DeSolve

Я создаю модель распределения населения с ODE и DDE в DeSolve с помощью метода lsoda. В этой модели я хочу иметь приток населения в определенное время (в определенный день). Очень простой пример: dn1dt=influx - mortality Приток (x) должен произойти …
26 янв '22 в 13:17