常微分方程数值解及其Matlab实现.pdfVIP

科技信息 博士·专家论坛 常微分方程数值积及其Matlab实坝 襄类学院数学系 武瑞蝉华东交通大学 邓黎 [摘要]本文对常豢分方程初值问道竹薮值解法作1进一步探讨，并运用Manab对其中较难以求解的疫进尤拉法和 四阶龙格库塔法进行编程实现，程序简洁、直观，求并速度快，方法实用性较强． [关键词]擞分方程数值解法M“lBb 常微分方程的求解是现代科学研究和工程技术中经常遇到的实际问题，然而，从实际问题中建立出来的微分方程往往具有 非常复杂的形式，有些解析式难以计算，有些则根本不能用解析式来表达，所以利用数值解法“1。求解实际问题就显得非常重 要。本文在总结常微分方程数值解法的基础上给出了Matlab“1求解常微分方程的多种方法．并对其中较为复杂的改进尤拉法和 四阶龙格库塔法给出了Matlab程序． 1一阶微分方程初值问题的数值解法 对于一阶方程的初值问题 』-y：2f(x，y)(口≤x≤6) ly(xo)2％ 从理沦上说，只要函数，(x，y)适当光滑，譬如关于，满足李氏条件 厂(x，y)一厂(x，歹)l≤三Iy一歹 就可以保证初值问题(1)的解y=y(x)存在且惟一． 求解问题(1)的数值解法又可称之为差分方法，因为常用的数值解法，如尤拉法，改进尤拉挂，龙格厍培珐等等本质上 均是将导数化为差商来做近似。 尤拉法的建立可以有多种，运用差商方法、数值积分方法和泰勒展式法等都可以得到如下关系式 』以“2以+矽(吒，儿) 。 (2) 【乩2y(‰) 尤拉法是其中展简单也摄容易分析的显式方法，求解精度为一阶“圳，它是其它数值解法推导的基础．由尤拉法的思想可 以相应推出隐式尤拉法-两步尤拉法以及梯形公式等-由于尤拉法是显式算法t计算量小·但精度低，而梯形公式虽然提高了 精度，但是一种隐式算法，计算量大．故将两者进行综合改进，得到改进尤拉法： j或+·2虬+6(¨一2_7) (3) ln+，=虬+Ⅱ(¨一2矗／％)+(只+。一2矗／只+．)]，2 定理改进尤拉法的精度是二阶的． 证明对(3)作如下变形 f％+。=％+向(^十如)／2 ． {毛=厂(_，儿) 【屯=厂(‰。M+饨) 基金项目，目孪自燕科学基仝啊目膏助(批准号 作者简介-武瑞蝉，士，讲师。主薹研宽期碱为复杂康垃、并行计算午 一b一 万方数据 科技信息 博士·专家论坛 在儿=y(_)前提下，将岛，心作泰勒屉开，则 尼。=／(％，％) ☆：=，(％+自，此+胁：，)=，(％，)+妩+^墨‘+o(^2) =y’(_)+咖”(“)+D(^2) 以+．=n+矗(y’(‘)+y’(矗)+砂”(％))／2+D(旷) ：，(矗)+吵-(矗)+嬖y”(‘)+。(矿) y(_。)：y(_)+砂t(％)+等y”(_)+。(向，) 故y(‘+．)一¨+，=D(矿)·故而此差分格式为二阶格式a证毕。 高阶龙格库塔方法是在尤拉法和改进尤拉法的基础上推导小来的，它的基本思想是通过增加计算厂(x，y)的次数来提高 问题的精度。由以上定理的证明思路可以得出以下四阶标准龙格库塔法的一般形式： 儿“=以+(向+2乞+2岛+缸)／6 屯=可(吒，此) 如=矿(％+五／2，咒+矗／2) 尼3=矿(■