你好,打扰您了,想知道对于一个二维的具有双时滞的捕食系统,这种分岔图用matlab怎么画

以下是我截取相关论文中的分岔图,分别对于时滞t1和时滞t2作出对应的分岔图,相关程序不清楚,希望您有空能看一下

img

因为具体论文我看不到,我基于截图信息发表下见解,这里他是基于τ2=0讨论的τ1。
实际上要画带时滞的图,只需要在原二维捕食系统上加一个滞后就好,就相当于t≥τ1的时候对应方程才起作用。
比如τ1=0.1,总时长t=1(举例),按照总共迭代100步,在第11步的时候加入含τ1的方程(10步以内含τ1方程不参与迭代)
(二维捕食系统里有经典的L-V系统,挺多二维论文模型都是基于那个改进的)