双曲型方程的差分法2.doc

1.3一阶线性常系数双曲型方程组 例：对波动方程的初值问题， 引进，， ， 其中 ，，，, 是线性双曲型方程组 。 对此模型问题，来构造差分格式。为方便起见，设 用不同的差商来代替，可得下列格式 显格式： （对t的偏导数都用向前差商，w对x在j点的偏导数用中心差商，v对x在点的偏导数用中心差商） 对初始条件：，， 由，，化为，，于是得：， 。 由初始条件，由第一个方程算出第1层的，再由第2个方程算出。依次进行，已知k层的和， 第一步： 第二步： (已知黑点，算出红点) 稳定性条件： 。 隐格式： （对t的偏导数都用向前差商，w对x在j点的偏导数用k 层和k+1层的中心差商的算术平均，v对x在点的偏导数用k 层和k+1层的中心差商的算术平均） (已知黑点，算出红点) 第一步： 第二步： 两步均为隐式。为减少工作量，可用第二个方程消去第1个方程中的， ，得： 是三对角方程组，先解出，再由第2个方程求，此时第2个方程成了显式，从而大大减少了工作量。 由Fourier方法可证明隐格式无条件稳定。 3、Lax-Friedrichs格式： 是显格式。全取整节点。截断误差为。当 时格式稳定。 4、Lax-Wendroff格式： 略去高阶项，设，得： 同理有： 是显格式。截断误差阶为。当时格式稳定。 作业：习题五 4，5。 对于一般的n个函数的一阶线性方程组的初值问题， ， 也可以建立一般的以向量和矩阵形式表示的差分方程。例如书上有Lax- Friedrichs格式表示为 当时就是上面的Lax- Friedrichs格式。记 上式也可改写为 同样，对于Lax-Wendroff格式，用向量和矩阵可表示为 其中当时就是上面的格式。 这两种格式的稳定性条件均为 另外书上还有一个与前面的一个方程时引进的跳蛙格式，以向量和矩阵形式可表示为 稳定性条件为 对于变系数双曲型方程的初值问题，即当矩阵A是x的函数时，也可建立相应的差分格式。如在书上P.136给出了如下的Lax-Wendroff格式： 一般的推导较繁，但对波动方程的初值问题，引进，，再按照上面的泰勒展开的推导思想易得上述公式。它的稳定性条件是 而对于一般的一阶方程组的初边值问题，要根据A的特征值的正负号给出相应的边界条件，才能使问题适定。A的若干特征值为正时，对应的方程的特征线斜率为正，初值向右无限传播，只能在x的变化区域的左边界上给出边界条件；反之，A的特征值若干为负时，对应的方程的特征线斜率为负，初值向左无限传播，只能在x的变化区域的右边界上给出边界条件。这里不进行深入的讨论。 §2 二阶线性双曲型方程的差分格式 2.1一维波动方程 1.二阶线性双曲型方程的最简模型：一维波动方程。 初值问题： ， ， 其中设 如上节所示，可化成一阶方程组来求解，已构造了一些格式。下面介绍直接对二阶方程的差分化方法。 2．显格式 设空间步长h，时间步长，作两族平行线。 用中心差商代替二阶偏导数，得 截断误差为。初始条件离散为 对后一个边值条件，也可用中心差商，误差为。但要消去。可令前面的差分格式中用代入，解出，代入边界条件的差分格式，可得： ， 其中为网比。于是差分格式也可化成 * K+1层 * * * K 层 * K-1层 由初始条件，得和。可算出…。计算层的一个点的值，用到k层上3个点和k-1层上一个点的值，是显格式。上述显格式也可用于解混合问题。 稳定性分析 上述显格式稳定的充分必要条件为（是稳定的必要条件，即稳定。反之不真） 的几何解释： (j,k) * K 层 * * * K-1 层 * * * * * K-2 层 * * * K-3 层 由差分格式，知依赖于前2层，， 和共4个值。这4个值又依赖于k-2层上5个值， k-3层上3个值，k-4层上一个值，等等。依次往下推，最终依赖于初始层层上的2k+1个值：，，…，，。 在中的网点为差分解的依存域。