[论文解读] Convergence of univariate non-stationary subdivision schemes via asymptotical similarity

Costanza Conti, Nira Dyn+2 more|arXiv (Cornell University)|Oct 10, 2014

Advanced Numerical Analysis Techniques参考文献 13被引用 8

一句话总结

本文引入渐近相似性——一种弱于渐近等价性的条件——作为保证单变量非平稳细分格式收敛的充分条件，该格式可再现常数。通过证明当平稳参考格式收敛时，若非平稳格式与之渐近相似，则其亦收敛，作者将收敛性分析扩展至更广泛的格式类别，包括生成指数多项式与样条的格式，在几何建模与等几何分析中具有应用价值。

ABSTRACT

A new equivalence notion between non-stationary subdivision schemes, termed asymptotical similarity, which is weaker than asymptotical equivalence, is introduced and studied. It is known that asymptotical equivalence between a non-stationary subdivision scheme and a convergent stationary scheme guarantees the convergence of the non-stationary scheme. We show that for non-stationary schemes reproducing constants, the condition of asymptotical equivalence can be relaxed to asymptotical similarity. This result applies to a wide class of non-stationary schemes of importance in theory and applications.

研究动机与目标

• 为解决非平稳细分格式收敛性分析的挑战，此类格式相较于平稳格式更为复杂。
• 通过引入一个更弱、更通用的条件，放松此前收敛性所需严格满足的渐近等价性条件。
• 建立渐近相似性与常数再现相结合可保证一大类非平稳格式收敛的理论结果。
• 为生成指数多项式与样条的格式提供一个理论框架，此类格式在几何建模与等几何分析中具有重要意义。
• 通过非平稳 de Rham 算法的变体，展示新准则的实际适用性。

提出的方法

• 引入非平稳细分格式之间渐近相似性的概念，定义为 $\lim_{k\to\infty}\|\mathbf{a}^{[k]} - \mathbf{a}^{*[k]}\| = 0$。
• 定义并分析差分格式，推导出可再现常数的非平稳格式收敛性的充分条件。
• 证明：若两个格式均可再现常数且渐近相似，且其中一个满足收敛的充分条件，则另一个亦满足。
• 使用符号表示 $a^{[k]}(z) = (1+z)q^{[k]}(z)$，其中 $q^{[k]}(1) = 1$，以刻画可再现常数的格式。
• 将理论框架应用于非平稳 de Rham 型格式，其掩码为 $\mathbf{a}^{[k]} = \{\frac{1}{2+\gamma_k}, \frac{1+\gamma_k}{2+\gamma_k}, \frac{1+\gamma_k}{2+\gamma_k}, \frac{1}{2+\gamma_k}\}$，其中 $\gamma_k = \gamma + \varepsilon_k$ 且 $\varepsilon_k \to 0$。
• 利用渐近相似性在参考格式收敛且格式可再现常数时保持收敛性的事实。

实验结果

研究问题

• RQ1 在弱于渐近等价性的条件下，是否可保证非平稳细分格式的收敛性？
• RQ2 当与常数再现结合时，渐近相似性（即掩码逐点收敛）是否足以保证收敛性？
• RQ3 在多大程度上，可通过其与收敛平稳格式的相似性来推断非平稳格式的收敛性？
• RQ4 新准则如何应用于实际格式，如非平稳 de Rham 算法或指数 B 样条格式？
• RQ5 当常数再现不成立时，渐近相似性存在哪些局限性？

主要发现

• 渐近相似性严格弱于渐近等价性，因其仅要求 $\lim_{k\to\infty}\|\mathbf{a}^{[k]} - \mathbf{a}^{*[k]}\| = 0$。
• 对于可再现常数的非平稳格式，若其与收敛的平稳格式渐近相似，则该非平稳格式亦收敛。
• 该结果适用于一大类格式，包括生成包含一个零指数的指数多项式空间的格式。
• 当 $\gamma_k = \gamma + \varepsilon_k$ 且 $\varepsilon_k \to 0$ 时，非平稳 de Rham 格式与经典平稳 de Rham 格式渐近相似，因此对所有 $\gamma > 0$ 均收敛。
• 由 $\mathbf{a}^{[k]} = \{\frac{1}{4} + \frac{1}{k}, \frac{3}{4} + \frac{1}{k}, \frac{3}{4} + \frac{1}{k}, \frac{1}{4} + \frac{1}{k}\}$ 定义的格式与 Chaikin 格式渐近相似，但不再现常数，因此不收敛，该例说明常数再现假设的必要性。
• 该理论框架具有鲁棒性，可扩展至非正则情形，如相关工作中所指出。