- 1、本文档共19页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
二阶常微分方程matlab的数值解和解析解分析总报告精选
Matlab解二阶常微分方程 方程:a*y(t)+b*y(t)+c=sin(wt) 求解:1.解析解 2.数值解(欧拉方法) 目的:1.比较两种求解方式的拟合情况 2.通解随w变化的规律 * 方程求解程序清单 a=-1,b=2,c=-1; w=1; m=2; n=1; h = 0.02; t=0:h:30; s1=dsolve(a*D2y+b*Dy+c*y=sin(w*t),y(0)=m,Dy(0)=n,t); s1_n = eval(s1); hold on plot(t,s1_n,ko); EulerOED(a,b,c,w,m,n,h); hold off function EulerOED(a,b,c,w,x0,x1,h) A = [x0;x1]; t=0:h:30; for i = 1:1:length(t)-1 A(:,i+1) = [1,h;(-(c/a)*h),(1-(b/a)*h)]*A(:,i) + [0;(h/a)]*sin(w*t(i)); end plot(t,A(1,:),r*); * * 对于二阶全微分方程a*y(t)+b*y(t)+c=sin(wt) ,不同的a,b,c,w取值和初始条件会求出不同的解,通解又是由齐次解和特解组成。其中,齐次解由特征方程决定,而特解的决定因素则比较复杂。 讨论思路 (1)通解随初始条件变化情况 (2)通解随a,b,c变化情况 b^2-4ac0(两个不同的实根) b^2-4ac=0(两个相同的重根) b^2-4ac0(两个不同的复数根) 1).b0 2).b=0 3).b0 (3)通解随w变化情况 * (1)通解随初始条件变化情况 Ex1: a=2,b=3,c=1,y(0)=0;y(0)=0,w=1 * Ex2: a=2,b=3,c=1,y(0)=2;y(0)=0,w=1 * Ex3: a=2,b=3,c=1,y(0)=2;y(0)=4,w=1 * (2)通解随a,b,c变化情况 Ex1: a=2,b=3,c=1,y(0)=0;y(0)=0,w=1 * Ex4: a=-2,b=3,c=1,y(0)=0;y(0)=0,w=1 * Ex5: a=2,b=-3,c=1,y(0)=0;y(0)=0,w=1 * Ex6: a=2,b=3,c=-1,y(0)=2y(0)=1,w=1 ? * * b^2-4ac=0情况 * EX: a=2 ,b=2*sqrt(2) ,c=1,y(0)=0;y(0)=0,w=1 * * b^2-4ac0情况 * (3).b^2-4ac0 EX:a=4,b=-1,c=2,y(0)=0;y(0)=0,w=1 * EX:a=4,b=1,c=2,y(0)=3,y(0)=0,w=1 * EX:a=4,b=0,c=1,y(0)=2;y(0)=0,w=1 * * (3)通解随w变化的规律 W属于(0,1)时,随w的增大在齐次解的旁边波动 w属于(1,+),随w的增大逐渐趋近于齐次解。 * * * * *
您可能关注的文档
- 事业集团经营分析报告模板-月度精选.doc
- 事业集团经营分析报告模板-月度final精选.doc
- 事故污水处理池施工组织设计方案精选.doc
- 事故的预防措施和应急预案精选.doc
- 事故污水处理池施工组织设计精选.doc
- 二、ID卡食堂售饭机系统介绍精选.doc
- 二。门诊多学科会诊中心制定精选.doc
- 二上写字校本教材精选.doc
- 二乙胺安全技术说明书精选.doc
- 事务所房地产建筑企业法律服务方案精选.pdf
- 分析let s单元56ago2卷纸zheng unit56.pdf
- 塑胶材料其它分类原料pa9t 12.pdf
- md16x16数字媒体切换器设备.pdf
- 者参考项目发起人学科类型单位序列承包商修订页代码顺序典型.pdf
- 届世界天然气大会阿姆斯特丹2006add10288.pdf
- 期测试记录表每周weekly g1g6 journeys tests level 6 lesson26.pdf
- modernize-whitepaper现代化您应用程序白皮书.pdf
- anybackup产品典型案例分析.pdf
- 约克金融工程课程tfeslide32.pdf
- 广州市妇女儿童医疗中心历份教学药历01tjy.pdf
最近下载
- 华为云平台运维复习试题附答案.doc
- 2022-2023学年湖北省荆州市荆州区八年级(下)期末物理试卷(含解析).pdf VIP
- 十年(2013-2023年)高考生物真题分项汇编专题解析版合集共20套.doc
- 矫健:天局(胜天半子!).pdf
- (完整版)高中英语3500单词表含音标.pdf
- 150吨履带吊性能表(参考)1.doc VIP
- 湖北省荆州市荆州区2022-2023学年八年级下学期期末考试数学试卷(含答案).pdf VIP
- 人教版(新教材)高中英语选修2Unit4Journey Across a Vast Land精品学案:Writing—An email about a journey.docx
- 2012-2021十年高考生物真题(含解析)共92套.docx
- 2023年湖南师范大学附属小学教师招聘考试真题及答案.docx
文档评论(0)