Rolling bearing fault diagnosis is crucial for ensuring the safe and reliable operation of mechanical equipment. Detecting faults directly from measurement signals is challenging due to severe noise and interference. Blind deconvolution (BD), as a preferred method, effectively recovers periodic pulses from the measured vibration signals of faulty bearings. This study introduces a simulated annealing-based BD approach to enhance the pulse signal components reflecting faults in vibration signals measured on rolling bearings. This method iteratively searches for the optimal coordinates in a high-dimensional orthogonal optimization space, where the optimal coordinates reflect the combination of the inverse filter coefficients. Compared to the generalized spherical optimization space used in the ‘Optimization-Blind Deconvolution’ method in previous works, the proposed finite high-dimensional optimization space helps overcome the problem of inverse filter coefficient convergence, allowing for the design of inverse filters without limit of its shape. To better accommodate the cyclostationarity characteristics of bearing signal measured in reality, the proposed method employs a target vector that allows for uncertainty in pulse occurrence instants, thus overcomes challenges introduced by pseudo-periodic phenomena resulting from bearing slippage. Numerical simulations and experimental results on real bearing vibration signals confirm that the proposed method can design more flexible filters to enhance pulse-like patterns in signals, effectively utilize limited filter resources. Its capacity to tolerate inaccurate fault period estimates, high background noise, and pulse randomness enables it to effectively address vibration measurement signals in real-world scenarios.