宇宙系统论的宇宙系统论——拉普拉斯方程

作者&投稿:荤亲 (若有异议请与网页底部的电邮联系)

以法国 P.-S. 拉普拉斯命名的二阶偏微分方程。在三维直角坐标系中,它的形式是:它的二次连续可微解称为调和函数,调和函数有极多的光滑性。拉普拉斯方程在物理吸广泛应用,因为它的解出现在电、磁、引力位势、稳态温度以及流体动力学各方面的问题中 。
拉普拉斯方程,又名调和方程,是一种偏微分方程。因为由法国数学家拉普拉斯首先提出而得名。求解拉普拉斯方程是电磁学、天文学和流体力学等领域经常遇到的一类重要的数学问题,因为这种方程以势函数的形式描写了电场、引力场和流场等物理对象(一般统称为“保守场”或“有势场”)的性质。
三维情况下,拉普拉斯方程可由下面的形式描述,问题归结为求解对实自变量x、y、z二阶可微的实函数φ :
<math>
{\partial^2 \varphi\over \partial x^2 } + {\partial^2 \varphi\over \partial y^2 } + {\partial^2 \varphi\over \partial z^2 } = 0. </math>
上面的方程常常简写作:
<math>
abla^2 \varphi = 0 </math>

<math>\operatorname{div}\,\operatorname{grad}\,\varphi = 0, </math>
其中div表示矢量场的散度(结果是一个标量场),grad表示标量场的梯度(结果是一个矢量场),或者简写作:
<math>\Delta \varphi = 0</math>
其中Δ称为拉普拉斯算子.
拉普拉斯方程的解称为调和函数。
如果等号右边是一个给定的函数f(x, y, z),即:
<math>\Delta \varphi = f</math>
则该方程称为泊松方程。拉普拉斯方程和泊松方程是最简单的椭圆形偏微分方程。偏微分算子<math>
abla^2</math>或<math>\Delta</math>(可以在任意维空间中定义这样的算子)称为拉普拉斯算子,英文是 Laplace operator 或简称作 Laplacian。
拉普拉斯方程的狄里克雷问题可归结为求解在区域<math>D</math>内定义的函数φ,使得<math>\varphi</math>在<math>D</math>的边界上等于某给定的函数。为方便叙述,以下采用拉普拉斯算子应用的其中一个例子——热传导问题作为背景进行介绍:固定区域边界上的温度(是边界上各点位置坐标的函数),直到区域内部热传导使温度分布达到稳定,这个温度分布场就是相应的狄里克雷问题的解。
拉普拉斯方程的诺依曼型边界条件不直接给出区域<math>D</math>边界处的温度函数φ本身,而是φ沿<math>D</math>的边界法向的导数。从物理的角度看,这种边界条件给出的是矢量场的势分布在区域边界处的已知效果(对热传导问题而言,这种效果便是边界热流密度)。
拉普拉斯方程的解称为调和函数,此函数在方程成立的区域内是解析的。任意两个函数,如果它们都满足拉普拉斯方程(或任意线性微分方程),这两个函数之和(或任意形式的线性组合)同样满足前述方程。这种非常有用的性质称为叠加原理。可以根据该原理将复杂问题的已知简单特解组合起来,构造适用面更广的通解。 两个自变量的拉普拉斯方程具有以下形式:
<math>\varphi_{xx} + \varphi_{yy} = 0.\,</math> 解析函数的实部和虚部均满足拉普拉斯方程。换言之,若z = x + iy,并且
<math>f(z) = u(x,y) + iv(x,y),\,</math>
那么f(z)是解析函数的充要条件是它满足下列柯西-黎曼方程:
<math>u_x = v_y, \quad v_x = -u_y.\,</math>
上述方程继续求导就得到
<math>u_{yy} = (-v_x)_y = -(v_y)_x = -(u_x)_x.\,</math>
所以u 满足拉普拉斯方程。类似的计算可推得v 同样满足拉普拉斯方程。
反之,给定一个由解析函数(或至少在某点及其邻域内解析的函数)f(z)的实部确定的调和函数,若写成下列形式:
<math>f(z) = \varphi(x,y) + i \psi(x,y),\,</math>
则等式
<math>\psi_x = -\varphi_y, \quad \psi_y = \varphi_x.\,</math>
成立就可使得柯西-黎曼方程得到满足。上述关系无法确定ψ,只能得到它的微增量表达式:
<math>d \psi = -\varphi_y\, dx + \varphi_x\, dy.\,</math>
φ满足拉普拉斯方程意味着ψ满足可积条件:
<math>\psi_{xy} = \psi_{yx},\,</math>
所以可以通过一个线积分来定义ψ。可积条件和斯托克斯定理的满足说明线积分的结果与积分经过的具体路径无关,仅由起点和终点决定。于是,我们便通过复变函数方法得到了φ和ψ这一对拉普拉斯方程的解。这样的解称为一对共轭调和函数。这种构造解的方法只在局部(复变函数f(z))的解析域内)有效,或者说,构造函数的积分路径不能围绕有f(z)的奇点。譬如,在极坐标平面(r,θ)上定义函数
<math>\varphi = \log r, \,</math>
那么相应的解析函数为
<math>f(z) = \log z = \log r + iheta. \,</math>
在这里需要注意的是,极角θ 仅在不包含原点的区域内才是单值的。
拉普拉斯方程与解析函数之间的紧密联系说明拉普拉斯方程的任何解都无穷阶可导(这是解析函数的一个性质),因此可以展开成幂级数形式,至少在不包含奇点的圆域内是如此。这与波动方程的解形成鲜明对照,后者包含任意函数,其中一些的可微分阶数是很小的。
幂级数和傅里叶级数之间存在着密切的关系。如果我们将函数f 在复平面上以原点为中心,R 为半径的圆域内展开成幂级数,即
<math>f(z) = \sum_{n=0}^\infty c_n z^n,\,</math>
将每一项系数适当地分离出实部和虚部
<math>c_n = a_n + i b_n.\,</math>
那么
<math>f(z) = \sum_{n=0}^\infty \left[ a_n r^n \cos n heta - b_n \sin n hetaight] + i \sum_{n=1}^\infty \left[ a_n \sin nheta + b_n \cos n hetaight],\,</math>
这便是f 的傅里叶级数。 设u、v 分别为满足定常、不可压缩和无旋条件的流体速度场的x 和y 方向分量(这里仅考虑二维流场),那么不可压缩条件为:
<math>u_x + v_y=0,\,</math>
无旋条件为:
<math>v_x - u_y =0. \,</math>
若定义一个标量函数ψ,使其微分满足:
<math>d \psi = v\, dx - u\, dy,\,</math>
那么不可压缩条件便是上述微分式的可积条件。积分的结果函数ψ称为流函数,因为它在同一条流线上各点的值是相同的。ψ的一阶偏导为:
<math>\psi_x = v, \quad \psi_y=-u, \,</math>
无旋条件即令 ψ 满足拉普拉斯方程。ψ的共轭调和函数φ称为速度势。 柯西-黎曼方程要求
<math>\varphi_x=-u, \quad \varphi_y=-v. \,</math>
所以每一个解析函数都对应着平面内的一个定常不可压缩无旋流场。解析函数的实部为速度势函数,虚部为流函数。
[编辑]在电磁学中的应用
根据麦克斯韦方程组,二维空间中不随时间变化的电场(u,v)满足:
<math>
abla imes (u,v) = v_x -u_y =0,\,</math>

