Skip to content

定理Theorem

谱簇子空间的残差与间隔界

Invariant subspace residual bound · Davis–Kahan sin theta theorem · Subspace angle residual certificate

把近似子空间的残差投到互补谱块,证明 Frobenius 与分隔区间下的二范数角度界,避免追踪重根内部任意基。

形式陈述 ​

一组方向怎样作为整体被认证 ​

设 A=A∗∈Cn×n,目标不变子空间 X 的维数为 m,1≤m<n。取其正交规范基 X,再取正交补基 X⊥,使 [X,X⊥] 酉。近似空间由 Q∈Cn×m 给出,要求 Q∗Q=I。

令

P=XX∗,H=Q∗AQ,R=AQ−QH,A⊥=X⊥∗AX⊥.

P 是目标空间的正交投影。用二范数和 Frobenius 范数定义两个基不变误差

d2=‖(I−P)Q‖2=‖X⊥∗Q‖2,dF=‖(I−P)Q‖F.

若

δ=minα∈σ(A⊥),β∈σ(H)|α−β|>0,

则

dF≤‖X⊥∗R‖Fδ≤‖R‖Fδ.

对二范数采用更明确的分隔条件:若 σ(H)⊂[a,b],而 σ(A⊥) 在开区间 (a−δ,b+δ) 外,δ>0,则

d2≤‖X⊥∗R‖2δ≤‖R‖2δ.

这里的谱都是按特征值组成的有限集合。第一条允许两组谱在实轴上交替,第二条要求整个目标区间和互补谱分开。不能只凭逐对最小距离,就把第一条中的所有 Frobenius 范数替成二范数。

子空间问题怎样落到 Sylvester 方程 ​

Hermitian 性使 X⊥ 也不变,因此左乘 X⊥∗ 得

A⊥Z−ZH=X⊥∗R,Z=X⊥∗Q.

这是Sylvester 方程。用谱定理分别把两侧 Hermitian 系数酉对角化,每个变换后元素满足 (αi−βj)zij=gij。逐项平方求和,便得 δ2‖Z‖F2≤‖G‖F2,证明 Frobenius 界。

二范数版本不能逐项平方,但分隔区间给另一条路。令 c=(a+b)/2、h=(b−a)/2。有 ‖H−cI‖2≤h,而 A⊥−cI 的最小奇异值至少为 h+δ。所以

‖A⊥Z−ZH‖2≥‖(A⊥−cI)Z‖2−‖Z(H−cI)‖2≥(h+δ)‖Z‖2−h‖Z‖2=δ‖Z‖2.

第二行的第一项可由 Z=(A⊥−cI)−1(A⊥−cI)Z 验证。证明完成后,未知目标基通常只留在理论解释中;可实际使用较粗的 ‖R‖/δ,前提是互补谱的下界已被另行认证。

直觉

一簇相近的特征值可能让簇内坐标轴大幅旋转,但只要整簇远离其余谱,所张成的平面仍可以稳定。R 测量 A 把近似平面推出自身的分量,δ 则测量把这部分误差解释为“错平面”时需要的放大。

把 Q 换为 QW,其中 W 酉,H 变成 W∗HW,R 变成 RW,各范数与误差不变。用SVD取 X∗Q 的奇异值为 cos⁡θj,则 Q∗(I−P)Q=I−(X∗Q)∗(X∗Q) 给出 d2=maxjsin⁡θj、dF2=∑jsin2⁡θj。这解释了“sin Θ”所度量的对象。

例子与边界

重根目标平面的完整残差计算 ​

取 A=diag(1,1,5),目标 X=[e1,e2],近似基

Q=[e1,(04/53/5)].

有

H=diag(1,61/25),R=[0,(0−144/125192/125)].

互补矩阵为 A⊥=(5),所以 δ=5−61/25=64/25。残差范数 ‖R‖2=‖R‖F=48/25,得到 d2≤3/4。实际 X⊥∗Q=(0,3/5),所以 d2=dF=3/5。

保留投影残差还能取到等号:‖X⊥∗R‖=192/125,除以 64/25 恰为 3/5。主角分别为 0 与 arcsin⁡(3/5),第一方向已经完全正确,误差只来自第二方向向外倾斜。

残差为零也未必选对目标簇 ​

对同一个 A,若取 Q=[e1,e3],则 R=0,但近似空间包含错误的高谱方向。此时 H 含 5,与目标空间的互补谱相交,δ=0,定理不适用。残差必须配合谱簇识别,不能凭 R=0 宣称已经找到指定簇。

另一方面,在正确平面内把 e1,e2 旋转 45∘,仍有 R=0、d2=0,尽管每列与原先选中的轴都不相同。目标重根没有唯一的单向量对应,但整个平面有明确意义。

推论与应用

若 Q 是 A+E 的某一不变空间正交基,那么关于 A 的压缩残差是 R=−(I−QQ∗)EQ,故 ‖R‖2≤‖E‖2。结合分隔条件就得到 Davis–Kahan 型扰动控制;选择对应簇及验证其维数和间隔仍是必要输入。

有限精度下必须先检查 Q∗Q 接近 I,再控制正交化和残差误差。若几列几乎相关,直接把未归一的 AQ−QH 当作子空间证书可能任意缩小数值。Ritz 压缩说明 H 的来源,单向量证书是 m=1 的对应机制。

参考资料
关系图谱17 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系