谱图与代数图论读书笔记2
Jun 22, 2020 00:00 · 413 words · 2 minute read
谱图理论第二章用于描述矩阵作为自然的优化方案存在的,并给出min max theorem,并通过谱定理来证明Courant-Fischer Theorem理论,之后通过该理论证明谱定理。
第三章用于描述拉普拉斯矩阵的一些性质以及通过拉普拉斯矩阵的特征值对图进行绘制。
第二章 特征值以及优化
Hermitian matrix(埃尔米特矩阵)
Hermitian matrix(埃尔米特矩阵)也被称为自伴随矩阵,矩阵的第i行第j列的元素都与第j行第i列的元素的复共轭。
由于本书并不是讨论线性代数的书,对于图所表示的矩阵,不会出现虚数,因此,我们只探求埃尔米特矩阵的特殊形式实对称矩阵,其具有以下性质:
- 实对称矩阵的特征值都是实数,并且不同特征值对应的特征向量相互正交
- 如果其特征值都是整数,称这个矩阵是正定的,若为非负,则称为半正定。
- 对于一个实对称矩阵,An与逆矩阵AT都是实对称矩阵
谱定理(T 1.3.1)
如果矩阵M是一个n*n并且是一个实对称矩阵,那么存在n个特征值,使得其对应的特征向量相互正交。
Courant-Fischer Theorem(T 2.0.1)
这里引入一个瑞丽商(Rayleigh quotient)的概念,对于一个向量x,关于M的瑞丽商定义为:
xTMxxTx
假设Mϕ=μϕ,很显然,其中的μ即M的特征值,ϕ为特征向量,对于x=ϕ我们可以知道
ϕTMϕϕTϕ=ϕTμϕϕTϕ=μ
CF定理所阐明的是当x是使瑞丽商最大化的向量时,其正好是M最大特征值的特征向量。
对于一个实对称矩阵,且其特征值μ1≥μ2≥…≥μn,我们可以得到:
μk=maxS⊆Rn&dim(S)=kminx∈S&x≠0xTMxxTx
通过上面的定理我们可以知道,对于一个矩阵的最大特征值μn=maxxTMxxTx ,该定理在接下来的证明中具有相当重要的作用。
首先我们为了证明2.0.1引入一个引理(L 2.1.1):
对于一个实对称矩阵M,具有特征值μ1,..,μn,以及一系列相互正交的特征向量ϕ1,…,ϕn,将x 通过这一组基向量进行表示:
x=n∑i=1ciϕi
这样我们可以将x带入瑞丽商中:
xTMx=(∑iciϕi)TM(∑iciϕi) =(∑iciϕi)T(∑iciμiϕi) =∑i,jcicjμjϕTiϕj
由于ϕ相互正交,则上式可以化简为: xTMx=∑ic2iμi
回到证明2.0.1上来,前置条件如上,设S是该向量空间的一部分,也就是说,对于S⊆Rn并且dim(S)=k, 我们可以扩展X得到:
x=k∑i=1ciϕi
带入瑞丽商得到xTMx=∑ki=1μic2i∑ki=1c2i≥∑ki=1μkc2i∑ki=1c2i=μk
所以 minx∈SxTMxxTx≥μk,我们可以得到在一个k维空间中,μk与瑞丽商之间的关系,但是我们很好的思考到,在整个n维空间中的任意k维子空间的关系:
于是,我们设T为ϕk,…,ϕn这一系列的特征向量组成的空间,因此,其维度为n−k+1对于所有的k维子向量S,很明显,两者之间会有一个维度是重叠的,因此我们尝试比较这两者的关系:
minx∈SxTMxxTx≤minx∈S∩TxTMxxTx≤maxx∈TxTMxxTx
(上式通过类似于上面瑞丽商的展开就可以得到,至于后面这一个max基本同理。)
又对于x∈T上进行求证:
xTMxxTx≤μk
当S与T只在μk对应的特征向量上相交时,minx∈SxTMxxTx=maxx∈TxTMxxTx=μk 从而我们可以推出:
μk=maxS⊆Rn&dim(S)=kminx∈S&x≠0xTMxxTx
定理 T 2.2.1
对于一个实对称矩阵M,使得x作为一个非0的向量,并使得瑞丽商可以得到最大,即:
x=argmaxx∈RxTMxxTx
那么x满足Mx=μ1x,其中,μ1是M最大的特征值。相反,瑞丽商最小时对应的是最小的特征值。
谱定理相关证明
首先对2.2.1进行证明:
为了求得瑞丽商的极值,我们计算其梯度:
∇xTx=2x
∇xTMx=2MX
我们通过使∇f(x)=0得到:
(xTx)Mx=(xTMx)x
Mx=xTMxxTxx
当瑞丽商最大时,瑞丽商是M最大的特征值。
由此,我们可以扩展出定理2.2.2:
对于M的特征向量,我们可以得到:
当μ1的时候,ϕ1∈argmax||x||=1xTMx(这里是由于x的模为1,那么分母就自然而然的等于1了)
对于 2≤μi≤n ,ϕi∈argmax||x||=1&xTϕi=0,forj>ixTMx
对于2.2.2的证明,这里就不赘述了,基本上是通过构造附加的矩阵来计算并且推导关于特征值与瑞丽商的关系。
更近一步,我们还能得到对于x1,x2,…,xn且满足
x∈argmax||x||=1&xTxj=0,forj<ixTMx
则每一个x都是M的特征向量。
第三章 拉普拉斯矩阵以及图像的绘制
第二章通过对于矩阵的几点定理来推进对于矩阵的理解,第三章开始将图与矩阵相结合,通过矩阵来得到图相关的性质以及相关的应用。
拉普拉斯矩阵
第一章曾经提到过有关拉普拉斯的定义以及证明,在这里我们对于通过二次型进行重新梳理一下,
对于一个加权图G=(E,V,w),w:E→R+,定义拉式矩阵的二次型为:
xTLGx=∑(a,b)∈Ewa,b(x(a)−x(b))2
对于上述二次型的展开,我们在第一章已经证明过了,在这里就不要再来一遍了。
我们接下来引入δa表示元素在方向a上面的基本单位向量,对于两个点的一张图,x(a)−x(b)=(δT1−δT2)x
xTLGx=∑(a,b)∈Ewa,b((δT1−δT2)x)2 =∑(a,b)∈Ewa,bxT(δ1−δ2)(δT1−δT2)x
因此,LG=∑(a,b)∈Ewa,b(δa−δb)(δa−δb)T=∑(a,b)∈Ewa,bLGa,b
对于每一个向量x都有(LGx)(a)=d(a)x(a)−∑(a,b)∈Ewa,bx(b)
在sagt书上写的还可以继续往下引申的式子,我认为只有在wa,b=1时才成立:
(LGx)(a)=∑(a,b∈E)wa,b(x(a)−x(b))
T 3.1.1 拉普拉斯矩阵对于联通图的性质
对于一个图G,其拉普拉斯矩阵的特征值为0=λ1≤λ2≤…≤λn,如果第二个特征值不为0,那么图是联通的。
虽然很显而易见,这里还是做一些说明,对于一个非联通图,我们可以将其拉普拉斯矩阵分为两个部分:
L=[LG1 00 LG2]
所以L必定有两个相互正交的特征向量[01] 以及[10] ,这个1代表的是全1矩阵,于是L一定会有两个特征值都等于0显然上述的特点是有效的。
关于拉普拉斯矩阵的特征值对于绘图的几点猜测
拉普拉斯矩阵根据其二次型可以看到,如果将x看作是点,那么(x(a)−x(b))2可以看作是在平面上两点之间的距离,为了使距离最短,我们可以对其进行化简,总之最后通过拉普拉斯最小的且不等于0的特征向量所画出的图通常都是比较易得的,另外在图着色问题上,矩阵也是会取得不错的效果,详细的在第四章可以得到证明。
我有一个不成熟的猜测,对于一个二维的图,我们可以得到关于其的特征向量以及特征值,对于特征向量,一个全联通图会有n-1个特征向量并且相互正交,我们在第一章可以看到,将整个特征向量绘制会呈现出跳跃式的数据迁移,有没有可能是每一个特征值对应的特征向量都是该矩阵在某个基向量的投影,也就是说当我们画图的时候,我们可以用最小的非零特征值的向量与第二小的来进行绘制,其绘制出的图形是非常易看的。当然这只是一个猜测,并没有进行证明,嘿嘿嘿,别信这里的话哦!!!