偏微分方程—matlab

《偏微分方程—matlab》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、基础知识偏微分方程的定解问题各种物理性质的定常（即不随时间变化）过程，都可用椭圆型方程来描述。其最典型、最简单的形式是泊松(Poisson)方程（1）特别地，当f(x,y)≡0时，即为拉普拉斯(Laplace)方程，又称为调和方程（2）带有稳定热源或内部无热源的稳定温度场的温度分布，不可压缩流体的稳定无旋流动及静电场的电势等均满足这类方程。Poisson方程的第一边值问题为（3）其中Ω为以Γ为边界的有界区域，Γ为分段光滑曲线，ΩUΓ称为定解区域，f(x,y),ϕ(x,y)分别为Ω,Γ上的已知连续函数

2、。第二类和第三类边界条件可统一表示成（4）其中n为边界Γ的外法线方向。当α=0时为第二类边界条件，α≠0时为第三类边界条件。在研究热传导过程，气体扩散现象及电磁场的传播等随时间变化的非定常物理问题时，常常会遇到抛物型方程。其最简单的形式为一维热传导方程（5）方程（5）可以有两种不同类型的定解问题：初值问题（也称为Cauchy问题）（6）初边值问题（7）其中ϕ为已知函数，且满足连接条件问题（7）中的边界条件称为第一类界条件。第二类和第三类边界条件为（8）其中。当时，为第二类边界条件，否则称为第三类边界

3、条件。双曲型方程的最简单形式为一阶双曲型方程（9）物理中常见的一维振动与波动问题可用二阶波动方程（10）描述，它是双曲型方程的典型形式。方程（10）的初值问题为（11）边界条件一般也有三类，最简单的初边值问题为如果偏微分方程定解问题的解存在，唯一且连续依赖于定解数据（即出现在方程和定解条件中的已知函数），则此定解问题是适定的。可以证明，上面所举各种定解问题都是适定的。§2偏微分方程的差分解法差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对

4、求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：（i）选取网格；（ii）对微分方程及定解条件选择差分近似，列出差分格式；（ii

5、i）求解差分格式；（iv）讨论差分格式解对于微分方程解的收敛性及误差估计。下面我们只对偏微分方程的差分解法作一简要的介绍。2.1椭圆型方程第一边值问题的差分解法以Poisson方程（1）为基本模型讨论第一边值问题的差分方法。考虑Poisson方程的第一边值问题（3）取h,τ分别为x方向和y方向的步长，以两族平行线将定解区域剖分成矩形网格。节点的全体记为为整数}。定解区域内部的节点称为内点，记内点集为。边界Γ与网格线的交点称为边界点，边界点全体记为Γhτ。与节点沿x方向或y方向只差一个步长的点和称为节

6、点的相邻节点。如果一个内点的四个相邻节点均属于ΩUΓ，称为正则内点，正则内点的全体记为Ω(1)，至少有一个相邻节点不属于ΩUΓ的内点称为非正则内点，非正则内点的全体记为Ω(2)。我们的问题是要求出问题（3）在全体内点上的数值解。为简便记，记。对正则内点，由二阶中心差商公式Poisson方程（1）在点处可表示为（12）在式（12）中略去，即得与方程（1）相近似的差分方程（13）式（13）中方程的个数等于正则内点的个数，而未知数,则除了包含正则内点处解的近似值，还包含一些非正则内点处的近似值，因而方程个

7、数少于未知数个数。在非正则内点处Poisson方程的差分近似不能按式（13）给出，需要利用边界条件得到。边界条件的处理可以有各种方案，下面介绍较简单的两种。(i)直接转移(ii)线性插值由式（13）所给出的差分格式称为五点菱形格式，实际计算时经常取h=τ，此时五点菱形格式可化为（14）简记为¸（15）其中¸。求解差分方程组最常用的方法是同步迭代法，同步迭代法是最简单的迭代方式。除边界节点外，区域内节点的初始值是任意取定的。例1用五点菱形格式求解Laplace方程第一边值问题其中。取。当时，利用点(k

8、,j),(k±1,j.1),(k±1,j+1)构造的差分格式（16）称为五点矩形格式，简记为£（17）其中£。2.2抛物型方程的差分解法以一维热传导方程（5）为基本模型讨论适用于抛物型方程定解问题的几种差分格式。首先对xt平面进行网格剖分。分别取h,τ为x方向与t方向的步长，用两族平行直线(k=0,±1,±2,…)，k(j=0,1,2,…)，将xt平面剖分成矩形网格，节点为(k=0,±1,±2,…,j=0,1,2,…)。为简便起见，记。2.2.1微分方程的差分近似在网