Aiming at the problem that the weak fault features of rolling bearings are difficult to extract, a bearing fault feature extraction method based on the combination of parameter adaptive optimization variable modal decomposition (VMD) and multi-point optimal minimum entropy deconvolution (MOMEDA) is proposed. Firstly, the VMD decomposition is performed on the rolling bearing time domain vibration signal, and then the best mode component (BIMF) is selected based on the principle of maximizing the index of impulse harmonic noise ratio (AIHN) of autocorrelation function and MOMEDA filtering is performed on it, and the fault characteristic frequency is obtained after envelope deconvolution, and finally the fault characteristic frequency can be clearly observed by applying the proposed method body to the numerical simulation signal 131.1Hz, which can be applied to the actual bearing fault signal to effectively identify the bearing fault characteristic frequency of 294.5Hz, which is closer to the theoretical fault characteristic frequency 294Hz compared with 311Hz extracted by the original envelope spectrum and 320Hz extracted by MCKD.