![MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程_第1页](http://file3.renrendoc.com/fileroot_temp3/2022-2/8/22d1c422-e75f-41d2-a25e-e679d868098c/22d1c422-e75f-41d2-a25e-e679d868098c1.gif)
![MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程_第2页](http://file3.renrendoc.com/fileroot_temp3/2022-2/8/22d1c422-e75f-41d2-a25e-e679d868098c/22d1c422-e75f-41d2-a25e-e679d868098c2.gif)
![MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程_第3页](http://file3.renrendoc.com/fileroot_temp3/2022-2/8/22d1c422-e75f-41d2-a25e-e679d868098c/22d1c422-e75f-41d2-a25e-e679d868098c3.gif)
![MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程_第4页](http://file3.renrendoc.com/fileroot_temp3/2022-2/8/22d1c422-e75f-41d2-a25e-e679d868098c/22d1c422-e75f-41d2-a25e-e679d868098c4.gif)
![MATLAB改进欧拉法与四阶龙格-库塔求解一阶常微分方程_第5页](http://file3.renrendoc.com/fileroot_temp3/2022-2/8/22d1c422-e75f-41d2-a25e-e679d868098c/22d1c422-e75f-41d2-a25e-e679d868098c5.gif)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、姓名:樊元君 学号:2012200902 日期:2012.11.06 一、实验目的 掌握MATLAB语言、C/C+语言编写计算程序的方法、掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题。掌握使用MATLAB程序求解常微分方程问题的方法。二、实验内容1、分别写出改进欧拉法与四阶龙格-库塔求解的算法,编写程序上机调试出结果,要求所编程序适用于任何一阶常微分方程的数值解问题,即能解决这一类问题,而不是某一个问题。 实验中以下列数据验证程序的正确性。 求,步长h=0.25。2、实验注意事项的精确解为,通过调整步长,观察结果的精度的变化三、程序流程图: 改进欧拉格式流程图: 四阶龙格库塔流程
2、图:四、源程序: 改进后欧拉格式程序源代码:function = GJOL(h,x0,y0,X,Y)format longh=input(h=);x0=input(x0=);y0=input(y0=);disp(输入的范围是:);X=input(X=);Y=input(Y=);n=round(Y-X)/h); i=1;x1=0;yp=0;yc=0;for i=1:1:n x1=x0+h;yp=y0+h*(-x0*(y0)2);%yp=y0+h*(y0-2*x0/y0);%yc=y0+h*(-x1*(yp)2);%yc=y0+h*(yp-2*x1/yp);%y1=(yp+yc)/2;x0=x1;
3、y0=y1;y=2/(1+x02);%y=sqrt(1+2*x0);%fprintf(结果=%.3f,%.8f,%.8fn,x1,y1,y);end end 四阶龙格库塔程序源代码:function = LGKT(h,x0,y0,X,Y)format longh=input(h=);x0=input(x0=);y0=input(y0=);disp(输入的范围是:);X=input(X=);Y=input(Y=);n=round(Y-X)/h); i=1;x1=0;k1=0;k2=0;k3=0;k4=0;for i=1:1:n x1=x0+h;k1=-x0*y02;%k1=y0-2*x0/y0;
4、%k2=(-(x0+h/2)*(y0+h/2*k1)2);%k2=(y0+h/2*k1)-2*(x0+h/2)/(y0+h/2*k1);%k3=(-(x0+h/2)*(y0+h/2*k2)2);%k3=(y0+h/2*k2)-2*(x0+h/2)/(y0+h/2*k2);%k4=(-(x1)*(y0+h*k3)2);%k4=(y0+h*k3)-2*(x1)/(y0+h*k3);%y1=y0+h/6*(k1+2*k2+2*k3+k4);%y1=y0+h/6*(k1+2*k2+2*k3+k4);%x0=x1;y0=y1;y=2/(1+x02);%y=sqrt(1+2*x0);%fprintf(结果
5、=%.3f,%.7f,%.7fn,x1,y1,y);end end五、运行结果:改进欧拉格式结果:四阶龙格库塔结果:步长分别为:0.25和0.1时,不同结果显示验证了步长减少,对于精度的提高起到很大作用,有效数字位数明显增加。六、实验小结: 通过这次实验学习,首先第一点对改进欧拉格式和四阶龙格库塔的原理推导有了深入的理解,改进欧拉格式采用(预报+校正)模式得到较精确的原函数数值解;而四阶龙格库塔则采用多预报几个点的斜率值,采用加权平均作为平均斜率的近似值的思想达到更高精度的数值解,二阶龙格库塔的特例就是改进后的欧拉格式。七、思考题: 如何对四阶龙格-库塔法进行改进,以保证结果的精度。答:可以通过计算结果的精度处理步长来保证结果的精度。(1)步长折半。对于给定精
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024版产品三方购销合同
- 2024版安装音响工程合同
- (高清版)GBT 7383-2020 非离子表面活性剂 羟值的测定
- (高清版)GBT 39812-2021 塑料 试样的机加工制备
- 职业技术学院学生社团指导老师聘任管理办法
- 关于建筑公司合同范本
- 投资建设意向合同范本
- 新版租店合同范本
- 地产中价合同范本
- 市政车位买卖协议合同范本
- 中传媒2024春离线(附答案)《电视艺术概论》24年6月作业考核
- 2024年1月浙江杭州海关技术中心招考聘用笔试历年典型考题及考点研判与答案解析
- 水工钢筋混凝土结构(本)(闭卷)国开机考答案
- 大学生诚信考试题目及答案
- 2022-2023学年江西省九江市八年级(下)期末数学试卷(含解析)
- 在线网课学习智慧知道《运动控制系统(山东科大)》单元测试考核答案
- 2024年能源集团所属辽宁能源煤电产业股份限公司招聘公开引进高层次人才和急需紧缺人才笔试参考题库(共500题)答案详解版
- 运动损伤与康复智慧树知到期末考试答案章节答案2024年同济大学
- 六年级语文下册第四单元综合性学习奋斗的历程-《心愿》公开课一等奖创新教学设计
- 湖北省2023年中考地理真题(含答案)
- 2024年江苏扬州市高邮市水务产业投资集团有限公司招聘笔试参考题库含答案解析
评论
0/150
提交评论