中文

球高斯-拉盖尔基函数的快速傅里叶变换

数值分析 2016-12-01 v3

摘要

球高斯-拉盖尔(SGL)基函数,即类型 Lnl1(l+1/2)(r2)rlYlm(ϑ,φ)L_{n-l-1}^{(l + 1/2)} (r^2) r^{l} Y_{lm}(\vartheta,\varphi) 的归一化函数,ml<nN|m| \leq l < n \in \mathbb{N}Lnl1(l+1/2)L_{n-l-1}^{(l + 1/2)} 为广义拉盖尔多项式,YlmY_{lm} 为球谐函数,构成了带高斯权重 exp(r2)\exp(-r^{2})R3\mathbb{R}^{3}L2L^{2} 空间的一组正交归一基。这些基函数被广泛使用,例如在生物分子动力学模拟中。然而,迄今为止还没有可靠的算法能快速计算函数关于SGL基函数的傅里叶系数。本文提出了此类广义FFT。我们从SGL采样定理出发,该定理允许精确计算带限函数的SGL傅里叶展开。通过分离变量方法并采用快速球傅里叶变换,我们揭示了一类通用的快速SGL傅里叶变换。所有这些算法的渐近复杂度为 O(B4)\mathcal{O}(B^{4})BB 为相应带宽,而 R3\mathbb{R}^{3} 上的采样点数量随 B3B^{3} 缩放。这明显改进了朴素边界 O(B7)\mathcal{O}(B^{7})。同时,我们的方法得到了与前向变换渐近复杂度相同的快速逆变换。我们在数值实验中证明了算法的实际适用性。值得注意的是,这是非紧域上首批广义FFT的实现之一。我们以讨论作结,包括真正 O(B3log2B)\mathcal{O}(B^{3} \log^{2} B) 快速SGL傅里叶变换及逆变换的布局,以及对未来发展的展望。

关键词

引用

@article{arxiv.1604.05140,
  title  = {Fast Fourier Transforms for Spherical Gauss-Laguerre Basis Functions},
  author = {Jürgen Prestin and Christian Wülker},
  journal= {arXiv preprint arXiv:1604.05140},
  year   = {2016}
}