Optimal Feedback Control of the Fokker-Planck Equation for an SIRS Epidemic Model
Abstract
We develop a Fokker--Planck feedback control framework for a stochastic susceptible--infected--recovered--susceptible (SIRS) epidemic model with three intervention mechanisms targeting disease transmission, loss of immunity, and recovery. The interventions are formulated as space- and time-dependent feedback controls that influence both the drift and diffusion coefficients of the underlying stochastic differential equations, allowing control of the probability distribution of epidemic states. We investigate the existence of optimal controls and derive first-order necessary optimality conditions involving a forward Fokker--Planck equation, a backward adjoint equation, and a Hamiltonian minimization condition. For numerical approximation, we combine a Chang--Cooper discretization of the diagonal drift--diffusion terms with a finite-difference treatment of cross-diffusion and a sequential quadratic Hamiltonian (SQH) algorithm for computing the optimal feedback policies. We also examine sufficient-decrease and convergence properties of the SQH iteration. Numerical experiments consider four epidemiological objectives involving cumulative infection burden, infection-threshold exceedance, terminal susceptible population, and a combination of running and terminal costs. The resulting feedback policies exhibit distinct spatial and temporal structures depending on the epidemiological objective. Monte Carlo simulations of the underlying stochastic system illustrate that the optimized interventions can substantially reduce infection growth and undesirable epidemic outcomes relative to uncontrolled dynamics. The proposed framework provides a distribution-based approach to stochastic epidemic control that integrates mathematical analysis, numerical optimization, and state-dependent intervention design.
Related
- Optimal Control of Traveler Testing in Epidemic Models
- Complex dynamics of a fractional-order epidemic system integrating with fractional-order optimal control and PRCC analysis
- Adaptive Grid Refinement for Optimal Feedback Control
- Optimal control of the Dyson equation
- Solution to a fully coupled McKean-Vlasov forward-backward stochastic difference equation and applications to optimal control with law-delay
BibTeX
@misc{parkinson2026optimal,
title = {{Optimal Feedback Control of the Fokker-Planck Equation for an SIRS Epidemic Model}},
author = {Christian Parkinson and Souvik Roy},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.08708},
eprint = {2611.08708},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.08708},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Optimal Feedback Control of the Fokker-Planck Equation for an SIRS Epidemic Model
AU - Christian Parkinson
AU - Souvik Roy
PY - 2026
DA - 2026/11/05
PB - arXiv
DO - 10.48550/arxiv.2611.08708
UR - https://arxiv.org/abs/2611.08708
AB - We develop a Fokker--Planck feedback control framework for a stochastic susceptible--infected--recovered--susceptible (SIRS) epidemic model with three intervention mechanisms targeting disease transmission, loss of immunity, and recovery. The interventions are formulated as space- and time-dependent feedback controls that influence both the drift and diffusion coefficients of the underlying stochastic differential equations, allowing control of the probability distribution of epidemic states. We investigate the existence of optimal controls and derive first-order necessary optimality conditions involving a forward Fokker--Planck equation, a backward adjoint equation, and a Hamiltonian minimization condition. For numerical approximation, we combine a Chang--Cooper discretization of the diagonal drift--diffusion terms with a finite-difference treatment of cross-diffusion and a sequential quadratic Hamiltonian (SQH) algorithm for computing the optimal feedback policies. We also examine sufficient-decrease and convergence properties of the SQH iteration. Numerical experiments consider four epidemiological objectives involving cumulative infection burden, infection-threshold exceedance, terminal susceptible population, and a combination of running and terminal costs. The resulting feedback policies exhibit distinct spatial and temporal structures depending on the epidemiological objective. Monte Carlo simulations of the underlying stochastic system illustrate that the optimized interventions can substantially reduce infection growth and undesirable epidemic outcomes relative to uncontrolled dynamics. The proposed framework provides a distribution-based approach to stochastic epidemic control that integrates mathematical analysis, numerical optimization, and state-dependent intervention design.
ER -