脓疱病（皮肤溃疡）确定性及随机建模与最优控制分析的混合人工智能框架

《Results in Chemistry》：Hybrid artificial intelligence framework for deterministic and stochastic modelling of skin sores disease with optimal control analysis

　　

脓疱病（impetigo）是偏远和资源有限社区儿童发病的主要原因，其中群体治疗而非疫苗接种是可操作的公共卫生杠杆。研究人员构建了一个四房室确定性模型，在已有公式基础上扩展了两个治疗参数：暴露个体的去定植（decolonisation）和活动性皮损的抗生素治疗。该模型具有正不变性和一致有界性，当基本再生数（basic reproduction number，R0）低于1时，无病平衡点在局部和全局均渐近稳定。敏感性分析确定活动性皮损治疗是最具影响的可控参数（弹性?0.691）。It?随机扩展允许唯一全局正解，灭绝阈值表明即使确定性再生数超过1，较大的环境波动也可根除感染。以疫苗接种和两种治疗率为时间相关控制，Pontryagin最大值原理（Pontryagin's Maximum Principle）结合前向–后向扫描在200天内将成本泛函降低93.3%、累积感染降低97.0%，最优投入强烈前置。前馈神经网络（feedforward neural networks）提供准确的确定性和分布替代模型，并从含噪每周病例计数中在确定性动态下恢复参数（R2≈0.89），但从单一随机实现中仅弱恢复（R2=0.36），表明单社区监测无法识别传播。活动性皮损治疗而非疫苗接种是主导杠杆。

研究背景与立项依据：脓疱病是由金黄色葡萄球菌（Staphylococcus aureus）和A组链球菌（Streptococcus pyogenes）引起的儿童最常见细菌性皮肤感染，在澳大利亚北部偏远土著社区儿童患病率可超40%，且是急性风湿热与风湿性心脏病的前驱环节。目前无上市疫苗，现场可执行杠杆仅为卫生促进、环境改善与抗生素治疗（含群体给药）。已有房室模型多把恢复视为被动生物过程，未把治疗作为可调干预，亦未整合随机波动与神经网络推断，难以回答“有限预算下努力如何随时间分配”这一动态问题。为此，研究人员在《Results in Chemistry》发表上述论文，建立确定性–随机–最优控制–前馈神经网络（feedforward neural network，FNN）统一框架，量化治疗与接种相对效果、环境波动对持续/灭绝的影响、资源最优配置及神经网络近似与参数反演能力。

主要关键技术方法：研究人员采用四房室（易感S、暴露/定植L、活动性皮损I、恢复R）标准发生率确定性常微分方程组，增设去定植率τ₁与抗生素率τ₂两个可控治疗参数及假设性接种率β；通过下一代矩阵法算R₀，用Castillo-Chavez分解证全局稳定，正规化前向敏感指数排序参数；随机部分改为It?漂移项加乘性白噪声，证唯一全局正解与波动灭绝阈值；最优控制以疫苗接种、τ₁(t)、τ₂(t)为时变控制，成本泛函含感染与投入项，用Pontryagin最大值原理配前向–后向扫掠数值求解；FNN分别以房室轨迹为输入输出做替代模型，并以含噪周病例计数反演确定性参数（R²≈0.89）与单条随机轨迹反演（R²=0.36）。参数基线取自远程社区纵向区间删失估计（传染期12–20天）。

研究结果

模型构建（2. Model formulation）：研究人员将总人口N=S+L+I+R，以κSI/N为标准发生率，L以ξ进展为I，以γ自然清除、τ₁去定植，I以σ自然愈、τ₂抗生素愈，S以λ招募、β接种、φ免疫衰退、μ自然死亡。设τ₁=τ₂=0退化为无治疗模型。

解的正性与有界性（2.1–2.2）：在非负初值下各边界通量非负，解恒留?⁴₊；总种群dN/dt=λ?μN，故N(t)→λ/μ，区域Θ={N≤λ/μ}正不变且吸引。

平衡点（2.3–2.4）：无病平衡点E₀=(λ(μ+φ)/(μ(β+μ+φ)),0,0,βλ/(μ(β+μ+φ)))；地方病平衡点E^*存在当且仅当R₀>1。

基本再生数（2.5）：下一代矩阵得R₀=κξ(μ+φ)/((β+μ+φ)(μ+σ+τ₂)(γ+μ+τ₁+ξ))。接种仅经(μ+φ)/(β+μ+φ)因子呈边际递减，τ₁、τ₂于分母加法无界下压R₀。

稳定性（2.6，3）：R₀<1时E₀ Jacobian感染块行列式正、Routh-Hurwitz满足，局部渐近稳定；用Castillo-Chavez M₁/M₂分解，B为Metzler矩阵、?≥0，证E₀在Θ内全局渐近稳定。

敏感性分析（4）：正规化敏感指数显示κ为+1.000，τ₂为?0.691（最强可控负向），σ为?0.309，β负向但边际递减；表明活动性皮损治疗是现场最具杠杆的可调参数。

随机扩展与灭绝：It?版存唯一全局正解；当环境波动强度超阈值，即便R₀>1感染仍随机灭绝，说明环境噪声本身可成控制因子。

最优控制：以β(t),τ₁(t),τ₂(t)为控制，200天 horizon 下Pontryagin前向–后向扫掠使成本泛函降93.3%、累积感染降97.0%，控制曲线强前置（front-loaded），即早期集中投入远优于均摊。

神经网络替代与反演：FNN对确定性轨迹与随机分布均高精度替代；用含噪周病例在确定性下反演参数R²≈0.89，但单条随机实现反演R²=0.36，证单社区噪声轨迹不可辨识传播参数。

讨论与结论翻译：现有脓疱病房室模型把恢复作被动过程、缺治疗作为可控项，且未把随机波动、最优时序与神经网络反演统合。本研究补入τ₁、τ₂两治疗控制，证R₀<1时无病平衡全局稳定，敏感性指出τ₂弹性?0.691为最杠杆项；It?扩展给出波动灭绝机制；Pontryagin解显示200天内成本与感染分别降93.3%与97.0%且投入前置；FNN可作替代但与单社区随机监测反演受限。结论：对无疫苗的脓疱病，活动性皮损抗生素处理（非疫苗接种）是主导公共卫生杠杆，资源应前期倾斜，环境波动可辅助灭绝，神经网络宜作确定性替代而非单社区随机辨识工具。