嗨,在床上,
最简单的方法来跟踪一个物种在一个模拟的最大可能是使用一个外部函数(见附件)与一个持久变量存储当前最大的值:
函数= getmax(时间、年代)
持续的马克斯
如果时间= = 0 | | isempty (max)
max =无穷;
结束
max = max (max, s);
=最大;
结束
然后,您可以使用此函数在一个重复的任务:
%建立模型。
m1 = sbiomodel (“测试”);
%添加舱
c1 = addcompartment (m1,“薪酬”1);
%添加物种间。
addspecies (c1,“年代”);
%规则添加到模型中。
addrule (m1,' s =罪(时间* 2 * pi + 2) + 2”,“repeatedAssignment”);
%添加虚拟歌唱模型
addspecies (c1,“s2”10);
addrule (m1,“s2 = s2 ',“速度”);
%将参数添加到模型中。
addparameter (m1,“马克斯”、常数= false);
addrule (m1,' max = getmax(时间、年代),“repeatedAssignment”);
%修改设置
c = getconfigset (m1);
cs.SolverOptions。MaxStep = 0.01;
%模拟
simfun = createSimFunction (m1, {}, {“马克斯”,“年代”},[][],AutoAccelerate = false);
stopTime = 1.5;
结果= simfun ([], stopTime);
sbioplot(结果);
我希望这可以帮助。
最好的问候,
杰里米
