Hybrid FSK-PSK Waveform Optimization for Radar Based on Alternating Direction Method of Multiplier (ADMM)

In this paper, a new radar signal modulated with a hybrid of the frequency shift keying (FSK) and the phase shift keying (PSK) signal—i.e., the FSK-PSK signal—is studied. Different phase encoding sequences are used to modulate the sub-pulses to obtain lower sidelobe levels and ensure signal orthogonality. In addition, to counter intra-pulse slice repeater jamming of specific length generated by the enemy jammer, an orthogonal waveform made of sub-pulses of equal length based on the FSK-PSK modulation scheme is designed. The simulation results show that the optimized discrete phase encoding sequence can significantly enhance the orthogonality of the sub-pulse in the FSK-PSK signal and effectively suppress the slice repeater jamming. Two algorithms are proposed: (1) the low sidelobe waveform optimization algorithm based on ADMM (LSW-ADMM); and (2) the anti-slice-repeater-jamming algorithm based on ADMM (ASRJ-ADMM). Both algorithms exhibit fast convergence speed and low computational complexity.


Introduction
In modern radar systems, the characteristics of radar waveform contributes greatly to various radar performance metrics, such as the detection probability, the resolution, the measurement accuracy, the probability of intercept, and anti-jamming capabilities. By optimizing the waveforms transmitted by radar, the radar detection and anti-jamming performance could be enhanced dramatically [1]. For example, the phased coded sequences set could be optimized for high peak-to-sidelobe level (PSL) and good quadrature performance [2,3].
Both phase shift keying (PSK) and frequency shift keying (FSK) are commonly used radar modulation schemes. Since the Doppler tolerance of the PSK signal is relatively low, it is not an optimum option for detecting target with high moving speed [4,5]. Compared with the PSK signal, the pulse compression ratio of the FSK signal is relatively high, which yields lower sidelobe levels. The hybrid modulation of FSK and PSK could be used to form a new class of modulated signal-i.e., FSK-PSK signal-which exhibits large time-bandwidth product and outstanding low probability-of-intercept (LPI) performance.
In recent years, extensive research has been carried out in the radar transmit waveform design area. In [6][7][8][9][10][11][12][13][14], in order to solve the problem of target masking caused by high sidelobes in the range-Doppler domain, different optimization algorithms have been designed for sidelobe suppression. Tao Fan [15] et al. designed a Doppler filter bank for non-uniform pulse repetition interval (PRI) radar to suppress clutters. For radar to coexist with other wireless systems, Khaled Alhujaili et al. optimized the transmit beam pattern for multiple-input multiple-output (MIMO) radar in [16]. Furthermore, Ehsan Raei et al. proposed a method for the joint design of the phase coding sequence and the Doppler filter at the receiver [17]. To minimize the side lobe level of the ambiguity function, the periodic phase sequence set is optimized in [18][19][20]. In [21], X. Wei et al. first established the slicing jamming model, then calculated the correlation coefficients to choose the appropriate location of the sampling window for jamming sampling, finally suppressed the slicing jamming with adaptive sidelobe cancellation. C. Zhou et al. designed a specific low-sidelobe phasecoded signal as the master signal, which was split into a number of sub-signals according to the modulation parameters of the repeater jamming. The mismatched filter was subsequently designed for the sub-signals via convex optimization to ensure the orthogonality between the different sub-signals [22]. In [23,24], an anti-jamming method was proposed, which was realized by generating the phase-coded (PC) signal and LFM-PC signal with high Doppler tolerance to improve the SAR imaging performance. S. Zhong [25] et al. designed an intra-pulse quadrature phase-encoded signal, and proposed a genetic algorithm based on the DNA coding to achieve low autocorrelation sidelobes and cross-correlations. In [26], S. Baher Safa Hanbali proposed an improved method based on both interrupted sampling repeater jamming (ISRJ) and frequency shifting modulation to protect the true target more effectively in the presence of amplitude mismatch between cancellation signal and true target echo. In [27], K. Zhou et al. proposed a method to suppress the ISRJ by jointly designing the radar waveform and mismatch filter.
The hybrid FSK-PSK modulated signal, which combines the advantages of the FSK signal and the PSK signal, exhibits good resolution, measurement accuracy, and antijamming performance. In recent years, it has been researched extensively in literature. In [28], Skinner et al. demonstrated the superiority of the Gaussian FSK signal modulated with pseudo-noise bipolar PSK sequence in the aspect of the range sidelobes. In [29], Donohoe et al. compared the ambiguity functions of various FSK signals, PSK signals and FSK-PSK signals modulated according to different coding schemes. The sidelobes of the ambiguity function associated with FSK-PSK signals are lower than those associated with the other two types of signals. In [30], each sub-pulse of the FSK-PSK signal is modulated with a pair of complementary Welti codes, while each sub-pulse performs frequency hopping according to the Costas array sequence. Compared with the Costas-coded FSK signal and the Welti-coded PSK signal, the FSK-PSK signal demonstrates lower ambiguity function sidelobes. In [31,32], Wong and Chung studied the characteristics of the FSK-PSK signal in many aspects-e.g., ambiguity function-and proposed a genetic algorithm (GA) to optimize the FSK-PSK signal. The simulation results show that the optimized FSK-PSK signal has lower autocorrelation function (AF) sidelobe levels and cross-correlations. In [33], J. Zhang et al. proposed a novel jamming detection and suppression method by exploiting the sub-pulse compression with hybrid FSK-PSK signals. Actually, limitations and conditions of the FSK-PSK signal exists in the aspects of ELINT (Electronic Intelligence) systems. [34] presents the procedure of radar source identification based on recognition of residually distinctive primal features using linear regression method (LRM) and concerns this problem.
The conventional FSK-PSK signals are usually modulated with intra-sub-pulse phasecoding and inter-sub-pulse frequency-coding-i.e., identical phase coding sequences are used for distinctive sub-pulses. To achieve low sidelobes and near-orthogonality, we propose a new FSK-PSK signal, which uses different phase-coding sequences to modulate each sub-pulse. To reduce the AF sidelobes while enhancing the anti-jamming capability of radar, two optimization problems are considered in this paper. The first is to optimize the autocorrelation function of the FSK-PSK signal so that low sidelobe levels are obtained. The second is to optimize the internal discrete phase encoding sequences of the sub-pulses so that the ultimate goal of anti-slice-repeater-jamming is realized.
To obtain an improved version of the FSK-PSK signal for enhanced detection performance and anti-jamming capability, this paper exploits the alternating direction method of multiplier (ADMM) algorithm [35]. The contribution of this work is summarized as follows: (1) The low sidelobe waveform (LSW) optimization algorithm based on ADMM, i.e., LSW-ADMM, is proposed. Distinctive phase encoding sequences are used for different subpulses, which are optimized with the LSW-ADMM algorithm for sidelobe reduction. Compared with the conventional FSK-PSK signals modulated by the classical code sequences, a higher-order of randomness is achieved, which guarantees enhanced anti-jamming capability of the radar system. (2) To encounter the intra-pulse slice repeater jamming, the ASRJ-ADMM algorithm is proposed. Each sub-pulse of the FSK-PSK signal, which is coded with a distinctive set of phases, performs frequency hopping according to the optimized sequence. By ensuring the orthogonality of the sub-pulses of the FSK-PSK signal, the intra-pulse slice repeater jamming is effectively suppressed.
Both optimization algorithms demonstrate fast convergence speed and low computational load. In numerical simulations, it is proved that the proposed algorithms outperform the widely adopted heuristic search algorithm, such as the GA algorithm.
The rest of this work is organized as follows. Section 2 presents the mathematical model of FSK-PSK signal, the AF of which is derived analytically. Section 3 proposes LSW-ADMM algorithm. Section 4 proposes the ASRJ-ADMM algorithm for countering intra-pulse slice repeater jamming. Section 5 verifies the outstanding performance of the proposed algorithm with extensive numerical simulations. Finally, concluding remarks and directions for future research are offered in Section 6.

