拉普拉斯算子

拉普拉斯算子是n维欧几里德空间中的一个二阶微分算子,定义为梯度(▽f)的散度(▽·f)。因此如果f是二阶可微的实函数,则f的拉普拉斯算子定义为:

f的拉普拉斯算子也是笛卡儿坐标系xi中的所有非混合二阶偏导数:

作为一个二阶微分算子,拉普拉斯算子把C函数映射到C函数,对于k ≥ 2。表达式(1)(或(2))定义了一个算子Δ : C(R) → C(R),或更一般地,定义了一个算子Δ : C(Ω) → C(Ω),对于任何开集Ω。

函数的拉普拉斯算子也是该函数的黑塞矩阵的迹:

坐标表示式

二维空间

其中x与y代表 x-y 平面上的笛卡儿坐标:

另外极坐标的表示法为:

三维空间

笛卡儿坐标系下的表示法

圆柱坐标系下的表示法

球坐标系下的表示法

N 维空间

在参数方程为(其中以及)的N维球坐标系中,拉普拉斯算子为:

其中是N − 1维球面上的拉普拉斯-贝尔特拉米算子。

恒等式

如果f和g是两个函数,则它们的乘积的拉普拉斯算子为: f是径向函数f(r)且g是球谐函数Ylm(θ,φ),是一个特殊情况。这个情况在许多物理模型中有所出现。f(r)的梯度是一个径向向量,而角函数的梯度与径向向量相切,因此: 球谐函数还是球坐标系中的拉普拉斯算子的角部分的特征函数:

因此:

推广

拉普拉斯算子可以用一定的方法推广到非欧几里德空间,这时它就有可能是椭圆型算子,双曲型算子,或超双曲型算子。

在闵可夫斯基空间中,拉普拉斯算子变为达朗贝尔算子:

达朗贝尔算子通常用来表达克莱因-高登方程以及四维波动方程。第四个项前面的符号是负号,而在欧几里德空间中则是正号。因子c是需要的,这是因为时间和空间通常用不同的单位来衡量;如果x方向用寸来衡量,y方向用厘米来衡量,也需要一个类似的因子。

拉普拉斯-贝尔特拉米算子

主条目:拉普拉斯-贝尔特拉米算子

拉普拉斯算子也可以推广为定义在黎曼流形上的椭圆型算子,称为拉普拉斯-贝尔特拉米算子。达朗贝尔算子则推广为伪黎曼流形上的双曲型算子。拉普拉斯-贝尔特拉米算子还可以推广为运行于张量场上的算子(也称为拉普拉斯-贝尔特拉米算子)。

另外一种把拉普拉斯算子推广到伪黎曼流形的方法,是通过拉普拉斯-德拉姆算子,它运行于微分形式。这便可以通过Weitzenböck恒等式来与拉普拉斯-贝尔特拉米算子联系起来。

拉普拉斯算子是n维欧几里德空间中的一个二阶微分算子,定义为梯度(▽f)的散度(▽·f)。因此如果f是二阶可微的实函数,则f的拉普拉斯算子定义为:

f的拉普拉斯算子也是笛卡儿坐标系xi中的所有非混合二阶偏导数:

作为一个二阶微分算子,拉普拉斯算子把C函数映射到C函数,对于k ≥ 2。表达式(1)(或(2))定义了一个算子Δ : C(R) → C(R),或更一般地,定义了一个算子Δ : C(Ω) → C(Ω),对于任何开集Ω。

函数的拉普拉斯算子也是该函数的黑塞矩阵的迹:

坐标表示式

二维空间

其中x与y代表 x-y 平面上的笛卡儿坐标:

另外极坐标的表示法为:

三维空间

笛卡儿坐标系下的表示法

圆柱坐标系下的表示法

球坐标系下的表示法

N 维空间

在参数方程为(其中以及)的N维球坐标系中,拉普拉斯算子为:

其中是N − 1维球面上的拉普拉斯-贝尔特拉米算子。

恒等式

如果f和g是两个函数,则它们的乘积的拉普拉斯算子为: f是径向函数f(r)且g是球谐函数Ylm(θ,φ),是一个特殊情况。这个情况在许多物理模型中有所出现。f(r)的梯度是一个径向向量,而角函数的梯度与径向向量相切,因此: 球谐函数还是球坐标系中的拉普拉斯算子的角部分的特征函数:

因此:

推广

拉普拉斯算子可以用一定的方法推广到非欧几里德空间,这时它就有可能是椭圆型算子,双曲型算子,或超双曲型算子。

在闵可夫斯基空间中,拉普拉斯算子变为达朗贝尔算子:

