Matlab语言中边值问题算法的改进-精品.pdfVIP

维普资讯 第 23卷 第 i期 武 汉 理 工 大 学 学 报 Vo1．g3 N0 1 2001年 1月 JOURNAL OFW UHAN UNIVERSrrY OFTECHNOLOGY Jan．2O01 Matlab语言中边值 问题算法的改进 陆永忠 廖道训 黄其柏 (华 中科技太学) 摘 要 指出在Matlab语言中波有现成计算边值问题的M 程序，并且在连用打靶法解边值问题的过程 中、往往 出现对初值要求 比较苛刘 。针对此缺陷本文对选代方程进行适当的延拓 ，并壤韵 了Matlab边值 问是 的程序 ，最后用实倒验证所缩写的程序 。 关键词： 边值； 算i}； 打靶击 中图法分类号： TP3¨ 众所周知，在当今 30多个数学类科技应用软件中，一般有数值计算 (NumberCrunching)型和数学分折 (MathAnalysis)型两类 。目前MATLAB已经占据了数值型软件市场 的主导地位 。它 由主包和可选工具包组 成，功能强大，已教广泛地应用于研究和解决各种具体工程问题。但也存在着不足，例如在解边值问题方面就 不是很令人满意，没有解边值问题的计算程序。为此，结合科研实践，对语言中关于边值问题的解法作了一些 发展。 在Matlab中，与微分方程的数值有关的函数有 ODE23，ODE45，ODE113，ODE]5S，ODE23S，ODE． SEF，ODEGET，ODEPIOT，ODEPHAS2，ODEPHAS3，ODEPRINT，ORBITODE，ORBT2ODE，RIGI IX)DE，VDPODE；与符号微分方程有关 的有 DSOIVE，s0IVE，SUBS。这些函数有丰富的晦涵 ，为解决实 际问题提供了极大的便利，但他们都是在初始条件已知的基础上进行的，很少涉及边值条件的问题。常微分 方程边值问题在应用科学与工程技术中是经常遇到的，由于边值问题解的存在性比初值问题复杂得多，至今 没有很好解决，因此求解边值问题无论在理论上或数值解法上都比初值问题田难，目前很多数值方法都是假 定边值问题解存在的前提讨论的．其 中很重要的一类方法是将边值问题转化为初值 问题求解 ，如打靶法 (Shooting)，连续法 (Continuation)，和不变嵌人法 (Invariantimbedding)均属这一类 ，有限差分法也是另一类 很重要的经典算法。本文是基于在非 自治系统中应用打靶法(Shooting)来进行讨论的。 l 打靶法 大多数非线性动力系统都是非 自治系统 ，不失一般性 ，将对如下的非 自治系统进行分析 ： ． 』l=f(x，，)，(，f，)∈(R，R，R) f_) If(x，t+T，)：f(x，，) 一 其中丁为系统的周期 。求这类非线性动力系统的周期解 ，就是求下列边值问题 ： f =f(x，，g )) Ir0)一 (0，)一 (丁，)= 0 … 其 中g()为方程中与 有关的部分。∈ 为待求参数，一旦 确定，初值x(O，)也就 已知 ，解初值问题 ，便可 求的z(，，，)。对于 的求解，采用Newton--Raphson遗代法，其迭代格式为： s 一 ～ D (5，一 r(s) 格就 Oz(T,s)aT 。‘ 一 一 一