Optimizing an Algorithm Designed for Sparse-Frequency Waveforms for use in Airborne Radars

: Low-frequency bands are an important way to realize stealth target detection for airborne radars. However, in a complex electromagnetic environment; when low-frequency airborne radar operates over land, it will inevitably encounter a lot of unintentional communication and intentional interference, while e ﬀ ective suppression of interference can not be achieved only through the adaptive processing of the receiver. To solve this problem, this paper proposes optimizing an algorithm designed for sparse-frequency waveforms for use in airborne radars. The algorithm establishes a joint objective function based on the criteria of minimizing waveform energy in the spectrum stop-band and minimizing the integrated sidelobe level of speci ﬁ ed range cells. The waveform is optimized by a cyclic iterative algorithm based on the Fast Fourier Transform (FFT) operation. It can ensure the frequency domain stopband constraint to realize the e ﬀ ective suppression of main-lobe interference while forming lower-range sidelobes at speci ﬁ ed range cells to improve the ability to detect dim targets. Theoretical analysis and simulation results have shown that the algorithm has good anti-interference performance.


Introduction
The appearance of stealth aircraft, such as the F-22 and F-35, has brought serious challenges to airborne early warning radar systems that take conventional targets as detection objects. Low-frequency bands are an internationally recognized effective means of stealth detection, yet the low-frequency band spectrum is very crowded, and there are a large number of radiation sources, which are especially dense over land [1], resulting in the widespread existence of main lobe interference in low-frequency-band airborne early warning radars. Thus, effective suppression of dense main lobe interference can not be achieved only through the adaptive processing of the receiver.
Cognitive radars can obtain information by sensing the battlefield environment and feeding it back to the transmitter, thereby changing the traditional one-way processing mode of the radar and forming a new architecture of a dynamic closed loop between the receiver, the transmitter, and the battlefield environment [2][3][4][5]. The processing information in the different stages flows into the closed loop, and finally achieves the goal of enhancing the performance of the radar system [6][7][8]. The cycle of information begins with the transmitter broadcasting signals to illuminate the surrounding environment. After the radar receives the signal and the signal enters the receiver, the scene analyzer extracts the target (SMIA). In the frequency domain, the objective function is established through an optimization criterion of minimum wave energy in the frequency stopband. Then, a new autocorrelation function sequence is formed by the relationship between the waveform autocorrelation function and power spectrum density, and the problem of integrated sidelobe level minimization of a specified range cells is converted into a power spectrum density fitting problem. Then, the weighting factor is employed to form a new objective function. Finally, the waveform is optimized by the circular iterative algorithm based on an FFT operation. Compared with traditional waveform-optimization design algorithms, the stopband depth of the proposed algorithm can be adjusted while keeping the stopband constraint in the frequency domain, and lower range sidelobes can be formed at specified range cells. Meanwhile, the algorithm is optimized based on an FFT operation, and has relatively lower computational complexity.

Materials
It is assumed that the discrete signal form of the radar-transmitted waveform after N-point sampling is The autocorrelation function of waveform s is defined as represents the sidelobe of the waveform autocorrelation function, and the autocorrelation function of Equation (2) refers to the aperiodic autocorrelation function.

Construction of Spectrum Stopband Constrained Objective Function
Assume that the radar operating frequency band is   For the construction of the objective function of spectrum stopband constraint, the main purpose is to minimize the waveform energy in the frequency stopband so as to effectively suppress the interference from the same frequency band, that is, to minimize the waveform energy in In order to be consistent with the sequence length when constructing the low-sidelobe-constrained objective function of specified range cells in Section 2.2 below, here, let , that is, a 2N point PSD of the transmitted waveform s is calculated by a 2N point DFT.    (5) where F is a 2N × 2N FFT transformation matrix.
To obtain the power spectrum density with the minimum waveform energy in the frequency stopband, the objective function can be constructed by introducing the weighting coefficient, where the weighted vector w can be obtained by prior information [36].

Construction of Low-Sidelobe-Constrained Objective Function of Specified-Range Cells
There is another mathematical model of the minimum waveform energy in the spectrum stopband: Since the autocorrelation function sequence of waveform and waveform PSD are Fourier transform pairs, a new autocorrelation function sequence   m s can be constructed [37] (See Appendix A for details).
 m s s s s s s s (11) Assuming that the ideal autocorrelation function sequence is s is the ideal waveform, and the minimum mean square error criterion is used to make the final optimized waveform s to the ideal waveform  s , then the objective function is It is assumed that the relevant information of the ideal waveform  s can be obtained from prior information, that is, the low sidelobe information of specified-range cells is known. Then, the new autocorrelation function is processed in the way of 0-1 weighting, 0 k z  for the region that needs to be suppressed, and 1 k z  for the region that does not need to be suppressed, where k z represents the value of the k-th cell.
