基于MATLAB的微分方程组的数值计算.pdf

臣翟圈絮㈣。洲乱…N，渊… 学术论坛 基于MATLAB的微分方程组的数值计算 刘双杰 (沈阳理工大学 沈阳 1 10168) 摘要：传统的解簟分方程诅的方法有近似分析解法，表解法和圈解法。这些方法都有局限挂．电子计算机墙码的出现及其应用，不仅 代瞢了繁重的人工求并．而且改变了传统的研完方法。 MATLAB是一种基于矩阵的数学软件包，蕾软停包包括了一个复值程序扩展库， 并且有高衄稿程格式。应用四阶五氲龙格岸譬法薯翻Matlsb程序对截分方程诅进行求并．蛄果表明无论是由茂最是特珠点与试验实测 值一致性郜比较好 关键词：|t分方程组 t值计算 仿真 MATLAB 中图分类号。0241 文献标识码：A 文章编号：1672—379l(2009)02(c)一0238一ol 传统的解微分方程组的方法有近似 分析解法．表解法和图解法。这些方 ，(血∞ 法有一定的局限性。MATLAB是一种 基于矩阵的数学软件包，该软件包包括 了一个数值程序扩展库，并且有高级编 程格式。应用MATLAB工具箱中自带 的四阶五级的龙格库塔法(ode45)子程 序对微分方程组进行求解，简化了繁琐 的过程。 求解常微分方程组的数值方法不止 龙格一库塔法一种，还有亚当姆斯法 图3 p-t曲线 图4 v--I曲线 等。龙格一库塔法虽然工作量较大，但 是它在计算值时，只用到，而不直接依 赖于和等。也就是说．在初值确定后， 就可以依次计算和等等。直至整个过 程结束，不存在计算起步的问题。另 外，这种方法没有规定后一步的步长与 前一步步长必须满足的某种关系，可以 任意改变计算步长。龙格一库塔方法每 一步须要4次计算，看起来计算过程烦 琐，实际上在MATLAB运用程序进行 计算非常快速准确。 1微分方程组的程序框图 炳|) 编程思路的逻辑框图为图1和图2。 圈5 Z-t曲线 2微分方程组的仿真计算 由某实例受力及运动分析，建立起如 下方程组： ’I，=zZ(I+)cZ十肛z。)【0≤Zs” (1) 囤2 RK子程序图 些：堕矿 (2) 联立，利用四阶五级龙格一库塔法。借助 dt e．‘ 参考文献 MATLAB软件编制微分方程组的求解程 7． 坐：。 f31 【l】楼顺天，姚若玉．沈俊霞．MATLAB d， 序，对该数学模型进行仿真计算，得出p， x程序设计语言【M】．西安电子科技大学 印=‘P…idv (4) y，1，z关于t的图像，计算结果如图3、图 出版社，2007． pS(／+‘)=细一兰币∥ …4和图5。 京：北京航空航天大学出版社，2003． 3结语 仿真计算结果与实测的p-t、v-t、Z- 1 0．35 0≤Z≤0．06 t曲线一致性较好．表明该微分方程组的计 o．06z≤o．56 UI(z)={o．796Z+0．3022 算正确