FSK-PSK Hybrid Modulated Signal Model
The conventional FSK-PSK signal is modulated by intra-sub-pulse phase-coding and inter-sub-pulse frequency agility. A wide pulse is divided into several phase-coded subpulses that perform frequency-hopping according to a specific frequency-coding sequence.

Ambiguity Function of the FSK-PSK Hybrid Modulated Signal
To improve the waveform properties of the FSK-PSK signal and to ensure the nearorthogonality between the sub-pulses, we propose to use distinctive phase encoding sequences to modulate each sub-pulse. Figure 1 shows the partial structure of the FSK-PSK signal. In Figure 1, f m represents the frequency hopping code of the m-th sub-pulse. If the length of the phase-coding sequence of the sub-pulse is N, then the phase-coding sequence for the m-th sub-pulse is c m (n) = e j2πk m (n)/K , n = 0, 1, · · · , N − 1, where K represents the number of discrete phases.
Suppose that the total pulse width of the FSK-PSK hybrid modulated signal composed of M sub-pulses is T, it follows that the sub-pulse width is T c = T/M. Hence, the FSK-PSK hybrid modulated signal can be expressed as where the m-th sub-pulse signal v m (t) can be written as v m (t) = e j2π f m t In (2), f m = σ m /T P , m = 0, 1, · · · , M − 1 is the frequency hopping code sequence, σ m ∈ N + is the frequency hopping coefficient, c = c 0 T , c 1 T , · · · c M−1 T T represents the coding sequence of M sub-pulses, and p n (t) is the ideal rectangular pulse with a width of T P = T c /N. The spectrum of the phase-coded signal is obtained as and the spectrum of the FSK-PSK signal is given by The FSK-PSK signal bandwidth B FSK−PSK is approximately equal to If the frequency hopping coefficients σ m , m = 0, 1, · · · , M − 1, are a series of continuous positive integers, then we have The ambiguity function of the FSK-PSK signal can be obtained as: where m and l represent the sub-pulses corresponding to the m-th and the l-th frequency point, and (l − m − 1)T c ≤ τ ≤ (l − m + 1)T c . χ m,l (τ, ζ) in (7) could be further written as: where n and r represent the n-th and the r-th code, and χ m,n,l,r (τ, ζ) is given by, Given that |τ + (n − r)T P + (m − l)T c | ≤ T P , (9) can be further simplified as The amplitude of the main peak of the FSK-PSK signal in the Doppler dimension is Therefore, the Doppler tolerance of the FSK-PSK hybrid modulated signal is the same as that of FSK signal and PSK signal, which is expressed as

