Abstract When a large amount of full life-cycle data are available, similarity-based methods are the preferred method for remaining useful life (RUL) prediction due to their reliability and accuracy. Traditional similarity-based RUL prediction methods use a single model and single-scale degradation features, which are incapable of fully capturing the degradation behavior of the system. Additionally, the similarity of spatial orientation is neglected in the similarity-matching process. To fill these research gaps, a novel method is developed based on multimodal degradation features and adjusted cosine similarity (ACS) to tackle complex-system RUL prediction in this paper. Complete ensemble empirical mode decomposition with adaptive noise is employed to decouple global degradation and random fluctuations in run-to-failure sensor data. Slow feature analysis is utilized to obtain local degradation features, and residual terms are used as global degradation features. Then, multimodal degradation features are transformed into one-dimensional health degradation indicators by bidirectional gated recurrent unit autoencoder. An ACS is developed to estimate the matching similarity between the test degradation curve and the training degradation curve. The proposed scheme captures the time-varying multimodal degradation behavior and provides libraries of health curves with multiple degradation patterns. The designed scheme is evaluated on the C-MAPSS dataset and the results illustrate the competitiveness and effectiveness of the proposed method.