什么是球谐函数计算器？
球谐函数 \(Y_{n}^{m}(\theta,\phi)\) 是拉普拉斯方程在球坐标系下解的角向部分。它们在球面上构成一组完备的标准正交基，广泛出现在物理与工程的各个领域：量子力学中的原子轨道、电磁学与引力理论中的多极展开、宇宙学中的宇宙微波背景辐射、地球物理中的地磁场模型，以及计算机图形学中的光照建模等。本计算器可针对你指定的阶数 \(n\)、次数 \(m\)、极角（天顶角）\(\theta\) 和方位角 \(\phi\)，求出对应的复数球谐函数值。
使用方法
首先选择一种约定：A 型（Wolfram/Mathematica）采用相位因子 \(e^{i m \phi}\)；B 型（Maple）采用 \(e^{i m (\phi+\pi)}\)，两者相差一个 \((-1)^{m}\) 因子。接着输入满足 \(n \ge 0\) 的整数阶 \(n\)，以及满足 \(-n \le m \le n\) 的整数次 \(m\)。再按所选角度单位（度或弧度——两个角度使用同一单位）填入天顶角 \(\theta\) 和方位角 \(\phi\)。最后设定要显示的有效数字位数，即可读出结果的实部、虚部，以及模和相位。
计算公式
令 \(x = \cos\theta\)。归一化系数为 $$N = \sqrt{\frac{2n+1}{4\pi}\cdot\frac{(n-m)!}{(n+m)!}}.$$ 连带勒让德函数 \(P_{n}^{m}(x)\) 由递推关系构造，初始项为 \(P_{m}^{m}(x) = (-1)^{m}(2m-1)!!(1-x^{2})^{m/2}\)。于是 A 型为：$$Y = N \cdot P_{n}^{m}(x) \cdot \left(\cos(m\phi) + i\cdot\sin(m\phi)\right).$$ 负次数则使用 \(P_{n}^{-m}(x) = (-1)^{m}\frac{(n-m)!}{(n+m)!} P_{n}^{m}(x)\)。在极点处（\(x = \pm 1\)），只有 \(m = 0\) 时结果才非零。
计算示例
采用默认值与 A 型约定：\(n=2\)，\(m=1\)，\(\theta=45^\circ\)，\(\phi=0\)。此时 \(x = \cos 45^\circ = 0.7071068\)，\(\frac{2n+1}{4\pi} = \frac{5}{4\pi} = 0.3978874\)，\(\frac{(n-m)!}{(n+m)!} = \frac{1}{6}\)，故 \(N = \sqrt{0.0663146} = 0.257516\)。$$P_{2}^{1}(x) = -3x\sqrt{1-x^{2}} = -3\cdot 0.5 = -1.5.$$ 由于 \(\phi=0\)，相位因子为 1，因此 $$Y = 0.257516 \cdot (-1.5) = -0.386274 + 0\,i,$$ 模为 0.386274。
常见问题
结果为什么是复数？ 因子 \(e^{i m \phi}\) 使 \(Y\) 一般为复数；只有当 \(m\phi\) 是 \(\pi\) 的整数倍（例如 \(\phi=0\)）时，虚部才会消失。
为什么我的参考书里显示的是 +0.386274？ 有些资料省略了 Condon–Shortley 相位因子 \((-1)^{m}\)。而 Wolfram/A 型约定将该因子包含在 \(P_{n}^{m}\) 之内，因此得到的是负值。
极点处会发生什么？ 当 \(\theta=0\) 或 \(180^\circ\) 时，因子 \((1-x^{2})^{m/2}\) 为零，所以对任意 \(m \ne 0\) 都有 \(Y = 0\)；只有 \(m=0\) 的项得以保留。