Autocorrelation Function of FSK-PSK Hybrid Modulated Signal
For the convenience of derivation, the FSK-PSK signal in (1)-(2) is rewritten as where f n = σ n /T P , n = 0, 1, · · · , N M − 1 represents the n-th hopping frequency point and is expressed as (13) represents the n-th phase code. The phasecoding sequence could be expressed as vector Suppose that N S points are sampled from u(t) in (13) with sampling period T S = T P /N S , and the total number of sampling points in a pulse is N = N M · N S . Equation (13) could be rewritten as where f n = σ n /T P , n = 0, 1, · · · , N − 1 represents the n-th hopping frequency point after sampling, and is expressed as (15) is the sampled phase-coding sequence given by Finally, p n (t) in (15) represents the ideal rectangular pulse with width T S . It follows that the autocorrelation function of the FSK-PSK signal is given by With some manipulation, (18) can be recast as Since T S is a constant, χ kT S is written as χ(k) in the following text for simplicity. We further define where J k is the shift matrix, and Finally, the compact form of the autocorrelation function of the FSK-PSK signal is obtained as

ADMM Algorithm
The ADMM algorithm is an effective way to solve separable convex optimization problems with multiple variables. By decomposing the original optimization problem into multiple easier-to-solve sub-problems and alternately solving them, the scale of the computational complexity can be reduced dramatically [35]. The effectiveness of this method in simplifying the complicated problem of radar waveform optimization has been demonstrated in [36][37][38][39].
The general ADMM framework is expressed as Equation (23) can be written as an augmented Lagrange function as where ρ is the coefficient for the penalty term. The updating equations for the ADMM algorithm are given by We set a small positive number δ as the convergence threshold, i.e., the algorithm

