织构分析：织构数学推导
球谐展开从头推：一张极图怎样变成一组系数
从圆周傅里叶级数走到球面上的勒让德方程，亲手求出最低几阶、系数、归一化和截断误差。
先记住这三点
- 球谐函数是球面上的正交基；极密度不必自身是调和函数。
- 必须把积分面积元、基函数归一化与系数定义配套使用。
- 展开阶数控制可表示的角尺度，不能把某个历史常用阶数视为普遍精度保证。
先从绕圆周的一个变量理解展开
如果沿某个纬圈测得极密度随方位角 ψ 的变化，因为 ψ 与 ψ+2π 是同一个方向，可以用 1、cosψ、sinψ、cos2ψ、sin2ψ……叠加。每一项描述一种起伏频率，系数描述它在实际信号中有多强。球谐展开保留这个思路，但还必须处理从北极到南极的倾角 θ。
单位方向写成 y=(sinθcosψ,sinθsinψ,cosθ)。球面面积元为 dΩ=sinθdθdψ，不能按等 θ 网格作普通平均。下面先在完整球面上推导；采用半球和正负极合并时，应先说明反极对称，再把积分范围与因子一致地改写。
球面微分算子怎样把两个角分开
球面 Laplace 算子为 ΔS=(1/sinθ)∂θ(sinθ∂θ)+(1/sin²θ)∂²ψ。寻找它的本征函数：ΔS Y=−λY。令 Y(θ,ψ)=Θ(θ)Ψ(ψ)，除以 ΘΨ 并分离变量，方位角部分得到 Ψ″+m²Ψ=0。周期条件 Ψ(ψ+2π)=Ψ(ψ) 要求整数 m，解可以选 e^(imψ) 或相应的正弦、余弦。
倾角部分于是成为 (1/sinθ)d/dθ(sinθ dΘ/dθ)+[λ−m²/sin²θ]Θ=0。令 x=cosθ，则 d/dθ=−sinθ d/dx，整理后得到 (1−x²)Θ″−2xΘ′+[λ−m²/(1−x²)]Θ=0。这样，一个看起来复杂的球面问题已经化成一条常微分方程。
勒让德多项式与连带勒让德函数
要求解在两极都正常，得到 λ=l(l+1)，l 为非负整数，|m|≤l。先取 m=0，解是勒让德多项式 Pₗ(x)。由 Rodrigues 公式 Pₗ(x)=[1/(2ˡl!)]dˡ(x²−1)ˡ/dxˡ，可以算出 P₀=1、P₁=x、P₂=(3x²−1)/2、P₃=(5x³−3x)/2。这里的 l 是阶数，不是晶面指数。
对非负 m，采用含 Condon–Shortley 相位的约定，Pₗᵐ(x)=(−1)ᵐ(1−x²)^(m/2)dᵐPₗ/dxᵐ。例如 P₁¹=−√(1−x²)=−sinθ，P₂¹=−3x√(1−x²)。负 m 可由标准比例关系生成。不同书有时把 (−1)ᵐ 放到别处；系数符号随之变化，物理极图不会因此改变。
把“彼此独立”写成可以计算的正交关系
复球谐函数取 Yₗᵐ=NₗₘPₗᵐ(cosθ)e^(imψ)，Nₗₘ=√[(2l+1)(l−m)!/(4π(l+m)!) ]。方位角积分 ∫₀²πe^[i(m−n)ψ]dψ 在 m≠n 时为零；倾角部分的正交性再区分不同 l。组合起来得到 ∫YₗᵐYₖⁿ* dΩ=δₗₖδₘₙ，星号表示复共轭。
这也解释了归一系数的作用：让每个基函数的平方积分都是 1，方便直接提取系数。若使用实数表示，可用 Yₗ⁰、√2ReYₗᵐ、√2ImYₗᵐ 组成同样正交的基 Kₗₐ；a 在同一 l 下枚举 2l+1 个实基函数。后面的取向矩阵展开采用这个实基，免去共轭符号造成的干扰。
从完整极图真正求出系数
写 P(y)=ΣₗΣₘAₗₘYₗᵐ(y)。两边乘 Yₖⁿ*，再对整个球面积分，所有不同下标的项都因为正交而消失，只剩 Aₖₙ=∫P Yₖⁿ* dΩ。这不是猜测拟合公式，而是逐项投影得到的结果。若 P 是 m.r.d. 密度，∫P dΩ=4π，因此 A₀₀=√(4π)。
离散测量时，Aₗₘ≈ΣᵢPᵢYₗᵐ*(yᵢ)ΔΩᵢ；一个倾角、方位角网格单元的准确面积为 ΔΩᵢ=(cosθlow−cosθhigh)Δψ。极点处不要把 sin0=0 误解为整个第一圈没有面积。若只测了部分球面，以上“正交提系数”已不再成立，需要后面的带掩膜最小二乘。
纸面例子：只包含一个二阶起伏的方向密度
设教学分布 P(θ)=1+aP₂(cosθ)，它只用于演算，不指定某种晶体对称。因为 ∫₋₁¹P₂(x)dx=0，它自动满足平均值为 1。P₂ 的取值范围是 −1/2 至 1，要使 P 处处非负，必须有 −1≤a≤2。取 a=1，极点密度为 2、赤道为 0.5。
由于 Y₂⁰=√(5/4π)P₂，非零系数只有 A₀₀=√(4π)、A₂₀=a√(4π/5)。以完整球面为总体，算北极 θ≤30° 的极份额：F=(1/2)∫cos30°¹[1+a(3x²−1)/2]dx=(1−c)/2+a(c−c³)/4，c=cos30°。a=1 时 F≈0.121。若把相反极合并到半球，同一轴向锥的份额应合并南北两帽，为约 24.2%。两个数字都来自积分，不能把峰顶 2 m.r.d. 改写为 2%。
为什么展开可以有负系数，却不能有负极密度
基函数本身会正负起伏，负系数只表示某种模式的反向贡献，并不等于负体积。物理要求作用于最后的总和 P，而不是要求每个 A 都非负。反过来，即使原始 P 非负，截断后的有限级数也可能在峰旁振荡到负值。
若所有项都保留且函数平方可积，正交基可以按均方意义逼近它。只保留 l≤L 时，平方误差等于被舍掉系数的平方和：∫|P−Pᴸ|²dΩ=Σₗ>LΣₘ|Aₗₘ|²。这是 Parseval 关系的直接结果，告诉我们误差来自哪些角频率，而不是保证任意材料在 L=22 时都足够精确。
纠正一个容易被旧教材表述带偏的地方
毛书印刷 pp.18–20 用 Laplace 方程引出球谐函数，便于建立数学联系。但不能据此说任意实测极密度都满足 ΔS P=0。在完整球面上，光滑且 ΔS P=0 的函数只能是常数：对等式乘 P 并积分，可得 ∫|∇S P|²dΩ=0。明显有峰的极图不满足这个结论。
正确理解是：不同 l 的球谐函数分别是 ΔS 的本征函数，具有不同本征值；它们共同构成展开一般球面函数的基。对于三维调和函数 rˡYₗᵐ，径向变化恰好抵消角向项，但这不要求极密度的球面分布本身无起伏。
把这套数学用到真实研究的顺序
毛书的工业纯铝轧制案例先测多个不完整极图，再计算 ODF，最后比较 β 纤维与组分体积分数。球谐系数在这里是连接测量与取向分布的中间量，并不是最后要解释的材料机制。初学者可先用本章二阶教学分布手算系数，再理解案例为什么不能只比较单个衍射峰。
处理实际数据时，先校正背景和几何响应，确认对称，再选择基与截断阶；最后回算已测极图并比较残差。若峰旁出现负波纹，检查截断、覆盖和误差；不要把每条等高线都解释成一个新织构组分。真实案例的材料结论仍应由取样、工艺和完整分析共同支持。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理，不复制受版权保护的标准正文。