当前位置:   article > 正文

瑞利商(Rayleigh quotient)与广义瑞利商(genralized Rayleigh quotient)_广义瑞利商的最优解

广义瑞利商的最优解

我们首先来看看瑞利商的定义。瑞利商是指这样的函数 R ( A , x ) R(A,x) R(A,x):
R ( A , x ) = x H A x x H x R(A,x )=\frac{x^HAx}{x^Hx} R(A,x)=xHxxHAx

其中 x x x为非零向量,而 A A A为n×n的Hermitan矩阵。所谓的Hermitan矩阵就是满足共轭转置矩阵和自己相等的矩阵,即 A H = A A^H=A AH=A。如果我们的矩阵A是实矩阵,则满足 A T = A A^T=A AT=A的矩阵即为Hermitan矩阵。

瑞利商 R ( A , x ) R(A,x) R(A,x)有一个非常重要的性质,即它的最大值等于矩阵 A A A最大的特征值,而最小值等于矩阵 A A A的最小的特征值,也就是满足:
λ m i n ≤ x H A x x H x ≤ λ m a x λmin ≤\frac{x^HAx}{x^Hx} ≤λmax λminxHxxHAxλmax
具体的证明这里就不给出了。当向量 x x x是标准正交基时,即满足 x H x = 1 x^Hx=1 xHx=1时,瑞利商退化为: R ( A , x ) = x H A x R(A,x)=x^HAx R(A,x)=xHAx,这个形式在谱聚类和PCA中都有出现。

以上就是瑞利商的内容,现在我们再看看广义瑞利商。广义瑞利商是指这样的函数 R ( A , B , x ) R(A,B,x) R(A,B,x) :
R ( A , x ) = x H A x x H B x R(A,x)=\frac{x^HAx}{x^HBx} R(A,x)=xHBxxHAx
其中 x x x为非零向量,而 A , B A,B A,B为n×n的Hermitan矩阵。B为正定矩阵。它的最大值和最小值是什么呢?其实我们只要通过将其通过标准化就可以转化为瑞利商的格式。我们令 x = B − 1 2 x ′ x=B^{−\frac{1}{2}} x ′ x=B21x,则分母转化为:
x H B x = x ′ H ( B − 1 2 ) H B B − 1 2 x ′ = x ′ H B − 1 2 B B − 1 2 x ′ = x ′ H x ′ x^HBx=x ′^H (B^{−\frac{1}{2}})^H BB^{−\frac{1}{2}}x′ =x ′^HB^{−\frac{1}{2}} BB ^{−\frac{1}{2}} x ′ =x′^Hx′ xHBx=xH(B21)HBB21x=xHB21BB21x=xHx
而分子转化为:
x H A x = x ′ H B − 1 2 A B − 1 2 x ′ x^HAx=x′^HB^{−\frac{1}{2}}AB^{−\frac{1}{2}}x′ xHAx=xHB21AB21x

此时我们的 R ( A , B , x ) R(A,B,x) R(A,B,x)转化为 R ( A , B , x ′ ) R(A,B,x ′) R(A,B,x) :
R ( A , B , x ′ ) = x ′ H B − 1 2 A B − 1 2 x ′ x ′ H x ′ R(A,B,x ′ )=\frac{x′^HB^{−\frac{1}{2}} AB^{−\frac{1}{2}}x′}{x′^Hx′} R(A,B,x)=xHxxHB21AB21x

利用前面的瑞利商的性质,我们可以很快的知道, R ( A , B , x ′ ) R(A,B,x ′ ) R(A,B,x)的最大值为矩阵 B − 1 2 A B − 1 2 B^{−\frac{1}{2}} AB^{−\frac{1}{2}} B21AB21 的最大特征值,或者说矩阵 B − 1 A B^{−1}A B1A的最大特征值,而最小值为矩阵 B − 1 A B^{−1}A B1A的最小特征值。即对矩阵进行标准化。

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/你好赵伟/article/detail/83765
推荐阅读
相关标签
  

闽ICP备14008679号