LSW-ADMM Optimization Method
In this subsection, the phase-coding sequence is optimized to minimize the AF sidelobes of the FSK-PSK signal. This integrated sidelobe level (ISL) is defined as To minimize the ISL in (26), the ADMM and the quasi-Newton method are employed to solve the convex optimization problem. The optimization problem is expressed as Introducing an auxiliary variable c = z to the constraint, we can obtain The augmented Lagrange equation corresponding to (28) is expressed as where λ r , λ i and ρ are Lagrange coefficients. Define u = (λ r + jλ i )/ρ, (29) can be recast as Define c (t) ,z (t) , u (t) as the values of c, z and u after t iterations. Given the initial value c (0) , z (0) , u (0) , then the solution of the problem can be carried out according to the following steps: (1) Given z (t) , u (t) , update c (t+1) as This problem is an unconstrained optimization problem, which can be solved by the quasi-Newton method [40]. First, we define Then c H H k c is expressed as (2) Given c (t+1) , u (t) , update z as Equation (33) could be recast as Since (34) is a linear equation and z(n), n = 0, 1, · · · , N − 1 are independent of each other, (34) can be decomposed into k(n) sub-problems as With some manipulation, the objective function in (35) could be written as Ignoring the constant coefficient, Equation (36) is simplified as To minimize the angle between the vectors z(n) = e j2πk(n)/K and c (t+1) (n) + u (t) (n) in the complex plane, the integer k(n)(0 ≤ k(n) < K) is obtained as where · stands for the floor function. The solution to (34) is obtained as (3) Given u (t) , c (t+1) , z (t+1) , update u (t+1) as The procedures of LSW-ADMM algorithm are as follows. To begin, initialize c,z, and u as c (0) , z (0) , u (0) . Next, repeat step 1) to step 3) until the number of iterations reaches the upper limit-i.e., ISL (t+1) − ISL (t) ≤ ε-with ε being the convergence threshold (or until the iteration stop condition is satisfied). The implementation steps of the LSW-ADMM algorithm are given in Table 1.

LSW-ADMM
(a) Set t as the number of iterations. Initialize the algorithm with t = 0 and set c (0) , z (0) , u (0) ; (b) Using quasi Newton method to solve the problem (e) t ← t + 1; (f) If the convergence condition is satisfied, the algorithm is completed; otherwise, returns to (b)

Intra-Pulse Slice Repeater Model
The signal structure of the intra-pulse slice repeater jamming generated by the enemy jammer is shown in Figure 2. The jammer intercepts part of the radar signal in the time domain according to the sampling period to generate the intercepted signal, which can be regarded as multiplying the radar signal with an ideal rectangular pulse with a sampling period of T s and a width of τ 0 to obtain an intra-pulse slice signal. The slice signal is processed by power amplification and delayed-forwarding to generate an intra-pulse slice repeater jamming signal. After the radar receives the intra-pulse slice repeater jamming, multiple deception targets are formed due to the coherence of the repeater jamming.

LSW-ADMM
(a) Set t as the number of iterations. Initialize the algorithm with 0 t = and set (b) Using quasi Newton method to solve the problem If the convergence condition is satisfied, the algorithm is completed; otherw turns to (b)

Intra-Pulse Slice Repeater Model
The signal structure of the intra-pulse slice repeater jamming generated by th jammer is shown in Figure 2. The jammer intercepts part of the radar signal in domain according to the sampling period to generate the intercepted signal, whic regarded as multiplying the radar signal with an ideal rectangular pulse with a s period of s T and a width of 0  to obtain an intra-pulse slice signal. The slice processed by power amplification and delayed-forwarding to generate an intra-pu repeater jamming signal. After the radar receives the intra-pulse slice repeater ja multiple deception targets are formed due to the coherence of the repeater jamm For the specific slice repeater jamming whose slice signal width is equal to th of the FSK-PSK intra-sub-pulse signal, the intra-sub-pulse orthogonal waveform used to suppress the repeater jamming [41][42][43].  For the specific slice repeater jamming whose slice signal width is equal to the width of the FSK-PSK intra-sub-pulse signal, the intra-sub-pulse orthogonal waveform can be used to suppress the repeater jamming [41][42][43].

