[论文解读] Regular polynomial interpolation and approximation of global solutions of linear partial differential equations

Joerg Kampen|ArXiv.org|Jul 9, 2008

Numerical Methods and Algorithms参考文献 5被引用 6

一句话总结

本文提出了一种递归的正则多项式插值方法，用于近似求解线性偏微分方程（PDE）的全局解，该方法利用椭圆型、双曲型和抛物型系统的事先估计。该方法构造的多项式族在插值点处保持PDE的结构特性，并实现了对解正则性的敏感误差估计，从而可与高维求解器（如稀疏网格和蒙特卡洛方法）兼容。

ABSTRACT

We consider regular polynomial interpolation algorithms on recursively defined sets of interpolation points which approximate global solutions of arbitrary well-posed systems of linear partial differential equations. Convergence of the 'limit' of the recursively constructed family of polynomials to the solution and error estimates are obtained from a priori estimates for some standard classes of linear partial differential equations, i.e. elliptic and hyperbolic equations. Another variation of the algorithm allows to construct polynomial interpolations which preserve systems of linear partial differential equations at the interpolation points. We show how this can be applied in order to compute higher order terms of WKB-approximations of fundamental solutions of a large class of linear parabolic equations. The error estimates are sensitive to the regularity of the solution. Our method is compatible with recent developments for solution of higher dimensional partial differential equations, i.e. (adaptive) sparse grids, and weighted Monte-Carlo, and has obvious applications to mathematical finance and physics.

研究动机与目标

• 开发一种稳定且递归的多项式插值算法，用于近似求解线性PDE的全局解，且不施加对插值点配置的限制。
• 确保插值多项式在离散插值点处满足相同的线性PDE系统，从而保持结构特性。
• 基于椭圆型、双曲型和抛物型PDE的事先估计，推导出误差估计，使方法对解的正则性敏感。
• 通过与现代数值技术（如稀疏网格和加权蒙特卡洛）集成，实现对高维问题的应用。
• 将方法扩展至通过精确的导数近似，计算抛物型方程基本解的高阶WKB近似。

提出的方法

• 该方法采用递归构造多元多项式，将每个先前的多项式扩展至在额外一组插值点处满足PDE系统。
• 将牛顿型插值推广至多变量情形，具备单调性与稳定性，且不依赖于插值点集合的约束。
• 该算法被适配于向量值函数，并确保若目标函数满足线性PDE系统，则其插值多项式在插值点处也满足相同的PDE系统。
• 误差估计基于标准PDE类（如椭圆型与双曲型系统）的事先估计，将近似误差与解的光滑性相关联。
• 对于抛物型方程，该方法可通过准确近似高达10阶的导数，实现高阶WKB展开的计算。
• 该方法支持部分并行化，并可将定义在互不相交插值集合上的局部多项式近似拼接为全局解。

实验结果

研究问题

• RQ1 能否设计一种递归多项式插值方案，用于近似求解线性PDE的全局解，且不施加对插值点选择的限制？
• RQ2 如何在离散点处保持插值过程的底层PDE结构，以确保与控制系统的协调一致？
• RQ3 可为多项式近似推导出何种误差估计，使其能反映真实解的正则性？
• RQ4 该方法在多大程度上可与高维数值求解器（如稀疏网格或加权蒙特卡洛）集成？
• RQ5 该方法能否扩展至计算抛物型PDE基本解的WKB展开中的高阶项？

主要发现

• 当结合椭圆型与双曲型方程的事先估计时，递归多项式插值收敛于适定线性PDE系统的全局解。
• 显式推导出误差界，其大小依赖于解的正则性，光滑性越高，误差越小。
• 对于模型问题 f(x) = 1/(1+x)，插值系数 a_i 的幅值随 i 增大而减小，表明具有数值稳定性。
• 通过近似高达10阶的导数，该方法可准确计算一大类线性抛物型方程基本解的WKB展开至5阶。
• 该算法支持部分并行化与局部近似的拼接，可在合并后的插值点集合上构建全局插值多项式。
• 该方法与自适应稀疏网格和加权蒙特卡洛等先进数值技术兼容，显著提升了其在高维场景下的实用性。