The invention relates to the field of
computer simulation and scientific calculation, and provides a method for carrying out numerical solution on a dynamic
system described by a generalized
differential equation. According to the method, a
system is divided into two subsystems which interact through a port function, a zero-input independent solution and an
impulse response of the first subsystem are obtained in a computer,
time domain convolution is carried out on the zero-input independent solution and the
impulse response of the zero-input independent solution and the
impulse response of the zero-input independent solution and the impulse response of the first subsystem and port input of a second subsystem, and a cross-subsystem effect is decomposed into an independent solution item, a historical effect item and an instant effect item; and forming an equivalent equation of the second subsystem and executing numerical integration to obtain a
system time domain solution. The method is suitable for electrical, mechanical, hydraulic, control, structural, thermal, electromagnetic, fluid and other systems represented by generalized differential equations, and can be popularized to other
time evolution models. According to the method, repeated solution of a linear part is replaced by impulse response, and the matrix
operand is reduced while the spectral characteristics and precision are reserved, so that the
simulation efficiency is remarkably improved, and the stability is improved.