Duffing system was linearized by equivalent linearization method based on energy difference

then

the random response of the linear system was solved by combining the pseudo excitation method and four-order symplectic geometry algorithm.It is concluded that an approach to combine the pseudo excitation method

symplectic geometry algorithm and equivalent linearization method is accurate and efficient.