![Matlab作业龙格库塔欧拉方法解二阶微分方程.ppt_第1页](http://file.renrendoc.com/FileRoot1/2019-11/19/a5095ac7-6bb1-40c8-ab94-0116cc87ea91/a5095ac7-6bb1-40c8-ab94-0116cc87ea911.gif)
![Matlab作业龙格库塔欧拉方法解二阶微分方程.ppt_第2页](http://file.renrendoc.com/FileRoot1/2019-11/19/a5095ac7-6bb1-40c8-ab94-0116cc87ea91/a5095ac7-6bb1-40c8-ab94-0116cc87ea912.gif)
![Matlab作业龙格库塔欧拉方法解二阶微分方程.ppt_第3页](http://file.renrendoc.com/FileRoot1/2019-11/19/a5095ac7-6bb1-40c8-ab94-0116cc87ea91/a5095ac7-6bb1-40c8-ab94-0116cc87ea913.gif)
![Matlab作业龙格库塔欧拉方法解二阶微分方程.ppt_第4页](http://file.renrendoc.com/FileRoot1/2019-11/19/a5095ac7-6bb1-40c8-ab94-0116cc87ea91/a5095ac7-6bb1-40c8-ab94-0116cc87ea914.gif)
![Matlab作业龙格库塔欧拉方法解二阶微分方程.ppt_第5页](http://file.renrendoc.com/FileRoot1/2019-11/19/a5095ac7-6bb1-40c8-ab94-0116cc87ea91/a5095ac7-6bb1-40c8-ab94-0116cc87ea915.gif)
已阅读5页,还剩2页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
Matlab应用使用Euler和Rungkutta方法解臂状摆的能量方程,1.背景,单摆求解单摆的运动一般使用角动量定理化简得到这样在小于5度的时候容易简化为,这样比较容易解。实际上这是一个解二阶常微分方程的问题。,2.问题,现在求解的是一个类似的问题,在这里的单摆是一种特别的单摆,具有均匀的质量M分布在长为2的臂状摆上。使用能量法(动能定理)建立方程化简得到(重力加速度取9.80665m/s2),计算,边值条件y(0)=0,y(0)=0.1.使用Euler方法精度随着h的减小而更高,因为向前欧拉方法的整体截断误差与h同阶,(因为用了泰勒公式)所以欧拉方法的稳定区域并不大。通过减小h增加了稳定性。,h=0.0001,h=0.01,计算,2.RK4-四阶龙格库塔方法使用四级四阶经典显式Rungkutta公式,误差很小:RK4法是四阶方法,每步的误差是h5阶,而总积累误差为h4阶。所以在同样步长h时候比欧拉方法准确。接下来进行对比,计算,运行第三个程序:在一幅图中显示欧拉法和RK4法,随着截断误差的积累,欧拉法产生了较大的误差h=0.01h=0.0001,总结,通过这两种方法计算出角度峰值y=3.141593,周期是1.777510。Euler方法结构简单,但是由于截断误差,使误差较大。RK4是很好的方法,很稳定,由于到五阶的时候精度并没有相应提升,所以四阶是很常用
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 【正版授权】 IEC 60512-20-3:2010 EN-FR Connectors for electronic equipment - Tests and measurements - Part 20-3: Fire hazard tests - Test 20c: Flammability,glow-wire
- 【正版授权】 IEC 60424-3:1999 EN-FR Ferrite cores - Guide on the limits of surface irregularities - Part 3: ETD-cores and E-cores
- 【正版授权】 IEC 60383-1:2023 EN-FR Insulators for overhead lines with a nominal voltage above 1000 V - Part 1: Ceramic or glass insulator units for a.c. systems - Definitions,test methods and ac
- 【正版授权】 IEC 60371-3-8:1995 EN-FR Insulating materials based on mica - Part 3: Specifications for individual materials - Sheet 8: Mica paper tapes for flame-resistant security cables
- 【正版授权】 IEC 60368-1:2000/AMD1:2004 EN-D Amendment 1 - Piezoelectric filters of assessed quality - Part 1: Generic specification
- 【正版授权】 IEC 60364-1:2005 EN-FR Low-voltage electrical installations - Part 1: Fundamental principles,assessment of general characteristics,definitions
- 【正版授权】 IEC 60357:2002/AMD1:2006 EN-FR Amendment 1 - Tungsten halogen lamps (non-vehicle) - Performance specifications
- 【正版授权】 IEC 60335-2-98:2023 RLV EN Household and similar electrical appliances - Safety - Part 2-98: Particular requirements for humidifiers
- 买房资金赠与合同范本
- 外贸中合同范本
- 科研伦理与学术规范期末考试1题库
- 2024年四川省自贡市中考数学试卷真题(含答案逐题解析)
- 2024年上海市普通高中学业水平等级性考试地理试卷(含答案解析)
- 酿酒工厂设计概论智慧树知到期末考试答案章节答案2024年齐鲁工业大学
- 2024年度保密教育线上培训考试题库及答案
- 中高端酒店销售方案(2篇)
- 物理与文化智慧树知到期末考试答案章节答案2024年山东大学
- 2024年江苏省南通张謇文化旅游景区招聘3人历年公开引进高层次人才和急需紧缺人才笔试参考题库(共500题)答案详解版
- 4.1 走向生态文明(教学设计)高二地理同步备课系列(人教版2019选择性必修3)
- 湖南省长沙市宁乡市2022-2023学年五年级下学期期末素养检测英语试卷
- 【公文基础】公文写作培训公文写作格式以及要点
评论
0/150
提交评论