[论文解读] Mathieu functions computational toolbox implemented in Matlab

E. Cojocaru|ArXiv.org|Nov 12, 2008

Advanced Theoretical and Applied Studies in Material Sciences and Geometry被引用 7

一句话总结

本文提出了一套自包含的Matlab计算工具箱，用于计算Mathieu函数，这些函数在椭圆柱坐标系中求解波动方程时至关重要。该工具箱提供了显式公式、归一化符号、展开系数的递推关系，以及基于贝塞尔函数的级数表达式，适用于角向和径向Mathieu函数，使科学领域内无需深入了解专业文献即可实现准确且便捷的计算。

ABSTRACT

The Mathieu functions are used to solve analytically some problems in elliptical cylinder coordinates. A computational toolbox was implemented in Matlab. Since the notation and normalization for Mathieu functions vary in the literature, we have included sufficient material to make this presentation self contained. Thus, all formulas required to get the Mathieu functions are given explicitly. Following the outlines in this presentation, the Mathieu functions could be readily implemented in other computer programs and used in different domains. Tables of numerical values are provided.

研究动机与目标

• 解决科学软件中计算Mathieu函数缺乏可访问、标准化算法的问题。
• 简化并统一文献中用于Mathieu函数的不一致符号和归一化方式。
• 提供一个完整、自包含的Matlab实现，使角向和径向Mathieu函数的计算更加准确，适用于波动传播和散射等应用。
• 通过清晰列出所有必要的数学公式和递推关系，促进该方法在其他编程环境中的应用。

提出的方法

• 通过变量分离法，从椭圆柱坐标系中的标量波动方程推导出角向和径向Mathieu方程。
• 利用递推关系计算四种对称类型（偶-偶、偶-奇、奇-偶、奇-奇）的角向Mathieu函数的展开系数。
• 通过在特定点处的角函数值导出的连接因子，将第一类径向Mathieu函数表示为贝塞尔函数的组合。
• 在u=0和v=π/2处应用归一化条件，将角函数的值和导数与径向函数系数关联起来。
• 采用包含贝塞尔函数J_j(v₁)和J_j(v₂)的级数展开，其中v₁ = √q exp(−u) 且 v₂ = √q exp(u)，以实现快速收敛的数值计算。
• 基于其傅里叶级数表示，提供角向和径向Mathieu函数导数的显式公式。

实验结果

研究问题

• RQ1 如何在不同科学应用中，通过统一、自包含的符号体系，实现Mathieu函数的精确且高效计算？
• RQ2 计算角向和径向Mathieu函数的展开系数及连接因子，需要哪些递推关系和归一化条件？
• RQ3 如何将第一类径向Mathieu函数表示为贝塞尔函数的组合，以确保快速且稳定的收敛？
• RQ4 角向和径向Mathieu函数在全部四种对称类中的导数，其正确的解析表达式是什么？
• RQ5 该方法在多大程度上可以推广并实现在Matlab以外的编程语言中？

主要发现

• 该工具箱为所有四种对称类型的角向Mathieu函数的展开系数A_{pm}^{(m)}(q,n)提供了显式递推关系。
• 通过从角函数展开中导出的系数，利用贝塞尔函数级数，准确计算了第一类径向Mathieu函数。
• 连接因子g_{pm}(q,n)由角函数在v=π/2处的值和导数导出，实现了不同函数类型间的一致归一化。
• 角向和径向Mathieu函数的导数以傅里叶级数形式表达，其系数与原始展开系数成正比。
• 该方法在|u|较小时确保快速收敛，并为波散射和椭圆柱问题中的Mathieu函数计算提供了数值稳定的框架。
• 整个公式体系自包含，所有必要的公式和归一化惯例均明确列出，可直接在任何计算环境中实现。