数值分析常微分方程初值问题数值解法.docVIP

数学与计算科学学院 《数值分析》实验报告 目： 常微分方程初值问题数值解法 业： 信息与计算科学 学 号： 姓 名： 指导教师： 成 绩： 二零一五年六月五日 数值分析》实验报告 (7) 学号： 姓名： 班级： 成绩: 实验名称：常微分方程初值问题数值解法 实验地点：化工楼317 所使用的工具软件及环境：Matlab 一、实验任务与目的 用MATLAB实现欧拉法 用改进的欧拉法计算常微分方程的数值解 应用Runge ? Kutta方法解常微分方程 进一步熟悉MATLAB矩阵的应用。 二、实验内容或问题描述 通过给出儿个节点和函数值去求解最逼近原函数的方法，并且尽可能让误差达到最小。 三、解题思路与方法，或解题详细步骤 欧拉法 改进的欧拉法 四阶R-K法 四、程序设计代码 %% function y=oulafa(x) % y’=f(x，y) =y—2*x/y; y(0)=1; 0xl h=0. 05; n:l/h+l; x=0:h:1; y=zeros(l, n): %%应该给定y的大小，否则出错。 y ⑴=1; %欧拉法 for i=2:n y(i)=y(i-l)+h*(y(i-l)-2*x(i-l)/y(i-l)); end; size(x) size (y) % plot (x, y); % hold on; % plot (x, sqrt (1+2氺x)，’ r’)； %四阶龙格库塔法 z=zeros (1, n); z(l)=l; for i=2:n kl=y(i-l)-2*x(i-l)/y(i-l); k2=y(i-l)+h/2*kl-(2*x(i-1)+h)/(y(i-1) +h/2*kl); k3=y (i-l)+h/2*k2-(2*x(i-1)+h)/(y(i-1)+h/2*k2); k4=y(i-1)+h*k3-(2*x(i-l)+h)/(y(i-1)+h*k3); z(i)=z(i-1)+h/6*(kl+2*k2+2*k3+k4): end; % hold on; % plot (x, z, ’ in’)； yl=sqrt(l+2*x); plot (x, y, x, yl, ’ r’，x, z,’m’); legend C欧拉法’，’解析解’，’四阶龙格库塔法’)； [x;y;z;yl] 五、运行结果 ans = Columns 1 through 8 0 0.0500 0.1000 0.1500 0.2000 0.2500 0.3000 0.3500 1.0000 1.0500 1.0977 1.1435 1.1876 1.2301 1.2713 1.3113 1.0000 1.0492 1.0963 1.1415 1.1851 1.2273 1.2682 1.3079 1.0000 1.0488 1.0954 1.1402 1.1832 1.2247 1.2649 1.3038 Columns 9 through 16 0.40000.4500 0.4000 0.4500 0.5000 0.5500 1.3501 1.3880 1.4250 1.4612 1.3465 1.3842 1.4211 1.4571 1.3416 1.3784 1.4142 1.4491 0.6000 0.6500 0.7000 0.7500 1.4966 1.5313 1.5654 1.5990 1.4925 1.5272 1.5612 1.5948 1.4832 1.5166 1.5492 1.5811 Columns 17 through 21 0.8000 0.8500 0.9000 0.9500 1.0000 1.6320 1.6646 1.6968 1.7286 1.7600 1.6278 1.6605 1.6927 1.7245 1.7561 1.6125 1.6432 1.6733 1.7029 1.7321  六、心得体会 除常系数线性微分方程可用特征根法求解，少数特殊方程可用初等积分法求解外，大部 分微分方程无限lit界，应用中主要依靠数值解法。原欧拉法只有一阶，改进的欧拉法是两阶， 原欧拉法做出的图像明显可以看出顶点显著地偏离丫原來的积分曲线，可见欧拉法是相当粗 糙的，改进的欧拉法明显改善了精度。 数学与计算科学学院上机实验成绩评定表 组长：姓名，班级学号 专业 实验项目： 实验报告评阅、答辩记录： 成绩评定依据： 评定项目 自评成绩 评分成绩 1.实验项目合理、目的明确（10分） 2.实验方案正确，具有可行性、创新性（20分） 3.实验结果（例如：硬件成果、软件程序）（40分） 4.态度认真、学习刻苦、遵守纪律（10分） 5.报告 、总体整体结构及概貌（10分） 6.头验报告的规氾