求助，时滞微分方程组，不求程序，只求分析和计算的方法

flashimba

也就是说只需指导思想，特别要指点计算步骤，并不需要具体代码或者算式。多谢了。

如下的共12个的时滞微分方程组（4个一组，公式不好打就直接贴程序里的）。参数为ylag1和g，我想分析Y1 、Y5 、Y9。

function [dydt]=ddex3de(t,y,Z)
r1=(0.1*y(1)+4)/(1-exp(-0.1*y(1)-4));
r2=4*exp((y(1)+65)/-18);
r3=0.07*exp((y(1)+65)/-20);
r4=1/(1+exp(-0.1*y(1)-3.5));
r5=(0.01*y(1)+0.55)/(1-exp(-0.1*y(1)-5.5));
r6=0.125*exp((y(1)+65)/-80);
%------------------------------------------------------
q1=(0.1*y(5)+4)/(1-exp(-0.1*y(5)-4));
q2=4*exp((y(5)+65)/-18);
q3=0.07*exp((y(5)+65)/-20);
q4=1/(1+exp(-0.1*y(5)-3.5));
q5=(0.01*y(5)+0.55)/(1-exp(-0.1*y(5)-5.5));
q6=0.125*exp((y(5)+65)/-80);
%------------------------------------------------------
s1=(0.1*y(9)+4)/(1-exp(-0.1*y(9)-4));
s2=4*exp((y(9)+65)/-18);
s3=0.07*exp((y(9)+65)/-20);
s4=1/(1+exp(-0.1*y(9)-3.5));
s5=(0.01*y(9)+0.55)/(1-exp(-0.1*y(9)-5.5));
s6=0.125*exp((y(9)+65)/-80);
%------------------------------------------------------
ylag1=Z(:,1);
%延迟-----------------------------------------------------
%主方程-----------------------------------------------------
dydt=[20-120*(y(2)^3)*y(3)*(y(1)-50)-36*(y(4)^4)*(y(1)+77)-0.3*(y(1)+54.4)+g*(ylag1(5)-y(1))+g*(ylag1(9)-y(1))
r1*(1-y(2))-r2*y(2)
r3*(1-y(3))-r4*y(3)
r5*(1-y(4))-r6*y(4)
20-120*(y(6)^3)*y(7)*(y(5)-50)-36*(y(8)^4)*(y(5)+77)-0.3*(y(5)+54.4)+g*(ylag1(1)-y(5))+g*(ylag1(9)-y(5))
q1*(1-y(6))-q2*y(6)
q3*(1-y(7))-q4*y(7)
q5*(1-y(8))-q6*y(8)
20-120*(y(10)^3)*y(11)*(y(9)-50)-36*(y(12)^4)*(y(9)+77)-0.3*(y(9)+54.4)+g*(ylag1(1)-y(9))+g*(ylag1(5)-y(9))
s1*(1-y(10))-s2*y(10)
s3*(1-y(11))-s4*y(11)
s5*(1-y(12))-s6*y(12)];

我的问题是：
1 原文有 “周期TP=10.43，算TP=10.43处的解算子的特征根即floquet乘子”这句话不明白。
那么我想要分析这个方程组的周期解的稳定性（主要分析Y1 、Y5 、Y9），那么这个方程的floquet乘子怎么算，大致步骤。
2   还有就是庞加莱截面的做法 ，也是主要分析Y1 、Y5 、Y9.
3  如果能说分岔图怎么画的大致思想。

pjfly

用DDE23试试