The operation of rolling bearings is directly related to the reliability of the whole rotating machinery. If the rolling bearing failure cannot be diagnosed and repaired in time, it will probably lead to equipment shutdown. The feature extraction process is an important part of bearing fault diagnosis. Some existing feature extraction methods are sensitive to noise and interference, and cannot fully explore and characterize useful information in nonlinear and non-stationary signals in complex scenes, and sometimes even lose important information contained in the data, resulting in low diagnostic accuracy. Therefore, a new method combining Maximum Second-order Cyclostationarity Blind Deconvolution (CYCBD) and Complementary Ensemble Empirical Mode Decomposition (CEEMD) was proposed for feature extraction of rolling bearing faults. Firstly, the cyclic frequency is set according to the failure frequency, and the original signal is filtered by CYCBD. Then, the filtered signal is decomposed into multiple Intrinsic Mode Function (IMF) by CEEMD, and the effective IMF components are selected by kurtosis criterion for reconstruction. Finally, the Teager energy operator is used to enhance the transient impact of the reconstructed signals. The results of simulation and comparison experiments show that the proposed method can extract bearing characteristics in different fault states more effectively, and the accuracy of network model diagnosis is improved compared with the traditional methods.