<math>
abla \cdot (u,v) = ho,\,</math>
其中ρ为电荷密度。第一个麦克斯韦方程便是下列微分式的可积条件:
<math>d \varphi = -u\, dx -v\, dy,\,</math>
所以可以构造电势函数φ使其满足
<math>\varphi_x = -u, \quad \varphi_y = -v.\,</math>
第二个麦克斯韦方程即:
<math>\varphi_{xx} + \varphi_{yy} = -ho,\,</math>
这是一个泊松方程。 基本解
拉普拉斯方程的基本解满足
<math>
abla \cdot
abla u = u_{xx} + u_{yy} + u_{zz} = -\delta(x-x',y-y',z-z'), \,</math>
其中的三维δ函数代表位于<math> (x',\, y', \, z')</math>的一个点源。由基本解的定义,若对u 作用拉普拉斯算子,再把结果在包含点源的任意体积内积分,那么
<math> \iiint_V
abla \cdot
abla u dV =-1. \,</math>
由于坐标轴旋转不改变拉普拉斯方程的形式,所以基本解必然包含在那些仅与到点源距离r 相关的解中。如果我们选取包含点源、半径为a 的球形域作为积分域,那么根据高斯散度定理
<math> -1= \iiint_V
abla \cdot
abla u \, dV = \iint_S u_r dS = 4\pi a^2 u_r(a).\,</math>
求得在以点源为中心,半径为r 的球面上有 <math> u_r(r) = -\frac{1}{4\pi r^2},\,</math>
所以 <math> u = \frac{1}{4\pi r}.\,</math> 经过类似的推导同样可求得二维形式的解 <math> u = \frac{-\log r}{2\pi}. \,</math>



拉普拉斯方程是什么啊~


拉普拉斯方程为:△u=d^2u/dx^2+d^2u/dy^2=0,其中△为拉普拉斯算子,这里的拉普拉斯方程为二阶偏微分方程。

Pierre Simon Laplace的生平介绍?
答:拉普拉斯(Laplace,Pierre-Simon,marquisde),法国著名数学家和天文学家,拉普拉斯是天体力学的主要奠基人,是天体演化学的创立者之一,是分析概率论的创始人,是应用数学的先躯。拉普拉斯用数学方法证明了行星的轨道大小只有周期性变化,这就是著名拉普拉斯的定理。他发表的天文学、数学和物理学的论文有270...