《第9章微分方程求解》由会员分享,可在线阅读,更多相关《第9章微分方程求解(39页珍藏版)》请在金锄头文库上搜索。
1、第第9章章 微分方程求解微分方程求解9 Solution of Differential Equations1浙江大学机械电子工程9.1 Introduction to Differential Equation2浙江大学机械电子工程3浙江大学机械电子工程热传导原理?4浙江大学机械电子工程5浙江大学机械电子工程简单的一阶常微分方程初值问题简单的一阶常微分方程初值问题ab6浙江大学机械电子工程7浙江大学机械电子工程8浙江大学机械电子工程欧拉方法9浙江大学机械电子工程10浙江大学机械电子工程欧拉公式的截断误差与精度分析11浙江大学机械电子工程12浙江大学机械电子工程后退(隐式)欧拉法后退(隐式)欧
2、拉法13浙江大学机械电子工程改进欧拉法改进欧拉法(休恩方法)休恩方法)14浙江大学机械电子工程15浙江大学机械电子工程16浙江大学机械电子工程17浙江大学机械电子工程18浙江大学机械电子工程9.4 Taylor Series Method19浙江大学机械电子工程20浙江大学机械电子工程21浙江大学机械电子工程22浙江大学机械电子工程23浙江大学机械电子工程24浙江大学机械电子工程龙格龙格-库塔方法库塔方法这就是欧拉公式25浙江大学机械电子工程26浙江大学机械电子工程27浙江大学机械电子工程28浙江大学机械电子工程29浙江大学机械电子工程30浙江大学机械电子工程经典龙格经典龙格-库塔法库塔法31
3、浙江大学机械电子工程32浙江大学机械电子工程33浙江大学机械电子工程34浙江大学机械电子工程通过MATLAB解常微分方程lRunge-Kutta 法: ODE解函数: ode23、 ode45、 ode113、 ode15s、 ode23s 参数选择函数: odeset、 odeget 输出函数: odeplot、odephas2 、 odephas3、 odeprint35浙江大学机械电子工程例 y=-2y+2x2+x其中(0X0.5),y(0)=1fun.mfunction f=fun(x,y)f=-2*y+2*x2+2*x;x,y=ode23(fun, 0,0.5,1);x,y=ode2
4、3(fun, 0:0.1:1.0,1)36浙江大学机械电子工程符号解法ldsolve(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) ldsolve(D2y=-a2*y,y(0)=1,Dy(pi/a)=0)ans =cos(a*t)37浙江大学机械电子工程思考题思考题1、采用最简洁格式的、采用最简洁格式的ODE文件和解算文件和解算指令指令,研究围绕地球旋转的卫星轨道。,研究围绕地球旋转的卫星轨道。 2、38浙江大学机械电子工程This Class Is Over, Thanks for Your Attention!39浙江大学机械电子工程