非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程v=Hv+f(v,t),结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的vk+i/s(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算eHt,而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。