抛物型方程周期问题的高精度精细积分法.pdfVIP

内蒙古工业大学学报 OFINNERMONGOLlA JOURNAL 第23卷第4期 UNIVERSITYOFTECHNOLoGY V01．23No．42004 文章编号：1001—5167(2004)04一0245一05 抛物型方程周期问题的高精度精细积分法+ 王玉兰1’2，杜 红2～，王玉荣4 (1．内蒙古工业大学理学院．呼和浩特010062；2．哈尔滨工业大学数学系，哈尔滨150001 3．黑龙江科技学院数力系，哈尔滨150027；4．呼和浩特第九中学．呼和浩特010038) 摘要：对于抛物型方程(詈一暴)2摊一窘一2焘+翥一。的初始值与周期 边值问题，提出了局部截断误差阶为O(△≯)的精细时程积分法．由于这种方法是 半解析方法，在时间域上可以精确计算．所以本方法不仅精度高，而且还绝对稳定． 文末的数值算例进一步表明．精细积分法对大的时间步长和长时间计算均有效，因 此是一种很实用的方法． 关键词：精细积分法；高阶抛物型方程；无条件稳定；局部截断误差；长时间计算 中图分类号：()241．8文献标识码：A o 引 言 求解偏微分方程有有限元法、有限差分法等，当计算域较大时，如果用有限元法求解，因空间坐标离 散后未知数多，势必导致计算机耗时量大；而用有限差分法求解，尽管能减少计算机存储量，但差分法固 有缺陷会限制其应用范围．因为差分法是从整体上考虑用一组差分方程来代替微分方程，对于复杂的结 构形式及边界条件难以处理．钟万勰院士提出的一种求解偏微分方程的新数值方法埔§细积分法n’21，该 方法属于半解析法，具有数值精度高及带宽小的优点，本文利用四阶中心差分格式，将抛物型方程 【詈～砉)。“墨擎一2磊叁+器一。化为常微分方程组，采用了时程精细积分法，放弃了传统的差分格 式，对有限的时间步长又进一步作了2Ⅳ个(Ⅳ为正整数)细化，使计算精度得到充分提高．目前，精细积 分法在许多领域得到广泛应用，多数是化为关于时间变量f的一阶偏导数，但对于关于时间变量f的二 ，1 吨、2 12．． 阶偏导数的高阶抛物型方程用精细积分法求解的文献尚未见到．本文求解了方程卜兰一笔l雒兰笔～2 呻．． 一。． 嘉丧+等=o的周期问题．数值计算表明，精细积分法对大的时间步长和长时间计算均有效，因此是一 种很实用的方法． 。 1 “。。，“。一．的四阶差分构造 记“(z。。￡)=“。，其中z。=棚^，将‰+1，“。一l在‰作Taylor级数展开，有i ·收稿日期：2004一lo一24 基金项目：国家自然科学基金、哈尔滨工业大学(威海)校基金(2002～15，14)资助 作者简介：王玉兰．女。讲师． 万方数据 246 内蒙古工业大学学报 zD04年 (2) “。一。一甜。一^(“，)。+筹(“。)。一筹(“一)。+筹(比一)。一等(“一)。+o(^6) (1)+(2)得 ’ (3) “。+。+“。一。=2‰+2筹(距。)。+2筹(M一。)。+o(危6) 同理可得 (4)