[论文解读] Optimal control and numerical software: an overview

Helena Sofia Rodrigues, M. Teresa T. Monteiro+1 more|arXiv (Cornell University)|Jan 28, 2014

Advanced Control Systems Optimization参考文献 24被引用 27

一句话总结

本文全面概述了用于求解最优控制（OC）问题的数值软件工具，通过麻疹流行病学模型对比了直接法与间接法。评估了五种软件包——MATLAB的ode45、OC-ODE、DOTcvp、IPOPT和SNOPT，展示了它们在使用直接转录和非线性规划技术求解复杂OC问题时的性能、精度和鲁棒性。

ABSTRACT

Optimal Control (OC) is the process of determining control and state trajectories for a dynamic system, over a period of time, in order to optimize a given performance index. With the increasing of variables and complexity, OC problems can no longer be solved analytically and, consequently, numerical methods are required. For this purpose, direct and indirect methods are used. Direct methods consist in the discretization of the OC problem, reducing it to a nonlinear constrained optimization problem. Indirect methods are based on the Pontryagin Maximum Principle, which in turn reduces to a boundary value problem. In order to have a more reliable solution, one can solve the same problem through different approaches. Here, as an illustrative example, an epidemiological application related to the rubella disease is solved using several software packages, such as the routine ode45 of Matlab, OC-ODE, DOTcvp toolbox, IPOPT and Snopt, showing the state of the art of numerical software for OC.

研究动机与目标

• 评估用于求解复杂现实应用中最优控制问题的最先进数值软件。
• 通过一个实际的流行病学案例研究，对比直接法与间接法在性能、精度和鲁棒性方面的表现。
• 通过基准麻疹传播模型，展示现代非线性规划求解器在求解最优控制问题中的有效性。
• 为研究人员提供一份实用指南，基于问题结构、计算成本和收敛行为选择合适的软件工具。

提出的方法

• 本研究采用直接法，通过使用欧拉法和均匀时间划分，将最优控制问题离散化为非线性规划（NLP）问题。
• 麻疹传播模型以Bolza形式表述，状态变量代表易感者、感染者、康复者和接种者人群。
• 性能指标在疾病控制（通过最小化感染者数量）与干预成本（通过控制努力惩罚）之间进行权衡。
• 使用五种软件包——MATLAB的ode45、OC-ODE、DOTcvp、IPOPT和SNOPT——求解离散化问题，AMPE作为NEOS服务器上的建模语言。
• IPOPT和SNOPT分别采用内点法和序列二次规划（SQP）方法，利用稀疏雅可比和海塞矩阵结构以提高效率。
• 求解过程结合自动微分和稀疏性利用，以提升收敛速度和数值稳定性。

实验结果

研究问题

• RQ1 不同数值软件包在求解具有现实世界复杂性的最优控制问题时表现如何比较？
• RQ2 在非线性流行病学模型上应用时，直接法与间接法在相对精度和鲁棒性方面有何差异？
• RQ3 内点法与基于SQP的求解器在最优控制问题的收敛速度和解质量方面表现如何？
• RQ4 标准非线性规划工具在处理具有生物动力学的复杂、约束性最优控制问题方面的能力如何？
• RQ5 稀疏性利用和导数信息在提升最优控制软件效率方面起到何种作用？

主要发现

• 当与现代NLP求解器如IPOPT和SNOPT结合时，直接法即使在初始猜测较差的情况下，也能提供高度准确和鲁棒的解。
• IPOPT和SNOPT由于有效利用了雅可比矩阵和海塞矩阵中的稀疏性，表现出优越的收敛行为和计算效率。
• 在AMPL支持的建模环境中使用自动微分显著提高了导数计算的可靠性和速度。
• 使用IPOPT和SNOPT求解的解达到了最小的目标函数值，表明疾病控制与干预成本之间实现了有效权衡。
• MATLAB的ode45在最优控制方面表现较差，因其缺乏优化能力，凸显了专用NLP求解器的必要性。
• 本研究证实，直接法在复杂、高维最优控制问题中更具实用性，尤其是在涉及约束和非线性动力学时。