Due to the variable working environment of bearings, the collected data often follow different probability distributions. It is hard to directly use the trained models to identify the bearing fault with different operating con-ditions. In addition, it is a high cost to label the samples for every work condition. To solve these problems, a multi-scale adversarial subdomain adap-tation bearing fault diagnosis method is proposed, which is based on Contin-uous Wavelet Transform (CWT) and our constructed Multi-scale Adversarial SubDomain Adaptation Network (MASDAN). Firstly, to extract the features of non-stationary signals with different frequency bands, CWT is used to con-vert continuous vibration signals into two-dimensional time-frequency images. Secondly, to enhance the correlation of the features across frequency bands, a multi-scale ConvNeXt is proposed, which adds a multi-scale module based on ConvNeXt to obtain features with different scales. Finally, to reduce the distribution discrepancy between the source domain and the target domain and to avoid feature confusion in different domains, a domain adaptive align-ment network introducing domain information is constructed. Two modules are included: the Domain Alignment Classification Network Module (DACNM) based on Multi-kernel Local Maximum Mean Discrepancy (MK-LMMD) and the Domain Adversarial Network Module (DANM) based on domain discrim-ination. Thus, the MASDAN consists of the multi-scale ConvNeXt module, DACNM, and DANM, which can realize the multi-scale feature extraction and adaptively align the fault features at different domains. The experimental re-sults on the Qingdao University of Technology (QUT) bearing dataset and the Case Western Reserve University (CWRU) bearing dataset demonstrate the proposed method can effectively diagnose bearing faults under different operating conditions.