数值计算与法设计课程设计.docVIP

课程设计说明书 课程名称: 数值计算与算法设计课程设计 题 目: 导弹追踪微分方程模型的数值解法 院 系： 理学院 ＿ 专业班级：＿应用数学200200513794 ＿ 学生姓名：＿ ＿储素霞＿＿ 指导教师：＿＿ 许 峰＿＿＿ 2008 年 月 日 学 号 2005姓名 专业班级 应用数学200-2 设计题目 任务起至时间 2008年7月日至2008年7月1日 设 计 要 求 及 任 务 总 述 。 工 作 计 划 及 安 排 月日 分析 月日月日 算法设计月日月日 月日 参 考 资 料 . [2] 李庆扬, 王能超, 易大义. 数值分析, 华中科大出版社, 武汉, 2005. [3] 王高雄, 周之铭. 常微分方程, 高等教育出版社, 北京, 2006. [4] 何青, 王丽芬. Maple教程, 科学出版社, 北京, 2006. 指导教师签字 系主任签字 2008年月日 学生姓名： 储素霞 学号： 200513794 专业班级： 应用数学200-2 课程设计题目： 导弹追踪微分方程模型的数值解法 指导教师评语： 成绩： 指导教师： 2008年月1 日 一、问题及分析 1 二、模型的建立 1 三、算法的设计 3 1、解析方法 3 2、数值方法 4 3、仿真方法 5 四、程序及结果分析 8 1、解析方法 8 2、数值方法 8 3、仿真方法 10 设计总结 13 参考文献 14 导弹追踪问题及其求解 一、问题及分析 某军导弹基地发现正北方向120km处海面上有敌舰一艘以90km/h的速度向正东方向行驶。该导弹基地立即发射一枚导弹跟踪追击敌舰，导弹速度为450km/h，自动导航系统使导弹在任一时刻都能对准敌舰。试问导弹在何时何处击中敌舰？ 二、模型的建立 如图建立坐标系，取导弹基地为原点，轴指向正东方，轴指向正北方。 当时，导弹位于点，敌舰位于点，其中。设导弹在时刻的位置为，由题意 (1) 其中。 在时刻，敌舰位于处，其中。由于导弹轨迹的切线方向必须指向敌舰，即直线PM的方向就是导弹轨迹上点P的切线方向，故有 或 (2) 方程(1)、(2)连同初值条件 (3) 构成了一个关于时间变量的一阶常微分方程组的初值问题。 为了获得与的关系，要设法消去变量，由(2)式得 两边对求导得 ，即 将上式与(1)式合并，再加上初值条件，则得如下初值问题 这就是导弹轨迹的数学模型。 三、算法的设计 1、解析方法 模型中的二阶方程可以降阶。令 ，记， 则方程可降为一阶可分离变量方程 即 易得 由初值条件，得，从而 注意到上式可改写为 于是有 这样我们又得到一个可分离变量方程 积分得 利用，得，从而导弹轨迹方程为 设导弹击中敌舰于，以代入上式，得 击中敌舰的时刻为 代入具体数据得 2、数值方法 可以用数值分析中介绍的Euler公式、改进的Euler公式和四阶Runge-Kutta公式来求解上述初值问题。 (1) Euler公式 将在上积分，，得，用数值积分法求。 (1)，得 Euler公式 (2)，得 后退的Euler公式 (3)， 得 梯形公式（隐式） (2) 改进的Euler公式 Euler公式计算简便，但精度差，梯形公式为隐式，计算较复杂，但精度较高，可将两者结合。 ， 称为改进的Euler公式，上式也可写为 。 (3) 四阶Runge-Kutta方法 三阶Runge-Kutta方法较少使用，仿二阶Runge-Kutta方法，可得四阶Runge-Kutta公式，经典的四阶Runge-Kutta公式为 。 特点：①单步、自开始；②精度高，误差为，四阶；③数值稳定；④要计算四次函数值；⑤对解的光滑性要求高。 3、仿真方法 如果建立微分方程很困难，或者微分方程很复杂而难以做出数值处理，常常可以用仿真方法。 所谓仿真方法，顾名思义，指的是模仿真实事件行为和过程的方法。在这个具体问题中，就是一步步地模拟导弹追踪敌舰的实际过程。而计算机仿真，则是在计算机上通过相应的程序和软件来实现对事件运行的