原创力文档- 1、本文档共39页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
第9章微方程求解
浙江大学机械电子工程 第9章 微分方程求解 9.1 Introduction to Differential Equation 简单的一阶常微分方程初值问题 欧拉方法 后退（隐式）欧拉法 改进欧拉法(休恩方法） 9.4 Taylor Series Method 龙格-库塔方法 经典龙格-库塔法 通过MATLAB解常微分方程 Runge-Kutta 法: ODE解函数: ode23、 ode45、 ode113、 ode15s、 ode23s 参数选择函数： odeset、 odeget 输出函数： odeplot、odephas2 、 odephas3、 odeprint 例 y’=-2y+2x2+x其中（0X0.5),y(0)=1 fun.m function f=fun(x,y) f=-2*y+2*x^2+2*x; [x,y]=ode23(‘fun’, [0,0.5],1); [x,y]=ode23(fun, [0:0.1:1.0],1) 符号解法 dsolve(Df=f+sin(t),f(pi/2)=0) ans =-1/2*cos(t)-1/2*sin(t)+1/2*exp(t)/(cosh(1/2*pi)+sinh(1/2*pi)) dsolve(D2y=-a^2*y,y(0)=1,Dy(pi/a)=0) ans =cos(a*t) 思考题 1、采用最简洁格式的ODE文件和解算指令，研究围绕地球旋转的卫星轨道。 2、 This Class Is Over, Thanks for Your Attention! * * 9 Solution of Differential Equations 热传导原理？ a b 欧拉公式的截断误差与精度分析 这就是欧拉公式
文档评论（0）