四阶抛物型方程的样条子域精细积分配置法.pdfVIP

维普资讯 2007年 l0月 重庆师范大学学报(自然科学版) Oct．2007 第24卷 第4期 JournalofChongqingNormalUniversity(NaturalScienceEdition) Vo1．24 No．4 四阶抛物型方程的样条子域精细积分配置法 刘桂牙0，刘牙1斌 (广西民族大学 数学与计算机科学学院，南宁530006) 摘 要：对于四阶抛物模型方程周期初值 问题，可用有限差分方法进行求解。通常的有限差分方法在使用过程中受 到精度和稳定性的限制。本文首先将四阶抛物型方程转化为一个二阶的偏微分方程组，然后对时间项采用子域精 细积分的方法、空间项采用三次样条基本公式进行离散，得到了一个含参数 0(ot~h)的无条件稳定的差分格式， 所得到的差分方程为五点、两层隐格式，它的局部截断误差为0(r+盯 +h)。r，h分别为时间及空间步长，最后 的数值实验表明，本文的方法具有很好的数值精度和 良好的实用性。 关键词：四阶抛物型方程；子域精细积分；配置法 中图分类号：0241．8 文献标识码：A 文章编号：1672—6693(2007】04—0033—40 近年来对高阶抛物型方程和方程组的研究逐渐 二阶偏微分方程组 增多。文献[1]中构造了两个格式，一个格式的截 ， aM a0 q 2 ＼ llaf 一～, (2) 断误差为0( +h)，另一个为 0( +h)；文献 … { 2 [2]构造了一族三层 (特殊情况下为两层)、含双参 ，0 1，0 T (3) 数、绝对稳定、高精度、五对角线型的隐式差分格式， 其局部截断误差为0( +h)；文献[34]提出了一 取时间步长 ，空间步长h=L／M(M为正整数)，函 类含有多个参数的差分格式，其局部截断误差阶能 数U(，t)在节点 (，t)= ( ，n)处的数值解记 为U，函数 (，t)在节点( ，t)=( ，n)处的数 值解记为 n。令p=一 ，为了构造精度高而且稳定 性好 的子域精细积分格式，引入 一个 附加 项 flot+O ：0，一。。 。。，0≤f≤T x ’ O／U( 0，O／《h是参数)，加到(2)式的两端，得 1M(+ ，f)=M(，f)，一。。 。。，0≤f≤T dui _—d- +十O／Ui：q (4斗) tu(，0)=／ )，一。。 。。 t 其中q=P+O／U。如果q取某一个常数，则利用常数 首先将方程(1)式转化为二阶的偏微分方程组，然 变易法可解得(4)式的通解 后利用钟万勰提出的子域精细积分思想 刮和三次 M(f)：Ce一 +q (5) 样条函数逼近 ，构造 了一个