大连理工大学高等数值分析偏微分方程数值解(双曲方程书稿).docx

双曲型方程的有限差分法 线性双曲型方程定解问题： （a）一阶线性双曲型方程 ?u  ? a?x??u ? 0 ?t ?x （b）一阶常系数线性双曲型方程组 ?u ? A ?u ? 0 ?t ?x 其中 A ， s 阶常数方程方阵， u 为未知向量函数。 二阶线性双曲型方程（波动方程） ? ??2u ? ? ? ? ? a?x??u ? ? 0 a?x?为非负函数 ?t 2 ?x ? ?x ? 二维，三维空间变量的波动方程 ???2u ? ?2u ?2u ? ? ? ?t 2 ? ? ?x2 ? ? 0 ?y2 ? ??2u ? ?2u ?2u ?2u ? ? ?t 2 ? ? ?x2 ? ?y2 ? ?z2 ? ? 0 ? §1 波动方程的差分逼近 波动方程及其特征 线性双曲型偏微方程的最简单模型是一维波动方程： （1.1） 其中a ? 0 是常数。 ?2u ?t 2  ? a2 ?2u ?x2 （1.1）可表示为： ?2u ?t 2  a2 ?2u ?x2 ? 0 ，进一步有 ? ? ? a ? ? ? ? ? ? a ? ?u ? 0 ? ?t ?x ? ? ?t ?x ? ? ? ? ? 由 于 ? ? a ? 当 dx ? ?a 时 为 u?x, t ? 的 全 导 数 ?t ?x dt （ du ? ?u ? ?u ? dx ? ?u ? a ?u ）,故由此定出两个方向 dt （1.3） ?t ?x dt ?t ?x dt ? ? 1 dx a 解常微分方程（1.3）得到两族直线 （1．4） x ? a ? t ? C 1 和 x ? a ? t ? C 2 称其为特征。 特征在研究波动方程的各种定解问题时，起着非常重要的作用。比如，我们可通过特征给出（1.1）的通解。（行波法、特征线法） 将（1.4）视为(x, t) 与(C , C 1 2 ) 之间的变量替换。由复合函数的微分 法则 ?u ?u ?C ?u ?C ?u ?u 12?x ? ?C ? ? ? ?C ? ? ? ?C ? ?C 1 2 x x 1 2 1 2 ????2u ? ? ?u ?u ? ?C ? ? ?u ?u ? ?C ? ? ? ? ? ? ? 1 ? ? ? ? 2 ??x2 ? ?C ?C 1 1 ?C ?x ?C ?C 2 2 1 ?C ?x 2 ? ?2u ? ?2u ? ?2u ? ?2u ?C 2 1 ?C ?C 1 2 ?C ?C 2 1 ?C 2 2 ? ?2u ? 2 ?2u ? ?2u ?C 2 ?C ?C ?C 2 同理可得 1 1 2 2 12?C ?t ， ?C 1 2 ?t ? ?a ?t ? ?a ?t2 ? a ?u ?u ?C ?u ?C ? ?u ?u ? 1?t ? 1 ?C ? ?t 1 ? ?C ? ?t 2 ? a? ?C ?1 ? ? ?C ? ?2 ? ?2u ? ? ? ?u ?u ?? ?C ?u ? ? ?u ?u ?? ?C 12? ? a ? ? 1 2 ? ?? 1 ? ? a ? ? ? ?? 2 1221?t 2 1 2 2 1 ?C ? ? ?C ?C ?? ?t ?C ? ? ?C ?C ?? ?t ? ?2u ?2u ? ? ?2u ?2u ? ? ?a2 ? ? ? ? a2 ? ? ? ??C ?C ?C 2 ? ??C 2 ?C ?C ? 2 1 2 2 1 2 ? ?2u ?2u ?2u ? ? a2 ? ? 2 ? ? ???C 2 ? ? 1 将 ?2u 和 ?2u 代入（1.1）可得： ?C ?C ?C 2 1 2 2 ?x2 ?t 2 ?? ?2u ?2u ?2u ? ? ?2u ?  ?2u ?2u ? a2 ? ? 2 ? ? ? a2 ? ? 2 ? ? ????C 2 ? ? ? 1 即有 ?C ?C ?C 2 1 2 2 ?C 2 1 ?2u ? 0 ?C ?C ?C 2 1 2 2 ?C ?C 1 2 求其对C 的积分得： ?u ? f ?C ? 其中 f ?C ?是C 的任意可微函数。 2 ?C 1 1 1 1 再求其对C 1 的积分得： （1.5） u?x, t ?? ? f ?C 1 ?dC 1 ? f ?C 1 1 ?? f ?C 2 2 ?? f 1 ?x ? at ?? f 2 ?x ? at ? 其中 f ???和 f ???均为任意的二次连续可微函数。 1 2 （1.5）为（1.1）的通解，即包含两个任意函数的解。 为了确定函数 f 1 ?x ? at ?和 f