计算方法7课件.pptVIP

* 电磁场中的计算方法 第七章 偏微分方程数值解法 7.1 差分方法的基本概念 1.补充知识 “高数”中接触了一些简单偏微分，也接触了简单偏微分方程，如： (1). 其中 (2). 满足： ( 3). 2sin(x+2y-3z)=x+2y-3z 满足： (4). 满足： (5). 满足： (6). 满足： 上面是已知函数: ，验证满足等式，反过来，将等式视为方程，则是求解方程，得到解函数。 因此偏微分方程： (1). 含偏微分的等式， (2). 求解偏微分方程、求含多个自变量的函数 (3). 带有初值、边界条件。 常微分方程的求解已很困难，通过分门别类研究，能求得一些特殊类型方程的解（只含一个变量），即便是一阶方程，也很难求出解析解表达式，也因此，在上一章我们研究了一阶微分方程的数值解法。 要求解偏微分方程比求解常微分方程更难，因此寻求偏微分方程的数值解更显重要，实际上，绝大部分偏微分方程不可能求到解析函数解，基本上都是数值解法。 2.偏微分方程定解问题 一般来说，偏微分方程从实际问题抽出后，多是下列几种类型: （1）泊松方程（Poisson），又称为椭圆型方程： ?：自变量的变化区域，有界区域。 ?：? 的边界，分段光滑曲线。 当 称为拉普拉斯方程（Laplace）或调和方程， 例如 满足： 相应第一边值条件： 第二、第三边值条件： 为边界?的外法线方向， 为第二边界条件， 为第三边界条件。 各种物理性质的定长问题（不随时间变化过程），都可用椭圆型方程描述。如带有稳定热源或内部无热源的稳定场的温度分布，不可压缩流体的稳定克旋流动及静电场的电热等均满足上述方程。 （2）热传导方程（抛物型） 相应有：柯西（Cauchy）初值条件： 初边值条件为： 第一边值条件： 第二 边值 条件： 第三边值条件为： 在热传导过程的研究中，气体的扩散现象及电磁场的传播等随时间变化的非定常物理问题，都可用上述方程来描述。 （3） 波动方程（双曲型） 最简单形式为线性双曲方程： 其初边值 条件为： 边值条件同热传导方程。 物理中常见的一维振动及各类波动问题，均可用波动方程描述。 3 有限差分法的基本概念 被称为函数 的差分(一阶差分)。 差分与微分是不同的，因为差分是有限量的差，故通常也被称为有限差分。但是，只要增量 很小，差分 与微分 之间的差异将很小。 根据差分的定义，在差分运算中还经常用到一阶中心差分： 由一阶差分的差分，得到 ，称为原函数 二阶差分，同样，当 很小时，二阶差分 很接近于二阶微分 。 设有一函数 ，当其独立变量 有一很小的增量 时，相应地函数 地增量为： 一阶差分 除以增量 的商称为一阶差商，即： 当增量 很小时，它很接近于一阶导数 , 即： 上述三式分别称为一阶向前、向后和中心差商，它们的截断误差可由泰勒公式的展开式得知： 由一阶中心差分得到二阶差分为： 所以，一阶向前差商公式的截断误差为 。 所以，一阶向后差商公式的截断误差为 。 所以，一阶中心差商公式的截断误差为： 。 同理，由泰勒公式可知： 同样，当增量 很小时，二阶差商很接近于二阶导数 ，即： 由此可知，上述三种差商表达式中，以中心差商的截断误差最小。 所以，二阶中心差商公式的截断误差为： 。 同理，偏导数也可以近似地用相应的差商来表达，若设定函数为为 ，当其独立变量 得到一个很小的增量 时，则 方向的一阶偏导数可以近似表达为： 同样，相应的二阶偏导数也可近似地表达为： 有限差分法是用来求解由偏微分方程