织构分析：取向代数深入
从一条转轴到 Rodrigues 空间：旋转公式、奇点和基本区逐步推导
把取向矩阵背后的几何拆成平行与垂直分量，亲自推出轴角公式，再理解 Rodrigues 坐标为什么有直线纤维、非均匀体积和边界跳跃。
先记住这三点
- 坐标有奇点不意味着晶体发生突变，必须回到旋转本身判断。
- 转轴分量依赖所用参考系，错取向角相同不代表晶界相同。
- Rodrigues 空间的欧氏距离和体积不能直接冒充旋转距离和概率。
从已经熟悉的矩阵出发
仍规定 R 将晶体正交向量转到试样系，采用右手主动旋转。暂时把晶体对称性放在一边：一个旋转由单位轴 n 和角度 θ 描述，0≤θ≤π。n 的三个分量满足 n·n=1，故只有两个独立数，加上 θ 正好是三个自由度。
设任意向量 v，先沿转轴拆开：v∥=(n·v)n，v⊥=v−v∥。转动时平行部分不动，垂直部分在一个平面内转动。这幅几何图已经包含完整公式，下面只需要把平面三角函数写出来。
第一步：亲自构造 Rodrigues 旋转公式
在垂直平面内，v⊥ 和 n×v 是互相垂直且等长的两根基。右手旋转 θ 后，v⊥ 变成 v⊥cosθ+(n×v)sinθ。加回不动部分，Rv=v cosθ+(n×v)sinθ+n(n·v)(1−cosθ)。
定义叉乘矩阵 [n]×，使 [n]×v=n×v；其三行依次为 (0,−n₃,n₂)、(n₃,0,−n₁)、(−n₂,n₁,0)。于是 R=I cosθ+(1−cosθ)nnᵀ+sinθ[n]×。这称为 Rodrigues 旋转公式；不要把它与稍后由三个数构成的 Rodrigues 向量混为一谈。
第二步：从矩阵反求角度和转轴
取矩阵迹：trI=3，tr(nnᵀ)=n·n=1，反对称矩阵迹为零，因此 trR=1+2cosθ，θ=acos[(trR−1)/2]。再相减得 R−Rᵀ=2sinθ[n]×，所以 n=(R₃₂−R₂₃,R₁₃−R₃₁,R₂₁−R₁₂)/(2sinθ)。
这套符号对应本课的主动旋转 R。若采用书中的被动矩阵 g=Rᵀ，反对称项变号，轴也随之变号。程序中还应把 acos 输入因舍入造成的微小越界截到 [−1,1]；但矩阵本身严重不正交时，应先查坐标或拟合错误，不能靠截断掩盖。
零度与一百八十度怎样处理
θ=0 时 R=I，任意轴都给出同一个不转动操作，转轴根本没有唯一答案。θ=π 时 R=2nnᵀ−I，因此 (R+I)/2=nnᵀ；对角元只能给出 nᵢ²，非对角元给出 nᵢnⱼ，必须一起确定分量相对符号。
实作时可求 R 的特征值 1 对应的单位特征向量；接近 π 时也可从 R+I 中最长的一列归一化。近零角时轴对噪声极敏感，所以微小错取向角虽然可报告，其轴方向未必有可信精度。不能把所有零角点强行赋成 [111] 后统计轴偏好。
一组可以手算的完整例子
令 n=(1,1,1)/√3、θ=120°。cosθ=−1/2、sinθ=√3/2，代入得到 R 的三行为 (0,0,1)、(1,0,0)、(0,1,0)。所以 e₁→e₂、e₂→e₃、e₃→e₁。迹为零，反求 θ=acos(−1/2)=120°；反对称差除以 √3 得到原来的 n。
这个矩阵对带标签的三根轴是非平凡旋转，却是立方晶体的一种对称操作。因此两个取向相差这个矩阵，若已考虑立方晶体对称性，物理错取向可为零。这里第一次清楚看到：旋转矩阵的角度与对称约化后的晶体错取向是两个层次。
第三步：为什么 Rodrigues 向量使用半角正切
定义 r=n tan(θ/2)，其长度记录角度，方向记录转轴。利用 cosθ=(1−t²)/(1+t²)、sinθ=2t/(1+t²)，其中 t=|r|，代入前面的旋转式得到 R=[(1−r·r)I+2rrᵀ+2[r]×]/(1+r·r)。
这样三个实数即可直接生成旋转，不用先求三角函数。θ→0 时 r≈nθ/2；θ→π 时 |r|→∞。后一现象是这个坐标图无法在有限位置表示所有半周旋转，不是物理取向变成无限大。需要跨过这里时可改用单位四元数。
对称基本区怎样从旋转距离选出来
同一立方取向可写为 RC，C 遍历晶体的 24 个正旋转对称操作。对每个代表求相对参考取向的旋转角，选角度最小者；因 tan(θ/2) 单调，这等价于选最短 Rodrigues 向量。所有这样选出的点组成基本区，表面上可能有多个并列代表。
立方基本区沿坐标轴的截距是 tan22.5°=√2−1≈0.4142，因为绕该轴 45° 已走到两个相差 90° 的对称代表的中间。沿 [111] 的边界对应 60°，给 r=(1/3,1/3,1/3)。用对称操作直接计算基本区，比脱离所选点群死记一个多面体表可靠。
为什么画成直线，仍不能用尺子量错取向
绕固定轴、从单位取向出发的旋转有 r=t n，显然在过原点的直线上。一般固定轴旋转轨迹在这类投影中也可表示为直线段，但 t 与物理角 θ 的关系非线性。两点 r₁、r₂ 的距离 |r₂−r₁| 并不等于错取向角。
概率同样不能忽略坐标变换。在未作晶体对称约化的全 Rodrigues 空间，Haar 测度为 dμ=dr₁dr₂dr₃/[π²(1+|r|²)²]。它可从轴角测度 sin²(θ/2)dθ dΩ/(2π²) 及 t=tan(θ/2) 推得。等体积小立方格越远离原点，代表的随机取向概率越小，必须加权。
放回晶界图与取向梯度
若 R₁、R₂ 都映晶体到试样，Δₛ=R₂R₁ᵀ 是在试样系中将晶粒1旋到晶粒2的旋转；Δ𝚌=R₁ᵀR₂ 给出对应的晶体参考表达。它们由相似变换联系，角度相同，轴分量不同。先说明采用哪一个，才能判断轴是否接近平行表面或某条晶向。
Engler 等第2.7–2.9节正是沿此路径连接晶界错取向与计算表示。应用到真实 EBSD 数据时，基本区边界两侧的点可能在 Rodrigues 图上相隔很远，实际却很近。应先作对称最小距离和局部展开，再求梯度；直接差分三个 r 分量会在边界处制造假高梯度。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理，不复制受版权保护的标准正文。