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.