Mechanomyography (MMG) signals have extensive applications in muscle function assessment and human intention recognition. However, during signal acquisition, MMG signals are easily contaminated by noise and artifacts, which seriously affects the recognition of their characteristics. To address these issues, a novel noise suppression and artifact removal method based on recursive least square (RLS), improved Gray Wolf Optimizer-optimized variable mode decomposition (IGWO-VMD), and complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) is proposed. In this paper, the RLS algorithm is first applied to adaptively filter out the power line interference (PLI). Then, IGWO is designed to select the appropriate VMD parameters and use the VMD to decompose the noisy signal into band-limited intrinsic mode functions (BLIMFs). In addition, the BLIMFs are classified into the low-frequency part and high-frequency part according to the given correlation coefficient (CC) threshold value. The effective components of the low-frequency part are identified by the center frequency. Meanwhile, the high-frequency part is decomposed by CEEMDAN, and its effective components are obtained according to the proposed sample entropy threshold range. Finally, the effective components of the low and high-frequency parts are reconstructed to obtain the denoised signal to realize the extraction of useful signals. Simulation experiment results demonstrate that the proposed method outperforms the classical methods and the designed IGWO-VMD method in terms of denoising performance. The effectiveness of the proposed method is verified through the measured MMG signal experiments. The proposed method not only effectively suppresses noise and artifacts but also overcomes the limitations of VMD and CCEMDAN.