Abstract: In order to ensure that the vibrations of a system under test are in accordance with a reference spectral density Z(w), the system is driven by a signal x.sub.n (t) provided by an inverse discrete Fourier transform (IDFT) with random phases introduced by a function R.sub.n (w) and the response spectral density Y.sub.n (w) provided by a DFT is brought into agreement with Z(w) by a feedback loop. Moreover, in order to ensure that this is done in one iteration, X'.sub.n (w) is derived from x.sub.n (t) by a DFT and divided by Y.sub.n (w) and the quotient is multiplied with Z(w) to provide the new input x.sub.n (w) for the IDFT. The division cancels out the random factor which is present in both X'.sub.n (w) and Y.sub.n (w).The system thus equalizes both rapidly and stably.