奇异微分方程边值问题数值解法.pdfVIP

靖ll卷鼻}l期 计 算 物 理 v。1．11，N。．1 1994年 3月 CHINESE JOURNALOF COM PUTATIoNAL PHYSICS Mar．，1994 奇异微分方程边值 问题的数值解法 ，7 A摘要… 出求核 究~…61化0041…) … … 关键词 奇异边值问题 正则化方法 差分方法 引 言 豸 { 柽／ 耘簋咩 在整个 区域上直接应用差分方法求解奇异边值 问题 ，将遇到奇点处理 以及跨奇点差分近 似的困难 配置法、Ritz—Garlerkin方法解这类 问题 ，误差估计阶数低，收敛速度差．本 文的方法避开奇点处直接 分近似，从而有效地克服了微分方程奇异性和解的奇异行为对数 值计算带来 的困难，准确地刻画出解的奇异行为 考虑奇异边值 问题 +L-,() +．． ( ) (1_1) 系数 ()中至少有一个在求解 区域上一点 = 存在奇异性 ，即函数值无界 ，除奇点 ，外 ()连续可微 ，F(x)为连续 函数 微分方程 的奇点可 以出现在 内部点和边界点上 ，由于奇 点的存在，使得解 中可能存在着奇异的因素，因而对内部奇点和边界奇点的边值问题解的要求 有一定的差异。奇点出现在边界点上，边界条件对解加 以限制，则要求 “和直到 阶导数存 在而且有界 ，P为边界条件 中奇点处导数 的最高阶数。奇点出现在 内部点上 ，边界点非解 的 奇点不必作此要求 。 线性独立的边界条件为 ： 三1[AIj!_)+。 ]=岛，，=l，2，…， (1．2) 其 中 ≤n，解的奇异行为可能使边界奇点 问题 的边界条件个数 比n来得少 设在奇点的某邻域 内 (11)的通解可写为 “()：∑ R()+R什t() ， ≤n (1_3) 其 中R．()， 1，2，…，s是齐次微分方程 Lu=0的线性独立解，R ()是 (1．1)的特解。 R()由具有正则形式的项m()∑口 组成，m()是初等函数，R()的一般推导参见 1992年 7月 16几 到原稿，1993年 1月 19 口收到修改稿 l8 计 算 物 理 第 l1卷 X~ [II。并假定边界奇点处的导数 【 ] 存在。 2 奇异边值 问题的正则化方法 2．I 内部奇点边值 问曩 设 ．、 为常数 ，口 ，2b，且 。Xs ：。在奇点 的邻域 【 ，叫 内解展开成 』l_3)， =n，余下 区间 【口，。]和 【 ，b】上存在两个正则边值 问题， 通过 ． 以及端点处的边界条件耦合在一起，应用差分近似将两个正则问题离散后，构成 一 个完整的线性代数系统 。下面推导 6．6处 的边界条件 ，微分 (1．3)得 r ．( )g+…+ ( )g= “( )一R+．(．) J ： ： (2．1．1) LI ’ ’ ：” )g．+…+-qd~-‘(．) =”“_。( )一 ”( ) 和 fR．(6)g．+ ·-·+R (6) 一 “(6)一R+。【6：) (2．1．2) “ “ ”