常微分方程的数值求解及其应用程序.doc

数值分析设计报告 常微分方程的数值求解及其应用 一、设计目的 2 二、任务简介 2 三、理论基础 3 3.1二阶龙格-库塔方法 3 3.2三阶龙格-库塔方法 4 3.3四阶龙格-库塔方法 5 3.4二阶亚当姆斯格式 6 3.5三阶亚当姆斯格式 7 3.6四阶亚当姆斯格式 8 3.7亚当姆斯预报-校正系统 9 四、案例的运算结果 10 4.1案例一运算结果 10 4.2案例二运算结果 11 五、数值分析设计的GUI界面 12 六、结果分析 13 七、设计心得 13 八、附录 13 一、设计目的 在matlab环境下熟悉用熟悉的计算机语言编程上机完成在区间上满足条件：时，的特解，画出数值解的图像，进行比较。 案例二：放射性废物的处理 美国原子能委员会以往处理浓缩的放射性废料的方法，一直是把它们装入密封的圆桶里，然后扔到水深为90多米的海底。生态学家和科学家们表示担心，怕圆桶下沉到海底时与海底碰撞而发生破裂，从而造成核污染。原子能委员会分辨说这是不可能的。为此工程师们进行了碰撞实验，发现当圆桶下沉速度超过12.2m/s与海底相撞时，圆桶就可能发生碰裂。这样为避免圆桶碰裂，需要计算一下圆桶沉到海底时速度是多少？这时已知圆桶重量为239.46kg，体积为0.2058m3， 海水密度为1035.71 kg /m3。如果圆桶速度小于12.2m/s，就说明这种方法是安全可靠的，否则就要禁止用这种方法来处理放射性废料。假设水的阻力与速度大小成正比例，其正比例常数k=0.6。 （1） 建立解决上述问题的微分方程数学模型。 （2） 用数值和解析两种方法求解微分方程，并回答是否要禁止用这种方法来处理放射性废料。 三、理论基础 3.1二阶龙格-库塔方法 只要成立，满足这一条件的格式统称为二阶龙格-库塔格式。在本题中 当时，该式为改进的欧拉格式。 当时，该式为变形的欧拉格式。 二阶龙格库塔-算法流程图： = = 图3.1 二阶经典的龙格-库塔方法 3.2三阶龙格-库塔方法 类似前面改进的欧拉方法公式的推导，将在处作Taylor展开，然后再将在处作Taylor展开，只要将两个展开式前四项相同便有。于是得到三阶龙格-库塔公式为： 三阶龙格库塔-算法流程图： = 图3.2 三阶经典的龙格-库塔方法 3.3四阶龙格-库塔方法 类似前面三阶龙格-库塔的推导方法，如果每步计算四次函数f（x，y）的值，完全类似的，可以导出局部截断误差为的四阶龙格-库塔公式，其公式为： 四阶龙格库塔-算法流程图： = 图3.3四阶经典的龙格-库塔方法 3.4二阶亚当姆斯格式 为使计算格式具有二阶精度，须取，这样可推出二阶亚当慕斯格式为： 类似可推出三阶亚当姆斯格式和四阶亚当姆斯格式。 二阶亚当姆斯方法流程图： = 图3.4 二阶亚当姆斯预方法流程图 3.5三阶亚当姆斯格式 三阶亚当姆斯方法流程图： = 图3.5 三阶亚当姆斯方法流程图 3.6四阶亚当姆斯格式 四阶亚当姆斯方法流程图： = 图3.6 四阶亚当姆斯方法程图 3.7亚当姆斯预报-校正系统 预报： 校正： 亚当姆斯预报-校正系统流程图： = 图3.7 亚当姆斯预报校正系统流程图 四、案例的运算结果 4.1案例一运算结果 二、三、四阶龙格-库塔方法求解方程组所得图