织构分析：织构数学推导
从 ODF 推到 PF 再反演：把积分、系数关系与鬼峰逐步接起来
推导极图的弱形式和球谐系数关系，用一个可手算的零空间反例区分数据约束、正值约束与模型补全。
先记住这三点
- PF 对绕极方向的自由旋转取平均，因而是 ODF 的投影。
- 普通反极对称极图消去奇阶信息，但实际完整 ODF 不必只有偶阶。
- 非负与组分模型能增加约束，却不会自动把未测信息变成唯一观测结果。
先给极图一个无需特殊函数的定义
固定晶体单位法向 h。每个取向 R 把它送到 y=Rh。若旋转空间中取向概率为 f(R)dμ，则这些方向在球面上的分布就是 PF。这里定义 Pₕ 为 m.r.d. 密度，因此方向落进小面积 dΩ 的概率是 Pₕ(y)dΩ/(4π)。
用任意观察函数 q(y) 求平均，两种描述必须相同：(1/4π)∫Pₕ(y)q(y)dΩ=∫f(R)q(Rh)dμ。这叫弱形式，但含义很朴素：先把所有晶粒投到球面再平均，与先在每个晶粒上算方向再平均，结果应该一致。后面的系数关系直接由它推出。
把自由旋转写出来就得到极图积分
设某个 R₀ 满足 R₀h=y。所有产生同一个方向的取向都可写成 R(γ)=Rot(y,γ)R₀，0≤γ<2π。它们组成一条绕固定极的旋转轨迹。沿这条轨迹的条件平均给 Pₕ(y)=(1/2π)∫₀²πf(R(γ))dγ。
若改变起点 R₀，只相当于给 γ 换一个零点，积分不变。若 f=1，则 Pₕ=1，归一化也正确。真实面族测量需对所有等价 h 求平均；常规反射再将 h 与 −h 的贡献合并，不能把单个有向法向的公式当成已经处理完面族与 Friedel 对称。
不用记忆常数，直接推导系数关系
采用上一章的实正交球谐 Kₗₐ 和旋转元素 Dˡₐᵦ，满足 Kₗₐ(Rh)=ΣᵦDˡₐᵦ(R)Kₗᵦ(h)。定义 PF 系数 Aₗₐ(h)=∫PₕKₗₐdΩ。在弱形式中取 q=Kₗₐ，得到 Aₗₐ=4πΣᵦKₗᵦ(h)∫fDˡₐᵦdμ。
再用 Cˡₐᵦ=(2l+1)∫fDˡₐᵦdμ，最终得到 Aₗₐ(h)=[4π/(2l+1)]ΣᵦCˡₐᵦKₗᵦ(h)。这就是毛书式(2-39)的实基版本。4π 来自球面 m.r.d. 定义，2l+1 来自旋转基正交性；换归一化后两个常数会随定义变化，不能孤立抄一个公式。
跟着推导 · 信息怎样丢失
测得的曲线相同，原分布可以不同
先用圆周上的简单模型理解“零空间”。仪器若只记录相反两个方向的平均，能否分辨隐藏的方向偏好？改变 a，看真实分布与测量曲线各怎样变化。
f(θ) = 1 + a cosθ，|a| ≤ 1 保证非负；∫₀²π f(θ)dθ/(2π) = 1
P(θ) = [f(θ)+f(θ+π)]/2 = 1，因为 cos(θ+π) = −cosθ
此处 f(60°) = 1.300，反向密度 = 0.700，平均仍为 1.000
增加采样点，或要求分布非负，能恢复 a 吗？
都不能。无论测多少角度，P 始终为 1；而区间 −1≤a≤1 内每个解都已非负。把 a 选成 0 可以是某种平滑或最小范数选择，但不能说测量证明真实 a=0。需要能区分相反方向的新观测，才能增加这项约束。
这是圆周平均算子的严格教学例，不是某晶体的完整 ODF。真实 PF–ODF 变换还含绕极点积分、SO(3) 基函数与晶体/试样对称；请接着读球谐与极图反演章节。
反演为什么变成一组小线性方程
测量法向 h₁,…,hJ 的 PF 后，可对每个 l、a 写 bⱼ=ΣᵦHⱼᵦcᵦ，其中 bⱼ=Aₗₐ(hⱼ)，Hⱼᵦ=[4π/(2l+1)]Kₗᵦ(hⱼ)，cᵦ=Cˡₐᵦ。若已采用晶体对称基，列数从 2l+1 减为 M(l)。
只有 H 满列秩，才可能在该截断模型下唯一解出这一组系数。无噪声、完整有向 PF 的一般情形至少需要相应数量的独立方向；对称性可减少未知数。但真实测量还有反极合并、未测区和误差，所以书中“两个完整极图可展开到 l=22”的计数结论不能当作任意实验只测两张就足够。
Friedel 合并怎样把奇阶项消掉
普通反极等价处理是 P̄ₕ(y)=[Pₕ(y)+Pₕ(−y)]/2。球谐有 Kₗₐ(−y)=(−1)ˡKₗₐ(y)，所以合并后的系数 Āₗₐ=[1+(−1)ˡ]Aₗₐ/2。偶阶保留，奇阶严格为零。这里零的是测量对这些分量的响应，不是证明材料里相关取向信息本来为零。
Friedel 近似来自通常非共振运动学衍射中实散射因子的结构因子关系 F(−H)=F(H)*，进而 |F(−H)|²=|F(H)|²。异常散射、极性测定及动力学效应的适用条件需要另议。因此，不宜像旧式简写那样把“任何衍射永远区分不了正负极”作为普遍定律。
一个严格的反例：极图相同，ODF 却不同
取不施加晶体对称的教学分布 f₊(R)=1+aR₃₃、f₋(R)=1−aR₃₃，0<|a|≤1。两者都非负、归一，而且不同。固定 h，条件 Rh=y 时，对剩余自转平均可得 ⟨R⟩=yhᵀ，因此 P±ₕ(y)=1±a y₃h₃。
把 y 和 −y 合并，两个分布的所有 P̄ₕ 都等于 1。也就是说，哪怕完整测量任意多张这种反极合并 PF，也分不开这两个分布。本例展示一般反演的零空间；立方或其他晶体对称会删去其中某些阶，但不能据此把普通极图反演的所有非唯一性都忽略。
鬼峰、截断振荡与噪声不是一个问题
零空间意味着有些不同 ODF 在测量上完全等价；只保留偶阶信息得到的传统简化 ODF 可能在不真实的位置产生峰，这是毛书 §3-5 重点讨论的鬼峰。另一方面，有限阶展开一个很窄的峰会产生旁瓣和负振荡，即使测量没有同类零空间也可能发生。
噪声、背景误差、峰重叠、未测角域又会造成第三类伪结构。判断时分别检查：改变展开阶数后旁瓣怎样变；增加独立极图或完整取向数据后是否仍存在；更换可接受的先验模型后组分是否稳定。不能把所有小峰都叫鬼峰，也不能把某个孪晶位置的小峰直接当作发生了孪生。
正值和 Gaussian 模型提供了什么信息
如果限定 f≥0，可排除大量负密度解；若再限定由若干 Gaussian 组分和随机底组成，未知数减少为组分中心、宽度和份额。毛书 §3-8 用这种模型匹配可观测偶阶系数，再计算模型对应的奇阶部分以补全 ODF。
这一步得到的是“符合所选组分模型并与数据相容”的完整分布。模型之外的非负解未必被排除，拟合也可能存在多组中心和宽度。应把新增奇阶信息标记为先验支持的推断，不能说是衍射独立测出的。现代反演同样需要区分数据残差小和解唯一这两个判断。
真实应用：为什么铝板组分不能只靠峰高
毛书 pp.96–102 的纯铝冷轧案例把 ODF、纤维密度、组分体积和散布宽度并列分析；其中 S 组分的密度峰并不必然最高，却可拥有较大体积分量。这个例子提醒我们，极图反演的目标是得到可解释的取向分布及其不确定性，而不是找最高红点。
Hielscher 与 Schaeben 的 MTEX 反演论文提供了另一种实现路线。使用任何程序时都应留存测量掩膜、反射列表、权重与回算 PF。一个未参与拟合的可靠反射，或独立 EBSD 取向分布，可以提供额外约束；但比较时仍需保证相、取样位置和统计权重一致。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理，不复制受版权保护的标准正文。