织构分析：织构数学推导
广义球函数不再跳步：从旋转球谐基到 Bunge 的 ODF 展开
把广义球函数理解成旋转作用于球谐基时的矩阵元素，从而推出乘法、正交、系数提取和傅里叶计算。
先记住这三点
- 球面函数处理两个自由度，取向分布需要 SO(3) 上三个自由度的基。
- 广义球函数可以由“旋转一个球谐函数会混合出哪些同阶球谐函数”来构造。
- 本章采用实正交基的旋转矩阵元素，避免把不同文献的相位和归一系数混用。
为什么不能直接给极图再添一个普通傅里叶变量
完整取向 R 同时决定晶体轴的朝向和绕轴自转。它不是三个相互独立的平移坐标；旋转复合的顺序也不能交换。若在三个欧拉角上直接作普通长方体傅里叶级数，会忽略 sinΦ 权重、坐标边界识别与旋转对称。
因此要建立与旋转本身相容的基。约定 R 把晶体正交坐标中的向量变到试样坐标，Haar 测度为 dμ=sinΦdφ₁dΦdφ₂/(8π²)，完整空间积分为 1。毛书的 g 与本站 R 的映射方向可能相反；下面先用明确的 R 和实基从头建立公式，而不直接复制书中的复相位。
第一步：旋转球面上的一个函数
取同一阶 l 的 2l+1 个实正交球谐函数 Kₗ₁,…,Kₗ,₂ₗ₊₁，满足 ∫KₗₐKₖᵦdΩ=δₗₖδₐᵦ。空间旋转不会改变球面 Laplace 本征值，所以 Kₗₐ(Rh) 仍能由同一阶函数展开：Kₗₐ(Rh)=ΣᵦDˡₐᵦ(R)Kₗᵦ(h)。
两边乘 Kₗ𝚌(h) 并积分，立即得到 Dˡₐ𝚌(R)=∫Kₗₐ(Rh)Kₗ𝚌(h)dΩ。右侧每一项都已知，所以这也是一条可执行的定义：给定旋转，旋转多项式的自变量，再向原来的基投影。Dˡ 的每个矩阵元素都是 R 的函数，也就是一种广义球函数。
最小实例：l=0 与 l=1
l=0 只有 K₀₁=1/√(4π)，旋转它不变，因此 D⁰(R)=1。l=1 可取 K₁=√(3/4π)(x,y,z)ᵀ；旋转后 K₁(Rh)=R K₁(h)，所以 D¹(R) 就是普通 3×3 旋转矩阵。广义球函数并不是凭空出现的新对象，最低非平凡阶就是熟悉的方向余弦。
以绕 z 轴转 α 为例，D¹=[[cosα,−sinα,0],[sinα,cosα,0],[0,0,1]]。这已经包含常数、正弦和余弦。到 l=2，可选 x²−y²、2xy、xz、yz、3z²−1 的适当归一化组合，旋转代入这些多项式，就能逐项算出 5×5 矩阵。
跟着二阶例子看见双倍角
令 q₁=x²−y²、q₂=2xy。绕 z 转 α 后 x′=xcosα−ysinα，y′=xsinα+ycosα。展开平方并合并，得到 q₁′=q₁cos2α−q₂sin2α、q₂′=q₁sin2α+q₂cos2α。这说明二阶旋转基里自然出现频率 2α，而不是人为添上一个二倍频参数。
同一旋转对 xz、yz 仍产生 α 频率，对 3z²−1 保持不变。这样一个五维基恰好包含 0、±1、±2 的角频率。一般 l 阶对应 −l 至 l 的频率，这就是计算广义球函数时有限傅里叶展开的来源。
加法定理其实是矩阵乘法
连续施加 R₂、R₁，有 K(R₁R₂h)=D(R₁)D(R₂)K(h)。另一方面，按照定义它也等于 D(R₁R₂)K(h)。比较各个独立基的系数，得到 Dˡ(R₁R₂)=Dˡ(R₁)Dˡ(R₂)，即 Dˡₐᵦ(R₁R₂)=Σ𝚌Dˡₐ𝚌(R₁)Dˡ𝚌ᵦ(R₂)。
毛书式(2-34)的广义球函数加法定理正是这种关系在另一套基下的写法。它让“把一个完整旋转拆成几步”成为处理极图积分与对称操作的工具。单位旋转给 D(I)=I，反向旋转给 D(R⁻¹)=D(R)ᵀ，都是同一乘法关系的结果。
欧拉角表达与广义连带勒让德函数
先固定记号：复球谐采用 Yₗᵐ∝e^(imψ)，并继续按 Y(Rh)=Dʸ(R)Y(h) 定义旋转矩阵，则主动绕 z 旋转 α 给 Dʸₘₙ(Rz)=δₘₙe^(imα)。令 dʸ(β)=Dʸ(Ry(β))，在 z–y–z 分解下便有 Dʸₘₙ=e^(imα)dʸₘₙ(β)e^(inγ)。这就是两个简单复指数乘一个广义连带勒让德函数的结构。常见 Wigner 记号另取 Dstd=(Dʸ)*，写作 e^(−imα)dstdₘₙ(β)e^(−inγ)；负号属于共轭约定，不能直接与前式混用。
Bunge 常用 z–x–z。对本站主动旋转，Rx(Φ)=Rz(−π/2)Ry(Φ)Rz(π/2)，所以可令 α=φ₁−π/2、β=Φ、γ=φ₂+π/2 后转换到 z–y–z 计算。不能把同一个三元数组同时当成两套欧拉角。毛书式(2-20)的 T 也需要按其基与映射方向一起使用；前文实基 D 和这里复基 Dʸ 则由基变换连接。
广义连带函数为什么也可作有限傅里叶计算
任何固定轴旋转都可写 Q Rz(β)Q⁻¹，其中 Q 把 z 轴送到目标轴；作用于向量时从右向左，先用 Q⁻¹ 把目标轴变到 z，再转 β，最后变回去。代入同一复基约定的加法定理，得到 Dʸₘₙ(QRz(β)Q⁻¹)=ΣₖDʸₘₖ(Q)e^(ikβ)Dʸₖₙ(Q⁻¹)，k=−l,…,l。Q 固定时，前后两个矩阵的元素都是常数，中间只剩简单指数。
这就是毛书 §2-5 将广义连带勒让德函数展开为有限傅里叶和的核心理由。实数基中把 ±k 配对便得到 sin(kβ)、cos(kβ)。计算时可预先保存 Q 对应的系数，再为许多取向重复调用；附录中的系数表服务于这个步骤，而不是额外的物理参数。
从旋转基到 ODF 展开系数
旋转矩阵元素的正交关系为 ∫Dˡₐᵦ(R)Dᵏ𝚌𝚍(R)dμ=δₗₖδₐ𝚌δᵦ𝚍/(2l+1)。这条关系是旋转群正交定理；Haar 测度保证对所有方向一视同仁，同阶矩阵维数 2l+1 给出分母。于是写 f(R)=ΣₗₐᵦCˡₐᵦDˡₐᵦ(R)，乘相应基并积分，就得 Cˡₐᵦ=(2l+1)∫fDˡₐᵦdμ。
l=0 时 C⁰=∫f dμ=1。若已有 EBSD 的完整取向 Rᵢ 和权重 wᵢ，可直接估计 Cˡₐᵦ≈(2l+1)ΣwᵢDˡₐᵦ(Rᵢ)/Σwᵢ。这一步不需要先把数据压缩成极图，再从极图反演回来；但仍需处理有限样本、核平滑与晶体对称。
一个完全可以手算的取向分布
忽略材料对称，取教学 ODF f(R)=1+aR₃₃，|a|≤1。因为随机旋转中 R₃₃ 的均值为零，分布已归一；R₃₃∈[−1,1] 保证非负。在 l=1 的上述基下，C¹₃₃=a，其余非零常数项是 C⁰=1。
验证系数：∫R₃₃²dμ=1/3，因此按提系数公式得到 3∫(1+aR₃₃)R₃₃dμ=a。它表示有方向的晶体 z 轴偏向试样 +z 的程度。普通反极等价的晶面极图会丢失这个奇阶信息；这个例子将在下一章用于精确展示零空间，不能冒充立方金属的实际 ODF。
系数、织构指数和分辨率怎样对应
用正交关系计算 J=∫f²dμ，可得 J=Σₗₐᵦ|Cˡₐᵦ|²/(2l+1)。常数项贡献 1，其余项衡量偏离随机分布的程度。上面的例子 J=1+a²/3。它给出可以逐项核算的“集中程度”，并不识别峰的位置与力学效应。
实际计算应记录基函数约定、L、平滑策略和对称处理。高阶项可以表达更细的角变化，也更容易放大测量不足和噪声。毛书的 L=22 是其历史计算方案的取舍；对尖锐单晶近邻峰、细窄纤维或高角分辨 EBSD，不能直接沿用其精度结论。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理，不复制受版权保护的标准正文。