这个例子展示了如何用两种物质(a和B)和一个反应构建一个简单的模型。的反应是A - B >
,它遵循质量作用动力学和正向速率参数k
.因此变化率是
.
创建一个名为SimBiology的模型simpleModel
.
m1 = sbiomodel (“simpleModel”);
加入一个包含两种物质的反应一个
和B
,在那里一个
被转换为B
.
r1 = addreaction (m1,“A - B >”);
SimBiology会自动添加物种一个
和B
到模型中。
m1.species
ans = SimBiology Species Array Index:隔间:名称:值:单位:1未命名的A 0 2未命名的B 0
设置第一个物种的初始数量(一个
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);
sbiomodel
|addreaction(模型)
|kineticlaw addparameter(模型)
|addkineticlaw(反应)
|sbiosimulate