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

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