Risk Analysis and Genomic Epidemiology Unit, Istituto Zooprofilattico Sperimentale della Lombardia e dell'Emilia Romagna, Via dei Mercati 13, 43126, Parma, Italy.
Department of Mathematical, Physical and Computer Sciences, University of Parma, Parco Area delle Scienze 53/A, 43124, Parma, Italy.
J Math Biol. 2021 Sep 22;83(4):36. doi: 10.1007/s00285-021-01668-1.
Mathematical models are formal and simplified representations of the knowledge related to a phenomenon. In classical epidemic models, a major simplification consists in assuming that the infectious period is exponentially distributed, then implying that the chance of recovery is independent on the time since infection. Here, we first attempt to investigate the consequences of relaxing this assumption on the performances of time-variant disease control strategies by using optimal control theory. In the framework of a basic susceptible-infected-removed (SIR) model, an Erlang distribution of the infectious period is considered and optimal isolation strategies are searched for. The objective functional to be minimized takes into account the cost of the isolation efforts per time unit and the sanitary costs due to the incidence of the epidemic outbreak. Applying the Pontryagin's minimum principle, we prove that the optimal control problem admits only bang-bang solutions with at most two switches. In particular, the optimal strategy could be postponing the starting intervention time with respect to the beginning of the outbreak. Finally, by means of numerical simulations, we show how the shape of the optimal solutions is affected by the different distributions of the infectious period, by the relative weight of the two cost components, and by the initial conditions.
数学模型是与现象相关的知识的形式化和简化表示。在经典的传染病模型中,一个主要的简化是假设传染期呈指数分布,从而意味着恢复的机会与感染后的时间无关。在这里,我们首先尝试通过使用最优控制理论来研究放宽这一假设对时变疾病控制策略性能的影响。在基本的易感-感染-清除(SIR)模型框架中,考虑了传染期的爱尔朗分布,并搜索了最优隔离策略。要最小化的目标函数考虑了每个时间单位的隔离工作成本和因疫情爆发而产生的卫生成本。通过应用庞特里亚金最小原理,我们证明最优控制问题仅存在至多两个开关的 bang-bang 解。特别是,最优策略可能会推迟干预时间相对于疫情爆发的开始。最后,通过数值模拟,我们展示了不同的传染期分布、两个成本组成部分的相对权重以及初始条件如何影响最优解的形状。