什么是连带勒让德多项式表计算器?
本工具用于计算连带勒让德函数 \(P_n^m(x)\)(阶数 \(n\)、次数 \(m\))在指定 \(x\) 区间内的取值表,并绘制对应曲线。它属于纯数学范畴,在任何地区结果完全一致,不涉及任何单位或与国家相关的设定。连带勒让德多项式在物理学和应用数学中无处不在:球谐函数、球坐标系下拉普拉斯方程的求解、多极展开,以及角动量的量子力学等都会用到它。
使用方法
输入整数阶数 \(n\)(0、1、2……)和整数次数 \(m\),要求满足 \(-n \le m \le n\)。再设定 \(x\) 的起始值(介于 \(-1\) 与 \(1\) 之间)、步长增量以及行数。默认参数为 \(n = 2\)、\(m = 1\)、起始值 \(= -1\)、步长 \(= 0.02\)、共 101 行,即让 \(x\) 从 \(-1\) 扫到 \(+1\)(含端点)。还可选择 Type A(Wolfram 约定)或 Type B(Maple 约定);对于区间 \((-1, 1)\) 内的实数 \(x\),两者数值大小相同,仅在前置系数的符号/相位上有差异。
公式详解
对于整数 \(n\) 且 \(0 \le m \le n\),我们采用 $$P_n^m(x) = (-1)^m(1-x^2)^{m/2}\,\frac{d^m}{dx^m}P_n(x)$$ 并用数值稳定的递推关系来求值: $$P_m^m = (-1)^m(2m-1)!!\,(1-x^2)^{m/2}$$ $$P_{m+1}^m = x(2m+1)P_m^m$$ 随后 $$(l-m)P_l^m = (2l-1)x\,P_{l-1}^m - (l+m-1)P_{l-2}^m$$ 对于负的 \(m\),则有 $$P_n^{-m} = (-1)^m\frac{(n-m)!}{(n+m)!}P_n^m$$ 采用这种封闭式递推,可以避免对正整数 \(m\) 直接套用 \({}_2F_1\) 形式时出现的 Gamma 函数发散问题。
实例演算
取 \(n = 2\)、\(m = 1\) 时,函数为 $$P_2^1(x) = -3x\sqrt{1-x^2}$$ 在 \(x = 0\) 处其值为 0;在 \(x = 0.5\) 处为 \(-3(0.5)(0.866025) = -1.299038\);在 \(x = -0.5\) 处为 \(+1.299038\)。曲线从 \(x = -1\) 处的 0 出发,在 \(x = -0.577\) 附近上升到约 \(+1.1547\),在 \(x = 0\) 处穿越零点,在 \(x = +0.577\) 附近下探到约 \(-1.1547\),最后在 \(x = +1\) 处回到 0。
闭式伴随勒让德函数 P_n^m(x)
整数次数 \(n\) 和阶数 \(0\le m\le n\) 的伴随勒让德函数 \(P_n^m(x)\) 由 \(P_n^m(x)=(-1)^m(1-x^2)^{m/2}\dfrac{d^m}{dx^m}P_n(x)\) 定义。因子 \((-1)^m\) 是 A型约定中包含的Condon–Shortley相位(与Wolfram相匹配);B型约定(Maple)省略了它,因此其奇数 \(m\) 项仅在符号上有所不同。下表列出了A型下的显式形式。
| \(n\) | \(m\) | \(P_n^m(x)\)(A型,带符号) |
|---|---|---|
| 0 | 0 | \(1\) |
| 1 | 0 | \(x\) |
| 1 | 1 | \(-\sqrt{1-x^2}\) |
| 2 | 0 | \(\tfrac{1}{2}(3x^2-1)\) |
| 2 | 1 | \(-3x\sqrt{1-x^2}\) |
| 2 | 2 | \(3(1-x^2)\) |
| 3 | 0 | \(\tfrac{1}{2}(5x^3-3x)\) |
| 3 | 1 | \(-\tfrac{3}{2}(5x^2-1)\sqrt{1-x^2}\) |
| 3 | 2 | \(15x(1-x^2)\) |
| 3 | 3 | \(-15(1-x^2)^{3/2}\) |
作为一个完成的检查,在 \(x=0.5\) 处,条目 \(P_2^1\) 给出 \(-3(0.5)\sqrt{1-0.25}=-1.5\sqrt{0.75}=\) -1.299038。\(m=0\) 列再现了普通勒让德多项式 \(P_n(x)\),例如 \(P_3^0(x)=\tfrac12(5x^3-3x)\),可以用勒让德多项式表计算器制成表格。
关键术语和变量
- 次数 \(n\)
-
一个非负整数(
degreeN),设置底层勒让德多项式 \(P_n(x)\) 的阶数,该多项式是次数为 \(n\) 的多项式。 - 阶数 \(m\)
-
一个整数(
orderM),控制取多少个导数。对于在 \((-1,1)\) 上的实值结果,通常使用 \(0\le m\le n\);当 \(m>n\) 时,函数恒等于零,因为次数为 \(n\) 的多项式的 \(m\) 阶导数消失。 - 自变量 \(x\)
-
计算点(
initialX加 \(i\cdot\)stepX)。这些函数在 \(-1\le x\le 1\) 上为实数;在物理中 \(x=\cos\theta\)。 - A型(Wolfram / Condon–Shortley)
-
包含相位因子 \((-1)^m\)。这是Wolfram的
LegendreP和标准量子力学教科书使用的约定。 - B型(Maple)
- 省略 \((-1)^m\) 相位,因此 \(P_n^m\)(B型)\(=(-1)^m\,P_n^m\)(A型)。大小相同;仅奇数 \(m\) 项的符号不同。
- 双阶乘 \((2m-1)!!\)
- 奇数整数的乘积 \((2m-1)(2m-3)\cdots 3\cdot 1\),其中 \((-1)!!=1\)。它出现在首项系数 \(P_m^m(x)=(-1)^m(2m-1)!!\,(1-x^2)^{m/2}\) 中;例如 \(P_3^3\) 使用 \(5!!=15\)。查看双阶乘计算器以获得这些值。
- 负阶关系
- 对于 \(m>0\),\(P_n^{-m}(x)=(-1)^m\dfrac{(n-m)!}{(n+m)!}\,P_n^{m}(x)\),通过阶乘连接正负阶。
解释表格和图形
几个结构性质让你可以理智地检查制表值和绘制的曲线:
- 奇偶性。 \(P_n^m(-x)=(-1)^{n+m}P_n^m(x)\)。当 \(n+m\) 为偶数时,图形关于 \(x=0\) 对称;当 \(n+m\) 为奇数时,它是反对称的(因此通过原点)。
- 内部零点。在开区间 \((-1,1)\) 上,\(P_n^m(x)\) 恰好有 \(n-m\) 个简单零点。例如 \(P_3^1\) 有两个内部零点,而 \(P_n^n\) 没有。
- 端点行为。由于因子 \((1-x^2)^{m/2}\),每个具有 \(m>0\) 的函数在 \(x=\pm 1\) 处消失。对于 \(m=0\),值为 \(P_n(1)=1\) 和 \(P_n(-1)=(-1)^n\)。
- 靠近边缘的大小。对于较高的 \(m\),\((1-x^2)^{m/2}\) 因子在 \(x\to\pm1\) 时急剧抑制曲线,因此最大偏差出现在范围的中间附近。
这些函数是球面谐波 \(Y_n^m(\theta,\phi)\) 的 \(\theta\) 依赖部分:写作 \(x=\cos\theta\),我们有 \(Y_n^m\propto P_n^m(\cos\theta)\,e^{im\phi}\)。内部零点变成纬度的球面结点,\(m>0\) 的端点消失对应于在两极倾向于零的谐波。因此相同的 \(P_n^m\) 值直接馈入在选定的 \(\theta\) 和 \(\phi\) 处的球面谐波计算。
常见问题
为什么 n 和 m 必须是整数?多项式能够终止(有限项)的形式要求 \(n\) 为非负整数;而递推关系和 \((n\pm m)!\) 因子则要求 \(m\) 为整数,并满足 \(-n \le m \le n\)。
顶部显示的样本值是什么?顶部信息框会显示取值表中间一行(中位索引)对应的 \(x\) 与 \(P_n^m(x)\),可用于快速核对曲线。
m = 0 时得到什么?此时 \(P_n^0(x)\) 就是普通的勒让德多项式 \(P_n(x)\)。