Intra-Pulse Slice Repeater Jamming Suppression Model
Suppose N s points are sampled with a sampling period of T s = T p /N s , and the total number of sampling points is N = N · N s for a FSK-PSK signal made of a single sub-pulse. Equation (2) is rewritten as where p n (t) represents an ideal rectangular pulse with a width of T S , and c m represents the phase-coding sequence of the m-th sampled sub-pulse signal, which is given by According to the definition of the ambiguity function, it follows from (41) that By setting τ = kT S , −(N − 1) ≤ k ≤ N − 1 , (43) can be recast as Since T p is a constant, χ ml kT p is simplified as χ ml (k) in the following text. We further define H ml,k = ξ ml,k · J k · D ml (45) The cross-correlation of the sub-pulses of the FSK-PSK signal is then obtained in the discrete form as when m = l, (48) represents the discrete-form of the AF of the sub-pulse signal To ensure the near-orthogonality between the sub-pulses of the FSK-PSK signal, optimization problem can be formulated to minimize the AF sidelobes and the weighted sum of the cross-correlations between the sub-pulses. Assume that the frequency-hopping coding sequence is fixed, we seek to optimize the discrete phase-coding sequence. The objective function of the optimization problem is formulated as , and α ∈ [0, 1] is the weighting coefficient. It could be easily observed from (50) that when α = 1, the sole purpose of (50) is to minimize the AF sidelobes. While for α= 0, the sole purpose of (50) switches to the minimization of the cross-correlations.

ASRJ-ADMM Optimization Method
The ADMM architecture and the quasi-Newton method are used to minimize the weighted sum of the AF sidelobes and the cross-correlations of the sub-pulses of the FSK-PSK signal. The optimization problem is formulated as c m (n) = e j2πk m (n)/K n = 0, 1, · · · , N − 1; Equation (51) can be recast as Introducing an auxiliary variable c = z to the constraint, we can obtain Note that the problem in (52) is equivalent to that in (53). The augmented Lagrange equation corresponding to (53) is expressed as where λ r , λ i , and ρ are Lagrange coefficients. Let u = (λ r + jλ i )/ρ, then (54) can be written as The implementation procedures of the ASRJ-ADMM algorithm are summarized in Table 2. To begin, initialize c,z, and u as c (0) , z (0) , u (0) . Next, repeat steps b) to e) until the number of iterations reaches the upper limit-i.e., δ (t+1) − δ (t) ≤ ε-with ε being a small positive number (or until the iteration stop condition is satisfied).

ASRJ-ADMM
(a) Set t as the number of iterations. Initialize the algorithm with t = 0 and set c (0) , z (0) , u (0) ; (b) Using quasi-Newton method to solve the problem (e) t ← t + 1; (f) If the convergence condition is satisfied, the algorithm is completed. Otherwise, returns to (b)

Computational Complexity Analysis
In the ADMM algorithm, (31) and (53) dominate the main computation. We need to calculate the gradient of (31) and (53). Considering the subpulse length N , we can suppose Ψ contains αN 2 points, where 0 < α < 1. Then, the computation of gradient contains times of real number multiplications, and times of real number additions. Therefore, the computational complexity of ADMM algorithm is on the order of o(N DT N 5 ).
We note that GA algorithm can also be successfully applied in discrete-phase sequence design problems. However, with increased variables in FSK-PSK optimization problem, the efficiency of GA algorithm is not satisfactory.
In each iteration of GA algorithm, the number of variables in c is N , and for each variable c n , we need to calculate the number αN 2 locations in the DAF. The problem of c n in [31], has the computational complexity of o(N DT N 3 ).The total computational complexity of GA algorithm is o(N DT N 6 ). Therefore, the computational complexity of GA algorithm is higher than that of the ADMM algorithm.

