[论文解读] Fast, exact (but unstable) spin spherical harmonic transforms

Jason D. McEwen|arXiv (Cornell University)|Jul 28, 2008

Geophysics and Gravity Measurements参考文献 27被引用 7

一句话总结

该论文通过将问题重新表述为在二维环面（two-torus）上的傅里叶变换，提出了一种快速且精确的球面上正向与反向旋量球谐变换算法，利用香农采样定理实现理论上的精确性。尽管实现了 O(L³) 的复杂度和精确性，但由于线性系统的病态性，当 L > 32 时正向变换在数值上变得不稳定；而反向变换虽预期稳定，但因缺乏稳定正向实现而无法验证。

ABSTRACT

In many applications data are measured or defined on a spherical manifold; spherical harmonic transforms are then required to access the frequency content of the data. We derive algorithms to perform forward and inverse spin spherical harmonic transforms for functions of arbitrary spin number. These algorithms involve recasting the spin transform on the two-sphere S^2 as a Fourier transform on the two-torus T^2. Fast Fourier transforms are then used to compute Fourier coefficients, which are related to spherical harmonic coefficients through a linear transform. By recasting the problem as a Fourier transform on the torus we appeal to the usual Shannon sampling theorem to develop spherical harmonic transforms that are theoretically exact for band-limited functions, thereby providing an alternative sampling theorem on the sphere. The computational complexity of our forward and inverse spin spherical harmonic transforms scale as O(L^3) for any arbitrary spin number, where L is the harmonic band-limit of the spin function on the sphere. Numerical experiments are performed and unfortunately the forward transform is found to be unstable for band-limits above L~32. The instability is due to the poorly conditioned linear system relating Fourier and spherical harmonic coefficients. The inverse transform is expected to be stable, although it is not possible to verify this hypothesis.

研究动机与目标

• 开发一种适用于任意旋量数的球面上快速、理论上精确的旋量球谐变换算法。
• 通过球面上新型采样定理的引入，实现计算复杂度为 O(L³) 的同时保持精确性。
• 提供一种与旋量数无关、且无需带宽限制为 2 的幂次的通用算法。
• 研究所提算法的数值稳定性，特别是针对高阶谐波带宽限制的情况。

提出的方法

• 通过变量替换，将球面上的旋量球谐变换重新表述为在二维环面 T² 上的傅里叶变换。
• 利用快速傅里叶变换（FFTs）计算环面上的傅里叶系数，再通过线性变换将其与球谐系数关联。
• 利用环面上的香农采样定理，确保带限函数的理论精确性。
• 采用由 θₜ 和 φₚ 定义的基于网格的像素化方法对函数进行采样，从而支持 FFT 的使用。
• 使用线性系统将傅里叶系数映射到旋量球谐系数，该系统在 L > 32 时变得病态。
• 探索其他公式化方法，包括归一化和几何代数，以提升数值稳定性。

实验结果

研究问题

• RQ1 能否推导出一种快速、精确的旋量球谐变换，使其对任意旋量数均保持 O(L³) 的复杂度？
• RQ2 将球面变换重述为环面上的傅里叶变换，是否能在球面上实现理论上精确的采样定理？
• RQ3 为何当谐波带宽限制超过 L ≈ 32 时，正向变换会变得不稳定？
• RQ4 尽管正向变换不稳定，反向变换是否仍具有数值稳定性？
• RQ5 是否可通过正则化或重新表述病态线性系统来恢复数值稳定性？

主要发现

• 通过依赖环面上的香农采样定理，所提算法对带限函数实现了理论精确性。
• 正向与反向变换的复杂度均为 O(L³)，使该方法在任意旋量数下计算高效。
• 数值实验验证了计算时间的 O(L³) 缩放特性，与理论复杂度一致。
• 正向变换在 L > 32 时被发现不稳定，L=50、m₀=38 时条件数高达 6×10¹⁴，表明存在严重病态性。
• 反向变换预期稳定，但因缺乏稳定正向变换以实现精确恢复，无法验证其稳定性。
• 目前尝试通过数值方法稳定系统尚未成功，正探索如几何代数等替代公式化方法以期提升稳定性。