Therefore, the ideal autocorrelation function sequence    m s can be expressed as where s w is the weight vector in the same order as the new autocorrelation function se- Convert Equation (10) to matrix form as follows: a y a y a y a y a y a y y y (15) where F is a 2N × 2N FFT transformation matrix, and y is where B is According to Equation (15) where Fy is the spectrum of the actual waveform,  g is the spectrum of the ideal waveform, and the objective function of Equation (18) can be expressed as the spectrum fitting function, namely Introducing phase assist vector Equation (21) can be converted to is the amplitude of the ideal spectrum,  represents the square root of the absolute value of each element in the vector, and Combining Equation (9) and Equation (23), a weighting factor  is introduced to form a new objective function: It should be noted that matrix F , matrix B , matrix T , vector  w , and vector h are known, and the objective function contains s , u , v , three variables, and the optimal solution of the third variable can be obtained by fixing two variables and then iteratively solving it continuously.

Optimization Algorithm under Constant Modulus Constraint
In order to maximize the utilization of transmitter power and avoid nonlinear distortion of the output waveform, constant modulus constraints are usually added in the waveform design process [38], where the objective function becomes The main idea of solving the objective function is to assume that s and u are known and minimize Equation (25) to obtain the optimal solution of v ; then, assuming that s and v are known, the optimal solution of u is obtained. Then, u and v are updated, the optimal solution of s is obtained, and the iterative cycle continues until the preset stopping condition ( ) s is the sequence obtained by the i-th iteration and  is the preset threshold.
When s and u are known, the optimal solution of v is When s and v are known, the optimal solution of u is When u and v are given, the objective function can be written as where const is the constant term independent of s , 1 c is the first N elements of Then, the optimal solution of s is Continuously iterate until convergence.
The specific steps of the SMIA algorithm under constant modulus constraint are provided as follows (Algorithm 1).

Algorithm 1: SMIA algorithm under constant modulus constraint
Input: randomly initialize the sequence s ; Step 1: For the current sequence s and u , calculate the optimal solution of v (see Equation (26)); Step 2: For the current sequence s and v , calculate the optimal solution of u (see Equation (27)); Step 3: For the current sequence u and v , calculate the optimal solution of s (see Equation (29)); Step 4: Iterate through Step 1, Step 2, and Step 3 until the preset stop condition

Optimization Algorithm under PAR Constraint
From the point of view of maximization of transmitter power, the ideal transmission waveform can be obtained under constant modulus constraint, but the constraint is too harsh. The Peak-to-Average-power Ratio (PAR) constraint is a more general constraint than the constant modulus constraint, and its expression is The transmitted waveform meets the energy constraint, and Similarly, the SMIA algorithm is used to solve the problem. When s and u are known, the optimal solution of v is When s and v are known, the optimal solution of u is When u and v are known, according to reference [39], the Tropp alternate projection method is used to solve the problem. For Step 1: For the current sequence s and u , calculate the optimal solution of v (see Equation (32)); Step 2: For the current sequence s and v , calculate the optimal solution of u (see Equation (33)); Step 3: For the current sequence u and v , calculate the optimal solution of s by the alternate projection method; Step 4: Iterate through Step 1, Step 2, and Step 3 until the preset stop condition Output: the final optimized waveform.

Results
In this section, it is assumed that the waveform length N = 250; without loss of generality, the frequency is normalized, and its range is 0~1, the stop threshold of iterative convergence is set to 10 −6 , and the maximum number of iterations is 10 4 . If no special description is made, the spectrum uses a 0-1 weighted vector. In this section, the proposed SMIA algorithm is compared with the traditional WeSCAN algorithm under the constant modulus constraint and PAR constraint so as to verify the effectiveness of the proposed SMIA algorithm in reducing the sidelobes of specified-range cells and imposing spectrum constraint.

