В этом примере показано, как создать простую модель с двумя видами (A и B) и реакцией. Реакция следующая A -> B
, который следует кинетике действующих масс с параметром прямой скорости k
. Следовательно, скорость изменения .
Создайте модель SimBiology с именем simpleModel
.
m1 = sbiomodel('simpleModel');
Добавьте реакцию, которая включает два вида A
и B
, где A
преобразуется в B
.
r1 = addreaction(m1,'A -> B');
SimBiology автоматически добавляет виды A
и B
в модель.
m1.species
ans = SimBiology Species Array Index: Compartment: Name: Value: Units: 1 unnamed A 0 2 unnamed B 0
Установите начальное количество первого вида (A
) к 10.
m1.species(1).InitialAmount = 10;
Определите кинетический закон реакции, чтобы следовать кинетике действующих масс. Вы можете добиться этого, добавив объект кинетического закона к реакции r1
.
kineticLaw = addkineticlaw(r1,'MassAction');
Добавьте параметр константы скорости в кинетический закон действия масс. Вы должны задать ParameterVariableNames
свойство объекта кинетического закона соответствует имени параметра 'k'
так, чтобы можно было определить скорость реакции.
p1 = addparameter(kineticLaw,'k',0.5); kineticLaw.ParameterVariableNames = 'k';
Симулируйте модель.
sd = sbiosimulate(m1);
Постройте график результатов симуляции.
sbioplot(sd);