[论文解读] Trivariate polynomial approximation on Lissajous curves

Len Bos, Stefano De Marchı+1 more|arXiv (Cornell University)|Feb 13, 2015

Advanced Numerical Analysis Techniques参考文献 27被引用 5

一句话总结

本论文提出了一种基于三维Lissajous曲线的三元多项式逼近框架，用于在3-立方体中构造超插值和离散极值点集（Fekete与Leja型），适用于秩-1切比雪夫格栅。通过利用一种新型的三维Lissajous曲线，其参数由满足特定正交性条件的整数三元组定义，该方法可通过单次一维快速切比雪夫变换实现高效计算，获得较低的Lebesgue常数和高逼近精度，具有在磁性粒子成像（MPI）中应用的潜力。

ABSTRACT

We study Lissajous curves in the 3-cube, that generate algebraic cubature formulas on a special family of rank-1 Chebyshev lattices. These formulas are used to construct trivariate hyperinterpolation polynomials via a single 1-d Fast Chebyshev Transform (by the Chebfun package), and to compute discrete extremal sets of Fekete and Leja type for trivariate polynomial interpolation. Applications could arise in the framework of Lissajous sampling for MPI (Magnetic Particle Imaging).

研究动机与目标

• 将Lissajous曲线框架从二元推广至3-立方体中的三元多项式逼近。
• 利用位于三维Lissajous曲线上的一组节点，构造乘积切比雪夫测度的代数求积公式。
• 通过沿曲线的单次一维快速切比雪夫变换，实现高效的三元超插值。
• 从所得切比雪夫格栅中提取离散极值点集（Fekete与Leja型），以实现稳定多项式插值。
• 提供一种计算高效、数值稳定的框架，适用于新兴应用如磁性粒子成像（MPI）。

提出的方法

• 通过基于n奇偶性的分段整数公式导出的频率三元组(aₙ, bₙ, cₙ)，定义三维Lissajous曲线，确保总次数不超过2n的单项式之间正交。
• 证明：任意总次数≤2n的三元多项式在3-立方体上关于乘积切比雪夫测度的积分，等于其沿Lissajous曲线积分的π²倍。
• 通过沿曲线采样函数值并应用单次一维快速切比雪夫变换，构造三元超插值，实现高效正交投影。
• 通过求解具有Vandermonde型结构的线性系统并选取解中非零分量对应的索引，提取近似Fekete点（AFP）。
• 通过在相同矩阵上进行带行主元的LU分解，计算离散Leja点（DLP），确保嵌套性与单变量性，适用于逐步插值。
• 通过Matlab工具包进行数值实验验证，结果表明Lebesgue常数较低，对测试函数的插值精度高。

实验结果

研究问题

• RQ1 能否构造一条三维Lissajous曲线，使得任意总次数≤2n的多项式在3-立方体上的积分，等于其沿曲线积分的常数倍？
• RQ2 所得到的曲线是否能生成支持通过一维变换实现稳定且高精度多项式超插值的秩-1切比雪夫格栅？
• RQ3 能否从格栅节点中提取离散极值点集（Fekete与Leja型），以实现低Lebesgue常数的高精度插值？
• RQ4 所提取点集的Lebesgue常数与理论界及连续Fekete点的常数相比如何？
• RQ5 该框架是否可在采样受限场景（如磁性粒子成像MPI）中实际应用，其中数据沿单一轨迹采集？

主要发现

• 所提出的三维Lissajous曲线确保了所有总次数≤2n的多项式在3-立方体上的切比雪夫积分与沿曲线的一维积分完全相等。
• 该方法仅通过一次一维快速切比雪夫变换即可实现三元超插值，与多维方法相比显著降低了计算成本。
• 近似Fekete点与离散Leja点的Lebesgue常数远低于理论上的上界dim(ℙₙ³)，且随n近似呈二次增长。
• 离散Leja点构成嵌套序列，允许对任意阶数r≤n实现逐步插值，且每一级均保持单变量性。
• 数值结果表明，对测试函数f₁与f₂的插值具有高精度，相对误差随阶数增加迅速减小，证实了该方法的鲁棒性。
• 该框架已实现于Matlab工具包中，并提供可下载的节点坐标，可直接用于MPI等应用，其中采样受限于单一轨迹。