八阶奇异边值问题精确解的表达形式.pdfVIP

数学物理学报 http：／／actams．wipm．ac．cn 八阶奇异边值问题精确解的表达形式 ， 。么焕 民 林迎珍 f 哈尔滨Z-业大学数学系 山东威海 264209； 。哈尔滨师范大学数学科学学院 哈尔滨 150025) 摘要：该文在再生核空间 0【，1]中给出了求解八阶奇异边值问题的新算法．方程的精确解 以级数形式给出．算例及数值结果验证了方法的实用性和有效性． 关键词：八阶奇异边值问题；常微分方程；再生核空间． MR(2000)主题分类：34B05；34K28；47B32 中图分类号：0241．81 文献标识码：A 文章编号：1003—3998(2010)01—103—11 1 引言 奇异常微分方程出现在气体动力学，牛顿流体力学，流体力学，流体动力学，弹性过程， 反应扩散过程，化学动力学和应用数学的其它一些分支．近几年，奇异边值理论已经成为非 常重要的研究领域 (见文献 『卜51)． 在文献 f61中，作者介绍了高阶方程最大值理论并发展了在下解和上解情况下的单调方 法 文献 [7】给出奇异边值问题解存在的一个充分条件．文献 [81讨论了超线性奇异边值问 题．文献 9【】和 [10】研究了2n阶奇异次线性和超线性m 点边值问题正解的存在性． 然而，很少有文献研究八阶边值和特征值 问题的数值解．文献 f111利用 Adomian分解 方法讨论了八阶边值问题的解．文献 1『21指出：如果从下面加热流体的无限横向层并进行旋 转，将发生不稳定．当这种不稳定过度发生时，将以八阶常微分方程为模型．文献 [13]给出 八阶边值 问题解存在和唯一的条件定理．文献 f141发展了八阶边值问题的有限差分方法． 文献 f151进一步研究了在热不稳定中出现的八阶，十阶，十二阶特征值问题的数值解．文献 f16]用octic样条法给出了八阶边值问题的解．在文献 [17]中，作者用nonic样条法讨论了 以下特殊情形的线性八阶边值问题的数值解 f()()+f(x)y(x)=夕()，∈ 6]， {(n)=o，(6)=Ql，yra)=加，yl(b)= (1．1) I a)= o， (6)：l，yma)=，0，ym(扫)= l， 其中OL，， 和 ui(i=0，1)是常数，_厂()，g(x)是 [n，b]上的连续函数．运用导数的连续性 得到了在结点处样条函数和其导数值之间的一致性关系．他们还用非多项式样条法研究了特 收稿 日期：2008—12—25；修订 日期：2009—10—09 E—mail：liliy55@163．corn；yhm hsd@126．com 基金项 目：国家自然科学基金 、黑龙江省自然科学基金 (A2007—11)和黑龙江省骨干教师计划项 目 (1151G019)资助 104 数 学 物 理 学 报 Vo1．30A 殊情形的十阶线性边值问题的数值解并发展了文献 1【8]中的quintic样条法，得到的结果同 文献 [19—20】中的结果进行了比较． 本文在再生核空间 [o．1]中给出了求解以下八阶边值问题的有效方法 “ (8’()+血()(6()+b()(4(z)+c(z)(2(22)+d()“()=，()．X∈(0，1)， (1．2) Ju(0)=。，u(1)=l?u(0)=。-u()：1一 f1．3) Iu，，(0)= o，u(1)=一y1，u，，(o)=o，1“，，(1)= 1． 其中 ．， 和 5i(i=0，1)是常数．a(z)．b()，e(z)，d(x)和 f(x)是 (0，1)上的连续函数． n()，6()，c(z)和 d(z)可能在 =0或