《表3 经ode15s、精细逐块积分方法、RK4法和Newmark方法求得的q1值Tab.3 Numerical results of q1obtained by ode15s, PBIM, RK4a

《表3 经ode15s、精细逐块积分方法、RK4法和Newmark方法求得的q1值Tab.3 Numerical results of q1obtained by ode15s, PBIM, RK4a   提示:宽带有限、当前游客访问压缩模式
本系列图表出处文件名:随高清版一同展现
《非线性动力学系统的精细逐块积分求解方法》


  1. 获取 高清版本忘记账户?点击这里登录
  1. 下载图表忘记账户?点击这里登录

令ω=50,用精细逐块积分方法求解式(24)。设置绝对误差和相对误差均为10-14,将Matlab求解器ode15s的求解结果作为参考解。表3列出了经精细逐块积分方法、四阶Runge-Kutta法(RK4)和Newmark方法解出的q1值。可以看出,精细逐块积分方法在步长设置为0.001s时具有9或10位有效数字的精度;当步长增大为0.005s时,精细逐块积分方法仍然具有7位有效数字的精度。然而,在设置步长为0.001s时,RK4方法只有6位有效数字的精度;Newmark方法直到步长小至0.0001s时才有所收敛,并且只有1或2位有效数字的精度。