Создайте объект модели, и затем добавьте объект реакции.

modelObj = sbiomodel('my_model');
reactionObj = addreaction(modelObj, 'a -> c + d');

Установите Reversible свойство для reactionObj к true и проверьте эту установку.

set (reactionObj, 'Reversible', true)
get (reactionObj, 'Reversible')

MATLAB^® возвращается:

ans =

     1

MATLAB возвращает 1 для true и 0 для false.

На следующих шагах пример иллюстрирует, как уравнение скорости реакции присвоено для обратимых реакций.

Создайте кинетический объект закона для объекта реакции типа 'MassAction'.

kineticlawObj = addkineticlaw(reactionObj, 'MassAction');

reactionObj KineticLaw свойство сконфигурировано к kineticlawObj.

'MassAction' кинетический закон для обратимых реакций имеет две переменные параметра ('Forward Rate Parameter' и 'Reverse Rate Parameter') то, что необходимо установить. Переменные разновидностей для MassAction автоматически определяются. Чтобы установить переменные параметра, сначала создайте переменные параметра как объекты параметра (parameterObj1, parameterObj2) названный Kf и Kr и присвойте объект kineticlawObj.

parameterObj1 = addparameter(kineticlawObj, 'Kf');
parameterObj2 = addparameter(kineticlawObj, 'Kr');

Установите имена переменных для кинетического объекта закона.

set(kineticlawObj,'ParameterVariableNames', {'Kf' 'Kr'});

Проверьте, что скорость реакции описывается правильно в объекте ReactionRate реакции свойство.

get (reactionObj, 'ReactionRate')

MATLAB возвращается:

ans =

Kf*a*b - Kr*c*d