版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 食品业合同范例
- 小吃类转让合同范例
- 规范性合同范例
- 广场商铺购买合同范例
- 大量物资采购合同范例
- 人防产品合同范例
- 进口原料招商加盟合同范例
- 人事用工合同范例
- 经纪中介服务合同范例
- 求购采砂船转让合同范例
- 2023年全国统一高考化学试卷(辽宁卷)含答案与解析
- 临床医学检验:抗微生物药物和敏感性试药试题及答案真题
- 2024年中国铁道科学研究院集团有限公司招聘笔试参考题库附带答案详解
- 天津市和平区2023-2024学年高一年级上册1月期末英语试题(解析版)
- 组合机床动力滑台液压系统设计
- 子宫腺肌瘤的个案护理
- 管道护理小组工作总结
- 零售商店人员培训:零售业的法律与合规问
- 《21世纪学生发展核心素养研究》
- 2023版道德与法治教案教学设计专题3第1讲 理想信念的内涵及重要性
- 山东省济南市济钢高级中学2024届高一上数学期末调研模拟试题含解析
评论
0/150
提交评论