求解一类三阶边值问题的Adomian分解方法.pdfVIP

第 7卷第 1期 杭州师范学院学报(自然科学版) Vo1．7 NO．1 2008年 1月 Journal of Hangzhou Teachers College(Natural Science Edition) JaD-．2008 文章编号：1008—9403(2008)01—0035—05 求解一类三阶边值问题的Adomian分解方法 金 怡 ，陈 靓 (1．杭州师范大学理学院，浙32-杭州 310036；2．杭州第一技师学院 园林系，浙江杭州 310005) 摘 要：利用 Adomian分解方法近似求解一类产生于物理问题中的三阶障碍边值问题．给出了实现该方法 的数值例子，并与其他方法比较，以验证该方法的有效性． 关键词：Adomian分解法；三阶边值问题；障碍问题；变分不等式 中图分类号：O175．8 MSC2000：34B10 文献标志码：A 0 引 言 考虑利用Adomian分解方法去获得如下一类三阶边值问题的近似解： r厂(z)， a≤ z≤ C， “ 一 J g(z)“(z)+厂(z)+r， c≤z≤d， (1) l厂(z)， d≤z≤b， 其中“(n)一a ，“ (n)一 ，“ (6)一 ，且“，“ 和“ 在c和d点连续．函数厂和g分别在[n，6]和[c， ] 上连续，系数r，a-， 和 是实常数．该系统产生于物理中的障碍、单侧、移动和自由边值问题，在纯理论科 学和应用科学的其他分支中有重要应用．2003年A1-Said和Noorc 使用三次样条方法解决了(1)的特殊形式： {0 0≤z≤丢， 1 “ 一 “ z)一1，丢≤z≤ 3， (2) 1 0 导≤z≤1， 其中“(0)一0,u／(0)一0，“ (1)一0，且“，“ 和“ 在{和 3点连续．此前，A1一said等 5 已给出了求解 问题(1)的一阶、二阶有限差分方法．最近，Siraj—ul—Islam等 。 又给出了求解该问题的非多项式样条方法． 该文的目标是要利用Adomian 和Wazwaz[ 提出的Adomian分解方法来求解问题(1)，得到高精度 的近似解．在第2部分，构造了求解三阶边值问题(1)的分解格式．在第3部分，用了一个典型的例子并与 数种方法比较以验证该方法的有效性． 1 Adomian分解方法 对于三阶边值问题： 收稿日期：2007—06—22 作者简介：金 怡(1983 )，女，浙江杭州人，应用数学硕士研究生，主要从事计算数学研究 36 杭州师范学院学报(自然科学版) 2008年 U =g( )“( )+ (z)，C≤ ≤ d， (3) 记微分算符L 为L 一参．将算符L 作用与方程(3)的两端得 L U==：g(x)u+，( )． (4) 定义逆算彳奇Lx_,为_个三重积分算符：L (·)一j JIJ(·)d d dx．将L 作用于方程(4)有。 ，，， 、 “( )一“(c)+“ (c)( —c)+ 嘉 ( —c) +L 一 (g( )“)+L 一 ( ( ))==： A+B(x—c)+ c( —c)。+L 叫(g( )“)+L 1(，( ))，