Hermite矩阵的特征值估计——courant-fischer定理

2023-12-01 05:15

本文主要是介绍Hermite矩阵的特征值估计——courant-fischer定理,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Hermite矩阵的特征值估计——courant-fischer定理

一、courant-fischer定理(min-max定理)

将hermite矩阵的特征值表示为一系列最优化问题的解。

  1. 一个函数 R ( x ) = x H A x x H x R(x)=\frac{x^HAx}{x^Hx} R(x)=xHxxHAx,称为Rayleigh商,A是hermite矩阵
  2. λ m i n = min ⁡ x ≠ 0 R ( x ) \lambda_{min} = \min_{x\ne0} R(x) λmin=minx=0R(x)
  3. λ m a x = max ⁡ x ≠ 0 R ( x ) \lambda_{max} = \max_{x\ne0} R(x) λmax=maxx=0R(x)
  4. λ k = m i n d i m ( U ) = k max ⁡ x ∈ U , x ≠ 0 R ( x ) \lambda_{k} = min_{dim(U)=k} \max_{x\in U,x\ne0} R(x) λk=mindim(U)=kmaxxU,x=0R(x)
  5. λ k = m a x d i m ( U ) = n − k + 1 min ⁡ x ∈ U , x ≠ 0 R ( x ) \lambda_{k} = max_{dim(U)=n-k+1} \min_{x\in U,x\ne0} R(x) λk=maxdim(U)=nk+1minxU,x=0R(x)

该定理的2、3是容易证明和理解的。如果不能深入理解子空间的含义,将难以理解4、5。

首先给出4、5的另外一种表述:
λ 1 ≤ λ 2 ≤ ⋯ ≤ λ n , k = n , n − 1 , ⋯ , 1 λ k = min ⁡ w 1 , w 2 , ⋯ , w n − k ∈ C n max ⁡ x ≠ 0 , x ∈ C n , x ⊥ w 1 , w 2 , ⋯ , w n − k R ( x ) λ k = max ⁡ w 1 , w 2 , ⋯ , w k − 1 ∈ C n min ⁡ x ≠ 0 , x ∈ C n , x ⊥ w 1 , w 2 , ⋯ , w k − 1 R ( x ) \lambda_1 \le \lambda_2 \le \cdots \le \lambda_n,k=n,n-1,\cdots,1\\ \lambda_k = \min_{w_1,w_2,\cdots,w_{n-k}\in C^n} \max_{x\ne0,x\in C^n,\atop x\perp w_1,w_2,\cdots,w_{n-k}} R(x)\\ \lambda_k = \max_{w_1,w_2,\cdots,w_{k-1}\in C^n} \min_{x\ne0,x\in C^n,\atop x\perp w_1,w_2,\cdots,w_{k-1}} R(x) λ1λ2λn,k=n,n1,,1λk=w1,w2,,wnkCnminxw1,w2,,wnkx=0,xCn,maxR(x)λk=w1,w2,,wk1Cnmaxxw1,w2,,wk1x=0,xCn,minR(x)
关注4。

1️⃣当k=n时,U是n维复空间 V ( C n ) V(C^n) V(Cn),当x旋转到和 λ n \lambda_n λn的特征向量共线的时候, λ n = λ m a x \lambda_n=\lambda_{max} λn=λmax

2️⃣当k=n-1时,U是n-1维复空间, U ⊂ V ( C n ) U \subset V(C^n) UV(Cn),这样的子空间有无穷多个(可以想象,三维空间中有无穷多个二维平面)。当这个空间包含了 λ n \lambda_n λn的特征子空间的时候, max ⁡ x ∈ U , x ≠ 0 R ( x ) = λ m a x = λ n \max_{x\in U,x\ne0} R(x)=\lambda_{max}=\lambda_n maxxU,x=0R(x)=λmax=λn。注意,本定理名叫min-max定理,此时还仅仅只考虑极大那部分。

对于 m i n d i m ( U ) = k min_{dim(U)=k} mindim(U)=k,该min符号的含义是,在无穷个n-1维复空间中,要找到一个U,使得以U为可行域的函数 max ⁡ x ∈ U , x ≠ 0 R ( x ) \max_{x\in U,x\ne0} R(x) maxxU,x=0R(x)取得最小值。U的取法就是排除 λ n \lambda_n λn的特征子空间,也就是说 λ n \lambda_n λn的特征向量不属于U,或者说和U的基向量正交。

此时的max部分,由于U不包含 λ n \lambda_n λn的特征子空间,那么 max ⁡ x ∈ U , x ≠ 0 R ( x ) = m a x ( λ ( A ) − { λ n } ) = λ n − 1 \max_{x\in U,x\ne0} R(x)=max(\lambda(A)-\{\lambda_n\})=\lambda_{n-1} maxxU,x=0R(x)=max(λ(A){λn})=λn1

3️⃣对于k=n-2,是同样的理解方式,只不过取min的时候,要排除掉 λ n 和 λ n − 1 \lambda_n和\lambda_{n-1} λnλn1

