Hidden Markov Models (HMMs) have become an immensely popular tool for health assessment and fault diagnosis of rolling element bearings. The advantages of an HMM include its simplicity, robustness, and interpretability, while the generalization capability of the model still needs to be enhanced. The Dempster-Shafer theory of evidence can be used to conduct a comprehensive evaluation, and Stacking provides a novel training strategy. Therefore, the HMM-based fusion method and ensemble learning method are proposed to increase the credibility of quantitative analysis and optimize classifiers respectively. Firstly, vibration signals captured from bearings are decomposed into intrinsic mode functions (IMFs) using ensemble empirical mode decomposition (EEMD), and then the Hilbert envelope spectra of main components are obtained; Secondly, multi-domain features are extracted as model input from preprocessed signals; Finally, HMM-based intelligent health assessment framework and fault diagnosis framework are established. In this work, the life cycle health assessment modeling is performed using a few training samples, the bearing degradation state is quantitatively evaluated, normal and abnormal samples are effectively distinguished, and the accuracy of fault diagnosis is significantly improved.