四_微分方程实验报告.docVIP

撰写人姓名： 撰写时间： 审查人姓名： 实验四 全 过 程 记 录 实验 名称 利用Matlab求解微分方程 时间 4学时 地点 数学实验室 姓 名 学 号 班 级 实验项目性质：普通实验 所属课程名称：数学建模 实验参考资料：详见电子课件； 实验目的： 1、利用MATLAB求解微分方程的解析解或初值问题的数值方法； 2、通过实际例题学习用微分方程模型解决简单的实际问题； 3、了解龙格-库塔方法的基本思想和计算公式，及稳定性等概念。 二、实验内容 1、掌握用MATLAB求解微分方程的解析解或初值问题的数值方法 2、用微分方程模型解决简单的实际问题。 三、实验报告要求 1、实验报告格式严格按哈尔滨理工大学有关规定要求； 2、应在理解的基础上简单扼要的书写实验原理,实验方法和步骤； （包括程序、运行结果、结果的解释） 3、就观察到的现象,变化的规律给出相应的解释； 4、对实验中存在的问题,进一步的想法等进行讨论。 四、实验用仪器设备及材料 软件需求： 操作系统：Windows 2000或更新的版本 实用数学软件：Matlab 硬件需求： Pentium Ⅲ 450以上的CPU处理器、64MB以上的内存、500MB的自由硬盘空间、 CD-ROM驱动器、打印机、打印纸等。 五、实验原理 常微分方程相关理论。 六、实验步骤 6.1 求解微分方程解析解与初值问题数值解的相关命令基本用法 1、求解微分方程； 2、求解微分方程； 3、求解线性微分方程组，其中矩阵 。 4、求解初值问题的解析解及数值解： 6.2 建模实例 田鼠与天敌数学模型 在生物的种群关系中，一种生物以另一种生物为食，称为捕食。一般，由于捕食关系，当捕食动物数量增长时，被捕食动物的数量即逐渐下降，捕食动物由于食物来源短缺，数量也随之下降，而被捕食动物数量却随之上升。这样周而复始，捕食动物与被捕食动物的数量随时间形成周期性的震荡。表1和表2分别给出了每隔两个月在田间调查的田鼠及其天敌种群数量的记录。 表1 田鼠种群数量 29.7 33.1 32.5 69.1 134.2 236.0 269.6 162.3 69.6 39.8 34.0 20.7 21.7 37.6 57.6 124.6 215.8 272.7 195.7 95.0 41.9 25.7 10.9 22.6 33.6 48.1 92.5 183.3 268.5 230 111.1 表2 田鼠的天敌种群数量 1.6 1.3 1.1 1.2 1.1 1.3 1.8 2.2 2.4 2.1 1.9 1.5 1.5 1.2 0.9 1.1 1.3 1.6 2.3 2.4 2.1 1.7 1.8 1.4 1.2 1.0 0.9 1.1 1.3 1.9 2.3 试用表1和表2分别给出的数据说明田鼠及其天敌的田间种群消长动态规律也是如此。 试用数学模型来概括这种现象，并总结出其数量变化的近似公式。 6.3 刚性方程数值解（ode23、ode45等与ode23s、ode15s等的区别） 求解上述方程组的解析解； 分别使用ode23与ode23s求解上述方程组的数值解，取计算区间为[0，1]； 分别使用ode23与ode23s求解上述方程组的数值解，取计算区间为[0，10]。 成绩评定： 指导教师：