形式陈述
一组方向怎样作为整体被认证
设 ,目标不变子空间理路不变子空间Invariant subspace在线性算子作用下保持封闭、可被单独限制研究的子空间。 的维数为 ,。取其正交规范基 ,再取正交补基 ,使 酉。近似空间由 给出,要求 。
令
是目标空间的正交投影理路正交投影Orthogonal projection把向量映到子空间上最近点并使误差与子空间正交的线性算子。。用二范数和 Frobenius 范数理路矩阵范数与诱导算子范数Matrix norm · Induced matrix norm · Operator norm of a matrix用诱导范数和常用可计算矩阵范数度量线性映射的放大能力,并区分算子范数、Frobenius 范数与谱半径。定义两个基不变误差
若
则
对二范数采用更明确的分隔条件:若 ,而 在开区间 外,,则
这里的谱都是按特征值理路特征值与特征向量Eigenvalue and eigenvector满足 Tv=λv 且 v 非零的标量 λ 与向量 v。组成的有限集合。第一条允许两组谱在实轴上交替,第二条要求整个目标区间和互补谱分开。不能只凭逐对最小距离,就把第一条中的所有 Frobenius 范数替成二范数。
子空间问题怎样落到 Sylvester 方程
Hermitian 性使 也不变,因此左乘 得
这是Sylvester 方程理路Sylvester 方程与谱分离Sylvester equation刻画 AX−XB=C 对所有右端唯一可解的充要条件,并用非正规矩阵说明特征值间距不能代替 sep 敏感性。。用谱定理理路有限维谱定理Finite-dimensional spectral theorem有限维复正规算子存在正交规范特征基;实数情形对应自伴算子。分别把两侧 Hermitian 系数酉对角化,每个变换后元素满足 。逐项平方求和,便得 ,证明 Frobenius 界。
二范数版本不能逐项平方,但分隔区间给另一条路。令 、。有 ,而 的最小奇异值至少为 。所以
第二行的第一项可由 验证。证明完成后,未知目标基通常只留在理论解释中;可实际使用较粗的 ,前提是互补谱的下界已被另行认证。
直觉
一簇相近的特征值可能让簇内坐标轴大幅旋转,但只要整簇远离其余谱,所张成的平面仍可以稳定。 测量 把近似平面推出自身的分量, 则测量把这部分误差解释为“错平面”时需要的放大。
把 换为 ,其中 酉, 变成 , 变成 ,各范数与误差不变。用SVD理路奇异值分解Singular value decomposition · SVD任意有限维线性映射都可在正交规范基下表示为非负对角伸缩。取 的奇异值为 ,则 给出 、。这解释了“sin Θ”所度量的对象。
例子与边界
重根目标平面的完整残差计算
取 ,目标 ,近似基
有
互补矩阵为 ,所以 。残差范数 ,得到 。实际 ,所以 。
保留投影残差还能取到等号:,除以 恰为 。主角分别为 与 ,第一方向已经完全正确,误差只来自第二方向向外倾斜。
残差为零也未必选对目标簇
对同一个 ,若取 ,则 ,但近似空间包含错误的高谱方向。此时 含 ,与目标空间的互补谱相交,,定理不适用。残差必须配合谱簇识别,不能凭 宣称已经找到指定簇。
另一方面,在正确平面内把 旋转 ,仍有 、,尽管每列与原先选中的轴都不相同。目标重根没有唯一的单向量对应,但整个平面有明确意义。
推论与应用
若 是 的某一不变空间正交基,那么关于 的压缩残差是 ,故 。结合分隔条件就得到 Davis–Kahan 型扰动控制;选择对应簇及验证其维数和间隔仍是必要输入。
有限精度下必须先检查 接近 ,再控制正交化和残差误差。若几列几乎相关,直接把未归一的 当作子空间证书可能任意缩小数值。Ritz 压缩理路Hermitian 压缩、Ritz 值与交错Cauchy interlacing theorem · Poincaré separation theorem · Hermitian Ritz values证明子空间压缩的有序谱交错,算出 Ritz 向量与原残差,并说明压缩特征值何时不是真特征值。说明 的来源,单向量证书理路Hermitian 特征对的残差与方向证书Hermitian eigenpair residual bound · Residual-gap eigenvector angle bound把可计算特征对残差分别转成谱距离、结构化后向误差与有谱隙条件的方向界,并核验三阶扰动。是 的对应机制。
参考资料