拉普拉斯矩陣定義爲:
L=D−A
其中,A表示鄰接矩陣,D表示度矩陣。
拉普拉斯矩陣的元素級別定義
Lij=⎩⎨⎧deg(vi)−10ifi=jifeij∈Eotherwise
拉普拉斯矩陣正則化形式(symmetric normalized
laplacian):Lsym=D−21LD−21
Lsym=[i,j]=⎩⎪⎨⎪⎧1deg(v1)deg(v2)−10ifi=jifeij∈Eotherwise
- 是對稱矩陣
- 有實特徵值
- 有實正交特徵矩陣(意味着VVT=E,即V−1=VT)
這意味着可以被正交對角化,即L=VΛV−1,又V是正交矩陣,則L=VΛVT。
拉普拉斯算子矩陣的定義來源於拉普拉斯算子,拉普拉斯算子是n維歐式空間中的一個二階微分算子:Δf=∑i=1n∂xi2∂2f。
把拉普拉斯算子用於圖像,就變成了邊緣檢測算子:⎣⎡˙0˙1˙01−41010˙⎦⎤
Δf===∂x2∂2f(x,y)+∂y2∂2f(x,y)[f(x+1,y)−f(x,y)−f(x−1,y)]+[(f(x,y+1)−f(x,y))−(f(x,y)−f(x,y−1)][f(x+1,y)+f(x−1,y)+f(x,y+1)+f(x,y−1)]−4f(x,y)
爲什麼拉普拉斯矩陣的定義爲:L=D−A?
思考,既然拉普拉斯算子是二階微分,那麼
我在b站看到這兩個問題的時候,對拉普拉斯矩陣的理解就突然有感覺了。
注:下面的解釋是從b站看的,我還沒看到其他的參考資料。所以這裏只做理解用。頻域、譜域以及拉普拉斯算子在圖上的應用 ,參考:https://b23.tv/av51204684/p7
第一種定義
- 圖上的函數是什麼?
對任一結點i,f(i)=結點的出度
- 圖上的函數梯度是什麼?
f(⋅)在結點i處沿着結點j方向的導數=f(i)−f(j)
f(⋅)在結點i處的梯度變化率=結點i出度方向上的導數和−結點i入度方向上的導數和
考慮下圖
可得
ff′f′變化率======[3010]T[e1e2e3e4]T//只考慮存在的邊⎣⎢⎢⎡f(1)−f(2)f(1)−f(3)f(1)−f(4)f(3)−f(4)⎦⎥⎥⎤//方向是邊的方向,完整的⎣⎢⎢⎡0−3−2−330102−10−13010⎦⎥⎥⎤[3231]T⎣⎢⎢⎡e1+e2+e3−e1e4−e2−e3−e4⎦⎥⎥⎤[8−3−1−4]T
第二種定義
- 圖上的函數是什麼?
對任一結點i,f(i)=結點的出度
- 圖上的函數梯度是什麼?
設K爲關聯矩陣
f(⋅)的梯度=KTf
梯度變化率=KKTf
fKKTfKKTf====[3010]T⎣⎢⎢⎡1−10010−10100−1001−1⎦⎥⎥⎤[3231]T[8−3−1−4]T
注意KKT=⎣⎢⎢⎡˙3˙−1˙−1˙−1−1100−102−1−10−12˙⎦⎥⎥⎤爲對應無向圖的拉普拉斯矩陣
驗證
L===D−Adiag([3,1,2,2])−⎣⎢⎢⎡0111100010011010⎦⎥⎥⎤⎣⎢⎢⎡3−1−1−1−1100−102−1−10−12⎦⎥⎥⎤