我想确定周期性延迟微分方程系统（季节性宿主寄生虫模型）的稳定性。我尝试实现本文引理 2.5 中描述的方法：

• 梁 X.、张 L. 和赵 XQ。周期性抽象泛函微分方程的基本再现比率（应用于莱姆病的空间模型）。J Dyn Diff Equat 31 , 1247–1278 (2019)。：

引理 2.5假设 (E, E ₊ ) 是一个有序 Banach 空间，其中 E ₊是正规的并且 Int(E ₊ )$\neq$∅，符合范数||·|| _E。 _ 令 L 为正有界线性算子。选择 v ₀ ∈ Int(E ₊ ) 并定义_n = ||Lv _n-1 || _E , v _n =$\frac{Lv_{n-1}}{ > a_{n}}$, ∀n ≥ 1。如果$\lim_{(n→+∞)}a_n$存在，则 r(L) = 。$\lim_{n→+∞}a_n$

但我不确定它是否正常工作（我可能在实现的某个地方犯了一个错误）。问题是它似乎并不总是对应于我在很长一段时间内简单地模拟模型时发现的平衡点。我研究了为延迟方程近似 Floquet 乘数/单数矩阵等的不同方法，以便我可以仔细检查我的工作，并且似乎有许多不同的方法（例如，半离散化、有限元分析、切比雪夫多项式等.)

对于比我更熟悉该文献的任何人，我的问题是：目前有首选方法吗？有没有容易快速实现的（我不是一个非常老练的程序员或数学家）？我会对带有已发布代码的示例特别感兴趣。一些建议将不胜感激。