二维波动方程数值解一种离散方法.pdfVIP

维普资讯 南京大学学报数学半年刊 第 13卷 第 1期 JOURNALOFNANJINGUNIVERSITY Vo1．13，No．1 1996年 5月 MATHEMATICALBIQUARTERLY May，1996 』。7厂j? 二维波动方程数值解的一种离散方法 02~／, 蔡 日增 (江苏石油化工 常州 213016) (无儡 214036) A 摘■ 本文蛤出了二维诚动方程的一种离散方法，文中蛤出了离散逼近的收敛性．鼓值计算结 果裹明这种方法收敛速度是相当快的． 二 分囊号 ．82 枨幽 咝 岛I 改a 9241 1． 引 言 波动方程由于其在地璋物理 振动力学等学科中的地位始终是一个引人注 目的问题 ． 由于一维波动方程的求解方法方便 有效，因此不少文章将二维波动方程化为一维波动方程 的耦台方程来求解叫[|]．离散方法是处理高维问题的常用方法啪_[|．本文先将一十变量离 散(半高散)，将二维波动方程化为一维波动方程组．再化为积分方程高教求解 ． 为简单起见，考虑如下半无界区域上二维波动方程定解问题：求 ，，)满足 一 (耄+雾)=o． 0，l0．01，(1l1) (，，0) 毋( ．)． 堡丛： 三 一 (，j，)． z≥0， 0≤ j，≤ 1 (1．2) u(O．，)；，1(，)， ≥O，0≤ ≤1． (1．3) · u(x，0，)； hl( ，)， #(，1，)= hl(，f)， ≥ 9， ≥ o n．4) 其中a ，)≥c0 t0表示波速 ． 本文先将变量 离散，将定解问题化为一维披动方程组，用求解一维被动方程的方法求 解波动方程组 ．在第二部分将给出求解半离散问题的选代方法及数值模拟结果，第三部分 给出收敛性结果 ． · 收翦 日期 199s一09—21． 奉{寰囊 由江苏省教妥优秀青年教师基金，￡助 维普资讯 ·110 · 南京大学掌掇数学半年刊 1995年5月 2． 求解方法 在区间[0，1]上离散变量 给定结点数 Ⅳ，步长^= ， ，； tj=o，l一。，N-记 ， )=(，yJ，})，i(，幻=嘉[“，)+“一，。，})一2uif·)]· 则定解问题 (1，1)-- (i．4)的半离散问题为求 ，J‘， j=1,2’．．’t～一1，满足 挚 ()[≥+V )]一。， o (2．1) (，O)= J(工)， ，f ，O)； ，()， ≥ 0 (2．2) 嘶(0．)= (f)， I≥0 (2．3) 口。妇 ，f)= ^1 ，})， “ ，})= ^2 ，)， ≥ 0，}≥ 0 (2，4) 这里a／x)~a(x，)， )= ，yJ)，％() 仇(，)， () ^(J，})以下设 =1， 2．… N 一 1． 作变换 = Gjc— (2．5) 由于aj(x)O，嚣而逆变换x=Gf