This　paper　is　aiming　to　synthesize　Non-stationary　Spatially　Variable　Ground　Motions　(NSVGMs).　After　analyzing　the　relationship　between　time-frequency　non-stationary　characteristics　and　phase　differences,　the　simulating　method　of　Non-stationary　Spatially　Variable　Ground　Motions　(NSVGMs)　is　proposed　in　this　paper.　In　the　simulation　progress,　amplitudes　and　related　phases　are　calculated　from　the　CHOLESKY　decomposition　of　the　power　spectrum　matrix,　and　the　initial　phase　is　obtained　from　the　phase　differences　spectrum.　Then　the　NSVGMs　are　synthesized　by　IFFT　method.　Based　on　the　proposed　method,　a　program　to　synthesize　NSVGMs　with　interactive　interface　is　compiled.　Finally,　a　series　NSVGMs　of　certain　spatial　sites　are　synthesized　by　this　program.　Results　of　the　computational　example　show　that　this　synthesis　program　can　effectively　provide　NSVGMs　for　multi-excitation　seismic　analysis.　©　2017,　Editorial　Department　of　Journal　of　Wuhan　University　of　Technology.　All　right　reserved.