色散方程的六点半显式差分格式.pdfVIP

簟 ¨誊 卑2辩 华 侨 大 学 学 报 自然 科 学 版 VoI．14 ．2 JOURNALOFHUAQIAO UNIVERSITY 1●’3年 ‘月 (NATURAL SCIENCE) Apr．1993 色散方程的六点半显式差分格式 王 子 丁l 一 一 — — — 一 — 《管理信息科学系1 铺要 奉文构造色散方程_l=4-～的一麦三层六点的差分格式．其截断误差为0(一十 +rth’． 格式是无条件穗定的．且用的同格点觳步 ．精度高 ．可以用显式进行计算．文束用效值例子说明r 格式耐定■问曩的应用． 关键调 !截·断最差一无条件穗定，墨 格 妻泛 _I_ 0 引言 在数值计算中，对于差分格式面言，稳定性好的显式格式一般耍比隐式格式经挤 {，时 于色散方程 一口#一 (n为常数 ，可正可负)．在文C2)中提出一类三层无条件稳定的半显式差 分格式 A ，其截断误差为O(r+^+ ／^’)，本文是对格式 ^．的推广 ．提出一类三层的无条件 稳定的半显式格式．此格式只用到六个网格点，其截断误差可达O(+ +一／^’)．它在计算的 一 层网格点仍不是牵涉劐相郎的网格点．而是牵涉到相隔一个节点的两个网格点．佳!与 (或 )不同的是第t堪只用到两个网格点．在计算过程 中，采用半显式计算 ，即先从左向右显式 地计算偶效节点的网格函数值，然后再显式地计算奇数节点的网格函数值．从而格式是无条件 稳定，又大大减少工作量．文朱用数值倒子说明了本格式对定解问题的应啊． 为了分析差分格式的稳定性，先介绍如下 Miller准则 “． 当l l—lCl时t复系数二次方程 A + B^+ C = O，(A ≠ O) 具有模为 l的不等复根 ，其充要条件为；B=盈_c，且 IBI2IAJ． 1 差分格式的构造 设 f为时问步长，^为空问步长，同域 由求解区域的点集( 、￡，)构成 (打t、n为整数)，在节 本文 lggl-12—28收到． 第 2期 色t方程的六点半曼式差分格式 149 点(如 、“)处甩 ．I=表示傲分方程的解。即嵋一_(如 ，t．)，甩 表示差分方程的解． 将下科的效值徽分公式 ( )：一 (1／2)C( )二-L-I-()：+L)-I-o(k)‘， (1) (_。)：一I一 (1／2f) ： 一 。 )+ 口(一)， (2) (嘶)：+I一 (1／2)『C 一_ )+ o(rt)。 (3) (搿一)：皇 (1／2^’)(q：+I一 2H：+I+ 2_：．I一誓：一-)-I-o(k’)， (4) “)：¨一 (112)(。： + _： )-I-o(rt)， (5) (-)：一l一 (t／2) 盘 -I-_： )+ 口(一)， (6) 代入方程 (_I)：一d(_一 )： 中，舍去余习lO(一+^l+fI，^I)，便得封以下l的差分格式 (1／t‘)(喘 一 口： -I- 一_： ) 一 (a／2k’)C口二+l_((口： -I-蠕 )+ ( -I- )一 一1)， (7) 令 r一一r，^’，I哼得到奉文的差分格式