Comparisons under Constant Modulus Constraint
The autocorrelation function and the power spectral density obtained by the SMIA algorithm when 0.9  = are shown in Figure 1; the frequency stopband is [0.2 0.3) ， and the region of reducing the sidelobes of specified range cells is The image obtained by the WeSCAN algorithm is shown in Figure 2, and other conditions are the same as those in Figure 1. By comparing Figures 1 and 2, the average PSD notch depth of the SMIA algorithm is −30.12 dB, and the average PSD notch depth of the WeSCAN algorithm is −32.23 dB. The mean sidelobe of the specified-range cells autocorrelation function of the SMIA algorithm is −61.32 dB, and the mean sidelobe of the WeSCAN algorithm is −43.25 dB. It can be found that, when 0.9  = , that is, the stopband constraint weight is large, both the SMIA algorithm and WeSCAN algorithm can form a notch at the stopband, and the notch depth is basically the same; but, for the autocorrelation function sidelobe level of specified range cells, the SMIA algorithm can generate a lower range sidelobe than WeSCAN algorithm, with a difference of about 18 dB.

Comparisons under PAR Constraint
The autocorrelation function and power spectrum density of the SMIA algorithm and WeSCAN algorithm under the PAR constraint are shown in Figures 3 and 4, respectively, where PAR ≤ 2 and other conditions are the same as in Section 4.1. It can be found that the results obtained by PAR constraint are consistent with those obtained by constant modulus constraint. However, due to the relaxation of the peak average ratio by PAR constraint, the degree of freedom of waveform design increases. Therefore, compared with the constant modulus constraint, the range sidelobe level obtained was relatively lower (about 20 dB lower), while the PSD spectrum constraint was less affected.

Anti-Interference Performance
In this experiment, the region of reducing the sidelobes of specified-range cells is (3, 250) k  and the weighting factor is 0.9   . In Figure 5, the transmitted waveform , the interference-tonoise ratio is 15 dB, and the PSDs of the transmitted waveform and the interference waveform are shown in Figure 5a,b, respectively. Figure 5c is a PSD diagram of the echo waveform, including the transmitted waveform, interference waveform, and noise. Figure 5d shows the output of matched filtering between the echo waveform and the transmission waveform. It can be seen that the interference frequency band does not fall within the frequency stopband of the transmitted waveform, so the interference is not effectively suppressed after matched filtering. The simulation conditions in Figure 6 are essentially the same as those in Figure 5. The difference is that the frequency band of the interference waveform is (0.22 0.28) ， ; this frequency band is in the frequency stopband of the transmitted waveform. Figure 6d shows the output after matched filtering. It can be found that the interference band falls within the frequency stopband of the transmitted waveform, and the interference is effectively suppressed after matched filtering because of the orthogonality of the frequency domain.

