偏微分方程数值解法2014.11.05.ppt

偏微分方程数值解法（第四讲） 2014年11月05日星期三 1.一阶线性常系数双曲型方程组 考虑一阶线性常系数方程组 其中 为常系 数矩阵。 定义： 如果A的特征值是实的，并存在非奇异矩阵S使得 其中 为A的特征值，称（1）为双曲型方程组，如果A是对称阵，则称（1）为对称双曲型方程组，如果A的特征值是实的且互不相同，则称（1）为严格双曲型方程组。 给定初始条件 与（1）构成初值问题。 4.3 一阶双曲型方程的差分格式 常用求解该线性方程组的格式有 a. 后退的Euler方法 此格式不能用。 b.Lax—Friedrichs 格式 c. 蛙跳格式 4.3.1 Lax-Friedrichs格式 2. 二阶线性双曲型方程 为推导简单起见，令 令 因此我们有 4.3.1 Lax-Friedrichs格式 这里 A的特征值为-1,1，因此这个方程组为双曲型方程组。 边界条件 所以 所以 初始条件 4.3.1 Lax-Friedrichs格式 也可直接求解： 截断误差为 稳定性条件为 3 二维问题 求解的格式有： a. 隐格式 该截断误差为 令 代入格式后可得， 格式是无条件稳定的。 b. Crank-Nicolson格式 该截断误差为 格式为无条件稳定。 用隐格式求解二维问题得到的线性方程组其系数矩阵为宽带状。因此求解 不是特别便利。 c. Lax-Friedrichs 格式 显格式，条件稳定 令 注意到 如果 ， 即 则格式稳定，若b=a,稳定性条件化为 d.局部一维格式 在x方向的差分把 推进到 在y方向上的差分把 推进到 其中 是关于x方向、y方向的差分算子，这样的二步法称为分数 步长法，也称为局部一维格式 格式1： 可改写为 这里， 为分析截断误差，消去 上述格式等价于 可得截断误差为 对应第一式的增长因子是 对应第二式的增长因子是 整体格式的增长因子为 格式为无条件稳定。 格式2：Lax-Wendroff 这里，x方向定义的算子有 y方向定义的算子有 格式2：Lax-Wendroff 将第一式代入第二式可得等价式子： 得到截断误差为 第一式当 时有 第二式当 时有 由 可得上述格式的稳定性条件为 e. Beam-Warming格式 C-N 格式 将C-N格式等价地写成 或写成 上式也可等价为 上式中最后一项是 的项，去掉后可得一个二阶格式 由上可得ADI格式为 增长因子， 该格式为二阶精度，无条件稳定。 另一个ADI格式： 在 式两边减去