常微分方程数值解法课程设计.docVIP

课程设计报告 课程设计题目： 常微分方程数值解法 学生姓名: 专 业： 班 级： 指导教师： 时 间： 题目：常微分方程数值解法 用欧拉方法、改进欧拉方法、3阶龙格－库塔法以及4阶龙格－库塔法求解常微分方程初值问题： 一、摘要 在matlab环境下熟悉用熟悉的计算机语言编程上机完成 将作Taylor展开，得 ，那么当h充分小时，略去误差项，用近似替代、近似替代，并注意到，便得 上述方法称为Euler方法。 2.改进Euler方法【3】 在应用梯形方法的迭代公式进行运算时，每迭代一次都要重新计算函数的值，且还要判断何时可以终止或转下一步计算。为了控制计算量和简化算法，通常只迭代一次就转入下一步计算。具体说，我们先用Euler公式求得一个初步的近似值，称之为预测值，然后用梯形方法的迭代公式作一次迭代得，即将校正一次，这样建立的预测-校正方法称之为改进的Euler方法： 预测：， 校正：. 3.三阶龙格-库塔方法【3】 类似前面改进的Euler方法公式的推导方法，将在处作Taylor展开，然后再将在处作Taylor展开，只要将两个展开式前四项相同便有。于是得到三阶龙格-库塔公式为：  4.四阶龙格-库塔方法【3】 类似前面三阶龙格-库塔的推导方法，如果每步计算四次函数f（x，y）的值，完全类似的，可以导出局部截断误差为的四阶龙格-库塔公式，其公式为： 三、程序代码及运算结果 用欧拉方法、改进欧拉方法、3阶龙格－库塔法以及4阶龙格－库塔法求解常微分方程初值问题： 求精确解 首先可以求得其精确解为：y=exp(-cos(x)) x=0:0.05:1; y=exp(-cos(x)); plot(x,y,b*-); Data=[x,y] Data = 0 0.3679 0.0500 0.3683 0.1000 0.3697 0.1500 0.3720 0.2000 0.3753 0.2500 0.3795 0.3000 0.3847 0.3500 0.3909 0.4000 0.3981 0.4500 0.4064 0.5000 0.4158 0.5500 0.4263 0.6000 0.4381 0.6500 0.4511 0.7000 0.4654 0.7500 0.4811 0.8000 0.4982 0.8500 0.5169 0.9000 0.5371 0.9500 0.5590 1.0000 0.5826 图 1精确解 2、 用欧拉法求解 程序如下： 建立函数文件cwfa1.m function [x,y]=cwfa1(fun,x_span,y0,h) x=x_span(1):h:x_span(2); y(1)=y0; for n=1:length(x)-1 y(n+1)=y(n)+h*feval(fun,x(n),y(n)); end x=x;y=y; 在MATLAB输入以下程序： clear all fun=inline( y*sinx ); [x,y]=cwfa1(fun,[0,1],exp(-1),0.05); [x,y] plot(x,y,r+-) ans = 0 0.3679 0.0500 0.3679 0.1000 0.3688 0.1500 0.3706 0.2000 0.3734 0.2500 0.3772 0.3000 0.3819 0.3500 0.3876 0.4000 0.3944 0.4500 0.4023 0.5000 0.4113 0.5500 0.4216 0.6000 0.4332 0.6500 0.4462 0.7000 0.4607 0.7500 0.4768 0.8000 0.4947 0.8500 0.5145 0.9000 0.5364 0.9500 0