关于该定理数学上的证明,可参考特征值的重要定理:Courant-Fischer min-max theorem 极大极小定理

二、韦尔定理(Weyl定理)

对于一个n阶Hermite矩阵A,受到一个n阶Hermite矩阵B的扰动,那么A+B的第k个特征值满足:
λ k ( A ) + λ m i n ( B ) ≤ λ k ( A + B ) ≤ λ k ( A ) + λ m a x ( B ) \lambda_k(A)+\lambda_{min}(B)\le\lambda_k(A+B)\le\lambda_k(A)+\lambda_{max}(B) λk(A)+λmin(B)λk(A+B)λk(A)+λmax(B)
证明,由min-max定理容易证得。

这篇关于Hermite矩阵的特征值估计——courant-fischer定理的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/439936

相关文章

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

hdu 4565 推倒公式+矩阵快速幂

题意 求下式的值: Sn=⌈ (a+b√)n⌉%m S_n = \lceil\ (a + \sqrt{b}) ^ n \rceil\% m 其中: 0<a,m<215 0< a, m < 2^{15} 0<b,n<231 0 < b, n < 2^{31} (a−1)2<b<a2 (a-1)^2< b < a^2 解析 令: An=(a+b√)n A_n = (a +

hdu 6198 dfs枚举找规律+矩阵乘法

number number number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description We define a sequence  F : ⋅   F0=0,F1=1 ; ⋅   Fn=Fn

线性代数|机器学习-P35距离矩阵和普鲁克问题

文章目录 1. 距离矩阵2. 正交普鲁克问题3. 实例说明 1. 距离矩阵 假设有三个点 x 1 , x 2 , x 3 x_1,x_2,x_3 x1​,x2​,x3​,三个点距离如下: ∣ ∣ x 1 − x 2 ∣ ∣ 2 = 1 , ∣ ∣ x 2 − x 3 ∣ ∣ 2 = 1 , ∣ ∣ x 1 − x 3 ∣ ∣ 2 = 6 \begin{equation} ||x

数据集 3DPW-开源户外三维人体建模-姿态估计-人体关键点-人体mesh建模 >> DataBall

3DPW 3DPW-开源户外三维人体建模数据集-姿态估计-人体关键点-人体mesh建模 开源户外三维人体数据集 @inproceedings{vonMarcard2018, title = {Recovering Accurate 3D Human Pose in The Wild Using IMUs and a Moving Camera}, author = {von Marc

【线性代数】正定矩阵,二次型函数

本文主要介绍正定矩阵,二次型函数,及其相关的解析证明过程和各个过程的可视化几何解释(深蓝色字体)。 非常喜欢清华大学张颢老师说过的一段话:如果你不能用可视化的方式看到事情的结果,那么你就很难对这个事情有认知,认知就是直觉,解析的东西可以让你理解,但未必能让你形成直觉,因为他太反直觉了。 正定矩阵 定义 给定一个大小为 n×n 的实对称矩阵 A ,若对于任意长度为 n 的非零向量 ,有 恒成

Java验证辛钦大数定理

本实验通过程序模拟采集大量的样本数据来验证辛钦大数定理。   实验环境: 本实验采用Java语言编程,开发环境为Eclipse,图像生成使用JFreeChart类。   一,验证辛钦大数定理 由辛钦大数定理描述为: 辛钦大数定理(弱大数定理)  设随机变量序列 X1, X2, … 相互独立,服从同一分布,具有数学期望E(Xi) = μ, i = 1, 2, …, 则对于任意正数ε ,

数据集 Ubody人体smplx三维建模mesh-姿态估计 >> DataBall

Ubody开源人体三维源数据集-smplx-三维建模-姿态估计 UBody:一个连接全身网格恢复和真实生活场景的上半身数据集,旨在拟合全身网格恢复任务与现实场景之间的差距。 UBody包含来自多人的现实场景的1051k张高质量图像,这些图像拥有2D全身关键点、3D SMPLX模型。 UBody由国际数字经济学院(IDEA)提供。 (UBody was used for mesh r

python科学计算:NumPy 线性代数与矩阵操作

1 NumPy 中的矩阵与数组 在 NumPy 中,矩阵实际上是一种特殊的二维数组,因此几乎所有数组的操作都可以应用到矩阵上。不过,矩阵运算与一般的数组运算存在一定的区别,尤其是在点积、乘法等操作中。 1.1 创建矩阵 矩阵可以通过 NumPy 的 array() 函数创建。矩阵的形状可以通过 shape 属性来访问。 import numpy as np# 创建一个 2x3 矩阵mat

特征值分解(EVD)和奇异值分解(SVD)—应用于图片压缩

特征值分解(EVD)和奇异值分解(SVD)—应用于图片压缩 目录 前言 一、特征值分解 二、应用特征值分解对图片进行压缩 三、矩阵的奇异值分解 四、应用奇异值分解对图片进行压缩 五、MATLAB仿真代码 前言         学习了特征值分解和奇异值分解相关知识,发现其可以用于图片压缩,但网上没有找到相应代码,本文在学习了之后编写出了图片压缩的代码,发现奇异值分