达朗贝尔算子通常用来表达克莱因-高登方程以及四维波动方程。第四个项前面的符号是负号,而在欧几里德空间中则是正号。因子c是需要的,这是因为时间和空间通常用不同的单位来衡量;如果x方向用寸来衡量,y方向用厘米来衡量,也需要一个类似的因子。

拉普拉斯-贝尔特拉米算子

主条目:拉普拉斯-贝尔特拉米算子

拉普拉斯算子也可以推广为定义在黎曼流形上的椭圆型算子,称为拉普拉斯-贝尔特拉米算子。达朗贝尔算子则推广为伪黎曼流形上的双曲型算子。拉普拉斯-贝尔特拉米算子还可以推广为运行于张量场上的算子(也称为拉普拉斯-贝尔特拉米算子)。

另外一种把拉普拉斯算子推广到伪黎曼流形的方法,是通过拉普拉斯-德拉姆算子,它运行于微分形式。这便可以通过Weitzenböck恒等式来与拉普拉斯-贝尔特拉米算子联系起来。


相关内容

  • log算子边缘检测6
  • 基于双边滤波的LOG边缘检测算法 摘要:传统LOG边缘检测算法采用高斯函数滤波时,尽管抑制了噪声,但同时也损坏了部分低强度边缘.针对这一问题,本文提出了结合双边滤波的LOG边缘检测算法,首先采用双边滤波来替代传统LOG算子中的高斯滤波,然后计算平滑后图像的拉普拉斯,最后提取运算后的零交叉点作为图像的 ...

  • 数字图像处理大作业
  • 数字图像处理实验报告 学 院: 信息学院 专 业: 电科1004班 姓 名: 学 号: 辅导老师: 完成日期: 2013年6月29日 空域图像增强 实验要求: (1) 选择若干图像(两幅以上),完成直方图均衡化. (2) 选择若干图像(两幅以上),对图像文件分别进行均值滤波.中值滤波和拉 普拉斯锐化 ...

  • 空间域滤波器(实验报告)
  • 数字图像处理作业 --空间域滤波器 摘要 在图像处理的过程中,消除图像的噪声干扰是一个非常重要的问题.本文利用matlab软件,采用空域滤波的方式,对图像进行平滑和锐化处理.平滑空间滤波器用于模糊处理和减小噪声,经常在图像的预处理中使用:锐化空间滤波器主要用于突出图像中的细节或者增强被模糊了的细节. ...

  • 几种典型的图像边缘检测算法的分析比较
  • ISSN1009-3044Computer与技术ComputerKnowledgeKnowledgeandandTechnologyTechnology电脑知识电脑知识与技术Vol.8,No.7,March2012. Tel:+86-551-[1**********]964 E-mail:eduf@ ...

  • -图像边缘提取方法研究
  • 图像边缘提取方法研究 摘 要 图像边缘检测一直以来都是图像处理与分析领域的研究热点.边缘提取是图像处理的基础工作,如何精确.有效地提取边缘是图像处理领域相关学者讨论的热点问题,由此产生的各种边缘检测算法层出不穷并且得到了广泛的应用.该文对传统的具有代表性的各种图像边缘提取方法进行了阐述.对比和分析了 ...

  • 边缘检测算法
  • 图像处理中各种边缘检测的微分算子简单比较 (Sobel,Robert, 不同图像灰度不同,边界处一般会有明显的边缘,利用此特征可以分割图像.需要说明的是:边缘和物体间的边界并不等同,边缘指的是图像中像素的值有突变的地方,而物体间的边界指的是现实场景中的存在于物体之间的边界.有可能有边缘的地方并非边界 ...

  • 图像分割方法
  • MATLAB的图像分割方法 蒋翠翠 湖北 李明 430079) 武汉 并借助MATLAB 对几种边缘检测算子进行了详细的分析比较,指出MATLAB 后再用新的区域重复这一过程,直到没有可接受的邻近点时生长过程终止:对于有多个目标区域的图像,可在之前已划分的区域外再选取一种子点,按同样过程生成新的区域 ...

  • Halcon学习之边缘检测函数
  • sobel_amp (Image : EdgeAmplitude : FilterType, Size : )根据图像的一次导数计算图像的边缘 close_edges (Edges, EdgeImage : RegionResult : MinAmplitude : ) close_edges_le ...

  • (二阶)基于新型二阶复合导数的边缘检测算法
  • 基于新型二阶复合导数的边缘检测算法:图像边缘是图像的重要特征,边缘检测是图像处理的重要步骤.本文将由分数阶次微摘要摘要: 分和分数阶次积分组合成的新型二阶复合导数应用于图像处理.对于二维图像的边缘检测,本文推导出一种全新的复数模板用于实现分数阶次微分.响应一起实现复合导数,在此基础上,形成了一种新的 ...