数值分析实验报告之常微分方程数值解

《数值分析实验报告之常微分方程数值解》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数学与计算科学学院实验报告实验项目名称常微分方程数值解所属课程名称数值方法B实验类型验证实验日期2013.11.11班级学号姓名成绩16一、实验概述：【实验目的】1.掌握求解常微分方程的欧拉法；2.掌握求解常微分方程的预估校正法；3.掌握求解常微分方程的经典的四阶龙格库塔法；4.能用C语言或MATLAB将上述三种算法用程序运行出来；5.将算法实例化，并得出三种算法的相关关系，如收敛性、精度等；6.附带书中例题的源程序见附录1。【实验原理】1．欧拉格式（1）显式欧拉格式：局部截断误差：（2）隐式欧拉格式：局部截断误差：2．预估校正法预估：校

2、正：统一格式：平均化格式：3．四阶龙格库塔方法的格式（经典格式）16【实验环境】1.硬件环境：HPMicrosoft76481-640-8834005-23929HPCorporationIntel(R)Core(TM)I5-2400CPU@3.10GHz3.09GHz,3.16GB的内存2.软件环境：MicrosoftWindowsXPProfessional版本2002ServicePack3二、实验内容：【实验方案】方案一：用欧拉法，预估校正法，经典的四阶龙格库塔方法求解下列ODE问题：例题：在区间【0，1】上以h=0.1用欧拉法，

3、预估校正法,经典的四阶龙格库塔法求解微分方程dy/dx=-y+x+1，初值y（0）=1；其精确解为y=x+exp(-x),且将计算结果与精确解进行比较，对三个算法的收敛性的进行分析比较。方案二：用欧拉法，预估校正法,经典的四阶龙格库塔方法求解初值问题dy/dx=，初值y（0）=1；将计算结果与精确解为比较在区间[0,1]上分别取步长h=0.1;0.05时进行计算。对三个算法的收敛性进行分析比较，16【实验过程】（实验步骤、记录、数据、分析）注：以下图形是通过Excel表格处理数据得出，并未通过MATLAB编程序所得！1、由题可知精确解为：

4、，当x=0时，y(x)=0。h=0.1表1h=0.1时三个方法与精确值的真值表步长Euler法预估校正法经典四阶库精确值0.11.0100001.0050001.0048381.2490800.21.0290001.0190251.0187311.0554550.31.0561001.0412181.0408181.0912170.41.0904901.0708021.0703201.1318030.51.1314411.1070761.1065311.1768510.61.1782971.1494041.1488121.2260250.

5、71.2304671.1972111.1965861.2790160.81.2874211.2499751.2493291.3355360.91.3486781.3072281.3065701.3953221.01.4138111.3685411.3678801.458127图1h=0.1时三个方法走势图16h=0.05（此时将源程序中i的范围进行扩大，即for(i=0;i<20;i++)）表2h=0.05时三个方法与精确值的真值表步长Euler法预估校正法经典四阶库精确值0.051.0025001.0012501.0012291.011

6、7210.101.0073751.0048771.0048371.0249080.151.0145061.0107641.0107081.0395040.201.0237811.0188021.0187311.0554550.251.0350921.0288851.0288011.0727100.301.0483371.0409151.0408181.0912170.351.0634211.0547951.0546881.1109310.401.0802501.0704361.0703201.1318010.451.0987371.087

7、7521.0876281.1537910.501.1188001.1066621.1065311.1768510.551.1403601.1270871.1269501.2009420.601.1633421.1489541.1488121.2260250.651.1876751.1721931.1720461.2520620.701.2132911.1967361.1965851.2790160.751.2401271.2225201.2223671.3068520.801.2681211.2494851.2493291.3355360

8、.851.2972151.2775721.2774151.3650370.901.3273541.3067281.3065701.3953220.951.3584861.3369001.336