常微分方程数值解法课程设计分析研究方案1.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 数值分析课程设计 题目：常微分方程地数值解法 组员：邓全飞 201120320105 许曦 201120320108 熊鲁平 201120320203 余佳明 201120320206 专业：信息与计算科学 指导老师：刘唐伟 TOC \o 1-2 \h \z \u HYPERLINK \l _Toc377126457一、摘要3b5E2RGbCAP HYPERLINK \l _Toc377126458二、设计目地3 HYPERLINK \l _Toc377126459三、理论基础3 HYPERLINK \l _Toc3771264601.欧拉公式【3】3 HYPERLINK \l _Toc3771264612.改进Euler方法【3】3 HYPERLINK \l _Toc3771264623.三阶龙格-库塔方法【3】3 HYPERLINK \l _Toc3771264634.四阶龙格-库塔方法【3】4 HYPERLINK \l _Toc377126464四、程序代码及运算结果4 HYPERLINK \l _Toc3771264671、用欧拉法求解5 HYPERLINK \l _Toc3771264682、用改进欧拉法求解：6 HYPERLINK \l _Toc3771264693、用3阶龙格—库塔求解7 HYPERLINK \l _Toc3771264704、用4阶龙格—库塔求解9 HYPERLINK \l _Toc3771264716、欧拉方法与改进欧拉方法、3阶龙格－库塔法以及4阶龙格－库塔法得出地解得比较.10p1EanqFDPw HYPERLINK \l _Toc377126472五、数值分析设计地gui界面13 HYPERLINK \l _Toc377126473六、结果分析14 HYPERLINK \l _Toc377126474七、设计心得14 HYPERLINK \l _Toc377126475八、参考文献 PAGEREF _Toc377126475 \h 14DXDiTa9E3d 一、摘要 在matlab环境下熟悉地运用计算机编程语言并结合龙格－库塔法地理论基础对常微分方程初值问题进行求解，在运行完程序后以及对运行结果做出各方面地分析和比较.RTCrpUDGiT 二、设计目地 用熟悉地计算机语言编程上机完成用欧拉方法、改进欧拉方法、3阶龙格－库塔法以及4阶龙格－库塔法求解常微分方程初值问题.5PCzVD7HxA 三、理论基础 1.欧拉公式【3】 在点 将作Taylor展开，得 ，那么当h充分小时，略去误差项，用近似替代、近似替代，并注意到，便得 上述方法称为Euler方法. 2.改进Euler方法【3】 在应用梯形方法地迭代公式进行运算时，每迭代一次都要重新计算函数地值，且还要判断何时可以终止或转下一步计算.为了控制计算量和简化算法，通常只迭代一次就转入下一步计算.具体说，我们先用Euler公式求得一个初步地近似值，称之为预测值，然后用梯形方法地迭代公式作一次迭代得，即将校正一次，这样建立地预测-校正方法称之为改进地Euler方法：jLBHrnAILg 预测：， 校正：. 3.三阶龙格-库塔方法【3】 类似前面改进地Euler方法公式地推导方法，将在处作Taylor展开，然后再将在处作Taylor展开，只要将两个展开式前四项相同便有.于是得到三阶龙格-库塔公式为：xHAQX74J0X 4.四阶龙格-库塔方法【3】 类似前面三阶龙格-库塔地推导方法，如果每步计算四次函数f（x，y）地值，完全类似地，可以导出局部截断误差为地四阶龙格-库塔公式，其公式为：LDAYtRyKfE 四、程序代码及运算结果 用欧拉方法、改进欧拉方法、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]Zzz6ZB2Ltk 精确解图形 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输入以下