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

数值分析课程设计 题目：常微分方程的数值解法 组员：邓全飞 5 许曦 8 熊鲁平 3 余佳明 6 专业：信息与计算科学 指导老师：刘唐伟 一、摘要 3 二、设计目的 3 三、理论基础 3 1.欧拉公式【3】 3 2.改进Euler方法【3】 3 3.三阶龙格-库塔方法【3】 3 4.四阶龙格-库塔方法【3】 4 四、程序代码及运算结果 4 1、用欧拉法求解 5 2、用改进欧拉法求解： 6 3、用3阶龙格—库塔求解 7 4、 用4阶龙格—库塔求解 9 6、欧拉方法与改进欧拉方法、3阶龙格－库塔法以及4阶龙格－库塔法得出的解得比较。 10 五、数值分析设计的gui界面 13 六、结果分析 14 七、设计心得 14 八、参考文献 14  一、摘要 在matlab环境下熟悉用熟悉的计算机语言编程上机完成 将作Taylor展开，得 ，那么当h充分小时，略去误差项，用近似替代、近似替代，并注意到，便得 上述方法称为Euler方法。 2.改进Euler方法【3】 在应用梯形方法的迭代公式进行运算时，每迭代一次都要重新计算函数的值，且还要判断何时可以终止或转下一步计算。为了控制计算量和简化算法，通常只迭代一次就转入下一步计算。具体说，我们先用Euler公式求得一个初步的近似值，称之为预测值，然后用梯形方法的迭代公式作一次迭代得，即将校正一次，这样建立的预测-校正方法称之为改进的Euler方法： 预测：， 校正：. 3.三阶龙格-库塔方法【3】 类似前面改进的Euler方法公式的推导方法，将在处作Taylor展开，然后再将在处作Taylor展开，只要将两个展开式前四项相同便有。于是得到三阶龙格-库塔公式为： 4.四阶龙格-库塔方法【3】 类似前面三阶龙格-库塔的推导方法，如果每步计算四次函数f（x，y）的值，完全类似的，可以导出局部截断误差为的四阶龙格-库塔公式，其公式为： 四、程序代码及运算结果 用欧拉方法、改进欧拉方法、3阶龙格－库塔法以及4阶龙格－库塔法求解常微分方程初值问题： 精确解为：y0=[ -2.0000 -1.9802 -1.9231 -1.8349 -1.7241 -1.6000 -1.4706 -1.3423 -1.2195 -1.1050 -1.0000] 精确解图形 1、用欧拉法求解 程序如下： 建立函数文件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^2*x ); [x,y]=cwfa1(fun,[0,1],-2,0.1); [x,y] plot(x,y,r+-) ans = 0 -2.0000 0.1000 -2.0000 0.2000 -1.9600 0.3000 -1.8832 0.4000 -1.7768 0.5000 -1.6505 0.6000 -1.5143 0.7000 -1.3767 0.8000 -1.2440 0.9000 -1.1202 1.0000 -1.0073 结果及其图象： 图1 欧拉法解 2、用改进欧拉法求解： 程序如下： 建立函数文件cwfa2.m function [x,y]=cwfa2(fun,x_span,y0,h) x=x_span(1):h:x_span(2); y(1)=y0; for n=1:length(x)-1 k1=feval(fun,x(n),y(n)); y(n+1)=y(n)+h*k1; k2=feval(fun,x(n+1),y(n+1)); y(n+1)=y(n)+h*(k1+k2)/2; end x=x;y=y; 在MATLAB输入以下程序： clear all fun=inline( y^2*x); [x,y