Фон

Предположим, что у вас есть модель с одним отсеком с разновидностью под названием drug это представляет общую сумму препарата в теле. Препарат удален из тела через устранение первого порядка, представленное реакцией drug -> null, с константой скорости устранения ke. Другими словами, концентрация препарата по сравнению с профилем времени следует за моноэкспоненциальным снижением $$C_t=C_0{e}^{-k_{e}t}$$, где $$C_t$$ концентрация препарата во время t, $$C_0$$ начальная концентрация и ke константа скорости устранения. В этом примере показано, как настроить такую модель с одним отсеком и добавить капельное внутривенное введение на постоянном уровне 10 мг/час для суммы суммарной дозы 250 мг.

Создайте модель с одним отсеком

Создайте модель SimBiology под названием onecomp.

m1 = sbiomodel('onecomp');

Задайте устранение препарата от системы путем добавления реакции drug -> null к модели.

r1 = addreaction(m1,'drug -> null');

Разновидности drug автоматически создается и добавляется к отсеку. null разновидность является зарезервированной разновидностью, которая действует как приемник в этой реакции.

Добавьте массовую акцию кинетический закон в реакцию. Этот кинетический закон задает устранение препарата, чтобы следовать за кинетикой первого порядка.

k1 = addkineticlaw(r1,'MassAction');

Задайте параметр уровня устранения ke и добавьте его в кинетический закон.

p1 = addparameter(k1,'ke','Value',1.0,'ValueUnits','1/hour');

Задайте параметр уровня ke как параметр форвардного курса реакции путем установки ParameterVariableNames свойство кинетического объекта k1 закона. Это позволяет SimBiology определять скорость реакции для drug -> null реакция.

k1.ParameterVariableNames = 'ke';

Настройте Капельное внутривенное введение

Добавьте объект дозы в модель с помощью adddose метод. Задайте сумму дозы (Amount), цель дозы (TargetName), и скорость введения (Rate). Также необходимо установить Active свойство дозы возражает против true так, чтобы доза была применена к модели в процессе моделирования.

d1 = adddose(m1,'InfusionDose');
d1.Amount = 250;
d1.TargetName = 'drug';
d1.Rate = 10;
d1.RateUnits = 'milligram/hour';
d1.Active = true;

Симулируйте модель

Измените время остановки симуляции в 48 часов, чтобы видеть полный курс времени.

cs = getconfigset(m1);
cs.StopTime = 48;
cs.TimeUnits = 'hour';
sd = sbiosimulate(m1);

Постройте результаты

Постройте концентрацию по сравнению с профилем времени препарата в системе.

sbioplot(sd);