matlab在常微分方程的应用,Matlab在常微分方程教学中的应用

1 引言 科学技术和工程中大量的问题都表达为常微分方程的形式， 特别是描述系统的动态演变时， 如 R L C 电路、 数学摆、 人口模型、 生态模型、 化学反应过程等都表达为以时间 t 为独立变量的常微分方程或方程组， 所以常微分方程非常重要．但在实际中， 除了一些简单或特殊的微分方程， 我们可以用公式或特定的方法求解外， 大部分方程要给出其解析解常常是非常困难的， 我们只能借助于数学软件给出其数值解．而现行的各种常微分方程教材， 对数值解法的介绍都涉及不多， 而且受限于常规的教学方式， 给学生讲解数值解法也比较困难．另外， 传统的常微分方程的教学， 往往偏重于理论学习， 如一阶初等微分方程的初等解法，我们往往偏重于求解微分方程的一系列方法， 而对于抽象的方程的解对应的积分曲线和积分曲线族， 以及一些与几何联。 系紧密的概念如奇解等往往很难给学 生以直观的解释， 致使学生很难理解相关概念和解的几何意义． M a tla b 是由美国的 C le v e M o le r 博士于 1 9 8 0 年提出的， 它以矩阵运算为基础， 把计算、 可视化、 程序设计融合到了一个简单易用的交互式工作环境中．可实现工程计算、 算法研究、 符号运算、 建模和仿真、 原型开发、 数据分析及可视化、 科学和工程绘图、 应用程序设计等功能．M a tla b 强大的运算功能和图形使其成为目前世界上应用最为广泛的科学计算软件之一．鉴于此， 我们可以借助于 M a tla b 来介绍常微分方程的数值解法， 并且利用它的绘图功能来演示像奇解这样有明显几何意义的概念的几何特征． 下面我们对本文结构作一简要介绍， 第二部分对常用的求解常微分方程的数值方法如欧拉方法、 改进欧拉方法和 R u n g e - K u tta 方法等作一简单介绍； 第三部分结合具体实例介绍一下 M a tla b 在常微分方程教学中的一些应用； 第四部分我们对借助 M a tla b 进行常微分方程教学的意义作一简要评述． 2 数值方法介绍 下面我们针对微分方程的初值问题 M a tla b 在常微分方程教学中的应用 闫金亮 (武夷学院 数学与计算机系， 福建 武夷山 3 5 4 3 0 0 ) 摘 要： 本文主要介绍了 M a tla b 在常微分方程教学中的一些应用， 如利用数值方法求解常微分方程的数值解、 利 用 M a tla b 来描绘常微分方程解曲线及方向场、 利用 M a tla b 描绘常微分方程奇解的几何意义等， 最后我们对在常微分方程教学中使用 M a tla b 的意义作一简单的评价。 关键词： 常微分方程； M a tla b ； 方向场； 奇解 中图分类号： O 1 7 5 ．1 文献标识码： A 文章编号： 1 6 7 4 － 2 1 0 9 (2 0 1 2 )0 1 － 0 0 9 5 － 0 5 收稿日期： 2 0 1 2 - 0 2 - 2 1 作者简介： 闫金亮(1 9 7 9 - )， 男， 汉族， 硕士， 助教， 主要研究方向： 微分方程数值解。 第 3 1 卷 第 2 期 武夷学院学报 Vol．3 1 No．2 2 0 1 2 年 4 月 JOURNALOFWUYIUNIVERSITY A PR .2 0 1 2 武夷学院学报 2012 年第 2 期 d y d x = f( x ， y ) ， y ( x 0) = y 0 分别介绍欧拉方法、 改进欧拉方法以及 R u n g e - K u tta 方法． 微积分基本定理(牛顿— 莱布尼茨公式) 设 f