论文部分内容阅读
基于钟万勰等提出的指数矩阵精细算法，对n维未知向量v的一阶微分方程v=Hv+f(v,t)进行求解，其中Hv和f(v,t)分别是右端顶的线性齐次部分和非线性部分。将非线性部分在所论时刻tk处展成t-tk=τ的Taylor级数形式，并通过指数矩阵e^Ht及其精细算法对状态方程直接积分，推导出状态方程的级数形闭合解，此解的精度易于控制。算法不需对矩阵[H]求逆，数值计算的稳定性均可确保，对大型问题计算更为有利，算例验证了本算法的有效性。