偏微分方程数值解法教材.ppt

偏微分方程数值解法（第二讲） 2014年10月22日星期三 2. Fourier变换 Fourier积分公式，是指定义在 上的函数 的一个关系式。 设 有关系式 Fourier积分公式， 令 则有 称为 的Fourier变换， 称为 的逆变换。 对流方程的初值问题 有差分格式 上述差分格式中的解及其初值仅在网格点上有定义，为了应用Fourier方 法分析稳定性，必须扩充这些函数的定义域，使得它们在整个实轴上有定 义。 令 于是 对上式两边用Fourier积分来表示，可以得到 重新记为 这里 称为增长因子。 推广到一般形式的差分格式（限于常系数情形） 由于 不依赖于时间层n，可以得到 设 的任意次幂一致有界，界为K，则有 由Parseval等式 由定义可知 结 论：常系数差分格式稳定的充分必要条件是存在常数 使得当 时有 判别准则（Von Neumann条件）： 差分格式 稳定的必要条件是当 时，对任何 有 充分必要条件:正规矩阵，实对称矩阵，酉矩阵，Hermite矩阵 作 业 已知扩散方程 的一种差分格式如下， 求截断误差并推导其修正方程 。 第四章 双曲型方程的差分方法 1. 常系数双曲型方程数值解法 典型方程： a 对流方程（Advection equation） 第四章 双曲型方程的差分方法 b.二阶线性波动方程（Second order linear wave equations) 第四章 双曲型方程的差分方法 d.Burgers方程（Burgers equation) 非线性方程 第四章 双曲型方程的差分方法 2.数值方法 a.Lax-Friedrichs 格式 b.迎风格式 c.蛙跳格式 d.Lax-Wenddroff 格式 e.C-N格式 f.Beam-warming 格式 1） 相容性： 方程趋近 当差分方程中 ，时间与空间步长均趋近于0 时，差分方程的截断误差也趋近于0，则称差分方程与原微分方程是相容的。 2）收敛性： 解趋近（更强） * 当时间与空间步长均趋近于0 时，差分方程的解趋近于微分方程的解，则称差分方程的解收敛于原微分方程的解。 注意！ 方程互相趋近 解互相趋近 （根据Lax等价定理，只有稳定性条件满足的情况下，方程趋近才能保证解趋近） 分别为差分方程和微分方程的解 1. 复习：相容性、收敛性，稳定性 定义：称差分方程的初值问题是稳定的，如果当 做够小时，存在于 无关的常数C1和C2使得: 含义： 在差分方程的求解过程中，如果引入的误差随时间的增长有界，则称差分方程是稳定的。 * 4） Lax 等价定理 如果微分方程的初边问题是适定的，差分方程是相容的，则差分方程解的收敛性与稳定性是等价的。 含义： 如果微分方程不出问题（适定），差分方程性质好（稳定），则方程逼近就可保证解逼近。 如果方程逼近就可以导致解逼近，则差分方程的性质肯定好（稳定） 3）稳定性 3. Fourier分析法 4. 差分方程的修正方程 修正方程—— 差分方程准确逼近（无误差逼近）的方程 差分方程 截断误差 * 微分方程=差分方程+截断误差 差分方程 = 微分方程-截断误差 ≡ 新的微分方程（修正方程） 等价于 修正方程 * 通常要求： 修正方程中不出现时间的高价导数项 （便于进行空间分析） 修正方程： 主导项： 1阶； 耗散型 * c.一阶线性双曲型方程组（Linear first order hyperbolic system) 这里，u 和f都是向量，A的特征值是实的并且A是可对角化的，则上述方 程组称为双曲型方程组