[论文解读] A PIE Representation of Coupled 2D PDEs and Stability Analysis using LPIs

Declan S. Jagt, Matthew M. Peet|arXiv (Cornell University)|Sep 14, 2021

Numerical methods for differential equations参考文献 15被引用 4

一句话总结

本文提出了一种用于耦合二维线性PDE的偏积分方程（PIE）表示方法，将先前仅适用于一维系统的PIE框架扩展至二维空间。通过在L₂[x,y]上构建偏积分（PI）算子的微分学，作者建立了PIE解与原始PDE解之间的双射映射，从而实现基于线性PI不等式（LPI）的凸优化稳定性分析，并在二维热方程和波动方程上验证了其极低的保守性。

ABSTRACT

We introduce a Partial Integration Equation (PIE) representation of Partial Differential Equations (PDEs) in two spatial variables. PIEs are an algebraic state-space representation of infinite-dimensional systems and have been used to model 1D PDEs and time-delay systems without continuity constraints or boundary conditions -- making these PIE representations amenable to stability analysis using convex optimization. To extend the PIE framework to 2D PDEs, we first construct an algebra of Partial Integral (PI) operators on the function space L_2[x,y], providing formulae for composition, adjoint, and inversion. We then extend this algebra to R^n x L_2[x] x L_2[y] x L_2[x,y] and demonstrate that, for any suitable coupled, linear PDE in 2 spatial variables, there exists an associated PIE whose solutions bijectively map to solutions of the original PDE -- providing conversion formulae between these representations. Next, we use positive matrices to parameterize the convex cone of 2D PI operators -- allowing us to optimize PI operators and solve Linear PI Inequality (LPI) feasibility problems. Finally, we use the 2D LPI framework to provide conditions for stability of 2D linear PDEs. We test these conditions on 2D heat and wave equations and demonstrate that the stability condition has little to no conservatism.

研究动机与目标

• 将PIE框架从先前用于一维PDE和时滞系统的应用扩展至二维空间。
• 在L₂[x,y]上建立偏积分（PI）算子的完整代数结构，涵盖复合、伴随和求逆运算。
• 建立二维PDE解与其对应PIE表示之间的双射对应关系。
• 通过正定矩阵参数化二维PI算子的凸锥，以支持优化分析。
• 推导并验证基于LPI的二维线性PDE稳定性条件，且保守性极低。

提出的方法

• 在L₂[x,y]函数空间上构建PI算子代数，定义复合、伴随和求逆等运算。
• 将PI算子代数扩展至Rⁿ × L₂[x] × L₂[y] × L₂[x,y]，以适应耦合的二维PDE。
• 证明每个合适的耦合线性二维PDE均存在等价的PIE表示形式，并给出明确的解转换公式。
• 利用正定矩阵参数化二维PI算子的凸锥，以支持凸优化。
• 构建线性PI不等式（LPI）可行性问题，以推导二维PDE的稳定性条件。
• 将LPI框架应用于二维热方程和波动方程的稳定性分析，采用凸优化技术。

实验结果

研究问题

• RQ1 PIE框架能否从一维系统成功扩展至具有多个空间变量的二维PDE？
• RQ2 在二维空间中，PI算子的复合与求逆运算受何种代数结构支配？
• RQ3 是否存在耦合二维PDE与其对应PIE表示之间的双射解映射？
• RQ4 如何对二维PI算子的凸锥进行参数化，以支持基于优化的分析？
• RQ5 基于LPI的稳定性条件在应用于二维线性PDE时，其保守性如何？

主要发现

• 本文成功构建了在L₂[x,y]上的完整PI算子代数，使得二维PDE的PIE形式可系统化处理。
• 建立了原始二维PDE与其PIE表示之间的双射解映射，确保了解的等价性。
• 通过正定矩阵参数化二维PI算子的凸锥，实现了稳定性分析的凸优化。
• 在二维热方程和波动方程上的测试表明，基于LPI的稳定性条件几乎无保守性。
• 该框架无需连续性约束或显式边界条件，即可实现稳定性分析，显著提升了适用性。