Implicit algorithms for numerical simulation of Hamiltonian systems with external noise, which are symplectic in the no-noise limit, are constructed. The methods were tested on the Kubo oscillator with white noise and for long times of integration gave better results than ordinary stochastic algorithms.