有限差分法求解拉普拉斯方程
收稿日期:2009-09-01
第一作者简介:贾新民(1956-),男,四川邻水人,新疆昌吉学院计算机工程系,副教授,研究方向:计算机程序设计及其
语言教学和理论物理研究。
有限差分法求解拉普拉斯方程
贾新民1 严
文2(1.昌吉学院计算机工程系新疆昌吉831100;2.昌吉学院物理系 新疆
昌吉831100)摘 要:以极板上具有半圆截面沟槽的电容器内的电势分布为例,介绍了综合应用计算机软件利用有限
差分法求解复杂边界的拉普拉斯方程数值解的方法。并利用数值解的结果讨论了沟槽表面的电场分布和电荷分布。
关键词:拉普拉斯方程;有限差分法;五点差分格式
中图分类号:O411.2 文献标识码:A 文章编号:1671-6469(2009)05-0105-05
1 引言
无源空间的引力场、静电场、稳定的温度分布等问题都满足拉普拉斯(Laplace )方程 2u (x ,y ,z )=0(1)但由于方程(1)是偏微分方程,只有在问题具有高度对称的情况下,才能求出解析解,而这种情形是极少
的。有些情形看上去很简单,但却求不出解析解。对于这些情况,只能寻求数值解。
2 计算机数值解法方案文献[1][2][3]给出了拉普拉斯方程数值解的方法———有限差分法。
有限差分法的思想是用差分Δu (x +Δx
,y )Δx ,Δu (x ,y +Δy )Δy 代替导数5u 5x ,5u 5y
,用网格将求解区域覆盖,对于平面拉普拉斯方程,第i 行第j 列小格的电势由Laplace 方程的五点差分格式给出。
u ij =14(u ij -1+u ij +1+u i -1j +u i +1j )(2)
考虑图1所示的具有半圆形截面的槽的电容器内部的电势和电场分布。
为了能够对坑(槽)内部的电场进行比较细致的观察,
应该将半径R 取的大些,为了满足无限远的条件,应该使
求解区域尽量大些。
我们选Excel 为计算工具,因为Excel 具有不用编写
程序和直观的优点。Excel 的一个单元格代表求解区域的
一个网格,单元格的值表示该网格处的电势。根据Excel
的能力,取R =28,极板宽度为16×28=448,两极板距离
d 取8×28=224。由于Excel 的最大列数为256(A -
IG ),最大行数为65536,为了能够用Excel 计算,需要将问题的行列置换,即将224+28=252行换成252列,而将448列换成448行。
5
01
初始化:初始化的目的是设置边界上的电势。问题的边界为:第1列(A列)为电容器上极板,电位为-u0(设为-223V);第224列(电容器下极板)电位为0;第1行和第448行表示离坑槽无限远处,电势应该是均匀电场的电势,这两行从第1列到第224列的电势初试化为依次为-223,-222…,0;半圆形坑(槽)的边界电势为0。区域内部电位的初始值可以设为任意值,假设为0。
3 跌代计算
选中区域内部的任一单元格,输入公式(2),如图2所示为选择了单元格HN191的情形,单元格编辑栏显示了输入的计算公式。
图 2
公式输入完成后,再选中刚才的单元格,复制,按住Ct rl键,用鼠标拖动的方法选中除边界的内部所有单元格,粘贴,这时将弹出一个“循环引用”的错误对话框,按照对话框的提示,打开“工具”菜单下的“选项”,在“选项”对话框中选择“重新计算”选项卡,将“跌代计算”选中,将“最多跌代次数”设置为30000,默认最大误差为0.001,确定(如图3),经过足够长的时间(跌代次数大约为25000),计算完成。
图 3
4 绘制等位线
利用Matlab的contour函数很容易画出等位线。具体作法是:先将Excel计算的数据表中没有数据的地方(如坑的外部)补0,使有数据的区域成为矩形区域,保存。打开Matlab,利用文件菜单下的“Import data…”命令将保存好的Excel数据导入Matlab工作空间。将数据转置(设转置后的矩阵为601
u ),在Matlab 命令窗口中输入:contour (u ,n );(其中n 是欲画等位线的数量),再将表示带坑槽的电容器极板也画在同一张图上,等位线就绘制完成了。如图4中实线所示,除了标注为-1V 的等位线外,图中每条等位线之间的间隔约为6.67V
。
图 4
5 绘制电场线
Matlab 的gradient 函数和quiver 函数可以计算并绘制出表示电场强度的小箭头族,但绘制的箭头太密,不便于观察。图5就是用[ex ,ey ]=gradient (u );quiver (ex ,ey )语句绘制的电场图在平面和圆相交附近的放大图。
图 5
根据电场强度和电位的关系:E =-
u 知,电场中某点的电场强度的方向(即
电场线的切线方向)就是由该点出发、电位
下降最快的方向。因此,只要知道电场线
的出发点或终止点,从该点出发,连接与该
点最临近的等位线上的点,便确定了电场
线上的第二个点,以此类推,就可以绘制出
整个电场线。
我们知道,电容器上极板附近的电场
可视为匀强电场,并且是电场线的终止点。
在上极板处等间隔地取一些点(例如取65
个),在Matlab 的命令窗口,用[c ,h]=con 2
tur (u ,n )绘制等位线,然后用get (h (k ),’
xdata ’);和get (h (k ),’ydata ’)获取第i 条等位线的横坐标和纵坐标。再用比较相邻等位线距离的方法确定出65条电场线的坐标,然后连接各组点。图4中的虚线就是用Matalb 从448条等位线(相邻等位线的间隔为0.5V )获得的数据,绘制出来的电场线(图4中显示的等位线间隔约为6.67V )。
6 结果分析
从图4可以看出,离电容器下极板1.43R 处(从槽底数第7条等位线处)坑槽对电场的影响就已经很小了。从迭代计算的数据也可以看出这一点,该处边缘的电势为-40V ,相同纵坐标处电势变化最大为-37.181V ,变化率仅为7%。
7
01
布情况,根据面电荷密度和电场强度的关
系,只要知道了坑(槽)壁附近的电场强度,
就可以知道面电荷分布。为了找出坑(槽)
壁附近电场强度的分布规律,将-1V等位
线的第200-248行数据用Excel绘制出
来,并用添加趋势线的方法进行回归分析,
选择趋势线类型为2次多项式,并要求显示
公式和R2(R2越接近1,说明趋势线和原始
数据之间的误差就越小)。回归结果如图6。
用几何画板绘制出坑底半圆:
y=-sqrt(28^2-(x-224)^2)
同时绘制出-1V等位线的回归线,如图7。
分别构造圆弧上的点B和等位线上的
B点,分别度量出BC的距离,用鼠标拖
动C点使度量的距离最短。注意两点
距离最短时C点并不一定在等角度线
上。
最后将这样得到的角度和距离输入
Excel,距离的倒数1/d正比于电场强度。
计算出cos(theta)和1/d,并绘制出它们
的关系(如图8)。从度量结果可以看出,
84°处的电场强度与坑底中心的电场强度
比值为2.802/0.379≈7.4。
说明两处的面电荷密度之比也约
为7.4。图8还显示了相对电场强度与
圆心角余弦之间的回归关系:
E∝co s-0.8213(t heta)
当然,这个结论的准确性还有待证
实。
从跌代计算的结果看,坑底中心附近的电势差为0.198,而远离坑的平坦处电势差为1,表明坑底处电荷面密度仅为平坦处的20%。
从图2(图2中的数据实际上是跌代计算完成后的数据)显示的计算结果看,如果将H P197单元格的电位(-1.699)看成是坑沿处电场强度的水平分量,而将HO195单元格的电势(-1.903)看成是坑沿处电场强度的垂直分量,则坑沿处的电场强度(相对单位)约为
1.6992+1.9032=
2.531
因此,坑沿的电场强度与坑底的电场强度之比约为2.531/0.198≈12.8。即两处的面电荷密度比约为12.8。
801
7 结束语
前已述及,本文计算的数据和绘制的图形都不是半球形坑电容器极板间的电场,而是具有半圆形截面的槽的电容器极板间的电场。至于半球坑问题应当如何解,还有待于找出轴对称情形下相邻点电势之间的数值关系。
关于为什么用Excel而不用Matlab计算电位的问题,前面已经说过,除了Excel具有直观和无需编程的好处外,还有一点需要提及的就是Matlab只能显示10000个数据,不能一下子显示出全部数据(约10万个),而Excel则可以将全部计算结果显示在屏幕上,
有限差分法和Laplace方程的五点差分格式不仅可以求解规则边界的Laplace方程,而且可以求解不规则二维Laplace方程。
既然是数值解法,就一定伴随着误差。用差分代替求导本身就具有近似性质,这种替代带来的误差可以参见本文参考文献。[4][5]另一方面的误差来自对无限远的理解。本文计算出的数据是建立在将离坑槽中心8倍半径处看成无限远的基础上的。当然,如果将不受坑槽缺陷影响的距离取得更远,计算准确性就更好些,但这受到计算机软硬件的限制。除此而外,边界的数字化也将导致计算误差。本文采用的是四舍五入取整法(使用Excel的Round函数)对边界进行数字化的。还可以使用Ceiling函数(将圆弧上的坐标值向上舍入)或Floor函数(将圆弧上的坐标值沿绝对值减小的方向向下舍入)对边界进行数字化处理。不同的数字化处理方法的结果可能有所不同。
最后需要指出的是,用连接相邻等位线距离最近处的方法绘制电场线的方法,实际上具有很大的近似性,图4中央处的垂线是人工添加上去的。理论上说等位线取的越密集,电场线上的点就越多,绘制的电场线应该越准确。但等位线取的太多,又使得电场线局部锯齿状加剧,为了使电场线光滑,需要对等位线进行插值处理。究竟等位线应该取多密,才能使误差最小,是一个值得研究的问题。
参考文献:
[1][4]张文生.科学计算中的偏微分方程有限差分法[M].北京:高等教育出版社,2006.
[2][5]陆金甫,关治偏微分方程数值解法(第2版)[M].北京:清华大学出版社,2005.
[3]梁昆淼.数学物理方法(第三版)[M].北京:高等教育出版社,1998:462.
(责任编辑:马海燕)
901
(完整版)大连理工大学高等数值分析抛物型方程有限差分法
抛物型方程有限差分法 1. 简单差分法 考虑一维模型热传导方程 (1.1) )(22x f x u a t u +??=??，T t ≤<0 其中a 为常数。)(x f 是给定的连续函数。（1.1）的定解问题分两类： 第一，初值问题（Cauchy 问题）：求足够光滑的函数()t x u ,，满足方程（1.1）和初始条件： （1.2） ()()x x u ?=0,, ∞<<∞-x 第二，初边值问题（也称混合问题）：求足够光滑的函数()t x u ,，满足方程（1.1）和初始条件： ()13.1 ()()x x u ?=0,, l x l <<- 及边值条件 ()23.1 ()()0,,0==t l u t u ， T t ≤≤0 假定()x f 和()x ?在相应的区域光滑，并且于()0,0，()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。
现在考虑边值问题（1.1），（1.3）的差分逼近 取 N l h = 为空间步长，M T = τ为时间步长，其中N ，M 是 自然数， jh x x j ==, ()N j ,,1,0Λ=； τ k y y k ==, ()M k ,,1,0Λ= 将矩形域G {}T t l x ≤≤≤≤=0;0分割成矩形网格。其中 ()j i y x ,表 示网格节点； h G 表示网格内点（位于开矩形G 中的网格节点）的集合； h G 表示位于闭矩形G 中的网格节点的集合； h Γ表示h G -h G 网格边界点的集合。 k j u 表示定义在网点()k i t x ,处的待求近似解，N j ≤≤0，M k ≤≤0。 注意到在节点()k i t x ,处的微商和差商之间的下列关系 （(,)k j k j u u x t t t ????≡ ? ????）： ()() ()ττ O t u t x u t x u k j k j k j +??? ????=-+,,1 ()() ()2112,,ττ O t u t x u t x u k j k j k j +??? ????=--+ ()()()h O x u h t x u t x u k j k j k j +??? ????=-+,,1 ()() ()h O x u h t x u t x u k j k j k j +??? ????=--,,1 ()() ()2112,,h O x u h t x u t x u k j k j k j +??? ????=--+ ()()() ()2 222 11,,2,h O x u h t x u t x u t x u k j k j k j k j +???? ????=+--+ 可得到以下几种最简差分格式
差分方程的解法
1、常系数线性差分方程的解 方程（ 8）其中为常数，称方程（8）为常系数线性方程。 又称方程（9） 为方程（8）对应的齐次方程。 如果（9）有形如的解，带入方程中可得： （10） 称方程（10）为方程（8）、（9）的特征方程。 显然，如果能求出（10）的根，则可以得到（9）的解。 基本结果如下： （1）若（10）有k个不同的实根，则（9）有通解： ， （2）若（10）有m重根，则通解中有构成项： （3）若（10）有一对单复根，令：，，则（9）的通解中有构成项： （4）若有m 重复根：，，则（9）的通项中有成项：
综上所述，由于方程（10）恰有k 个根，从而构成方程 （9）的通解中必有k个独立的任意常数。通解可记为： 如果能得到方程（8）的一个特解：，则（8）必有通解： + （11） （1）的特解可通过待定系数法来确定。 例如：如果为n 的多项式，则当b不是特征根时，可设成形如形式的特解，其中为m次多项式；如果b是r重根时，可设特解：，将其代入（8）中确定出系数即可。 2、差分方程的z变换解法 对差分方程两边关于取Z变换，利用的Z 变换F（z）来表示出的Z变换，然后通过解代数方程求出F（z），并把F(z)在z=0的解析圆环域中展开成洛朗级数，其系数就是所要求的 例1设差分方程，求 解：解法1：特征方程为，有根： 故：为方程的解。 由条件得： 解法2：设F（z）=Z(),方程两边取变换可得：
由条件得 由F（z）在中解析，有 所以， 3、二阶线性差分方程组 设，，形成向量方程组 （12）则 （13）（13）即为（12）的解。 为了具体求出解（13），需要求出，这可以用高等代数的方法计算。常用的方法有： （1）如果A为正规矩阵，则A必可相似于对角矩阵，对角线上的元素就是A的特征值，相似变换矩阵由A的特征向量构成：。 （2）将A 分解成为列向量，则有 从而，
抛物形扩散方程的有限差分法及数值实例
偏微分方程数值解 所在学院：数学与统计学院 课题名称：抛物形扩散方程的有限差分法及数值实例学生姓名：向聘
抛物形扩散方程的有限差分法及数值实例 1.1抛物型扩散方程 抛物型偏微分方程是一类重要的偏微分方程。考虑一维热传导方程： 22(),0u u a f x t T t x ??=+<≤?? (1.1.1) 其中a 是常数，()f x 是给定的连续函数。按照初边值条件的不同给法，可将(1.1.1)的定解分为两类： 第一，初值问题（Cauchy 问题）：求足够光滑的函数()t x u ,，满足方程(1.1.1)和初始条件： ()()x x u ?=0,, ∞<<∞-x (1.1.2) 第二，初边值问题（也称混合问题）：求足够光滑的函数()t x u ,，满足方程(1.1.1)和初始条件： ()()x x u ?=0,, 0x l << (1.1.3) 及边值条件 ()()0,,0==t l u t u ， T t ≤≤0 (1.1.4) 假定()x f 和()x ?在相应的区域光滑，并且于()0,0，()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。 1.2抛物线扩散方程的求解 下面考虑如下热传导方程 22()(0.)(,)0(,0)()u u a f x t x u t u L t u x x ????=+????? ==??=??? (1.2.1) 其中，0x l <<，T t ≤≤0,a (常数)是扩散系数。 取N l h = 为空间步长，M T =τ为时间步长，其中N ，M 是自然数，用两族
平行直线jh x x j ==, ()N j ,,1,0 =和k t t k τ ==, ()M k ,,1,0 =将矩形域 G {}T t l x ≤≤≤≤=0;0分割成矩形网格。其中 (),j k x t 表示网格节点；h G 表示 网格内点（位于开矩形G 中的网格节点）的集合；h G 表示位于闭矩形G 中的网格节点的集合；h Γ表示h G -h G 网格边界点的集合。 k j u 表示定义在网点(),j k x t 处的待求近似解，N j ≤≤0，M k ≤≤0。 现在对方程进行差分近似： (一) 向前差分格式 =-+τ k j k j u u 111 2 2(())k k k j j j j j j u u u a f f f x h +--++= (1.2.2) ()j j j x u ??==0， k u 0=k N u =0 (1.2.3) 计算后得： 111(12)k k k k j j j j j u ru r u ru f τ++-=+-++ (1.2.4) 其中，2 a r h τ = ，1,,1,0-=N j ，1,,1,0-=M k 。 显然，这是一个四点显示格式，每一层各个节点上的值是通过一个方程组求解到的。方程组如下： 1000 121011000 232121000 3432310001121(12)(12)(12)(12)N N N N N u ru r u ru f u ru r u ru f u ru r u ru f u ru r u ru f ττττ----?=+-++?=+-++??=+-++? ???=+-++? (1.2.5) 若记 () T k N k k k u u u 1 21,,,-= u ，()()()()T N x x x 121,,,-=???? ，()()()()T N x f x f x f 121,,,-=τττ f 则显格式(1.2.4)可写成向量形式 10 ,0,1,,1 k k k M φ +?=+=-?=? u Au f u (1.2.6) 其中
拉普拉斯方程
拉普拉斯方程 拉普拉斯方程 拉普拉斯方程，又名调和方程，是一砍。因为由法 国数学家首先提出而得名。求解拉普 拉斯方程栯、和等领域经常遇到的 一类重要的数学问領，因为这种方程以的形式描写了、和等物理对象（一般统称为“保守场”栖“有势场”）的性质。三维情况下＠拉普拉斯方程可由下面的形式描 述，闠题归结为求解对实自变量x、y、z二阶的实函数 φ ：: + + = 0. 上面的方程常常简写作：: \nabla^2 \varphi = 0 或: \operatorname\,\operatorname\,\varphi = 0, 其中div表示的（结果是一个），grad 表示标量场的（结果是一个矢量场），或者简写作 : \Delta \varphi = 0 其中Δ称为 . 拉普拉斯方 程的解称为。如果等号右边是一个给定的函数f( , y, z)，即：: \Delta \varphi = f 则该方程称为。拉 普拉斯方程和泊松方程是最简单?? 。偏微 分算子\nabla^2或\Delta（可以在任意维空间中定义这样的算 堐）称为，英文是Laplace operator 或简称作 Laplacian。拉普拉斯方程的可归结为求解在区 域D内定义的函数φ，使得\varphi在D的边界上等于某给定 的函数。为方便堙述，以下采用拉普拉斯算子应用的其?一个例
子——作为背景进行介绍：固定区域边界上砄温度（是边界上各点位置坐标的函数 ，直到区域内部热传导使温度 分布达堰稳定，这个温度分布场就是相应的狄頌克雷问题的解。拉普拉斯方程的不直接给出区域D边界处的温度函数φ本身，而是φ ??D的边界法向的。从物理 的角度看，这种边界条件给堺的是矢量场的势分布在区域边界处的堲知效果（对热传导问题而言，这种效栜便是边界热流密度）。拉普拉斯方稠的解称为，此函数在方程成立的区域内是。任意两个函数，如果它们都满足拉栮拉斯方程（或任意线性微分方程），蠙两个函数之和（或任意形式的线性组堈）同样满足前述方程。这种非常有用砄性质称为。可以根据该原理将复杂问题的已知砀单组合起来，构造适用面更广的。 二维拉普拉斯方程 两个自变量的拉普拉斯方程具有以丠形式：:\varphi_ + \varphi_ = 0.\, 解析函数 的实部和虚部均满足拉普拉斯方程。栢言之，若z = x + iy，并且:f(z) = u(x,y) + iv(x,y),\, 那么f(z)是解析函数的是它满足下列柯西-黎曼方程：:u_x = v_y, \quad v_x = -u_y.\, 上述方程继续求导就得到:u_ = (-v_x)_y = -(v_y)_x = -(u_x)_x.\, 所以u 满足拉普拉斯方程。类似的计算可推堗v 同样满足拉普拉斯方程。反之，给定一个由解析函数（或至尠在某点及其邻域内解析的函数）f(z)的堞部确定的调
差分方程的解法分析及MATLAB实现(程序)
差分方程的解法分析及MATLAB 实现（程序） 摘自：张登奇,彭仕玉.差分方程的解法分析及其MATLAB 实现[J]. 湖南理工学院学报.2014(03) 引言 线性常系数差分方程是描述线性时不变离散时间系统的数学模型,求解差分方程是分析离散时间系统的重要内容.在《信号与系统》课程中介绍的求解方法主要有迭代法、时域经典法、双零法和变换域 法[1]. 1 迭代法 例1 已知离散系统的差分方程为)1(3 1)()2(81)1(43)(-+=-+--n x n x n y n y n y ,激励信号为)()4 3()(n u n x n =,初始状态为21)2(4)1(=-=-y y ，.求系统响应. 根据激励信号和初始状态,手工依次迭代可算出24 59)1(,25)0(==y y . 利用MATLAB 中的filter 函数实现迭代过程的m 程序如下: clc;clear;format compact; a=[1,-3/4,1/8],b=[1,1/3,0], %输入差分方程系数向量,不足补0对齐 n=0:10;xn=(3/4).^n, %输入激励信号 zx=[0,0],zy=[4,12], %输入初始状态 zi=filtic(b,a,zy,zx),%计算等效初始条件 [yn,zf]=filter(b,a,xn,zi),%迭代计算输出和后段等效初始条件 2 时域经典法 用时域经典法求解差分方程:先求齐次解；再将激励信号代入方程右端化简得自由项,根据自由项形 式求特解；然后根据边界条件求完全解[3].用时域经典法求解例1的基本步骤如下. （1）求齐次解.特征方程为081432=+-αα,可算出4 1 , 2121==αα.高阶特征根可用MATLAB 的roots 函数计算.齐次解为. 0 , )4 1()21()(21≥+=n C C n y n n h （2）求方程的特解.将)()4 3()(n u n x n =代入差分方程右端得自由项为 ?????≥?==-?+-1,)4 3(9130 ,1)1()43(31)()43(1n n n u n u n n n 当1≥n 时,特解可设为n p D n y )4 3()(=,代入差分方程求得213=D . （3）利用边界条件求完全解.当n =0时迭代求出25)0(=y ,当n ≥1时,完全解的形式为 ,)4 3(213 )41()21()(21n n n C C n y ?++=选择求完全解系数的边界条件可参考文[4]选)1(),0(-y y .根据边界条件求得35,31721=-=C C .注意完全解的表达式只适于特解成立的n 取值范围,其他点要用 )(n δ及其延迟表示,如果其值符合表达式则可合并处理.差分方程的完全解为
拉普拉斯方程
拉普拉斯方程 拉普拉斯方程又名调和方程、位势方程，是一种偏微分方程。因为由法国数学家拉普拉斯首先提出而得名。求解拉普拉斯方程是电磁学、天文学和流体力学等领域经常遇到的一类重要的数学问题，因为这种方程以势函数的形式描写了电场、引力场和流场等物理对象（一般统称为“保守场”或“有势场”）的性质。 拉普拉斯方程表示液面曲率与液体压力之间的关系的公式。一个弯曲的表面称为曲面，通常用相应的两个曲率半径来描述曲面，即在曲面上某点作垂直于表面的直线，再通过此线作一平面，此平面与曲面的截线为曲线，在该点与曲线相重合的圆半径称为该曲线的曲率半径R1。 通过表面垂线并垂直于第一个平面再作第二个平面并与曲面相交，可得到第二条截线和它的曲率半径R2，用 R1与R2可表示出液体表面的弯曲情况。若液面是弯曲的，液体内部的压力p1与液体外的压力p2就会不同，在液面两边就会产生压力差△P= P1- P2，其数值与液面曲率大小有关，可表示为：▽p=γ（1/R1+1/R2）式中γ是液体表面张力。该公式成为拉普拉斯方程。 在数理方程中
拉普拉斯方程拉普拉斯方程为：Δ u=d^2u/dx^2+d^2u/dy^2=0，其中Δ为拉普拉斯算子，此处的拉普拉斯方程为二阶偏微分方程。三维情况下，拉普拉斯方程可由下面的形式描述，问题归结为求解对实自变量x、y、z二阶可微的实函数φ：其中Δ称为拉普拉斯算子. 拉普拉斯方程的解称为调和函数。 如果等号右边是一个给定的函数f(x, y, z)，即： 则该方程称为泊松方程。拉普拉斯方程和泊松方程是最简单的椭圆型偏微分方程。偏微分算子或Δ（可以在任意维空间中定义这样的算子）称为拉普拉斯算子，英文是 Laplace operator或简称作Laplacian。 狄利克雷问题 拉普拉斯方程的狄利克雷问题可归结为求解在区域D内定义的函数φ，使得在D的边界上等于某给定的函数。为方便叙述，以下采用拉普拉斯算子应用的其中一个例子——热传导问题作为背景进行介绍：固定区域边界上的温度（是边界上各点位置坐标的函数），直到区域内部热传导使温度分布达到稳定，这个温度分布场就是相应的狄利克雷问题的解。
时间序列分析讲义 第01章 差分方程
第一章 差分方程 差分方程是连续时间情形下微分方程的特例。差分方程及其求解是时间序列方法的基础，也是分析时间序列动态属性的基本方法。经济时间序列或者金融时间序列方法主要处理具有随机项的差分方程的求解问题，因此，确定性差分方程理论是我们首先需要了解的重要内容。 §1.1 一阶差分方程 假设利用变量t y 表示随着时间变量t 变化的某种事件的属性或者结构，则t y 便是在时间t 可以观测到的数据。假设t y 受到前期取值1-t y 和其他外生变量t w 的影响，并满足下述方程： t t t w y y ++=-110φφ (1.1) 在上述方程当中，由于t y 仅线性地依赖前一个时间间隔自身的取值1-t y ，因此称具有这种结构的方程为一阶线性差分方程。如果变量t w 是确定性变量，则此方程是确定性差分方程；如果变量t w 是随机变量，则此方程是随机差分方程。在下面的分析中，我们假设t w 是确定性变量。 例1.1 货币需求函数 假设实际货币余额、实际收入、银行储蓄利率和商业票据利率的对数变量分别表示为t m 、t I 、bt r 和ct r ，则可以估计出美国货币需求函数为： ct bt t t t r r I m m 019.0045.019.072.027.01--++=- 上述方程便是关于t m 的一阶线性差分方程。可以通过此方程的求解和结构分析，判断其他外生变量变化对货币需求的动态影响。 1.1.1 差分方程求解：递归替代法 差分方程求解就是将方程变量表示为外生变量及其初值的函数形式，可以通过以前的数据计算出方程变量的当前值。 由于方程结构对于每一个时间点都是成立的，因此可以将(1.1)表示为多个方程： 0=t ：01100w y y ++=-φφ 1=t ：10101w y y ++=φφ t t =：t t t w y y ++=-110φφ 依次进行叠代可以得到： 1011211010110101)()1()(w w y w w y y ++++=++++=--φφφφφφφφ 0111122113121102)1(w w w y y φφφφφφφ++++++=- i t i i t t i i t w y y ∑∑=-=++=0 111 1 0φφφφ (1.2) 上述表达式(1.2)便是差分方程(1.1)的解，可以通过代入方程进行验证。上述通过叠代将 t y 表示为前期变量和初始值的形式，从中可以看出t y 对这些变量取值的依赖性和动态变化 过程。 1.1. 2. 差分方程的动态分析：动态乘子(dynamic multiplier) 在差分方程的解当中，可以分析外生变量，例如0w 的变化对t 阶段以后的t y 的影响。假设初始值1-y 和t w w ,,1 不受到影响，则有：
拉普拉斯方程
拉普拉斯方程应该和泊松方程是同胞兄弟了，都是扩散方程，用来描述散度场的。只不过拉普拉斯方程是无源场，泊松方程是有源场。预备内容：梯度、旋度、散度和拉普拉斯算子在曲线坐标下的表达式： 如果在某个曲线坐标系内位移微元（其中是坐标），那么便有： 梯度：散度：旋度：拉普拉斯算符： 对于直角坐标系、球坐标系和柱坐标系来说，的值为： 于是，我们便可以轻松地默写球坐标下拉普拉斯算符的表达式\^o^/ 下面进入正题 1.直角坐标系 当出现金属平板之类的边界条件时，使用直角坐标系较为方便。 在直角坐标系下，拉普拉斯方程的表达式为： i）二维问题 假设沿z轴平移V保持不变，于是方程便简化为二维形式： 我们假设V可以写成两个函数相乘的形式： （乍看之下这不是一个很合理的假设。但是我们很快可以看到为什么可以这样做）
代入原方程并在两边除以V: 因为两部分之和为0，因此我们可以假设一个是正数另一部分是负数：（这里以含x的部分为正含y的部分为负为例） 很显然，这两个方程的解就是： 注记：这里决定哪一部分是正数哪一部分是负数要由边界条件来确定。比如说，沿x方向到达无限远时电势为零，x就应该含有指数衰减项，因此令含x的部分为正数。 于是，方程的一个解是 对所有可能的k求和，可以得到通解： 常数A,B,C,D的值需要由边界条件来确定。通常情况下，通过边界条件可以把k化成含有正整数的式子。将求和号改成对n求和，可以看到，第二个括号里的项便是傅里叶级数。狄利克雷定理保证了这个级数可以拟合任何边界条件。傅里叶系数可以由积分来确定。 ii)三维问题 三维问题的处理方法与二维的情形类似。 同样，假设是这种形式： 同样，代入方程并在两边同除以V:
拉普拉斯方程
拉普拉斯方程，也称为谐波方程和势方程，是一种偏微分方程，最早由法国数学家拉普拉斯提出。 拉普拉斯方程是液体表面曲率和液体表面压力之间关系的公式。 曲面称为曲面。通常，使用两个相应的曲率半径来描述表面，即在表面上的某个点处绘制垂直于该表面的直线，然后通过该线制作一个平面。平面和表面的截面是曲线，并且在该点与曲线相切的圆的半径称为曲线的曲率半径R1。第二剖面线及其曲率半径R2可以通过使第二平面垂直于第一平面并与表面相交来获得。液面的弯曲可以用R1和R2表示。如果液体表面弯曲，则液体P1内部的压力将与液体外部的压力P2不同，并且液体表面的两侧之间将存在压力差△P = P1-P2，这称为附加压力。压力。其值与液体表面的曲率有关，可以表示为：其中γ是液体的表面张力系数，称为拉普拉斯方程。 在数学公式中 拉普拉斯方程是：其中∥是拉普拉斯算子，而这里的拉普拉斯方程是二阶偏微分方程。在三维情况下，拉普拉斯方程可按以下形式描述。可以将问题简化为求解对于实变量X，y和Z可二阶微分的实函数φ ?2称为拉普拉斯算子。 拉普拉斯方程的解称为谐波函数。 如果在等号右边是给定的函数f（x，y，z），即： 然后将该方程称为泊松方程。拉普拉斯方程和泊松方程是最简单的椭圆偏微分方程。偏微分算子（可以在任何维空间中定义）称为拉
普拉斯算子。 方程解 它称为谐波函数，可以在建立方程的区域进行分析。如果任何两个函数满足拉普拉斯方程（或任何线性微分方程），则这两个函数的总和（或它们的任何线性组合）也满足上述方程。这种非常有用的特性称为叠加原理。根据这一原理，可以将已知的复杂问题的简单特殊解组合起来，以构建具有更广泛适用性的一般解。
泊松方程和拉普拉斯方程
拉普拉斯方程和泊松方程 摘要：拉普拉斯方程，又名调和方程、位势方程，是一种偏微分方程。因为由法国数学家拉普拉斯首先提出而得名。求解拉普拉斯方程是电磁学、天文学和流体力学等领域经常遇到的一类重要的数学问题，因为这种方程以势函数的形式描写了电场、引力场和流场等物理对象的性质。 关键词：分离变量电磁场拉普拉斯 简史 1777年，拉格朗日研究万有引力作用下的物体运动时指出：在引力体系中，每一质点的质量m k除以它们到任意观察点P的距离r k，并且把这些商加在一起，其总和 即P点的势函数，势函数对空间坐标的偏导数正比于在 P点的质点所受总引力的相应分力。1782年，P.S.M.拉普拉斯证明：引力场的势函数满足偏微分方程： ，叫做势方程，后来通称拉普拉斯方程。1813年，S.D.泊松撰文指出，如 果观察点P在充满引力物质的区域内部，则拉普拉斯方程应修改为，叫做泊松方程，式中ρ为引力物质的密度。文中要求重视势函数 V在电学理论中的应用，并指出导体表面为等热面。 静电场的泊松方程和拉普拉斯方程 若空间分区充满各向同性、线性、均匀的媒质，则从静电场强与电势梯度的关系E=-V高斯定理微分式，即可导出静电场的泊松方程： 式中ρ为自由电荷密度，纯数εr为各分区媒质的相对介电常数，真空介电常数εo=8.854×10-12法／米。在没有自由电荷的区域里，ρ＝0，泊松方程就简化为拉普拉斯方程。在各分区的公共界面上，V满足边值关系，
, 式中i ，j 指分界面两边的不同分区，σ 为界面上的自由电荷密度，n 表示边界面上的内法 线方向。 边界条件和解的唯一性 为了在给定区域内确定满足泊松方程以及边值关系的解，还需给定求解区域边界上的物 理情况，此情况叫做边界条件。有两类基本的边界条件：给定边界面上各点的电势，叫做狄 利克雷边界条件；给定边界面上各点的自由电荷 ，叫做诺埃曼边界条件。 静电场的唯一性定理： 设区域V 内给定自由电荷分布)(x ，在V 内电势满足泊松方程 或拉普拉斯方程,在V 的边界S 上给定电势 ,或V 边界上给定电势的法线方向偏导数 ,则V 内场（静电场）唯一确定。 除了静电场之外，在电学、磁学、力学、热学等领域还有许多服从拉普拉斯方程的势场。 各类物理本质完全不同的势场如果具有相似的边界条件，则因拉普拉斯方程解的唯一性，任 何一个势场的解，或该势场模型中实验测绘的等热面或流线图，经过对应物理量的换算之后，可以通用于其他的势场。 静磁场的泊松方程和拉普拉斯方程 在SI 制中，静磁场满足的方程为 ,式中j 为传导电流密度。第一式表明静磁 场可引入磁矢势r)描述： 。 在各向同性、线性、均匀的磁媒质中，传导电流密度j 0的区域里，磁矢势满足的方程 为 。 选用库仑规范，，则得磁矢势A 满足泊松方程 ，式中纯数μr 为媒质的相对磁导率， 真空磁导率μo =1.257×10-6亨／米。在传导电流密度j=0的区域里，上 式简化为拉普拉斯方程 。
差分方程求解
例题：已知差分方程51 (2)(1)()(+1)+0.5()66 x k x k x k r k r k +-++=，其中r (k )=1，k ≥0，x (0)=1， x (1)=2。 (1) 试由迭代法求其全解的前5项； (2) 分别由古典法求其零输入解、零状态解，以及全解； (3) 用Z 变换法求解差分方程。 解：注：解题过程中出现的下标“zi ”和“zs ”分别表示零输入条件和零状态条件。 1. 迭代法 题目中给出的条件仅仅是零输入初始条件，进行迭代求解时的初始条件应该是全解初始条件。 (1) 零输入初始条件 本题已给出零输入时的两个初始条件x zi (0)=1，x zi (1)=2。 (2) 零状态初始条件 取k =-2时，则51 (0)(1)(2)(1)0.5(2)66x x x r r --+-=-+-，得x zs (0)=0； 取k =-1 时，则51 (1)(0)(1)(0)0.5(1)66 x x x r r -+-=+-，求得x zs (1)=1。 (3) 全解初始条件 x (0)= x zi (0)+ x zs (0)=1； x (1)= x zi (1)+ x zs (1)=3。 (4) 根据求出的全解x (0)和x (1)，利用迭代法求解 取k =0时，则51(2)(1)(0)(1)0.5(0)66x x x r r -+=+，求得23(2)6x =； 取k =1时，则51(3)(2)(1)(2)0.5(1)66x x x r r -+=+，求得151 (3)36x =； 取k =2时，则51(4)(3)(2)(3)0.5(2)66x x x r r -+=+，求得941 (4)216 x =。 2. 古典法 (1) 零输入解 令输入为零，则得齐次方程 51 (2)(1)()066 x k x k x k +-++= (a) 根据差分方程定义的算子()()n d x k x k n =+，可得它的特征方程251 066 d d -+= 求得特征根为： 112d = ，21 3 d =
拉普拉斯方程
拉普拉斯方程 拉普拉斯方程（Laplace's equation）又称调和方程、位势方程，是一种偏微分方程，因由法国数学家拉普拉斯首先提出而得名。 [1] 拉普拉斯方程表示液面曲率与液体表面压强之间的关系的公式。 中文名 拉普拉斯方程 外文名 Laplace's equation 别称 调和方程、位势方程 提出者 拉普拉斯 关键词 微分方程、拉普拉斯定理 涉及领域 电磁学、天体物理学、力学、数学 目录 .1基本概述 .?在数理方程中 .?方程的解 .2二维方程 .3人物介绍
基本概述 一个弯曲的表面称为曲面，通常用相应的两个曲率半径来描述曲面，即在曲面上某点作垂直于表面的直线，再通过此线作一平面，此平面与曲面的截线为曲线，在该点与曲线相切的圆半径称为该曲线的曲率半径R1。通过表面垂线并垂直于第一个平面再作第二个平面并与曲面相交，可得到第二条截线和它的曲率半径R2，用R1与R2可表示出液体表面的弯曲情况。若液面是弯曲的，液体内部的压强p1与液体外的压强p2就会不同，在液面两边就会产生压强差△P= P1- P2，称附加压强，其数值与液面曲率大小有关，可表示为： ，式中γ是液体表面张力系数，该公式称为拉普拉斯方程。 在数理方程中 拉普拉斯方程为： ，其中?2为拉普拉斯算子，此处的拉普拉斯方程为二阶偏微分方程。三维情况下，拉普拉斯方程可由下面的形式描述，问题归结为求解对实自变量x、y、z二阶可微的实函数φ ： 其中?2称为拉普拉斯算子。 拉普拉斯方程的解称为调和函数。 如果等号右边是一个给定的函数f（x，y，z），即： 则该方程称为泊松方程。拉普拉斯方程和泊松方程是最简单的椭圆型偏微分方程。偏微分算子 （可以在任意维空间中定义这样的算子）称为拉普拉斯算子，英文是Laplace operator或简称作Laplacian。 方程的解 称为调和函数，此函数在方程成立的区域内是解析的。任意两个函数，如果它们都满足拉普拉斯方程（或任意线性微分方程），这两个函数之和（或任意形式的线性组合）同样满足前述方程。这种非常有用的性质称为叠加原理。可以根据该原理将复杂问题的已知简单特解组合起来，构造适用面更广的通解。 [2] 二维方程
迭代法
题目：Newton-Raphson 迭代法 （1）计算原理 （2）编出计算机程序 （3）给出算例（任意题型） (1)计算原理： 牛顿-拉夫森(Newton-Raphson)迭代法也称为牛顿迭代法，它是数值分析中最重要的方法之一，它不仅适用于方程或方程组的求解，还常用于微分方程和积分方程求解。 用迭代法解非线性方程时，如何构造迭代函数是非常重要的，那么怎样构造的迭代函数才能保证迭代法收敛呢？牛顿迭代法就是常用的方法之一，其迭代格式的来源大概有以下几种方式： 1设()[]2,f x C a b ∈，对()f x 在点[]0,x a b ∈，作泰勒展开： 略去二次项，得到()f x 的线性近似式：()()()()000f x f x f x x x '≈+- 由此得到方程()0f x =的近似根(假定()00f x '≠)，() () 000f x x x f x =-' 即可构造出迭代格式(假定()00f x '≠)：() () 1k k k k f x x x f x +=- ' 这就是牛顿迭代公式，若得到的序列{}k x 收敛于α，则α就是非线性方程的根。 2 牛顿迭代法 牛顿切线法，这是由于()f x 的线性化近似函数()()()()000l x f x f x x x '≈+-是曲线()y f x =过点()()00,x f x 的切线而得名的，求()f x 的零点代之以求() l x !2))((''))((')()(2 0000x x f x x x f x f x f -+ -+= ξ
的零点，即切线与x 轴交点的横坐标，如左图所示，这就是牛顿切线法的几何解释。实际上，牛顿迭代法也可以从几何意义上推出。利用牛顿迭代公式，由 k x 得到1k x +，从几何图形上看，就是过点()(),k k x f x 作函数()f x 的切线k l ，切线k l 与x 轴的交点就是1k x +，所以有()() 1 k k k k f x f x x x +'=-，整理后也能得出牛顿迭 代公式： 3 要保证迭代法收敛，不管非线性方程()0f x =的形式如何，总可以构造： 作为方程求解的迭代函数。因为： 而且 在根附近越小，其局部收敛速度越快，故可令： 若0(即根不是0的重根)，则由得： ， 因此可令 ，则也可以得出迭代公式： 。 4 迭代法的基本思想是将方程改写成等价的迭代形式，但随之而来的问题却是迭代公式不一定收敛，或者收敛的速度较慢。运用前述加速技巧，对于简单迭代过程 ，其加速公式具有形式： ，其中 记，上面两式可以合并写成： 这种迭代公式称作简单的牛顿公式，其相应的迭代函数是： 。 需要注意的是，由于是的估计值，若取，则实际上便是的估计值。假设，则可以用代替上式中的， 就可得到牛顿法的迭代公式： 。 )(')(1k k k k x f x f x x - =+)()()(x f x k x x x -==?)0)((≠x k )(')()()('1)('x f x k x f x k x --=?) ('x ?α0)('=α?≠)('αf α=)(x f 0)('=α?)('1 )(ααf k = )('1 )(x f x k = )(')(1k k k k x f x f x x - =+0)(=x f )(x x ?=)(1n n n x f x x +=+θθ?--= +1)(1n n n x x x ) (111n n n x x x --+=++θθ )(1 n n x x ?=+1-=θL L x f x x n n n )(1- =+L x f x x )()(- =?L )('x ?)()(x f x x +=?)('x ?)('x f 0)('≠x f )('x f L )(')(1n n n n x f x f x x - =+
泊松方程和拉普拉斯方程
拉普拉斯方程和泊松方程 摘要：拉普拉斯方程，又名调和方程、位势方程，是一种偏微分方程。因为由法国数学家拉普拉斯首先提出而得名。求解拉普拉斯方程是电磁学、天文学和流体力学等领域经常遇到的一类重要的数学问题，因为这种方程以势函数的形式描写了电场、引力场和流场等物理对象的性质。 关键词：分离变量 电磁场 拉普拉斯 简史 1777年，拉格朗日研究万有引力作用下的物体运动时指出：在引力体系中，每一质点的质量m k 除以它们到任意观察点P 的距离r k ，并且把这些商加在一起，其总和 m k r k n k=1 = V x ，y ，z 即P 点的势函数，势函数对空间坐标的偏导数正比于在 P 点的质点所受总引力的相应分力。1782年，P.S.M.拉普拉斯证明：引力场的势函数满足偏微分方程： ?2V ?x +?2V ?y +?2V ?z =0，叫做势方程，后来通称拉普拉斯方程。1813年，S.D.泊松撰文指出， 如果观察点P 在充满引力物质的区域内部，则拉普拉斯方程应修改为?2V ?x 2 + ?2V ?y 2 + ?2V ?z 2 =?4πρ， 叫做泊松方程，式中ρ为引力物质的密度。文中要求重视势函数 V 在电学理论中的应用，并指出导体表面为等热面。 静电场的泊松方程和拉普拉斯方程 若空间分区充满各向同性、线性、均匀的媒质，则从静电场强与电势梯度的关系E=-?V 高斯定理微分式??E =ρ/εr ε0，即可导出静电场的泊松方程：?2V ?x 2+?2V ?y 2+?2V ?z 2=?2V =?ρ/εr ε0 式中ρ为自由电荷密度，纯数 εr 为各分区媒质的相对介电常数，真空介电常数εo =8.854×10-12 法／米。在没有自由电荷的区域里，ρ＝0，泊松方程就简化为拉普拉斯方程?2V =0 。 在各分区的公共界面上，V 满足边值关系V i =V j ， ε0εri ?V ?n i ?ε0εrj ?V ?n j =??,
用matlab实现线性常系数差分方程的求解
数字信号处理课程设计 题目：试实现线性常系数差分方程的求解 学院： 专业： 班级： 学号： 组员： 指导教师：
题目：用Matlab 实现线性常系数差分方程求解 一． 设计要求 1． 掌握线性常系数差分方程的求解 2． 熟练掌握Matlab 基本操作和各类函数调用 3． 结合Matlab 实现线性常系数差分方程的求解 二．设计原理 1．差分与差分方程 与连续时间信号的微分及积分运算相对应，离散时间信号有差分及序列求和运算。设有序列f(k),则称…,f(k+2),f(k+1),…,f(k －1),f(k －2),…为f(k)的移位序列。序列的差分可以分为前向差分和后向差分。一阶前向差分定义为 ()(1)()f k f k f k ?=+- （3.1—1） 一阶后向差分定义为 ()()(1)f k f k f k ?=-- （3.1—2） 式中Δ和Δ称为差分算子。由式（3.1—1）和式（3.1—2）可见，前向差分与后向差分的关系为 ()(1)f k f k ?=?- （3.1—3） 二者仅移位不同，没有原则上的差别，因而它们的性质也相同。此处主要采用后向差分，并简称其为差分。 由查分的定义，若有序列1()f k 、2()f k 和常数1a ，2a 则 1122112211221112221122[()()][()()][(1)(1)][()(1)][()(1)]()() a f k a f k a f k a f k a f k a f k a f k f k a f k f k a f k a f k ?+=+--+-=--+--=?+? （3.1—4） 这表明差分运算具有线性性质。 二阶差分可定义为 2()[()][()(1)]()(1) ()2(1)(2) f k f k f k f k f k f k f k f k f k ?=??=?--=?-?-=--+- （3.1—5） 类似的，可定义三阶、四阶、…、n 阶差分。一般地，n 阶差分
差分方程方法
第四章 差分方程方法 在实际中，许多问题所研究的变量都是离散的形式，所建立的数学模型也是离散的，譬如，像政治、经济和社会等领域中的实际问题。有些时候，即使所建立的数学模型是连续形式，例如像常见的微分方程模型、积分方程模型等等，但是，往往都需要用计算机求数值解。这就需要将连续变量在一定条件下进行离散化，从而将连续型模型转化为离散型模型，因此，最后都归结为求解离散形式的差分方程解的问题。关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。 下面就不同类型的差分方程进行讨论。所谓的差分方程是指：对于一个数列{}n x ，把数列中的前1+n 项()n i x i ,2,1,0=关联起来所得到的方程。 4．1常系数线性差分方程 4.1.1 常系数线性齐次差分方程 常系数线性齐次差分方程的一般形式为 02211=+?+++---k n k n n n x a x a x a x (4.1) 其中k 为差分方程的阶数，()k i a i ,,2,1 =为差分方程的系数，且()n k a k ≤≠0。对应的代数方程 02 211=++++--k k k k a a a λ λλ （4.2） 称为差分方程的（4.1）的特征方程，其特征方程的根称为特征根。 常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。下面分别就特征根为单根、重根和复根的情况给出差分方程解的形式。 1. 特征根为单根 设差分方程（4.1）有k 个单特征根 k λλλλ,,,,321 ，则差分方程（4.1）的通解为 n k k n n n c c c x λλλ+++= 2211， 其中k c c c ,,,21 为任意常数，且当给定初始条件 () 0 i i x λ= ()k i ,,2,1 = (4.3) 时，可以唯一确定一个特解。 2. 特征根为重根 设差分方程（4.1）有l 个相异的特征根()k l l ≤≤1,,,,321λλλλ 重数分别为 l m m m ,,,21 且k m l i i =∑=1 则差分方程（4.1）的通解为
拉普拉斯方程
拉普拉斯方程 一、概念：一个弯曲的表面称为曲面，通常用相应的两个曲率半径来描述曲面，即在曲面上某点作垂直于表面的直线，再通过此线作一平面，此平面与曲面的截线为曲线，在该点与曲线相切的圆半径称为该曲线的曲率半径R1。通过表面垂线并垂直于第一个平面再作第二个平面并与曲面相交，可得到第二条截线和它的曲率半径R2，用R1与R2可表示出液体表面的弯曲情况。若液面是弯曲的，液体内部的压强p1与液体外的压强p2就会不同，在液面两边就会产生压强差△P= P1- P2，称附加压强，其数值与液面曲率大小有关，可表示为：，式中γ是液体表面张力系数，该公式称为拉普拉斯方程。 二、在数理方程中 拉普拉斯方程为：，其中?2为拉普拉斯算子，此处的拉普拉斯方程为二阶偏微分方程。三维情况下，拉普拉斯方程可由下面的形式描述，问题归结为求解对实自变量x、y、z二阶可微的实函数φ： 其中?2称为拉普拉斯算子。 拉普拉斯方程的解称为调和函数。 如果等号右边是一个给定的函数f（x，y，z），即： 则该方程称为泊松方程。拉普拉斯方程和泊松方程是最简单的椭圆型偏微分方程。偏微分算子（可以在任意维空间中定义这样的算子）称为拉普拉斯算子，英文是Laplace operator或简称作Laplacian。
三、方程的解 称为调和函数，此函数在方程成立的区域内是解析的。任意两个函数，如果它们都满足拉普拉斯方程（或任意线性微分方程），这两个函数之和（或任意形式的线性组合）同样满足前述方程。这种非常有用的性质称为叠加原理。可以根据该原理将复杂问题的已知简单特解组合起来，构造适用面更广的通解。 四、二维方程 两个自变量的拉普拉斯方程具有以下形式： Δu =δ2u/δu2+δ2u/δy2=0 解析函数的实部和虚部均满足拉普拉斯方程
chenpc_文件下载_数理方法_实验四、拉普拉斯方程与泊松方程的求解
实验四 拉普拉斯方程与泊松方程的求解 一、拉普拉斯方程的求解 例题：求解定界问题： ()()()()()00,030,0,,sin 3,00,,sin cos xx yy u u x a y b y u y u a y b x x u x u x b a a πμππμ??+=≤≤≤≤????==? ?????????==? ? ?????? 任意选取定界问题中参数的值，例如取1,1,1a b μ===。用偏微分方程工具箱来求解的步骤如下。 1、画求解区域 在指令窗口中，输入pdetool ，打开偏微分方程工具箱的界面， 图1 微分方程工具箱的界面 选择菜单Options/Axes Limits ，打开对话框如图2所示。 图2 设置坐标变化范围的对话框
在X-axis range 和Y-axis range 栏中都输入[-0.1 1.1]，单击按钮Apply 确认，再关闭对话框。 单击左上角画矩形框按钮，在pdetool 的窗口中画一个矩形，然后，在刚画出的灰色矩形区域内部双击鼠标左键，出现如图3所示的对话框，设置左边界（Left ）参数为0，下边界（Bottom ）参数为0，宽度（Width ）参数为1，高度（Hight ）参数为1，点击OK 按钮，画出一个边长为1的正方形区域01,01x y ≤≤≤≤，这个正方形被自动命名为R1，并显示在区域上方的公告栏（Set Formula ）中。 图3 确定正方形区域的边界位置和名称的对话框 2、设定方程类型 单击按钮，打开如图4所示的对话框。 图4 设置方程类型的对话框 在方程类型中选择椭圆型，这时方程的形式为 ()c u au f -???+= 取1,0,0c a f ===，设置好参数后，单击OK 即可。 3、设定边界条件 单击按钮，进入边界模式。这时区域由灰色变成白色，而边界变成红色。选择菜单Boundary/Show Edge Labels ，给四条边界标上序号1,2,3,4。根据题意，双击边界1，打