庞特里亚金最小值原理进行最优化控的 MATLAB 代码

我正在做一个含时滞的传染病模型。我想利用庞特里亚金最小值原理进行最优化控制？有没有相应的 MATLAB 代码