[论文解读] Programs in Mathematica relevant to Phase Integral Approximation for coupled ODEs of the Schrödinger type

A. A. Skorupski|ArXiv.org|Oct 29, 2007

Numerical methods for differential equations被引用 3

一句话总结

本文提出三个 Mathematica 程序——Y2n、bvm 和 Ymsvm——用于自动计算一个或两个耦合的薛定谔型常微分方程的相位积分近似（PIA）修正。这些工具通过直接求和递推关系，生成高阶修正量 $Y_{2n}(x)$、向量 $<\mathbf{b}_m(x)\u0026gt;$ 和辅助量的符号表达式，以确保准确性并最大限度减少编程错误。

ABSTRACT

Three programs in Mathematica are presented, which produce expressions for the lowest order and the higher order corrections of the Phase Integral Approximation. First program is pertinent to one ordinary differential equation of the Schrödinger type. The remaining two refer to a set of two such equations.

研究动机与目标

• 自动化计算薛定谔型常微分方程的相位积分近似（PIA）修正的符号表达式。
• 为标量系统和耦合系统中的高阶修正提供可靠且错误最少的代码。
• 支持一般和特定的输入情况，包括三角函数和复数势能。
• 通过一致性检查（如奇数阶修正为零）实现验证。
• 通过自定义输出格式（TeX、Fortran、OutputForm）支持数值计算和结果简化。

提出的方法

• 实现三个 Mathematica 程序：Y2n 用于标量方程，bvm 和 Ymsvm 用于耦合系统。
• 使用直接嵌套的 Do 循环计算来自递推关系（如式 40、54）的多重求和，以保持数学清晰性。
• 允许通过外部 .dat 文件或默认参数输入用户数据，并通过交互式对话框选择输出格式和简化选项。
• 支持符号运算，包括 Expand、Simplify、Together 和 Apart，对三角函数可选使用 Trig 展开。
• 通过避免过早简化积分（如式 105）来保持完整的符号完整性，优先保证正确性而非速度。
• 支持追加自定义脚本（如 ap.mat）以实现目标数值计算和结果验证。

实验结果

研究问题

• RQ1 如何以高可靠性自动化计算耦合常微分方程的相位积分修正的符号表达式？
• RQ2 多重求和递推关系在计算高阶修正量 $Y_{2n}(x)$ 及其相关向量中的作用是什么？
• RQ3 如何通过一致性检查（如奇数阶修正为零）验证实现的正确性？
• RQ4 程序如何确保同一问题的厄米特与非厄米特形式之间结果的一致性？
• RQ5 如何灵活地将用户自定义输入和输出格式集成到符号计算工作流中？

主要发现

• 程序成功计算出 $Y_{2n}(x)$、$<\mathbf{b}_m(x)\u0026gt;$ 和 $<\mathbf{s}_m(x)\u0026gt;$ 的符号表达式，结果可保存为多种格式。
• 通过大量项的相互抵消，确认了奇数阶修正量的消失（$Y_{2n-1}(x) \equiv 0$），验证了代码的正确性。
• 在厄米特情形下，简化与非简化形式的计算结果完全一致，确认了结果的一致性。
• 使用直接求和（如嵌套 Do 循环）最大限度减少了编程错误，并保持了数学保真度。
• 通过追加脚本（如 ap.mat）进行的数值计算，与参考论文表 I 中的结果一致，验证了结果的可靠性。
• 该框架支持复数势能和三角函数势能（如示例 C.4 中的 $\cos \phi \pm i \sin \phi$），具有广泛的应用潜力。