Numerical Experiments
In this section, we provide two simulation examples to demonstrate the performance of the two optimization algorithms proposed in Sections 3 and 4. The first example shows that the LSW-ADMM algorithm proposed in this paper has faster convergence speed and lower computational complexity compared with the GA algorithm. Moreover, the FSK-PSK signal optimized with the LSW-ADMM algorithm exhibits lower AF sidelobes, which in turn contributes to better detection performance. The second example shows that, compared with the conventional FSK-PSK signal modulated with the classic coding sequence, the FSK-PSK signal optimized with the ASRJ-ADMM algorithm possesses better anti-jamming capability due to a higher-order of randomness of the radar signal. Moreover, unlike the GA algorithm, which suppresses the AF sidelobes at the cost of higher crosscorrelations, the ASRJ-ADMM algorithm can realize joint optimization of the AF and the cross-correlation. The FSK-PSK signal optimized with the ASRJ-ADMM algorithm guarantees near-orthogonality between sub-pulses, hence has the potential to suppress the repeater jamming of specific intra-pulse slices effectively.
Assuming that T = 2us, T P = 0.1us, and σ m = {1, 3, 2, 4}, The phase encoding sequence is a random two-phase code. The PSK, FSK, and hybrid FSK-PSK are plotted in Figure 3. The FSK-PSK signal ambiguity function is shown in Figure 4. The FSK-PSK signal ambiguity function is a pushpin ambiguity function, and the FSK-PSK signal has no range-Doppler coupling problem. signal. Moreover, unlike the GA algorithm, which suppresses the AF sidelobes at the cost of higher cross-correlations, the ASRJ-ADMM algorithm can realize joint optimization of the AF and the cross-correlation. The FSK-PSK signal optimized with the ASRJ-ADMM algorithm guarantees near-orthogonality between sub-pulses, hence has the potential to suppress the repeater jamming of specific intra-pulse slices effectively.
Assuming that The phase encoding sequence is a random two-phase code. The PSK, FSK, and hybrid FSK-PSK are plotted in Figure 3. The FSK-PSK signal ambiguity function is shown in Figure 4. The FSK-PSK signal ambiguity function is a pushpin ambiguity function, and the FSK-PSK signal has no range-Doppler coupling problem.
(a) (b) (c)  Let τ = 0 and ξ = 0, the zero-delay cut and the zero-Doppler cut of the 2D ambiguity function shown in Figure 5 could be obtained, which are depicted in Figure 5a,b, respectively. The velocity ambiguity function, main peak amplitude, and Doppler tolerance of FSK-PSK signals are all equivalent to PSK signals.

Optimization of the AF of the FSK-PSK Signal
According to the optimization model in (21), we define its average integrated sidelobe level (AISL) as where Ψ = k N S ≤ |k| ≤ N − 1 represents the sidelobe range of the AF of the FSK-PSK signal.
Assume that the number of sub-pulses is M = 4, the frequency hopping coefficient is σ = {σ m , m = 0, · · · , M − 1} = {1, 2, 3, 4}, the sub-pulse signal phase-coding length is N = 50, and the discrete phase number is K = 2/4/8. We also assume that the phasecoding sequence of the PSK signal is identical to that of the FSK-PSK signal. The AISL of the FSK-PSK signal and the PSK signal are plotted in in Figure 6. It can be seen from Figure 6 that the AISL of the FSK-PSK signal is significantly lower than that of the PSK signal. The results are summarized in Table 3.  Next, we set the maximum number of the iterations of the optimization algorithm to 50, then use the LSW-ADMM and the GA algorithm to optimize the same initial FSK-PSK signal. The AFs before and after optimization when K takes on different values (K = 2, K = 4, K = 8) are shown in Figure 7. The black lines represent the initial waveforms before optimization, the red lines in Figure 7a,c,e represents the AFs of the FSK-PSK signals optimized with the LSW-ADDM, and the red lines in Figure 7b,d,f represents the AFs of the FSK-PSK signals optimized with the GA. It can be seen clearly from Figure 6 that although the AF sidelobes of the FSK-PSK signal become lower by applying both optimization algorithms, the LSW-ADMM algorithm yields much lower AF side-lobes than that of the GA algorithm. Table 4 summarizes the AISL before and after waveform optimization with the LSW-ADMM algorithm and the GA algorithm. As can be seen from Table 4, the LSW-ADMM algorithm offers an AISL 3-5 dB lower than that of the GA algorithm. Specifically, with the increase of K from 2 to 8, the advantage of the LSW-ADMM algorithm becomes more obvious. However, as the number of discrete phases K further increases, the optimization effect of LSW-ADMM algorithm may not improve accordingly due to the limitation set by the coefficients for the penalty term in the Lagrange expression (see Section 3.2, though a detailed analysis of this phenomenon is beyond the scope of this work).  Figure 8 shows the convergence speed of the LSW-ADMM algorithm and the GA algorithm. It could be seen that the LSW-ADMM algorithm converges in less than 10 iterations, while the GA algorithm takes much more iterations to converge. As the number of discrete phase increases, the number of iterations of the LSW-ADMM algorithm requires to converge will slightly increase. Hence, we could conclude that the optimization effect of the GA algorithm is strictly constrained by the number of optimization iterations that are allowed within a specified period of time.

