基于三次B样条求解四阶线性边值问题-evcarexpocom.PDFVIP

第30卷 第6期 黑 龙 江 大 学 自然 科 学 学报 Vo1．30 No．6 2013年 12月 JOURNALOFNATURALSCIENCEOFHEILONGJIANGUNIVERSITY December，2013 基于三次 B样条求解 四阶线性边值 问题 李明珠 一， 丁效华 ， 王怡东 (1．青岛理工大学 理学院，青岛266033；2．哈尔滨工业大学 (威海)数学系，威海 264209) 摘 要：fi-1用三次B样条求解四阶线性微分方程边值 问题的近似解。通过求解一个(n+3)维 线性方程组，得到对应三次B样条函数的待定系数，进而得到边值 问题的三次B样条近似解。利 用本方法同时能够得到方程精确解的导数近似值，算法收敛且有二阶精度。数值算例说明了方法 的有效性；与其他方法相比，本方法具有更高的精度 。 关键词 ：三次B样条 ；边值 问题 ；微分方程 ；线性方程组 中图分类号 ：O241．8 文献标志码 ：A 文章编号：1001—7011(2013)06—0711—05 0 引 言 带边值条件的四阶线性常微分方程作为一类重要数学模型，广泛地应用在弹性流、非弹性流、变形梁理 论 以及板弯 曲理论中 ]。这类方程解 的存在性和唯一性在理论上已经解决 j，但精确解却很难获得，因此 对其数值方法的研究有很大的理论意义和应用价值。 本文考虑如下方程 ： Y‘()+a(x)y()=g()， E[a，b]， (1) 及边界条件 ： Y(a)=o／。， Y(a)= ， Y(b)=3／。，Y(b)=3／。 (2) 其中 )和 g(x)为[a，b]上的连续函数，％，Ol， ，卢 是常数。 对于 四阶线性边值问题的研究主要包括打靶法、有限差分法、配置法 、Ritz．Galerkin法、Sinc．Galerkin法 等 j。近年，许多学者利用样条函数解决边值问题，得到大量有意义的结果。文献[6]最早提出使用样条 函数处理边值问题；文献 [7]首先利用三次样条函数解决四阶线性边值问题，收敛阶达到二阶；文献[8]利用 五次样条函数解决四阶线性边值问题，收敛阶达到六阶，但计算量 比较大；文献 [9—10]基于非多项式 {1，， ， ， sinkx，coskx}的样条函数解决四阶边值问题，收敛阶为二阶但计算复杂；文献 [11]提出利用五次 B样 条函数解决带障碍 四阶边值问题 。三次B样条是最重要样条函数之一，具有许多优 良的性质，应用十分广 泛。本文将利用三次 B样条对方程 (1)～(2)进行求解 。 1 三次B样条 1．1 三次 B样条定义 将区间，=[a，b]进行 凡等分，分割为 凡个小区间 = ，…](i=0，1，…，凡一1)，其 中节点为 =n+ 访 ( 0，1，…，n)， =b-a ， 则三次 B样条 Ⅳ ()( 2，…n一2)表达式定义如下 ： 收稿 日期 ：2013—05—14 基金项 目：国家 自然科学基金资助项 目；山东省 自然科学基金青年项 目(ZR2012AQo27) 作者简介 ：李明珠 (1982一)，男，讲师，博士研究生，主要研究方 向：微分方程数值解 ，E—mail：limingzh126．con 通讯作者：丁效华(1968一)，男，教授 ，博士，博士生导师，E-mail：mathdxh@126．COB 引文格式：李明珠，丁效华 ，王怡东．基于三次B样条求解四阶线性边值问题[J]．黑龙江大学自然科学学报，2013，30(6)：711—715 · 712 · 黑 龙 江 大 学 自 然 科 学 学 报 第30卷 戈∈[。一2，一1