采访中有人问我这个问题。
假设我们有一个形式为的相关矩阵
给定该相关矩阵,要求我查找gamma的值。
我以为我可以对特征值做些什么,因为它们都应该大于或等于0(矩阵应该是正半定数)-但我认为这种方法不会产生答案。我想念一个把戏。
您能否提供解决该问题的提示?
采访中有人问我这个问题。
假设我们有一个形式为的相关矩阵
给定该相关矩阵,要求我查找gamma的值。
我以为我可以对特征值做些什么,因为它们都应该大于或等于0(矩阵应该是正半定数)-但我认为这种方法不会产生答案。我想念一个把戏。
您能否提供解决该问题的提示?
Answers:
这是一个更简单(也许更直观)的解决方案:
将协方差视为抽象向量空间上的内积。然后,对于向量,,,相关矩阵中的条目为,其中尖括号表示和之间的夹角。v 1 v 2 v 3 ⟨ v我,v Ĵ ⟩ v我v Ĵ
不难想象被。因此,其余弦()的边界为。然后,基本三角给出。| ⟨ v 1,v 2 ⟩ ± ⟨ v 1,v 3 ⟩ | γ COS [ ⟨ v 1,v 2 ⟩ ± ⟨ v 1,v 3 ⟩ ] γ ∈ [ 0.6 × 0.8 - 0.6 × 0.8 ,0.6 × 0.8 + 0.6 ×
编辑:请注意,最后一行的确实是 - -由于,巧合出现了0.6和0.8的第二次出现。COS ⟨ v 1,v 2 ⟩ COS ⟨ v 1,v 3 ⟩ ∓ 罪⟨ v 1,v 3 ⟩ 罪⟨ v 1,v 2 ⟩ 0.6 2 + 0.8 2 = 1
这是我在最初对答案的评论中的意思,以及我认为@yangle可能在说什么(尽管我没有关注/检查它们的计算)。
“矩阵应为正半定数” 表示变量向量在欧几里得空间中是一堆。相关矩阵的情况比协方差矩阵更容易,因为三个向量的长度固定为1。想象3个单位向量XYZ并记住是角度的余弦。因此,和。的边界可能是什么?该关联可以采用Z围绕Y定义的任何值(保持角度):COS α = - [R X ý = 0.6 余弦β = - [R Ŷ Ž = 0.8 COS γ = - [R X ž ř Ŷ ž = 0.8
旋转时,两个位置都非常引人注目,即X极限,这两个都是Z落入平面XY时的位置。一个在X和Y之间,另一个在Y的另一侧。这些由蓝色和红色矢量表示。在这两个位置上,配置XYZ(相关矩阵)正好是奇异的。这些是Z可以达到Wrt X的最小和最大角度(因此相关)。
选取三角公式来计算平面上的角度之和或差,我们有:
作为界限。
这种几何视图只是以代数形式(次要等)表示的@rightskewed的另一种外观(在3D情况下为特定且简单的视图)。
与主要未成年人一起玩在3 x 3或4 x 4的问题上可能很好,但是在更高维度上没有足够的动力和数值稳定性。
对于像这样的单个“自由”参数问题,很容易看出,构成矩阵psd的所有值的集合将是单个间隔。因此,找到这样的最小值和最大值就足够了。这可以通过数值解决一对线性半定规划(SDP)问题来轻松实现:
例如,可以在MATLAB下使用YALMIP来公式化和数值解决这些问题。
快速,简便和可靠。
顺便说一句,如果聪明的裤子面试官问这个问题不知道,SemiDefinite编程可以很好地解决这个问题,它已经发展成熟并且具有复杂且易于使用的数值优化器来可靠地解决实际问题。告诉他/她这已经不是1870年了,现在是时候利用现代计算技术的发展了。
让我们考虑以下凸集
这是一个spectrahedron命名维elliptope。这是这个椭圆形的描绘
将此椭圆形与由和定义的平面相交,我们得到一个线段,其端点用黄色着色y = 0.8
椭圆形的边界是由
如果且,则上述三次方程式可归结为二次方程式
因此,椭圆形与两个平面的交点是由
每个正半定矩阵都是一个相关/协方差矩阵(反之亦然)。
要看到这一点,请从一个正半定矩阵并进行特征分解(由于是对称的,因此由频谱定理存在)其中是正交本征向量矩阵,是对角线在对角线上具有特征值的矩阵。然后,令,其中是对角矩阵,其对数值为平方根。甲甲= û d Ü Ť Ù d 乙= û d 1 / 2 ù Ť d 1 / 2
然后,采取与IID均值为零,方差为1级的条目,向量,并注意也具有零均值和协方差(和相关)矩阵。 B x A
现在,看到每个相关/协方差矩阵都是正半定的很简单:令为一个相关矩阵。然后,很容易看到,并且因此瑞利商对于任何非零都是非负的,因此是正半定的。- [R = - [R Ť 一个 Ť ř 一个 = ë [ (一个 Ť X )2 ] ≥ 0 一 ř
现在,注意到对称矩阵仅在特征值非负的情况下才是正半定的,我们看到您的原始方法是可行的:计算特征多项式,查看其根以查看它们是否为非负。请注意,使用西尔维斯特(Sylvester)的准则测试正定性是容易的(如另一个答案的评论中所述;当且仅当主要未成年人都具有正行列式时,矩阵才是正定);存在半定数的扩展(所有未成年人均具有非负行列式),但是在这种情况下,您必须检查未成年人,而对于正定数则仅需检查。 n