什么是艾里函数数值表计算器?
本工具可在一段实数 x 区间上计算两个艾里函数 Ai(x) 与 Bi(x),并可选择同时给出它们的导数 Ai'(x) 和 Bi'(x)。艾里函数是艾里微分方程 y'' − x·y = 0 的两个线性无关解,在物理学中应用极广:量子力学中它们描述经典转折点附近的波函数,此外还出现在光学、渐近分析以及彩虹理论等领域。
如何使用
输入起始 x 值、终止 x 值和步长,计算器会从 x 起点到 x 终点(含端点)按步长逐点生成一行数据。勾选导数选项即可同时列出 Ai'(x) 与 Bi'(x)。图像将 Ai(x) 和 Bi(x) 随 x 的变化绘制出来,让你直观看到:x 为正时 Ai 衰减,而 x 为负时两个函数都呈现振荡。
计算公式
采用在原点附近展开的级数,其中 alpha = Ai(0) = 0.3550280539,beta = −Ai'(0) = 0.2588194038:
Ai(x) = alpha·f(x) − beta·g(x),Bi(x) = sqrt(3)·(alpha·f(x) + beta·g(x)),其中 f(x) = 1 + x³/6 + x⁶/180 + …,g(x) = x + x⁴/12 + x⁷/504 + …。当 |x| 大于约 8 时,计算器会改用以 zeta = (2/3)|x|^(3/2) 表示的渐近公式,以避免相消误差。
计算示例
当 x = 0 时:f(0)=1,g(0)=0,因此 Ai(0) = alpha = 0.3550281,Bi(0) = sqrt(3)·alpha = 0.6149266。当 x = 1 时:f(1) ≈ 1.1722535,g(1) ≈ 1.0853407,于是 Ai(1) ≈ 0.1352924,Bi(1) ≈ 1.2074236,与标准函数表中的数值一致。
定义和术语表
- 第一类艾里函数,\(\text{Ai}(x)\)
- 艾里方程的解,当 \(x \to +\infty\) 时衰减到零。对于较大的正 \(x\),其衰减形式为 \(\dfrac{e^{-\zeta}}{2\sqrt{\pi}\,x^{1/4}}\);对于负 \(x\),它的振荡波长增长缓慢。
- 第二类艾里函数,\(\text{Bi}(x)\)
- 第二个线性无关的解。当 \(x \to +\infty\) 时,增长形式为 \(\dfrac{e^{\zeta}}{\sqrt{\pi}\,x^{1/4}}\),且与艾里函数相同,在 \(x<0\) 时振荡。
- 艾里微分方程,\(y'' - xy = 0\)
- 最简单的具有在原点处转折点的二阶线性常微分方程。其通解为 \(y(x) = c_1\,\text{Ai}(x) + c_2\,\text{Bi}(x)\)。它出现在光学、量子力学(线性势中的粒子)和波问题的 WKB 分析中。
- \(\zeta = \tfrac{2}{3}|x|^{3/2}\)
- 艾里函数的自然相位/衰减变量。它控制 \(x>0\) 时的指数增长和衰减,以及 \(x<0\) 时的振荡相位,在渐近展开式中出现。
- 转折点
- 使方程行为特性改变的 \(x\) 值。对于 \(y'' - xy = 0\),转折点在 \(x=0\):当 \(x<0\) 时解是振荡的(其中系数 \(-x\) 为正),当 \(x>0\) 时是指数型的(增长或衰减)。
- 渐近展开
- 一个关于 \(\zeta\)(或 \(x^{3/2}\))的倒数幂级数,可精确近似大 \(|x|\) 时的艾里函数。它可能不收敛,但少数几项对远离原点的地方提供了极好的精度,而公式选项卡的幂级数在那里收敛缓慢。
- 朗斯基行列式
- 行列式 \(W = \text{Ai}(x)\,\text{Bi}'(x) - \text{Ai}'(x)\,\text{Bi}(x)\)。非零常数朗斯基行列式(这里为 \(1/\pi\))确认艾里函数和伯麦函数线性无关,因此构成完整的解基。
常见问题
为什么 Bi(x) 会迅速发散? 当 x 取较大正值时,Bi(x) 按 exp(zeta) 的速度增长,在 x 超过约 230 时会超出双精度浮点数的表示范围而溢出。建议把上限设得小一些。
为什么 x 为负时函数会来回摆动? 当 x 趋向负无穷时,两个函数都会振荡,振幅按 |x|^(−1/4) 逐渐衰减。
使用什么单位? 没有单位——x 是纯实数,输出结果也都是无量纲的数值。