泊松方程解法

原文链接

泊松方程是数学中一个常见于静电学、机械工程和理论物理的偏微分方程。是从法国数学家、几何学家及物理学家泊松而得名的。 泊松方程为

在这里 Δ 代表的是拉普拉斯算子,而 f 和 φ 可以是在流形上的实数或复数值的方程。 当流形属于欧几里得空间,而拉普拉斯算子通常表示为 ,因此泊松方程通常写成

在三维直角坐标系,可以写成

如果没有f, 这个方程就会变成拉普拉斯方程

另外

在数学以及物理中, 拉普拉斯算子或是拉普拉斯算符(Laplace operator 或 Laplacian)是一个微分算子,通常写成 Δ 或 ;这是为了纪念皮埃尔-西蒙·拉普拉斯而命名的。

拉普拉斯算子有许多用途,此外也是椭圆型算子中的一个重要例子。

在物理中,常用于波方程的数学模型、热传导方程以及亥姆霍兹方程。

在静电学中,拉普拉斯方程和泊松方程的应用随处可见。在量子力学中,其代表薛定谔方程式中的动能项。

在数学中,经拉普拉斯算子运算为零的函数称为调和函数;拉普拉斯算子是霍奇理论的核心,并且是德拉姆上同调的结果。

[编辑] 定义

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

    (1) 

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

    (2) 

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

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

 

坐标表示式

二维空间

其中x与y代表 x-y 平面上的笛卡儿坐标
另外极坐标的表示法为:

三维空间

笛卡儿坐标系下的表示法

圆柱坐标系下的表示法

球坐标系下的表示法

时间: 2024-10-06 07:55:17

泊松方程解法的相关文章

从泊松方程的解法,聊到泊松图像融合

从泊松方程的解法,聊到泊松图像融合 成指导 ? 字节跳动 算法工程师 283 人赞同了该文章 2004 年 SIGGRAPH 上,Microsoft Research UK 有篇经典的图像融合文章<Poisson Image Editing>.先看看其惊人的融合结果(非论文配图,本人实验结果): 这篇文章的实现,无关目前算法领域大火的神经网络,而是基于泊松方程推导得出. 泊松方程是什么? 很多朋友比较熟悉概率论里面的泊松分布.泊松方程,也是同一个数学家泊松发明的.但却和泊松分布没有什么关系,是

一个不简洁的约瑟夫环解法

约瑟夫环类似模型:已知有n个人,每次间隔k个人剔除一个,求最后一个剩余的. 此解法为变种,k最初为k-2,之后每次都加1. 例:n=5,k=3.从1开始,第一次间隔k-2=1,将3剔除,第二次间隔k-1=2,将1剔除.依此类推,直至剩余最后一个元素. 核心思路:将原列表复制多份横向展开,每次根据间隔获取被剔除的元素,同时将此元素存入一个剔除列表中.若被剔除元素不存在于剔除列表,则将其加入,若已存在,则顺势后移至从未加入剔除列表的元素,并将其加入.如此重复n-1次.面试遇到的题,当时只写了思路,没

斐波那契数列的递归和非递归解法

//递归解法 function fib(n){ if(n < 1){ throw new Error('invalid arguments'); } if(n == 1 || n == 2){ return 1; } return fib(n - 1) + fib(n - 2); } //非递归解法 function fib(n){ if(n < 1){ throw new Error('invalid arguments'); } if(n == 1 || n == 2){ return 1

HDU1013 POJ1519 Digital Roots(解法三)

该问题的最佳解法是利用数论的9余数定理来计算数根.一个数的数根等于该数的9的余数,若余数为0则结果为9. 问题链接:HDU1013 POJ1519 Digital Roots.入门练习题,用C语言编写程序. 问题简述:输入若干正整数,求其数根,直到输入为0为止. 问题分析:数根是指整数的各个位的数字之和.如果其和为1位整数,则为结果:如果其和为多位整数,则再将各位数字相加,直到其和为1位数为止.这个问题的大陷阱是,没有指出整数是多少位的.即使使用unsignde long long类型,也可能会

ACM/ICPC 之 四道MST-Prim解法(POJ1258-POJ1751-POJ2349-POJ3026)

四道MST,适合Prim解法,也可以作为MST练习题. 题意包括在代码中. POJ1258-Agri Net 水题 1 //Prim-没什么好说的 2 //接受一个邻接矩阵,求MST 3 //Time:0Ms Memory:220K 4 #include<iostream> 5 #include<cstring> 6 #include<cstdio> 7 #include<algorithm> 8 using namespace std; 9 #define

hdu4521 小明系列问题——小明序列(LIS变种 (线段树+单点更新解法))

链接: huangjing 题目:中文题目 思路: 这个题目如果去掉那个距离大于d的条件,那么必然是一个普通的LIS,但是加上那个条件后就变得复杂了.用dp的解法没有看懂,我用的线段树的解法...就是采用延迟更新的做法,用为距离要大于d啊,所以我们在循环到第i的时候,就对(i-d-1)这个点进行更新,因为如果在(i-d-1)这个点更新了,会对后面的造成影响,然后线段树的tree[]数组存的是以i结尾的最长lis,那么每次询问的时候就找最大的tree[]就可以了... 代码: 小明系列问题--小明

最近点对问题的解法

回到二维的情况.参考一维退化版本的解法,首先进行分解和求解步骤. 分解:将所有点按照横坐标从中心分成两部分. 求解:递归求解两半部分的最近点对,左右两半部分依次为 p_1,p_2p?1??,p?2?? 和 q_1,q_2q?1??,q?2??. 合并:令 d = min\{|p_1-p_2|, |q_1-q_2|\}d=min{∣p?1??−p?2??∣,∣q?1??−q?2??∣}.如下左图所示,除了两半部分各自的最近点对之外,只需要计算横跨图中左右两半部分区域的点对的距离最小值,再和 dd 

kruskal算法求最小生成树(jungle roads的kruskal解法)

注意: 注意数组越界问题(提交出现runtimeError代表数组越界) 刚开始提交的时候,边集中边的数目和点集中点的数目用的同一个宏定义,但是宏定义是按照点的最大数定义的,所以提交的时候出现了数组越界问题,以后需要注意啦. Description The Head Elder of the tropical island of Lagrishan has a problem. A burst of foreign aid money was spent on extra roads betwe

实现杨辉三角的10种解法--体验Python之美

本文收集了使用python实现杨辉三角的多种解法,主要为网上收集,也有一些是自己写的.从中可以体会python编写一个算法的不同思想和Python语法的特点. 杨辉三角是什么?还是度娘吧,看起来像是这样的:                          1                          1   1                           1   2   1                         1   3   3   1