Rolling element bearings are frequently used in rotary machinery, but they are also fragile mechanical parts. Hence, exact condition monitoring and fault diagnosis for them plays an important role in ensuring machinery's reliable running. Timely diagnosis of early bearing faults is desirable, but the early fault signatures are easily submerged in noise.In this paper, Wigner–Ville spectrum based on cyclic spectral density (CSWVS for a brief notation) is studied, which is able to represent the cyclostationary signals while reducing the masking effect of additive stationary noise. Both simulations and experiments show that CSWVS is a noise resistant time frequency analysis technique for extracting bearing fault patterns, when bearing signals are under influences of random noise and gear vibrations. The 3-D feature of the CSWVS is proved useful in extracting bearing fault pattern from gearbox vibration signals, where bearing signals are affected by gear meshing vibration and noise. Besides, CSWVS utilizes the second order cyclostationary property of the vibration signals produced by bearing distributed fault, and clearly extracts its fault features, which cannot be extracted by envelope analysis.To quantitatively describe the extent of bearing fault, Renyi information encoded in the time frequency diagram of CSWVS is studied. It is shown to be a more sensitive index to reflect bearing performance degradation, compared with the spectral entropy (SE), squared envelope spectrum entropy (SESE) and Renyi informations for WVD, PWVD, especially when SNR is low.