[论文解读] Splinets -- efficient orthonormalization of the B-splines

Xijia Liu, Hiba Nassar+1 more|arXiv (Cornell University)|Oct 16, 2019

Advanced Numerical Analysis Techniques参考文献 7被引用 4

一句话总结

本文提出了一种高效正交化B样条基——splinets，其在保持局部性和对称性的同时，通过分层的二元正交化过程实现计算效率。对于k阶B样条，该方法将内积计算量减少至O(n)，总支撑大小增长为O(log n)，显著优于以往正交化基的O(n)规模。

ABSTRACT

A new efficient orthogonalization of the B-spline basis is proposed and contrasted with some previous orthogonalized methods. The resulting orthogonal basis of splines is best visualized as a net of functions rather than a sequence of them. For this reason, the basis is referred to as a splinet. The splinets feature clear advantages over other spline bases. They efficiently exploit 'near-orthogonalization' featured by the B-splines and gains are achieved at two levels: locality that is exhibited through small size of the total support of a splinet and computational efficiency that follows from a small number of orthogonalization procedures needed to be performed on the B-splines to achieve orthogonality. These efficiencies are formally proven by showing the asymptotic rates with respect to the number of elements in a splinet. The natural symmetry of the B-splines in the case of the equally spaced knots is preserved in the splinets, while quasi-symmetrical features are also seen for the case of arbitrarily spaced knots.

研究动机与目标

• 开发一种高效正交化B样条基的方法，同时保持其局部性和结构对称性。
• 通过最小化正交化过程中的内积计算次数，降低功能信号分解的计算成本。
• 构建一种支撑性质接近正交的基（即splinets），其支撑大小随节点数的对数增长。
• 提供一种通用算法，用于带状矩阵的稀疏对角化，可推广至样条以外的应用。
• 通过R包实现实际应用，充分利用splinets的数值与结构优势。

提出的方法

• 在二元节点结构中，对B样条实施分层的一侧正交化过程，从最细层级开始。
• 采用递归过程，每个样条仅与其在层次结构中的两个父样条正交化，利用不相交的支撑区域。
• 应用改进的Gram-Schmidt过程并保持对称性，即使在非均匀节点下也能确保准对称特征。
• 推导出变换矩阵P，使得P^T H P = I，其中H为B样条的Gram矩阵，从而实现高效的谱分解。
• 采用基于误差范数的停止规则（例如，迭代3次后‖ε‖² < 10⁻¹⁶），确保数值误差可忽略不计。
• 引入一种计算便捷的样条表示方法，使用正交化基，便于数值实现。

实验结果

研究问题

• RQ1 B样条基能否以保持其局部性和计算效率的方式实现正交化？
• RQ2 正交化B样条基的渐近支撑大小是多少？与现有方法相比如何？
• RQ3 与经典Gram-Schmidt方法相比，该方法实现正交化所需的内积数量是多少？
• RQ4 该正交化过程能否推广至任意具有对称带状结构的带状矩阵？
• RQ5 迭代正交化过程中产生的数值误差在多大程度上影响最终的正交性和稳定性？

主要发现

• 与以往正交化基的O(n)相比，splinets的总支撑大小增长为O(log n)，实现了近似局部性。
• 对于k阶B样条，正交化所需的内积数量为J²ₙ = (5k−1)/4 · n − (2k²/log 2) log(n+1) + …，远低于经典Gram-Schmidt的O(n²)。
• 正交化循环迭代三次后，误差范数降至10⁻¹⁶以下，表明已收敛至机器精度。
• 变换矩阵P为稀疏且近似带状，绝大多数非零元素可忽略，支持高效存储与计算。
• 该方法在等距节点情况下保持对称性，在非规则节点下表现出准对称性，维持了结构保真度。
• 该方法可推广至任意具有2k−1条对角线的正定带状矩阵，通过相同算法实现高效对角化。