织构分析：测量与反演
PF 怎样反演 ODF：积分、级数、鬼峰与正则化
把极图反演拆成可跟随的数学步骤，理解 Bunge 级数与直接法分别求什么，以及为何“重算极图吻合”还不够。
先记住这三点
- PF 是 ODF 沿一族取向的积分，反演必然丢失部分信息。
- 级数截断误差和由不可见分量造成的鬼峰不是同一件事。
- 非负性和光滑性帮助选解，但不能把未测信息变成实测事实。
先复习前向问题
取晶体单位法向 h 与试样方向 y。所有满足 Rh=y 的取向并不唯一：满足条件后，晶体仍可绕 y 旋转任意角 γ。因此 PF 在 y 处收集的是这一整条取向轨道上的 ODF。用归一化约定可写 Pₕ(y)=(1/2π)∫₀²π f(R₀ Rot(h,γ))dγ，再按晶体等价极、反极等价进行相应平均。
R₀ 是任意一个把 h 送到 y 的旋转。它的选择只改变积分起点，不改变完整一圈的结果。三维函数经积分成为二维函数；多测几族极图能增加约束，但每一张并不是 ODF 的普通二维切片。
从积分变成线性方程
将取向空间划成小单元，令 xⱼ 为第 j 单元的密度或质量，二者不可混用。每个极图测点 i 对哪些取向单元有贡献，由晶体几何决定。把这个贡献记成 Aᵢⱼ，便得到 pᵢ≈ΣⱼAᵢⱼxⱼ，即 p≈Ax。若 x 是密度，A 必须包含相应单元体积；若 x 是质量，体积已在 x 内。
例如只有两条独立观测，未知却有三个单元，即使没有噪声，也可能存在多组非负 x 给出相同 p。增加测量方向与反射，是增加独立约束；把网格画得更密，只是增加未知数，并不凭空增加信息。
Bunge 级数为什么有用
如同周期函数可展开为正弦余弦，SO(3) 上的取向函数可展开为广义球谐函数：f(R)=ΣₗΣₘΣₙ cˡₘₙ Tˡₘₙ(R)。l 控制角度变化的细致程度，m、n 区分同阶基函数。晶体和试样对称性会消去不允许的系数，或把基函数组合成对称形式。
将展开式代入 PF 积分，利用积分的线性性，把求未知函数变为求有限组系数。保留 l≤L 就得到有限维模型。L 小时平滑而易漏掉尖锐峰，L 大时能表达尖峰但更容易放大噪声与产生振荡。这里的截断误差来自表示能力有限。
跟着推导 · 信息怎样丢失
测得的曲线相同，原分布可以不同
先用圆周上的简单模型理解“零空间”。仪器若只记录相反两个方向的平均，能否分辨隐藏的方向偏好？改变 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 与 −h 无法区分的条件。在标准球谐表述下，某些 ODF 分量对这样的极图不可见；典型讨论是可求的偶数阶部分和不可直接由极图确定的奇数阶部分。必须结合具体晶体对称性理解哪些阶真正允许，不能说任何晶体都“丢了一半有效信息”。
若仅凭可见部分重建，某些本不存在的峰可能出现，真实峰也可能被削弱。这称为鬼峰问题。继续增加同一种信息的级数阶数，不能补回本来不可见的分量。EBSD 单点取向测量直接提供完整旋转，不受相同的 PF 反演歧义约束，但仍受采样、角分辨率和索引错误影响。
非负条件怎样约束多解
ODF 是概率密度，必须 f≥0，且 ∫f dμ=1。若某个级数解在取向空间出现负密度，说明它不能直接解释为材料体积分布。传统方法利用零区、正定修正或分量模型来约束缺失部分；直接法通常把非负条件放进求解过程。
可用统一形式理解：在 x≥0、wᵀx=1 下，最小化 ||W(Ax−p)||²+λ||Lx||²。W 按观测可靠程度加权，L 惩罚过度起伏，λ 决定平滑倾向。此式是教学性的正则化模型，并不宣称每个软件内部都实现同一目标函数。熵方法则用概率分布的熵选择较少附加结构的解。
WIMV 与核函数方法各在做什么
WIMV 一类直接法先在取向格点上给出非负初值，再由实验极图与重算极图的比值，调整对这些测点有贡献的取向密度。重算强度过高时降低相关密度，过低时增加，并反复归一化。它把几何投影关系直接用于迭代，而不是先显式求一套球谐系数。
MTEX 的极图反演文献使用取向空间上的局域函数和优化，同时处理非负密度与归一化等问题。这里的“从极图估计 ODF”与“从 EBSD 已知取向做核密度估计”是两个问题：前者反演投影，后者平滑离散样本。菜单都写 ODF，并不代表输入信息相同。
怎样判断解是否可用
第一层检查重算 PF 与实测 PF 的残差图，而非只比较两个最大值。成片的有方向残差可能提示坐标、样品对称性或遗漏组分错误。第二层检查非负性、归一化、边界和缺测区域。第三层在合理范围内改变阶数、带宽或正则化，观察主要组分是否稳定。
若多种合理设置都给出相同的大组分，但很弱的尖峰随参数消失，那应把大组分作为稳健结论，把尖峰作为待证线索。若可能，用不同测量位置、未参与反演的反射或 EBSD 局部统计交叉核对；不要用参与拟合的同一组极图作为唯一“独立验证”。
把求解写成读者能复做的教程
准备经修正的极图及未测掩膜，说明晶体群、试样坐标和极族；选定反演方式及其参数；保存 ODF；计算回投影 PF；画残差；报告组分积分而不仅是峰高。若原始信息不支持唯一解，报告参数敏感性和可识别的部分。
毛卫民、张新明书中的 Bunge 展开与鬼峰章节适合追踪数学结构；Engler 等第 5 章可用于比较级数法、WIMV 与较新的核函数方法。两者结合阅读，能理解软件选择背后的数学代价。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理，不复制受版权保护的标准正文。