[论文解读] Fast algorithms for optimal control, anisotropic front propagation and multiple arrivals

James A. Sethian|ArXiv.org|May 8, 2003

Numerical methods for differential equations参考文献 14被引用 5

一句话总结

本文提出了一种快速、单次遍历的算法，用于求解最优控制、各向异性前缘传播以及多波前到达问题中出现的静态Hamilton-Jacobi方程。通过使用Ordered Upwind Methods求解粘性解，并采用带有逃逸方程的相空间公式求解非粘性解，该方法利用信息传播中的因果性，实现了O(N log N)的时间复杂度，从而在连续空间中高效、稳定且精确地计算出第一波前到达和多波前到达的解。

ABSTRACT

We review some recent work in fast, efficient and accurate methods to compute viscosity solutions and non-viscosity solutions to static Hamilton-Jacobi equations which arise in optimal control, anisotropic front propagation, and multiple arrivals in wave propagation. For viscosity solutions, the class of algorithms are known as ``Ordered Upwind Methods'', and rely on a systematic ordering inherent in the characteristic flow of information. For non-viscosity multiple arrivals, the techniques hinge on a static boundary value phase-space formulation which again can be solved through a systematic ordering.

研究动机与目标

• 开发高效的数值方法，求解最优控制、各向异性前缘传播及多波前到达问题中出现的静态Hamilton-Jacobi方程。
• 通过在求解过程中保持因果性和信息流顺序，将Dijkstra最短路径算法推广至连续空间。
• 通过将问题转化为相空间中的两点边值问题，高效计算非粘性解（多波前到达）。
• 实现O(N log N)的计算复杂度，适用于粘性和非粘性解，与离散最短路径算法的效率相当。
• 消除时间步长限制，确保Hamilton-Jacobi方程求解的无条件稳定性。

提出的方法

• 通过Ordered Upwind Methods将Dijkstra算法推广至连续空间，根据信息从已知值向未知值的传播顺序对网格点进行排序。
• 使用优先队列维护并更新“待处理”点的集合，确保每个点仅被处理一次，且计算开销最小。
• 采用相空间表示（x, p），其中x为位置，p为动量，对多波前问题进行建模，以追踪从源点出发的所有可能射线。
• 基于特征线和Liouville方程，推导相空间中的逃逸方程，描述边界处的首次到达时间、位置及梯度。
• 通过系统性地从边界反向推进，利用逃逸方程在相空间中构建单次遍历算法，确保无点被重复访问。
• 通过求解相空间边值问题，无需迭代重算，即可计算出完整的波前到达集合（首次及后续到达）。

实验结果

研究问题

• RQ1 如何将Dijkstra的单次遍历最短路径算法适配于由Hamilton-Jacobi方程控制的连续最优控制问题？
• RQ2 在波传播与前缘演化问题中，如何正确表述非粘性解（多波前到达）的计算方法？
• RQ3 能否在不进行迭代重算或时间步长推进的前提下，为多波前问题构建单次遍历算法？
• RQ4 如何在物理空间或相空间中通过系统性排序，同时在粘性和非粘性解中保持信息传播的因果性？
• RQ5 与传统时间推进格式相比，这些方法的计算复杂度与稳定性如何？

主要发现

• Ordered Upwind Method通过在连续空间中模拟Dijkstra算法，实现了粘性解的O(N log N)复杂度，确保单次遍历收敛。
• 采用相空间公式与逃逸方程，可同时计算所有多波前到达，避免时间步长推进或迭代重处理。
• 该方法无条件稳定，无需时间步长限制，对具有复杂各向异性和慢速区域的问题具有强鲁棒性。
• 数值结果表明，该算法能准确捕捉因波前自相交而产生的尖点与吞没结构，如向内传播的方形波前。
• 该方法在地球物理与光学应用中成功计算出第一波前及多波前到达，且模拟中正确反映了各波前的振幅。
• 计算成本随问题规模高效增长，且单次遍历特性确保相空间中每个网格点最多被访问一次。