The Time-Frequency Distribution Characteristics of the Optimized Waveform
In order to clarify the time-frequency distribution characteristics of the waveform more clearly, this experiment assumes that the sequence length N = 1000. The frequency stopband is [0.08, 0. 16 Figure 7. It can be seen from the figure that the spectrum power of the final optimized waveform is evenly distributed in the passband, and corresponding notches are formed in the stopband, so the optimized waveform has a good lowintercept performance.

Comparisons between Sparse-Frequency Waveform and Non-sparse-Frequency Waveform
The autocorrelation function of the sparse-frequency waveform and the non-sparsefrequency waveform is shown in Figure 8, where the red line is the sparse-frequency waveform; the frequency stopband is [0.15, 0.25) U [0.4, 0.55) U [0.7, 0.8), which means the passband is 0.65 and the stopband is 0.35. The blue line is a non-sparse-frequency waveform, and its bandwidth is consistent with the passband size of the sparse-frequency waveform. Combined with Figure 8a,b, it can be found that, compared with the nonsparse-frequency waveform, the sparse-frequency waveform has a narrower main-lobe width, but the range sidelobes of the sparse-frequency waveform are higher. Therefore, it is necessary to comprehensively consider the anti-interference performance, range resolution, and the influence of range-strong isolated point clutter when optimizing the transmitted waveform design under the background of airborne radar.

Weighting Factor
In this section, we analyze the influence of the weighting factor  on the SMIA algorithm through two performance indices: peak stopband power and peak autocorrelation sidelobe level. Assuming sequence length N = 100, the region of reducing the sidelobes of specified range cells is (2,50) k  and the frequency stopband is [0.2 0.3) ， . The peak autocorrelation sidelobe level is defined as Peak stopband power is defined as (35) is obtained from the optimal sequence   1 The graph of corr P and stop P is changing with  from 0.1 to 1, as shown in Figure   9. As can be seen from the figure, as  increases, more weight is given to the stopband constraint penalty function, so stop P goes down and corr P goes up. It should be noted that, since the SMIA algorithm takes a random sequence as the initial sequence, stop P is not monotonically decreasing with  , and corr P is not monotonically increasing.

The Number of Frequency Stopbands
The power spectrum density of the optimized waveform under the different number of frequency stopbands is shown in Figure 10. The frequency stopband of Figure 10a  It can be seen from the figure that, when there is only one frequency stopband, the waveform PSD notch depth is the lowest, which is −25 dB. With the increase in the number of stopbands, the notch depth of each frequency stopband becomes shallow, which will cause performance degradation of interference suppression. Therefore, in the process of waveform design, it is necessary to set the number of stopbands reasonably.

Frequency Stopband Weighting Effect
In practice, the interference intensity of each interference frequency band received by airborne radar is not the same. If the notch of the same depth is formed, the frequency band with stronger interference intensity will not be effectively suppressed, and the weaker frequency band will occupy more resources. Therefore, the frequency stopband resources can be reasonably allocated by changing the weight vector of the SMIA algorithm Equation (6). This section assumes that the sequence length is N = 100, the region of reducing the sidelobes of specified-range cells is

WISL and Merit Factor
The weighted integral sidelobe level is one of the criteria for measuring the quality of the sequence, which is defined as where the smaller the WISL value, the better the final sequence obtained by the optimization algorithm.
Minimizing the WISL is equivalent to maximizing the merit factor (MF), which is defined as The variation of the merit factor and normalized WISL value with the number of iterations under the constant modulus constraint and PAR constraint of different algorithms are shown in Figures 12 and 13, respectively. Here, the CAN algorithm in [26], the WeSCAN algorithm in [27], the SMIA-ISL algorithm in this paper under ISL criteria, and the SMIA-WISL algorithm in this paper under WISL criteria are mainly compared.
From Figure 12a, it can be found that the normalized WISL values of all algorithms showed a decreasing trend as the number of iterations increased. The curves tended to be horizontal and eventually converge as the number of iterations increased. The convergence rate of the SMIA algorithm was the fastest either under the ISL criterion or the WISL criterion. The curve tended to be horizontal at about 90 iterations, and its normalized WISL value eventually converged at −24.59 dB and −25.10 dB. Because the WeSCAN algorithm requires a large number of matrix operations and has high computational complexity, its convergence speed is slower than the SMIA algorithm, but it can finally obtain a smaller WISL value, so the WeSCAN algorithm is inferior to the SMIA algorithm but better than the CAN algorithm. Since the merit factor of Figure 12b is actually the reciprocal of the WISL value, the conclusion is basically consistent with Figure 12a, and the merit factor of the SMIA algorithm is the largest.
As shown in Figure 13a,b, the WISL value and merit factor of each algorithm under the PAR constraint are basically consistent with the constant modulus constraint, but the numerical value is much higher than the constant modulus constraint. This is because the PAR constraint relaxes the peak average ratio and increases the degree of freedom of waveform design, so a smaller WISL value and greater merit factor are obtained.

The Amount of Computation
The running time of the WeSCAN algorithm and SMIA algorithm was compared by a simulation shown in Figure 14. The simulation environment was MATLAB 2022a/Windows 10, Intel Core i7-8750H CPU@2.2GHz. In this experiment, the region of reducing the sidelobes of specified-range cells was

Conclusions
In this paper, an algorithm for designing sparse-frequency waveforms for airborne early warning radar based on the minimum wave energy in a stopband was proposed to solve the practical problem of dense interference encountered by the low-frequency-band airborne early warning radar working over land. Theoretical analysis and simulation experiments showed that the cognitive waveform optimization design algorithm proposed in this paper can effectively suppress dense interference and reduce the influence of strongly isolated point clutter and interference. It can also form low-range sidelobes at The number of iterations specified-range cells, which is beneficial for dim target detection. In addition, the waveform optimized by the algorithm in this paper has good low-intercept performance, and because the waveform is optimized by the circular iterative algorithm based on FFT operation, it has low computational complexity and is convenient for engineering implementation.
In practical engineering, the main lobe interference and side lobe interference exist simultaneously in complex electromagnetic environments; in the process of the optimization design of the transmitted waveform, the influence of main-lobe interference should be emphasized. At the same time, the bandwidth of the sparse waveform determines the clutter intensity of the airborne radar receiver. The parameters such as the stopband position and the number of stopbands can be obtained by detecting the spatial interference distribution in passive mode. The bandwidth of the transmitted waveform to some extent determines the clutter-to-noise ratio (CNR) of the airborne radar receiver, and its size is determined by the clutter suppression ability of the airborne radar system.