[论文解读] Calculation of oblate spheroidal wave functions with complex argument

A. L. Van Buren|arXiv (Cornell University)|Aug 31, 2020

Electromagnetic Scattering and Analysis参考文献 10被引用 10

一句话总结

本文提出了一种用于计算具有复参数的扁球函数的鲁棒数值方法，扩展了先前针对实参数的研究。通过修改基于矩阵的特征值计算方法，并采用渐近逼近方法处理类似长球形的特征值，该方法在广泛范围的复参数下实现了高精度结果，即使在具有显著虚部的情况下，也仅使用双精度浮点数运算，通过一个免费提供的 Fortran 程序 c o b l f c 实现。

ABSTRACT

A previous article showed that alternative expressions for calculating oblate spheroidal radial functions of both kinds can provide accurate values over very large parameter ranges using double precision arithmetic, even where the traditional expressions fail. The size parameter c was assumed real. This paper considers the case where c = cr + ici is complex with an imaginary part ci often used to represent losses in wave behavior. The methods for c real modified to complex arithmetic work reasonably well as long as ci is very small. This paper describes the substantial changes necessary to obtain useful results for larger values of ci. It shows that accurate eigenvalues can usually be obtained even though the matrix methods used to obtain them for c real provide increasingly inaccurate values, primarily for those with relatively small magnitude, as ci increases. It also shows that some of the eigenvalues can be prolate-like with values that are well approximated using asymptotic estimates for prolate eigenvalues where c is replaced with -ic. A method to order the eigenvalues is presented. The modifications necessary to compute accurately the radial and angular functions for complex c are discussed. A resulting Fortran program coblfcn provides useful function values for a reasonably wide range of c, the order m, and the radial coordinate when using double precision arithmetic. The results can be improved by using quadruple precision for the Bouwkamp procedure to ensure accurate double precision eigenvalues. Further improvement is obtained using full quadruple precision. Coblfcn is freely available at www.mathieuandspheroidalwavefunctions.com

研究动机与目标

• 将扁球函数的精确计算扩展至复参数，特别是当虚部代表波损耗时。
• 克服传统矩阵方法在大小参数虚部增大时特征值计算失败的问题。
• 开发一种可靠的方法，使用双精度浮点数运算对复参数下的径向函数和角向函数进行排序与计算。
• 提供一个免费、数值稳定的 Fortran 实现（c o b l f c），以供科学与工程应用中的实际使用。
• 评估四精度浮点数精炼在复杂参数区域中提升特征值精度的效果。

提出的方法

• 将原本专为实参数设计的基于矩阵的特征值计算技术适配至复数运算，通过修改以保持数值稳定性。
• 利用渐近估计方法，通过将复参数 c 替换为 -i c，对行为类似长球形模式的特征值进行近似。
• 在 Bouwkamp 方法中引入四精度浮点数运算，以提升双精度特征值计算的精度。
• 基于模长和渐近行为开发一致的特征值排序方案，确保在全参数范围内正确标记特征值。
• 对径向函数和角向函数的递推关系及归一化过程进行适配，以在不损失精度的情况下处理复参数。
• 将所有组件集成至一个可投入生产的 Fortran 程序 c o b l f c 中，包含输入验证和稳健的错误处理机制，以适应实际应用。

实验结果

研究问题

• RQ1 当大小参数 c 的虚部显著时，基于矩阵的扁球函数特征值计算方法是否仍能保持精度？
• RQ2 在复参数 c 的区域中，如何识别并近似那些表现出长球形特征值行为的特征值？
• RQ3 当 c 为复数且 |c i| 较大时，对现有算法需进行哪些修改以维持数值精度？
• RQ4 在复参数区域中，四精度浮点数运算能在多大程度上提升双精度特征值的精度？
• RQ5 一个单一、稳定的 Fortran 实现能否可靠地在广泛复数 c 值范围内计算特征值和函数？

主要发现

• 对于实 c 的矩阵方法在 c 的虚部较小时仍可使用，但随着 |c i| 增大，必须进行大量修改才能保持精度。
• 使用标准矩阵方法时，模长较小的特征值在 c i 增大时变得越来越不准确，但所提出的修改可恢复精度。
• 部分复参数 c 下的特征值表现出类似长球形的特性，可利用针对参数 -i c 的长球形函数渐近估计方法进行良好近似。
• 所提出的特征值排序方法成功区分并正确排序了全复参数范围内的特征值，即使特征值模长发生交叉。
• Fortran 程序 c o b l f c 使用双精度浮点数运算，可在广泛 m、c（包括复数）及径向坐标范围内产生高精度函数值。
• 在 Bouwkamp 方法中使用四精度浮点数运算可显著提升特征值精度，若完全使用四精度，精度提升更为明显。