Joint Optimization of Intra-Pulse Sub-Pulse Correlation Function
In this example, we set the number of FSK-PSK signal sub-pulses as M = 4, the frequency hopping coefficient as σ = {σ m , m = 0, · · · , M − 1} = {1, 2, 3, 4}, the phasecoding length of the sub-pulse as N = 50, the number of discrete phases K = 2, the weighting coefficient of the optimization model as α = {0, 0.1, · · · , 1}, and the upper limit of the optimization algorithm iteration number as 50. Identical initial waveforms are employed to test the waveform optimization performance of the ASRJ-ADMM algorithm and the GA algorithm. The auto-correlation average sidelobe level (AC-ASL) and cross-correlation average level (CC-AL) of the radar signals before and after optimization are shown in Figure 8. The AC-ASL and the CC-AL are defined, respectively, as It can be seen from Figure 9a that as the weighting coefficient increases, the AC-ASL of the FSK-PSK signals optimized with the ASRJ-ADMM algorithm and the GA algorithm decreases gradually, while the AC-ASL of the FSK-PSK signal untouched by the waveform optimization algorithms remains constant. It could also be seen that the AC-ASL offered by the GA algorithm is slightly lower than that of the ASRJ-ADMM algorithm by 2-3 dB. However, this marginal performance gain of the GA algorithm over the ASRJ-ADMM algorithm is at the cost of higher CC-AL. As shown in Figure 9b, although the CC-AL of the proposed ASRJ-ADMM algorithm increases with the weighting coefficient, its maximum value is approximately equal to that of the FSK-PSK signal untouched by either optimization algorithms, and is 5-6 dB lower than that of the FSK-PSK signal optimized with the GA algorithm.  Figures 10 and 11 show that, after the ADMM-based optimization process, the AF sidelobes and cross-correlations of the FSK-PSK signal are suppressed. Table 5 compares the AC-ASL and CC-AL before and after the ASRJ-ADMM optimization. It can be seen from Table 5 that both the AC-ASL and CC-AL are suppressed with the ASRJ-ADMM algorithm. However, based on the preliminary simulation results that are not presented in this work, as the number of discrete phases K further increases, the optimization effect may not improve accordingly.   By comparing the iteration convergence curve of the optimization algorithm in Figure 8 and the optimization algorithm time in Table 6, we can find that-compared with GA algorithm-LSW-ADMM and MSSD-ADMM algorithms have fewer iterations, faster convergence time, and better optimization effect.

Conclusions
In this work, two novel waveform optimization algorithms based on ADMM are proposed within the general FSK-PSK hybrid modulation framework for AF sidelobe suppression and waveform orthogonality enhancement. Specifically, the LSW-ADMM waveform optimization algorithm proposed in this work yields radar signals with AF sidelobes significantly lower than those of the conventional FSK-PSK signals. Meanwhile, the ASRJ-ADMM algorithm generates radar signals with outstanding anti-slice-repeaterjamming capability. In the follow-up works, we plan to further improve the detection performance and the anti-jamming capability of the FSK-PSK signal by jointly optimizing the frequency-hopping coding sequence and phase coding sequence.

Data Availability Statement:
The datasets generated during the current study are not publicly available but are available from the corresponding author on reasonable request.