Objective With the continuous increase in wireless communication traffic and the growing scarcity of spectrum resources, the overlapping frequency bands of communication and radar systems have led to mutual interference. Therefore, Integrated Sensing And Communication (ISAC) has emerged as a critical area of research. A key technology in this field is the design of ISAC waveforms, which is of significant research interest. The Doppler resilience and low sidelobe levels of ISAC waveforms are essential for effective target detection and information transmission in ISAC scenarios. However, designing waveforms that optimize both radar and communication performance presents substantial challenges. To address these challenges, a method for the joint design of ISAC waveforms and receiving filters is proposed. An ISAC model based on mismatched filtering is proposed, and an optimization problem is formulated. The Iterative Twisted appROXimation (ITROX) algorithm is presented to solve this nonconvex problem with guaranteed convergence. This approach enables the design of unimodular ISAC waveforms with Doppler resilience, achieving enhanced performance in both communication and radar functions. Methods To design ISAC waveforms that optimize radar and communication performance, the concept of mismatched filtering is introduced to formulate an optimization problem. The requirements for Doppler-resilient ISAC waveforms are first analyzed, followed by the proposal of a waveform model based on mismatched filtering. An optimization problem is then formulated, with the objective of minimizing the Weighted Integrated Sidelobe Level (WISL) and the Loss-in-Processing Gain (LPG). Constraints include the unimodular property of the transmitted waveform, the phase difference between the transmitted ISAC waveform and the communication data-modulated waveform, and the energy of the receiving filter. To solve this nonconvex optimization problem, the task is transformed into identifying a suitable Mismatched Filtering Sequences Pair (MFSP) under multiple constraints. An ISAC waveform design algorithm based on an improved ITROX framework is proposed to simplify the optimization process. The core concept of the ITROX algorithm is to iteratively search for the optimal projection of the matrix set, with the goal of maximizing the main lobe and minimizing the sidelobes within the region of interest. This approach minimizes WISL and LPG, satisfying the objective function requirements. Additionally, the combination of the three constraints ensures that the waveform meets both communication and radar sensing requirements. The SQUAREd Iterative Method (SQUAREM) is employed to improve the algorithm's convergence speed. The balance between WISL and LPG is controlled by adjusting the coefficients. Results and Discussions The ITROX-based ISAC waveform design method proposed in this paper effectively solves the formulated optimization problem, resulting in unimodular ISAC waveforms with Doppler resilience. Compared to existing ISAC waveform methods, the proposed ISAC waveform demonstrates a lower sidelobe level and Symbol Error Rate (SER) within the region of interest, with only a minor sacrifice in LPG. This leads to significant improvements in both radar sensing and communication performance. Simulation results validate the effectiveness of the proposed ISAC waveforms. These results show that the proposed method exhibits excellent convergence, with WISL rapidly converging to a stable value as iterations increase (Fig. 1). When the LPG coefficient is set to 0. 9, a low sidelobe level and SER are achieved, despite a processing gain loss of 0.91 dB (Fig. 2, Fig. 3). For the same phase difference threshold, the proposed ISAC waveform exhibits a lower SER than existing methods, indicating superior communication performance (Fig. 4). When comparing the ISAC waveform designed by this method to existing methods with the same time-delay interval width, the proposed waveform demonstrates a lower sidelobe level, with sidelobes nearly zero, approaching ideal correlation performance (Fig. 5, Fig. 6). This leads to significant improvements in target detection by the ISAC system. Furthermore, the proposed ISAC waveform exhibits excellent Doppler resilience, maintaining low sidelobe levels within the given Doppler interval (Fig. 6), which contributes to improved target detection performance. Conclusions This paper proposes a method for the joint design of ISAC waveforms and receiving filters based on Doppler resilience. By integrating the concept of mismatched filtering with the ISAC model, an optimization problem is formulated to minimize WISL and LPG without compromising communication quality. Additionally, an improved ITROX algorithm is proposed to effectively solve the formulated nonconvex optimization problem. The results demonstrate that the proposed scheme maintains near-ideal correlation performance within the region of interest under specified Doppler intervals, with only a minor sacrifice in LPG, and enables communication with a low SER. Compared to existing ISAC waveform methods, the proposed ISAC waveform exhibits a lower sidelobe level and SER, showing superior radar sensing and communication performance. Furthermore, low sidelobe levels can be achieved in one or more regions of interest to meet different requirements by appropriately adjusting the weighting coefficient. Future work could explore more efficient optimization algorithms to design ISAC waveforms with enhanced Doppler resilience.