matlab逐步积分,第17章 隐式逐步积分法.ppt

第十七章 MATLAB车辆工程应用实战 第17章 隐式逐步积分法 第十七章 MATLAB车辆工程应用实战 17.2 常用积分表达式 第十七章 MATLAB车辆工程应用实战 17.2 常用积分表达式 第十七章 MATLAB车辆工程应用实战 17.3 欧拉方法 欧拉公式如下： 欧拉方法的截断误差如下： 17.3.2 向后欧拉方法 向后欧拉公式又称为隐式欧拉公式，向后欧拉方法称为隐式欧拉方法，向后欧拉公式如下： for k=1:n X(k+1)=X(k)+h; Y(k+1)=Y(k)+h*(X(k)-Y(k)); k=k+1; end y=X-1+2*exp(-X); 第十七章 MATLAB车辆工程应用实战 17.3 欧拉方法 图17-1 数值计算结果 第十七章 MATLAB车辆工程应用实战 17.4 高阶微分方程(组)的数值计算 设有m阶常微分方程初值问题： 令 则初值问题就化为一阶常微分方程组： 代入相应的初值得： 第十七章 MATLAB车辆工程应用实战 17.4 高阶微分方程(组)的数值计算 function dz=dzdx3(x,z) dz(1)=z(2); dz(2)=5*(1- 2*z(1)^4)* z(2)-7*z(1); dz=[dz(1);dz(2)]; end 第十七章 MATLAB车辆工程应用实战 17.4 高阶微分方程(组)的数值计算 对于单摆运动求解过程，利用牛顿第二定律即得微分方程有： function xdot=danbai(t,x) g=9.8; l=25; xdot(1)=x(2); xdot(2)=-g/1*sin(x(1)); xdot=[xdot(1);xdot(2)]; end 第十七章 MATLAB车辆工程应用实战 17.4 高阶微分方程(组)的数值计算 图17-3 单摆运动数值计算 第十七章 MATLAB车辆工程应用实战 17.5 四阶龙格库塔法及亚当斯法 假设一个有阻尼的弹簧质量系统，动力学方程如下： 一阶微分方程为： %函数方程 function dy=lgkutta(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=-20*y(1)-4*y(2)-12*sin(2*t); 第十七章 MATLAB车辆工程应用实战 17.5 四阶龙格库塔法及亚当斯法 假设一个有阻尼的弹簧质量系统，动力学方程如下： 图17-4 弹簧质量系统输出响应图 第十七章 MATLAB车辆工程应用实战 17.6 隐式逐步积分法 17.6.1 线性加速度法 经过一个短暂的时间间隔 则上方程可以表示为： 图17-6 线性加速度示意图 第十七章 MATLAB车辆工程应用实战 17.6 隐式逐步积分法 17.6.1 线性加速度法 第十七章 MATLAB车辆工程应用实战 17.6 隐式逐步积分法 17.6.1 线性加速度法 对于一个振动系统来说，初始速度和初始位移一般是给定的，而初始时刻的加速度则根据动力学方程来得到。由上可见，线性加速度法是通过前一步的加速度、速度和位移，使动力学微分方程化为代数方程来求解的。 线性加速度法是一种有条件稳定的格式，如果积分步长 和系统的振动周期配合不理想，系统会发散。一般情况下，当系统的运动周期比积分步长大5倍以上时，才能取得比较精确的计算结果。 第十七章 MATLAB车辆工程应用实战 17.6 隐式逐步积分法 17.6.1 线性加速度法 例如一个系统的动力学方程为： 初始条件为 第十七章 MATLAB车辆工程应用实战 17.6 隐式逐步积分法 17.6.1 线性加速度法 第十七章 MATLAB车辆工程应用实战 17.6 隐式逐步积分法 17.6.2 威尔逊法 威尔逊 法通过建立在 时刻加速度的外插值 通过积分进 一步得到外插点的速度和位移，如图17-8所示，再通过 的动力方程求解出该时刻 的位移，进一步计算 时刻的位移、速度、加速度，该法是线性加速度法的一 种改进。 图17-8 加速度外插图 第十七章 MATLAB车辆工程应用实战 17.6 隐式逐步积分法 17.6.2 威尔逊法 得到标准节点的加速度、速度和位移的递推公式为： 时，威尔逊 法退化为线性加速度法，因