matlab第二章节常微分方程的数值解法.pptVIP

微分方程的数值解法;常微分方程（Ordinary differential equations, ODE);一.解ODE的基本机理:;例：著名的Van der Pol方程;3. 根据式(2.2)编写计算导数的M函数文件-ODE文件;4. 使编写好的ODE函数文件和初值 供微分方程解算指令（solver）调用;例题1：著名的Van der Pol方程;Van der Pol方程;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;解法指令;二. 四 阶 Runge-Kutta 法;初值问题的数值解法分为两大类;四阶Runge-Kutta公式;四 阶 Runge-Kutta 法计算流程图;三. Runge-Kutta 法解Van der Pol 方程的Matlab 程序结构主程序：RK_vanderpol.m 子程序：RK_sub.m（函数文件） ;解法2：采用Runge_Kutta法编程计算;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;子程序：RK_sub.m function ydot = vdpol (t, y) ydot=zeros(size(y)); ydot(1) = y(2); ydot(2) = -y(2)*(y(1)^2-1)-y(1);  或写为： ydot = [y(1) ;-y(2)*(y(1)^2-1)-y(1)];;四. Matlab对应命令：ode23,ode45;3月15日作业: 1.Van der Pol 方程的两种解法:1)采用ode45命令 2)Runge-Kutta方法 2.Duffing 方程的求解（Runge-Kutta方法，计算步长h=0.005，计算时间t0=0.0,tN=100） 要求：写出程序体，打印所绘图形，图形标题用个人的名字。;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;五. 动力学系统的求解;其中：;3. Matlab 程序(主程序：ZCX);Matlab 程序(子程序：ZCX_sub.m);例题2：三自由度质量弹簧系统;矩阵表示;动力学方程：;4阶龙格－库塔法的结果;例题3： 蹦极跳系统的动态仿真;整个蹦极系统的数学模型为：;令：; y0=[-30; 0]; % 初始位移和初始速度 [t,y]=ode45(‘bengji_sub’, [0:0.01:100], y0); x1=50. - y(:,1); % x1代表蹦极者与地面之间的距离 plot(t,x1); grid plot(t,y(:,1)); grid % y(:,1)代表位移;function ydot=f(t,y) m=70; k=20; a1=1; a2=1; g=10; x=y(1); % x代表蹦极者的位移 x_dot=y(2); % x_dot 代表 x 的速度 if x0 ydot=[0,1;-k/m,-a1/m-(a2/m)*abs(x_dot)]*y+[0;g]; else ydot=[0,1;0,-a1/m-(a2/m)*abs(x_dot)]*y+[0;g]; end;y(:,1);作业（书面作业，写出程序体）： （1）三自由度模型仿真 (自编Runge-Kutta 法） （2）蹦极跳模型仿真 (解算指令ode45)