浅谈常微分方程的数值解法及其应用文献综述范文参考

文献综述浅谈常微分方程的数值解法及其应用一、前言部分微分方程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解 .牛顿在建立微积分的同时， 对简单的微分方程用级数来求解 .后来瑞士数学家雅各布?贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论 .微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法 .微分方程也就成了最有生命力的数学分支.总之，力学、天文学、几何学等领域的许多问题都导致微分方程 .在当代，甚至许多社会科学的问题亦导致微分方程，如人口发展模型、交通流模型等 .因而微分方程的研究是与人类社会密切相关的 .[1]“常微分方程”是理学院数学系所有专业学生的重要专业基础课之一，也是工科、经济等专业必学内容之一.其重要性在于它是各种精确自然科学、社会科学中表述基本定律和各种问题的根本工具之一，换句话说，只要根据实际背景，列出了相应的微分方程，并且能（数值地或定性地）求出这种方程的解，人们就可以预见到，在已知条件下这种或那种“运动”过程将怎样进行，或者为了实现人们所希望的某种“运动”应该怎样设计必要的装置和条件等等.例如，我们要设计人造卫星轨道，首先，根据力学原理，建立卫星运动的微分方程，列出初始条件，然后求出解，即卫星运行轨道 .随着物理科学所研究的现象在广度和深度两方面的扩展，微分方程的应用范围更广泛.[2]从数学自身的角度看，微分方程的求解促使数学在函数论、变分法、级数展开、常微分方程、代数、微分几何等各方面进行发展.从这个角度说，微分方程变成了数学的中心[3].总之，微分方程从它诞生起即日益成为人类认识并进而改造自然、社会的有力工具，成为数学科学联系实际的主要途径之一.文章就常微分的数值解法以及应用展开简单的论述。二、主体部分微分方程概念介绍微分方程概况由一元函数得到的方程 .即：称含有自变量，未知函数及其导数的关系式F(x,y,dy,d2y2,...,dnny)0.（1）dxdxdx为常微分方程.其中出现的最高阶导数的阶数，叫做常微分方程的阶.例如dyx，dxdyd2g0是二阶常微分方程.设y(x)定义于y，是一阶常微分方程.2sindxdtp区间J上，有直到n阶的导数，将它代入（1），使（1）变成关于x的恒等式，即F(x,(x),d(x),...,dn(x))0,xJ.dxdxn就称y＝(x)为（1）的一个定义于J上的解，并称J为该解的定义区间.[4]如果一个微分方程中出现多元函数的偏导数，或者说如果未知函数和几个变量有关，而且方程中出现未知函数对几个变量的导数，那么这种微分方程就是偏微分方程.2.2微分方程产生的历史背景微分方程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解.牛顿在建立微积分的同时，对简单的微分方程用级数来求解。后来瑞士数学家雅各布?贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的.数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具.[5]牛顿研究天体力学和机械力学的时候，利用了微分方程这个工具，从理论上得到了行星运动规律.后来，法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置.这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量.微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法 .微分方程也就成了最有生命力的数学分支.总之，力学、天文学、几何学等领域的许多问题都导致微分方程 .在当代，甚至许多社会科学的问题亦导致微分方程，如人口发展模型、交通流模型等 .因而微分方程的研究是与人类社会密切相关的.[6]2.3微分方程发展现状及其基本功能在数学学科内部的许多分支中，微分方程是常用的重要工具之一，微分方程进一步发展的需要，有推动着其它数学分支的发展；相反，微分方程每一步进展都离不开其他数学分支的支援.数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对微分方程的发展产生了深刻的影响.当前计算机的发展更是为微分方程的应用及理论研究提供了非常有力的工具.时至今日，可以说微分方程在所有自然科学领域和众多社会科学领域都有着广泛的应用，如自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等.只要能够列出相应的微分方程，有了解方程的方法，利用它就可以精确地表述事物变化所遵循的基本规律.从微积分理论形成以来，人们一直用微分方程来描述、解释或预见各种自然现象，不断的取得了显著的成效[7].2.4常微分方程的数值求解方法2.4.1Euler法Euler法是最简单的数值方法，[a,b]为求解良态初值问题y'f(t,y)，y(a)y0的区间。实际上，下面的过程不是要找到满足该初值问题的可微函数，而是要生成点集{(tk,yk)}，并且将这些点作为近似解，即y(tk)yk。如何构造“近似满足微方程”的“点集”呢？首先为这些点选择横坐标，为方便起见，将区间[a,b]划分为M个等距子区间，并选择网络点tkak，hk=0，1，，Mba（1）其中h值h称为步长。然后近似解My'f(t,y)在[t0,tM]上，y(t0)y0（2）设y(t)，y'(t)和y''(t)连续，；；利用泰勒定理将y(t)在tt0处展开，对每个值t，存在一个t0和t之间的值c1，使得''2y(t)y(t0)y'(t0)(tt0)y(c1)(tt0)（3）2'f(t0,y(t0))和将代人等式（），得到1的表示：ht1t0)3y(ty(t1)y(t0)hf(t0,y(t0))y''(c1)h2（4）2word文档 可自由复制编辑如果步长h足够小，则可以忽略2次项（包含h2的项），得到y1y0hf(t0,y0)（5）这就是欧拉近似。重复该过程，就能得到近似解曲线 y y(t)的一个点序列。欧拉方法的一般步骤是tk1tkh，yk1ykhf(tk,yk)其中k=0，1，，M-1[8]（6）2.4.2泰勒级数法泰勒级数法有着广泛的应用，并且是比较求解初值问题的各种不同数值方法的标准，它可设计为任意指定的精度。下面首先将泰勒定理用新的公式表示，使之适合于求解微分方程。定理9.5（泰勒定理）设y(t)CN1[t0,b]，且y(t)在不动点ttk[t0,b]处有N次泰勒级数展开：y(tkh)y(tk)hTN(tk,y(tk))O(hN1)（1）其中，TN(tk,y(tk))Ny(j)(tk)hj1（2）j1j!y(j)(t)fj1(t,y(t))表示函数f关t的（j1）次全导数。求导公式可以递归地计算：y'(t)fy''(t)ftfyy'ftfyfy(3)(t)ftt2fyty'fyy''fyy(y')2ftt2fytffyyfy(4)fttt3fytty'3fyyt(y')23ftyy''fyy'''3fyyy'y''(fttt3fyttf3fyytf2fyyyf3)fy(ftt2fytffyyf3(ftfyf)(fytfyyf)fy2(ftfyf)并且一般有y(N)(t)P(N1)f(t,y(t))其中P为导数算子P(fy)t

2fyf)fy(ftfyyy(y')32)

（3）（4）区间[t0,tM]上的初值问题 y'(t) f(t,y)的近似数值解可由各子区间 [tk,tk1]上的公式（1）来推导。 N次泰勒方法的一般步骤为d2h2d3h3dNhN（5）yk1ykd1h2!3!...N!其中在各步k0,1,......,M1有djy(j)(tk),j1,2,......,N。N次泰勒方法的最终全局误差是O(hN1)阶的，因此可选择所需大小的N，使得误差足够小。如果N是固定，则理论上可以推导出步长h，使之满足任意的最终全局误差。然而在实际运算中，通常用h和h/2计算两个近似结果集，然后比较其结果[9]。2.4.3龙格—库塔方法泰勒方法的优点是最终全局误差的阶为 O(hN)，并且可以通过选择较大的 N来得到较小的误差。然而泰勒方法的缺点是，需要先确定 N，并且要计算高阶导数，它们可能十分复杂。每个龙格一库塔（ Runge-Kutta ）方法都由一个合适的泰勒方法推导而来，使得其最终全局误差为 O(hN)。一种折中方法是每步进行若干次函数求值，从而省去高阶导数计算。这种方法可构造任意 N阶精度的近似公式。最常用的是 N=4的龙格一库塔方法，它适用于一般的应用，因为它非常精确、稳定，且易于编程。许多专家声称，没有必要使用更高阶的方法，因为提高的精度与增加的计算量相抵消。 如果需要更高的精度，则应该使用更小的步长或某种自适应方法。阶龙格一库塔方法（RK4）可模拟N=4的泰勒方法的精度。它基于如下yk1，的计算：yk1ykw1k1w2k2w3k3w4k4（1）其中k1，k2，k3和k4形如k1hf(tk,yk)k2hf(tka1h,ykb1k1)（2）k3hf(tka2h,ykb2k1b3k2)k4hf(tka3h,ykb4k1b5k2b6k3)通过与N=4阶的泰勒级数方法的系数匹配，使得局部误差为O(h5)，龙格和库塔得出了如下方程组：word文档 可自由复制编辑b1a1b2b3a2b4b5b6a3w1w2w3w41（3）w2a1w3a2w4a3121w2a12w3a22w4a323w2a13w3a23w4a3314w3a1b3w4(a1b51a2b6)6w3a1a2b3w4a3(a1b5a2b6)18w3a12b3w4(a12b5a22b6)1121w4a1b3b624该方程组有11个方程和13个未知量，必须补充两个条件才可以求解。最有用的选择是a11，b20（4）2其余变量的解为a21,a31,b11,b31,b40,b50,b61222（5）1,w21,w31,w41w16336将式（4）和（5）中的值代入式（2）和式（1），得到标准的N4阶龙格—库塔方法，其描述如下。自初始点(t0,y0)开始，利用yk1ykh(f12f22f3f4)6[10]生成近似值序列，其中f1f(tk,yk)f2f(tkh,ykhf1)22f3f(tkh,ykhf2)22f4f(tkh,ykhf3)

（6）（7）预报—校正方法word文档 可自由复制编辑欧拉方法、休恩方法、泰勒方法以及龙格一库塔方法都称为单步长方法，因为它们只利用前一个点的信息来计算下一个点，即计算 (t1,y1)时只使用了初始点 (t0,y0)。一般地，只有yk用来yk1。当计算出若干个点之后，就可以利用几个已计算出的点来计算下一个点。以亚当斯一巴什福斯4步法的推导为例，计算yk1需要yk3，yk2，yk1和yk。该方法不是自启动的，要生成点{(tk,yk):k4}，必须先给出其4个初始点(t0,y0)，(t1,y1)，(t2,y2)，(t3,y3)（可用前面各节中的方法完成）。多步法的一个优点是，可以确定它的局部截断误差（localtruncationerror，简称L.T.E.），并可以包含一个校正项，用于在每一步计算中改善解的精确度。该方法还可以确定步长是否小到能得到yk1的精确值，同时又大到能够免除不必要的和费时的计算。使用预报子和校正子的组合在每一步只需要进行两次函数f(t,y)求值[11]。亚当斯一巴什福斯一莫尔顿方法亚当斯一巴什福斯一莫尔顿方法（Adams—Bashforth—Moulton）是由基本微积分定理推导出的多步法：tk1y(tk1)y(tk)f(t,y(t))dt(1)tk预报子使用基于点(tk3,fk3),(tk2,fk2),(tk1,fk1)和(tk,fk)的f(t,y(t))的拉格朗日多项式逼近值，并在区间[tk,tk1]上对式（1）积分，这个过程产生亚当斯一巴什福斯预报子：pk1ykh(9fk337fk259fk155fk)（2）24校正子的推导类似。这时可以实用刚刚计算出的值pk1。基于点(tk2,fk2)，(tk1,fk1)(tk,fk)和新的点(tk1,fk1) (tk1,f(tk1,pk1))构造f(t,y(t))的一个新的拉格朗日多项式逼近，然后在区间 [tk,tk1]上对该多项式积分，即可得到亚当斯一莫尔顿校正子[12]：yk1ykh(fk25fk119fk9fk1)（3）24米尔恩—辛普森方法word文档 可自由复制编辑米尔恩—辛普森方法是预报子基于区间[tk3,tk1]上的对f(t,y(t))的积分：tk1y(tk1)y(tk3)f(t,y(t))dt（4）tk3预报子使用f(t,y(t))基于(tk3,fk3),(tk2,fk2),(tk和(tk,fk)的拉格1,fk1)朗日多项式逼近，在区间[tk3,tk1]上对它积分，得到米尔恩预报子：pk1yk34h(2fk2fk12fk)(5)3校正子的推导类似。此时值pk1已知，基于点(tk1,fk1)，(tk,fk)和新点(tk1,fk1)(tk1,f(tk1,pk1构))造f(t,y(t))的新的拉格朗日多项式，然后在区间[tk1,tk1]上对该多项式积分，结果为大家所熟悉的辛普森公式[13]：h(fk14fkfk1)（6）yk1yk132.5常微分方程的应用微分方程在传染病模型中的应用随着卫生设施的改善、 医疗水平的提高以及人类文明的不断发展， 诸如霍乱、天花、等曾经肆虐全球的传染性疾病已经得到有效的控制。 但是一些新的、不断变异着的传染病毒却悄悄向人类袭来。 20世纪80年代四分险恶的艾滋病毒开始肆虐全球，至今仍在蔓延；2003年春来历不明的非典病毒突袭人间，给人们的生命财产带来极大的危害。长期以来，建立传染病的数学模型来描述传染病的传播过程， 分析受感染人数的变化规律，探索制止传染病蔓延的手段等，一直是各国有关专家和官员关注的问题。这里我们介绍的是按照一般传播机理建立的数学模型如下。dp1(t)p2(t)p3(t))dp1(t)p3(t)p1(t)p2(t)，b(p1(t)dtdp2(t)(t)(dd)p2(t)，（1）p1(t)p2dtdp3(t)p2(t)(d3)p3(t)，dtt0:p1p10,p2p20,p3p30,其中pi00(i1,2,3)对上述常微分方程组进行分析求解，就可以了解不同时刻传染病的动力学特征 (比如：[14]传染病病情的发展趋势，即各类人的人口数量 pi(t)的分布情况) .2.5.2常微分方程在捕鱼业的应用设某水域现有鱼量x，由于受资源限制所能容纳的最大鱼量xm，高自然增长率r，捕捞增长率k，按人口的逻辑模型建立微分方程。dxrx(1x)kxdtxm要保持鱼量平衡dx0，设平衡点为x0，解得x0rkxm设dxdtrf(x),考虑f(x)在x0的泰勒展式dtf(x)f(x0)(xx0)0(xx0)f(x0)kr当f(x0)＞0时f(x)与xx0同号x0为不稳定平衡点当f(x0)＜0时f(x)与xx0异号x0为稳定平衡点f(x0)＜0即r＞k设f1(x)rx(1x)f2(x)kx由于k＜rxm曲线f1(x)与f2(x)有交点，因f1(x)在原点切线为yrx解得，易知当x0xm时，取得最大捕捞量1r，21rx0rxmkf2(x)224最大捕捞量为r[15]4xm三、总结部分微分方程是在解决实际问题的过程中产生的， 微分方程的研究又促进实际问题的解决，同时也促进其他学科的发展 .微分方程在物理、 工程、 力学、天文学、生物学、 医学、 经济学等诸多领域都有重要作用 .目前，微分方程的实际背景广、应用性强的特点已受到广泛关注.由于同一类型的微分方程往往可以用来描述许多性质上颇不相同的自然现象，对一些重要的微分方程开展研究，可以有多方面的应用前景，并可望在新兴学科或边缘学科的开发中及时地发挥作用 .[16]微分方程在几百年的发展史上， 在每个时期最新的技术科学中建立了自已的落脚点，得到了实际的应用， 从而刺激它飞跃地发展， 在完善理论的同时继续深入其它新兴技术学科领域近几十年来，世界科学技术进入核能、火箭、人造卫星、数字时代，微分方程定性理论及方法不论在应用上、理论上均不断地扩展着自身的领域，显示出前所未有的强大生命力 .它的理论和方法，过去和现在都对力学、天文、物理、化学、生物、各种技术科学 (如自动控制、无线电电子学等)及若干社会科学(如人口理论、经济预测等)提供了有力的工具，后者反过来也不断地向它提出新的间题，刺激着它不断地向前发展.[17]四、参考文献张良勇,董晓芳.常微分方程的起源与发展[J].高等函授学报(自然科学版),2006(3):34-38.[2] 黄赛.常微分方程发展的主要历史沿革 [J]. 教育与职业.2006(2):36-41.黄焕福.常微分方程课程建设初探[J].教育与职业,2007(14):146-148.林建平.常微分方程早期发展概观[J].南京工程学院学报,2001(2):1-4.[5] 周仲旺.几类特殊的常微分方程 [J]. 潍坊学院学报 .2003(6):7.邵晓锋,徐卫卫,李龙星.试析常微分方程模型的归结方法[J].黄冈职业技术学院报,2007(1):37-39.张同斌,管军军,焦万堂.基本微分方程的应用研究[J].河南科学,2007(6):888-890.胡建伟,汤怀民.微分方程数值解方法[M].科学出版社,2003(4).[9] JohnH.Mathews,Kurtis D.Fink.Numerical MethodsUsingMATLABFourth Edition[M].电子工业大学出版社 .2005(12).王高雄,周之铭等.常微分方程[M].高等教育出版社,2008(1).[11] JohnH.Mathews,Kurtis D.Fink.Numerical MethodsUsingMATLABFourth Edition[M].电子工业大学出版社 .2005(12).J.W.Thomas.偏微分方程的数值解法[M].世界图书出版公司,1995(8).郭晓梅.常微分方程的数值解法[J].网络财富，2009(8).姜启源,谢金星,叶俊.数学模型[M].高等教育出版社.2003(8).[15]JamesStewart.Calculus[M].FifthEdition.Beijing: 高等教育出版社 .2004(7):1141-1161.欧阳瑞,孙要伟.常微分方程在数学建模中的应用[J].宿州教育学院学报,2008(2):[17]JamesStewart.Calculus[M].FifthEdition.Beijing: 高等教育出版社 .2004(7):1141-1161.卫生管理制度1 总则1.1 为了加强公司的环境卫生管理，创造一个整洁、文明、温馨的购物、办公环境，根据《公共场所卫生管理条例》的要求，特制定本制度。1.2 集团公司的卫生管理部门设在企管部，并负责将集团公司的卫生区域详细划分到各部室，各分公司所辖区域卫生由分公司客服部负责划分，确保无遗漏。2 卫生标准2.1 室内卫生标准2.1.1 地面、墙面：无灰尘、无纸屑、无痰迹、无泡泡糖等粘合物、无积水，墙角无灰吊、无蜘蛛网。2.1.2 门、窗、玻璃、镜子、柱子、电梯、楼梯、灯具等，做到