Laplace算子和Laplacian矩陣

1 Laplace算子的物理意義

Laplace算子的定義為梯度的散度。

在Cartesian坐標系下也可表示為：

或者，它是Hessian矩陣的跡：

以熱傳導方程為例，因為熱流與溫度的梯度成正比，那麼溫度的梯度的散度就是熱量的損失率。

由此可見，Laplace算子可用於表現由於物質分佈不均引起的物質輸送。

2 Laplace算子的數學意義

現在，在一維空間中簡單分析上面的式子：

也可以寫作：

把分子第一項和第二項分別按泰勒展開：

可以看出Laplace算子實際上是一個使函數取平均的算子。多維空間相似。

3 Laplace方程

若Laplace算子右邊為零，稱為Laplace方程。Laplace方程的解稱為調和函數。若右邊是一個函數，稱為泊松方程。

4 Laplace算子在圖像處理的運用

圖像處理是以像素作為基礎離散化，如下：

5 Laplacian 矩陣

是一種用于表示圖的矩陣。它的維度是 |V|-by-|V| ( |V| 是節點的數目 )。 James Demmel提供了一種由Incidence matrix轉化為Laplacian矩陣的方法。

In(G)是一個 |V|-by-|E| 矩陣( |E| 是邊的數目 ), 設邊e=(i,j)，這一列除了第i行（為+1）和第j行(為-1)外都為零。需要说明的是，根据这个定义，对于无向图 e=(i,j) 和 e=(j,i) 是等价的，看似会生成很多不同的In图（根据每条边不同的取向）。但是实际上可以证明，无论边的方向怎么取，由In图生成的L图都是唯一的。也就是说， e=(i,j) 和 e=(j,i) 怎么取是无关紧要的。如何使用In图生成L图：

可得知Laplacian矩陣的两个重要性质：一是为对称阵。二是存在一个为零的特征值（秩为|V|-1）。三是一個半正定矩陣。注意Laplace算子是负定的。

在求解含Laplacian矩陣的方程組時，常常要求為正定矩陣。觀察發現這是因為Laplacian矩陣每列相加等于零。這時只需要手動更改第一行和第一列（比如第一個元素設為1，其余設為零），破壞其結構，令秩等于|V|就可以了。

对于非正定矩阵，左乘个transpose of the matrix, 推导如下：

Ax – b = 0

最小化 ||Ax – b||^2，展開后對x求導數：

可轉化為正定方程組。

6 Laplace算子和Laplacian矩陣的關系

Laplace算子可以推廣到多維情況計算。Laplacian矩陣主要用于三維以下的圖形學計算，可以表現復雜的幾何結構。而Lapace方程使用了Laplace算子來表示Laplacian矩陣。

參考文獻：
丘成桐，幾何三十載，香港中文大學
James Demmel, related lectures, University of California, Berkeley

Laplace算子和Laplacian矩陣相关推荐

图像梯度——Sobel算子和Laplacian算子
一.Sobel算子 1.定义 Sobel算子是一种离散的微分算子,结合了高斯平滑和微分求导运算,利用局部拆分寻找边缘,计算所得的是一个梯度的近似值. Sobel算子=|左-右|/|下-上| Schar ...
12-图像梯度-Scharr算子和laplacian算子
Scharr算子 cv2.Scharr(img,cv2.CV_64F,1,0) 第一个参数:当前的图像对象名称第二个参数:当前图像的深度,通常情况下指定为-1,表示输出和输入的深度是一样的:cv2. ...
Opencv(C++)笔记--Sobel算子和laplacian算子
目录 1--Sobel算子 1-1--原理 1-2--Opencv API 1-3--代码实例 1-4--Scharr算子 1-5--扩展操作 2--laplacian算子 2-1--原理 2-2-- ...
python图像锐化增强边缘_[Python图像处理]十一.图像锐化与边缘检测之Roberts算子、Prewitt算子、Sobel算子和Laplacian算子,Schar算子...
Roberts算子 Roberts算子即为交叉微分算法,它是基于交叉差分的梯度算法,通过局部差分计算检测边缘线条.常用来处理具有陡峭的第噪声图像,当图像边缘接近于正45度或负45度时,该算法处理效果更 ...
matlab shortest函数,MATLAB函數graphallshortestpaths不返回對稱矩陣
我正在使用MATLAB函數graphallshortestpaths來計算無向網絡頂點之間的最短路徑.無向網絡作爲加權邊緣列表文件給出,您可以在其中找到here.MATLAB函數graphallsho ...
java兩個矩陣相乘_Day12 練習java-矩陣乘法
昨天已經挑戰完矩陣的加和減了,今天來挑戰稍難一點的矩陣乘法,再開始之前先來簡單介紹矩陣的乘法吧~~ 對於矩陣來說跟一般數字的四則運算有著很大的差異,就是加或減與乘跟除的不同,矩陣的乘除比較複雜,除法是 ...
[Python图像处理] 十八.图像锐化与边缘检测之Scharr算子、Canny算子和LOG算子
该系列文章是讲解Python OpenCV图像处理知识,前期主要讲解图像入门.OpenCV基础用法,中期讲解图像处理的各种算法,包括图像锐化算子.图像增强技术.图像分割等,后期结合深度学习研究图像识别 ...
数字图像处理---LOG算子和CANNY算子边缘提取(matlab)
LOG算子和CANNY算子边缘提取边缘的含义: 在数字图像中,边缘是指图像局部变化最显著的部分,边缘主要存在于目标与目标,目标与背景之间,是图像局部特性的不连续性,如灰度的突变.纹理结构的突变.颜色 ...
[R]如何利用rnorm生成亂數矩陣和時間序列亂數?
在數值計算軟體的應用上,生成各類符合要求的隨機亂數分布以進一步驗證模型,幾乎是基本需求.R語言中,可利用rnorm指令來生成符合某個平均值和某個標準差要求的常態亂數分布.另外,亦可善用rnorm來生成 ...

Laplace算子和Laplacian矩陣

Laplace算子和Laplacian矩陣相关推荐

最新文章

热门文章