第 46B 章 矩阵均值¶
前置：算子单调函数 (Ch46A) · 正定矩阵 (Ch16) · 矩阵分析 (Ch14)
本章脉络：均值概念的算子化 \(\to\) Kubo-Ando 公理（单调性、连续性、变比例性） \(\to\) 核心映射：均值与算子单调函数的一一对应 \(\to\) 经典算子均值：算术平均 (\(A \nabla B\))、几何平均 (\(A \# B\))、调和平均 (\(A ! B\)) \(\to\) 算子均值不等式 \(\to\) 矩阵几何平均值的显式公式 \(\to\) 应用：量子度量几何（Fisher 信息度量）、电气工程中的并联阻抗计算、信号处理中的协方差平均
延伸：矩阵均值理论是解决算子非交换性的“和事利器”；它通过算子单调函数构建了一套公理化的系统，证明了即便在非交换环境下，我们依然可以定义具有优良代数性质的“中间态”算子，是现代量子计量学的数学核心
在标量代数中， \((a+b)/2\) 和 \(\sqrt{ab}\) 是我们习以为常的概念。然而，当 \(A\) 和 \(B\) 是不交换的矩阵时，直接定义 \(\sqrt{AB}\) 会导致结果不再是正定的。矩阵均值（Matrix Means）理论通过 Kubo-Ando 公理，为正定矩阵提供了一套严谨的“平均”方案。本章将介绍如何通过算子单调函数生成各种具有物理意义的矩阵均值。
46B.1 Kubo-Ando 公理化定义¶
定义 46B.1 (算子均值 \(\sigma\))
映射 \((A, B) \mapsto A \sigma B\) 称为算子均值，如果满足： 1. 单调性：若 \(A \preceq A'\) 且 \(B \preceq B'\)，则 \(A \sigma B \preceq A' \sigma B'\)。 2. 变比例性：\(C(A \sigma B)C^* = (CAC^*) \sigma (CBC^*)\)。 3. 连续性：在单调序列下是连续的。 4. 单位性：\(I \sigma I = I\)。
46B.2 核心对应关系¶
定理 46B.1 (均值与函数的对应)
每一个算子均值 \(\sigma\) 都唯一对应一个正半轴上的正算子单调函数 \(f\)： $\(A \sigma B = A^{1/2} f(A^{-1/2} B A^{-1/2}) A^{1/2}\)$ - 算术平均：\(f(t) = (1+t)/2 \Rightarrow A \nabla B = (A+B)/2\)。 - 几何平均：\(f(t) = \sqrt{t} \Rightarrow A \# B = A^{1/2} (A^{-1/2} B A^{-1/2})^{1/2} A^{1/2}\)。 - 调和平均：\(f(t) = 2t/(1+t) \Rightarrow A ! B = 2(A^{-1} + B^{-1})^{-1}\)。
46B.3 算子均值不等式¶
定理 46B.2 (算子 A-G-H 不等式)
对于任何正定矩阵 \(A, B\)： $\(A ! B \preceq A \# B \preceq A \nabla B\)$ 意义：这一不等式链在算子层面保持了标量均值的顺序，确立了矩阵“中间值”的层级结构。
练习题¶
1. [基础] 计算 \(A = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}\) 与 \(B = \begin{pmatrix} 4 & 0 \\ 0 & 9 \end{pmatrix}\) 的算子几何平均值。
参考答案
由于 \(A, B\) 是对角阵（交换）： 1. 几何平均值简化为对角元的标量几何平均。 2. 分量 1：\(\sqrt{1 \cdot 4} = 2\)。 3. 分量 2：\(\sqrt{1 \cdot 9} = 3\)。 结论：\(A \# B = \begin{pmatrix} 2 & 0 \\ 0 & 3 \end{pmatrix}\)。
2. [公式验证] 计算 \(\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix} \# \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix}\)。
参考答案
计算步骤： 1. 这两个矩阵的乘积为零矩阵。 2. 虽然公式涉及 \(A^{-1/2}\)，但对于半正定情形，均值定义为极限。 3. 由于两个算子张成的子空间正交，它们的“重叠部分”为零。 结论：\(A \# B = O\)。
3. [性质] 证明：若 \(A \sigma B = A \# B\)，则 \(B \sigma A = A \# B\)。
参考答案
证明： 1. 算子几何平均具有对称性：\(A \# B = B \# A\)。 2. 验证：\(B \# A = B^{1/2}(B^{-1/2} A B^{-1/2})^{1/2} B^{1/2}\)。 3. 利用恒等式 \(X(X^{-1}Y)^{1/2} = (YX^{-1})^{1/2}X\)，可以证明这与 \(A \# B\) 的表达式等价。 结论：几何平均值是关于变量对称的。
4. [调和平均] 证明：若 \(A, B\) 为电阻矩阵，\(A ! B\) 对应于什么物理配置？
参考答案
物理逻辑： 1. 两个电阻 \(R_1, R_2\) 并联的总电阻为 \(R_{total} = (R_1^{-1} + R_2^{-1})^{-1}\)。 2. 算子调和平均 \(A ! B = 2(A^{-1} + B^{-1})^{-1}\)。 结论：除了系数 2，这正是描述多端口网络并联阻抗的数学模型。
5. [单调性] 若 \(A \preceq A'\)，证明 \(A \nabla B \preceq A' \nabla B\)。
参考答案
证明： 1. \(A \nabla B = (A+B)/2\)。 2. \((A'+B)/2 - (A+B)/2 = (A'-A)/2\)。 3. 由于 \(A \preceq A'\)，有 \(A'-A \succeq 0\)。 4. 故差值是半正定的。结论成立。
6. [计算] 计算 \(A = \begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}\) 与 \(B = \begin{pmatrix} 8 & 0 \\ 0 & 8 \end{pmatrix}\) 的算术-几何平均值之差。
参考答案
计算： 1. \(A \nabla B = (2+8)/2 \cdot I = 5I\)。 2. \(A \# B = \sqrt{2 \cdot 8} \cdot I = 4I\)。 3. 差值 \(5I - 4I = I \succeq 0\)。符合 A-G 不等式。
7. [唯一性] 为什么 \(A \# B\) 的定义中不能简单写成 \(\sqrt{AB}\)？
参考答案
理由： 1. 若 \(A, B\) 不交换，\(AB\) 的特征值虽然是正的，但 \(AB\) 本身通常不是对称矩阵（Hermite）。 2. 算子均值必须定义在 Hermite 空间内以保持偏序关系的物理意义。 3. \(A^{1/2}(A^{-1/2}BA^{-1/2})^{1/2}A^{1/2}\) 这种夹心结构通过合同变换确保了结果的对称正定性。
8. [应用] 简述矩阵几何平均值在图像处理中的作用。
参考答案
在处理扩散张量图像 (DTI) 时，每个像素是一个正定矩阵（张量）。直接用算术平均会导致“膨胀”效应。利用矩阵几何平均值可以在保持张量行列式（体积）特性的同时实现更自然的平滑插值。
9. [极限] 证明：若 \(A, B\) 交换，则 \(A \# B = \sqrt{AB}\)。
参考答案
证明： 1. 若交换，则 \(A, B\) 有共同特征基。 2. \(A^{1/2} (A^{-1/2} B A^{-1/2})^{1/2} A^{1/2} = A^{1/2} (A^{-1} B)^{1/2} A^{1/2}\)。 3. \(= A^{1/2} A^{-1/2} B^{1/2} A^{1/2} = B^{1/2} A^{1/2}\)。 4. 由于交换，\(B^{1/2} A^{1/2} = (BA)^{1/2} = (AB)^{1/2}\)。
10. [最大均值] 在所有算子均值中，哪一个是最大的？
参考答案
结论：算术平均值 \(A \nabla B\)。 理由：任何算子单调函数 \(f\) 若满足 \(f(1)=1\)，在 \((0, \infty)\) 上都满足 \(f(t) \le (1+t)/2\)。这保证了算术平均是所有符合公理的均值体系的上限。
本章小结¶
矩阵均值理论是处理算子相互作用的“中道代数”：
- 公理的平衡：Kubo-Ando 公理确立了均值的合法性边界，将直观的“平均”感受转化为单调性与变比例性等严谨的代数特征。
- 非交换的桥梁：通过精心设计的“夹心结构”，几何平均等概念成功跨越了非交换性的障碍，为正定算子空间引入了自然的测地线路径。
- 不等式的延伸：算子 A-G-H 不等式不仅是标量结论的推广，更是描述量子信息流转与物理阻抗叠加的深层能量规律。