微分方程数值解设计学案.doc

微分方程数值方法课程设计 Basic Theory of Ordinary Differential Equations Experiment Report 教 务 处 2014年 7 月 课 程 设 计 说 明 书 题目： 常微分方程数值解法课程设计 学院（系）： 理学院 年级专业： 计算科学11-1 学生姓名： 指导教师： 教师职称： 教授 燕山大学课程设计（论文）任务书 院（系）：理学院 教学单位： 学 号 XXX 学生姓名 XXX 专业（班级） 11计算数学 设计题目 常微分方程数值解法课程设计 设 计 内 容 比较Adams四阶PECE模式和PMECME模式。 求解贝塞尔方程并与精确解比较。 小型火箭初始重量为1400kg，其中包括1080kg 燃料。火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N 的推力，火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比于速度的平方，比例系数为0.4kg/m，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。 设 计 要 求 1、要通过课题背景的介绍，阐明选题的意义； 2、运用所学知识，建立问题的数学模型； 3、利用Matlab软件进行计算，并给出计算机程序框图及计算机程序； 4、运用所学知识，对计算结果进行分析； 工 作 量 1、程设计要在两周内完成； 2、根据已知题目找出解题方法，然后利用Matlab进行计算； 3、正文在3000字左右。 工 作 计 划 第一天：收集资料，阅读参考文献； 第二天：构建数学模型； 第三天：设计计算机流程图，编制计算机程序，上机运算，调试程序； 第四天：计算结果分析 第五天：撰写报告等 参 考 资 料 [1]胡建伟汤怀民。微分方程数值方法。科学出版社。1999年01月 [2] 陈渝周璐钱方等。数值方法（MATLAB版）。电子工业出版社。2002.6 [3]尹泽明，丁春丽等。精通MATLAB 6。清华大学出版社。2006.6 [4]Cleve B.Moler。MATLAB数值计算。机械工业出版社2006.6 指导教师签字 基层教学单位主任签字 说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。 年 月 日 燕山大学课程设计评审意见表 指导教师评语： 成绩： 指导教师： 年 月 日 答辩小组评语： 成绩： 评阅人： 年 月 日 课程设计总成绩： 答辩小组成员签字： 年 月 日 摘 要 本文对常微分方程初值问题现有的数值解法进行了综述研究。主要讨论了几种常用的数值解法：即欧拉法，改进欧拉法，龙格库塔方法，阿达姆斯PECE,PMECME格式等。文章最后结合常见数值解法，对较为典型的微分方程模型进行数值求解，探讨了上述数值算法在实际建模问题中的应用。 论文阐述的是常微分方程数值解法的几个问题，通过对以下问题的求解 一.比较Adams四阶PECE模式和PMECME模式。 二.求解贝塞尔方程并与精确解比较。 三.小型火箭初始重量为1400kg，其中包括1080kg 燃料。火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N 的推力，火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比于速度的平方，比例系数为0.4kg/m，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。 来加强对用数值解法解常微分方程实际问题的能力。 关键词：常微分方程数值解 MATLAB 目 录 摘 要 i 绪 论 1 问题解答…...………………………………………………………………2 总结……………………………………………………………….……...17 参考文献资料……………………………………..