Most processes in nature are coupled; however, extensive null models for generating such processes still lacks. We present a new method to generate two coupled Gaussian stochastic processes with arbitrary correlation functions. This method is developed by modifying the Fourier filtering method. The robustness of this method is proved by generating two coupled fractional Brownian motions and extending its range of application to Gaussian random fields.