The current multi-shaker random control methods concentrate on the simulation of stationary vibrations. This work presents an innovative technique for multi-shaker non-stationary random vibration environmental control test. Firstly, an amplitude modulation method is introduced to generate multi-input non-stationary random signals with time-varying root mean square values, whose kurtoses and moving RMS values are specified by amplitude modulation functions. Then, the generated signals are applied as excitation signals for the non-stationary random vibration environmental control test. The response signals of a linear vibration system under non-stationary excitations are derived in the frequency domain. The derived formulas for kurtoses and moving RMS values are also applicable to non-stationary response signals. Meanwhile, the stationary power spectral densities of response signals are put forward. When the amplitude modulation functions are zero mean Gaussian distributions, the power spectral densities of the response signals can be directly calculated. Ultimately, the control of kurtoses, moving RMS values and stationary power spectral densities of response non-stationary random signals is realized. Simulation and experimental examples both adopt a two-input two-output beam vibration system and the results prove the correctness and practicability of the proposed method. The presented control procedure is concise and efficient in simulating multi-shaker non-stationary random vibration environments.