Enhanced Doppler Resolution and Sidelobe Suppression Performance for Golay Complementary Waveforms

: An enhanced Doppler resolution and sidelobe suppression have long been practical issues for moving target detection using Golay complementary waveforms. In this paper, Golay complementary waveform radar returns are combined with a proposed processor, the pointwise thresholding processor (PTP). Compared to the pointwise minimization processor (PMP) illustrated in a previous work, which could only achieve a Doppler resolution comparable to existing methods, this approach essentially increases the Doppler resolution to a very high level in theory. This study also introduced a further ﬁltering process for the delay-Doppler map of the PTP, and simulations veriﬁed that the method results in a delay-Doppler map virtually free of range sidelobes.


Introduction
Due to complementarity, Golay complementary waveforms are effective at producing a satisfactory resolution range in a delay-Doppler map, as well as theoretical free-range side lobes at zero Doppler. Nevertheless, obvious range sidelobes are induced in nonzero Doppler intervals by their sensitivity to Doppler mismatch during matched filtering, and it is hard for conventional sidelobe suppression methods such as windowing to eliminate them.
A decade ago, Calderbank and Pezeshki et al. addressed the above problem by carefully designing Golay complementary waveforms in a specific transmitted order, named the Prouhet-Thue-Morse (PTM) design, which caused a satisfactory reduction in the range sidelobes in a narrow band around zero Doppler in the delay-Doppler map [1,2]. In a similar way, Suvorova et al. extended the idea of transmitted order design from the PTM sequence to the Reed-Müller codes, achieving a minimum range sidelobe level at a given Doppler bin in the delay-Doppler map [3]. This was further studied by Dang et al., who presented a binomial design (BD) algorithm that assigns weights to the matched filtering sequence of Golay complementary waveforms and significantly expands the sidelobe blanking area at the cost of an obvious decrease in the Doppler resolution [4]. Based on these works, Wu et al. employed semidefinite programming as a novel method to design complementary waveforms for improved sidelobe suppression as well as Doppler resolution [5,6].
However, the aforementioned methods only either preserve the Doppler resolution (none of them exceed the resolution of the conventional Golay pair) or enlarge the sidelobe blanking area, but they cannot achieve both at the same time. From another view, pointwise processing [7] (or cell-by-cell processing in some publications [8,9]) has been extensively researched in the processing of radar images to integrate the advantages of two figures and produce a further improvement in the signal-to-noise ratio (SNR). In our early work [10], a pointwise minimization processor (PMP) was proposed to combine the delay-Doppler maps of the BD algorithm with a weighted average Doppler (WD) algorithm under Golay waveforms, which maintained the Doppler resolution as well as a large sidelobe suppression area. However, as described before, the results of the PMP cannot exceed the original Doppler resolution of the waveform. Therefore, another designed pointwise thresholding processor (PTP) is proposed in this research to replace the PMP, which can further increase the Doppler resolution. A further filtering process is then applied based on the delay-Doppler map of the PTP, which almost eliminates the range sidelobes.
In the remainder of the manuscript, a brief introduction to Golay pairs and PMP is first given in Section 2, then the PTP is introduced. Section 3 presents the simulation results of the PMP and the PTP and compares the performance under fixed and randomized scenarios. The delay-Doppler map of a further filtering process after the PTP is also simulated to illustrate the improved sidelobe suppression performance. The conclusion and future directions are discussed in Section 4.

Golay Pairs
A Golay pair (complementary waveforms) consists of two length L sequences, x(l) and y(l), with several unimodular (±1) values/chips in each sequence [11]. The time width of each pair is LT c (T c for each chip). This waveform scheme is well known for its complementarity, i.e., the autocorrelation of the sequence pair is where C x (k) and C y (k) are the autocorrelation outputs of x(l) and y(l) at lag k, respectively, and δ(k) is the Kronecker delta function. The sequence pair cannot be transmitted in the time domain before modulating a baseband pulse Ω(t) with unit energy on each chip, which means the transmitted sequences are given as where Next, the transmission of either x(t) or y(t) is determined by a (P, Q) pulse train. Here, P = {p(n)} N−1 n=0 is a binary sequence and the transmitted pulses are presented as where P = {0, 1, 0, 1 . . .} denotes the standard transmission order and T represents the pulse repetition interval (PRI). On the other hand, Q = {q(n)} N−1 n=0 stands for the positive real number weights on the radar returns, where an all 1 sequence is the standard weighting. Next, the signal for matched filtering is written as Specifically, the BD algorithm [4] designs Q as a binomial sequence, i.e., Q = {C n N−1 } N−1 n=0 .
According to [12], we then calculate the delay-Doppler map of Golay complementary waveforms as follows: where "*" denotes the complex conjugation. Based on the Equations (4) and (5), the delay-Doppler map of the Golay pair is further expanded as The first item contains [C x (k) + C y (k)], which is an impulse function due to the complementarity; thus, the sidelobe in the delay-Doppler map is only influenced by the second item.
For transmission order design methods such as standard order and PTM design, q(n) is always 1; thus, the sub-item reaches 0 at θ = 2πnF D T = 2nπ N , which means it is free of a sidelobe along the θ-Doppler axes and the Doppler resolution of the target is 2π N . However, a significant sidelobe can be observed at other Doppler axes other than θ. Moreover, for the PTM, it is easy to calculate that Θ is negligible if 2πF D T is small.
For receiving weight design methods, i.e., the BD algorithm, p(n) alternates between 1 and 0, while q(n) is the coefficient of the binomial. Then, the sub-item Θ is further expressed as Obviously, the expression does not have zero points (which means the Doppler resolution is poor), while the value is an exponential function and its absolute value exponentially increases as 2πF D T increases. This is the reason why a large sidelobe blanking area is generated.
The delay-Doppler maps of previous works are plotted in Figure 1 for a better understanding. As is demonstrated, the delay-Doppler maps of standard order and PTM design are divided into "grids" by the zero points, while the BD algorithm obtains a large blanking area with an exponential-like sidelobe on the side. The delay resolutions of the aforementioned approaches are all 2T c , which is the width of the impulse function.  Figure 2 describes the procedure presented in our early work [10], where χ BD (t, F D ), χ WD (t, F D ) and χ(t, F D ), respectively, stand for the delay-Doppler maps of the BD algorithm, the WD algorithm and this procedure. Specifically the "Pointwise Processor" in [10], represents the pointwise minimization processor (PMP), whose result is denoted as The WD algorithm is demonstrated based on [3] to select the the optimal transmission order of Golay complementary waveforms (a standard weighted Q is set for this algorithm), which minimizes the sidelobes near a known Doppler value. Here, we employ the mean target Dopplerf d associated with their amplitudes to bring the sidelobe blanking area closer to the weak targets [13]:

Pointwise Minimization Procedure
whereÂ h andf d h are the normalized amplitude and Doppler of the hth target, respectively, and H is the number of targets in the delay-Doppler map. A tracker is usually used to estimate the target magnitude and Doppler from the past detections [14]. Already existing methods, such as [1][2][3][4], etc., are able to suppress the range sidelobes and improve the signal-to-noise ratio (SNR) near the targets to different extents. However, they cannot reduce the overall sidelobe magnitude in the underlying surveillance window. Therefore, a PMP was proposed in our previous paper as a nonlinear pointwise processor, which achieves sidelobe suppression involving sidelobe power reduction and does not cause target loss, as verified by technical simulations.
The advantage of the PMP is that it maintains the ideal large range sidelobe blanking region (in which the range sidelobes are less than −90 dB) provided by the BD algorithm and the acceptable Doppler resolution of the WD algorithm, based on the assumption that the targets are stable during the whole radar illumination. As described before, a drawback of the PMP is that it can only retain the improved Doppler resolution and the lower sidelobes of the two approaches. It still needs further enhancement when the Doppler resolution and sidelobe suppression performance fail to meet the requirements; thus, we propose the PTP in the following.

Pointwise Thresholding Procedure
The "Pointwise Processor" in Figure 2 is defined as the pointwise thresholding processor (PTP) and χ PTP (t, F D ) is defined as the output of the PTP. Then, the processor is expressed as where the threshold thr is artificially delimited considering the magnitude difference of targets and sidelobes. As the most important parameter in the processor, it directly influences the sidelobe level and the resolution (and also the performance of further filtering that will be discussed later). A small value will have no effect on the enhancement in Doppler resolution and sidelobe blanking performance, while a too large threshold may also blank the targets. Obviously thr ≥ 0 dB, and we commonly consider it nonsensical if thr > 10 dB, since the radar return fluctuation of two illuminations caused by target micromotion and other interferences normally cannot reach such a high level. In this paper, we choose thr = 2 dB as an example [15], but further research needs to be performed for better determination of the practical threshold. Under the same assumption, the PTP is expected to bring a further increase in the Doppler resolution and suppress the sidelobe magnitude compared to the PMP, which will be illustrated by the simulations in the next section.

Simulation and Further Discussion
The PTP is verified in simulations with the following global parameters when no other demonstration is presented. The targets in the simulations are set as Swerling II targets with 10% fluctuation in the radar cross-section (RCS).

Fixed Scenario
We first consider a fixed scenario with three targets (one weak and two strong targets), whose ground truth locations and magnitudes are listed in Table 1 and shown in Figure 3. Note that two strong targets can only be separated from the Doppler.   The delay-Doppler maps of the BD and WD algorithms and the outputs of the PMP and PTP with thr = 2 dB are given in Figure 4, and the comparison results of PMP and PTP in terms of the delay cross-section and Doppler cross-section, respectively, at all the targets location are then illustrated in Figures 5 and 6. Obviously, though PMP maintains the Doppler resolution of the WD algorithm, it performs much worse than the PTP. In addition, the overall SNR of the PTP in theory is also remarkably higher than the PMP (which may result in a higher performance during target detection). However, the processing times of the PMP and PTP will be twice those of separately using the BD or WD algorithm.  The PTP results at different thresholds are also compared in Figure 7. When thr = 1 dB, the weak target is nearly blanked by the PTP, while two false targets near the strong targets may be detected if thr is increased to 8 dB. This further verified the previous illustration of the thresholding.

Further Filtering for the PTP
Even though we have discussed the better target detection ability of the PTP than PMP, the targets are still not easily visually detected. Therefore, a further filtering for PTP is proposed in this subsection, which picks the targets through its displayed characters in the map and almost eliminates the range sidelobes after filtering.
According to the delay-Doppler maps we obtain, it is found that the targets are displayed typically (approximate to a rectangle, whose size is related to the threshold thr as well as the pulse number N) in the results of the PTP, whereas the range sidelobes are usually irregular and cannot be described by the common shape. This character gives us a chance to further filter the targets and suppress the sidelobes in the images. A particular rectangle with a similar size to the target can be employed to search the delay-Doppler map after the PTP, and the target is considered to be found when all the values in the rectangle are higher than −90 dB (since the original effect of the PTP suppresses some of the range sidelobes lower than −90 dB). Note that the size of the target should be evaluated first before the above operation.
Based on the previous illustration in [13] and our subsequent analysis, we learned that the delay resolution of the target after PTP mainly depends on the threshold thr, while the Doppler resolution is primarily influenced by the pulse number N. On the other hand, the original delay and Doppler resolution of target can be analytically calculated as 2T c and 2π N , as described in Section 2.1, which means they occupy 20 and 19 pixels in the delay and Doppler axes, respectively. Though it is still hard to analytically calculate the size of this particular rectangle, which is expressed as the number of pixels occupied in the row (delay axis, R r ) and column (Doppler axis, R c ) of the delay-Doppler map, a numerical fitting could be adopted to obtain an experimental Equation (13) with a certain calculable basis to illustrate the R r and R c of the rectangle. By observing Figures 4 and 7, we find that the delay and Doppler pixels of the target both shrink by about a half after the PTP, while the number of delay pixels needs to be further increased, which is nearly equal to the value of thr. Therefore, the experimental Equation (13) is written as follows, where the operator "round" means calculating to the nearest integer. (13) This experimental equation is only used for an explanation of the filtering process that is handled in this work. A more careful deduction of the calculation of this rectangle in practical studies will be the next step of our research.
The outputs of further filtering of the PTP at different thresholds and pulse numbers are demonstrated in Figure 8, which exhibit that further filtering makes the delay-Doppler map almost free of range sidelobes and the targets can be clearly visually recognized. Nevertheless, the lower threshold may lead to a higher probability of miss detection (Target 3 is lost under thr = 1 dB), and the increase in pulse number may generate more false targets (some false targets arise when N = 64).

Randomized Scenario
In this scenario, several cases of Swerling II targets with different numbers are uniformly distributed in the delay-Doppler map. Based on the previous global parameters, we consider the following four cases: (1) Target number: 2 (one strong and one weak); (2) Target number: 3 (one strong and two weak); (3) Target number: 4 (two strong and two weak); (4) Target number: 5 (three strong and two weak).
thr and N are fixed as 2 dB and 32, respectively, in all cases. For the sake of a clearer explanation, the target detection thresholds in these cases are set to the magnitude of the weakest target, which means that all the targets can be detected but false targets may also exist in the range sidelobes. Again, the proper setting of realistic detection thresholds requires further consideration.
A Monte Carlo simulation was operated for 1000 iterations for each case above, and the number of correct detections was calculated. A correct detection is counted when targets are all detected without any false targets. The correct detection occurrences of the PMP, the PTP and the PTP after further filtering are shown in Figure 9.  As is observed, the PTP provided in this paper outperforms the previous proposed PMP due to a higher overall SNR, as discussed before. The output of the PTP after further filtering has even more correct detection occurrences than the others since the improved range sidelobe effect is achieved.

Conclusions
In this paper, a signal processing method involving a PTP is proposed for Golay complementary waveforms to achieve an essentially enhanced Doppler resolution compared to the previously proposed PMP, which can only maintain the original Doppler resolution of this waveform scheme. To solve the visual recognition problem of targets in a delay-Doppler map, further filtering of the PTP by extracting the targets more precisely through a particular rectangle is also employed for a significant improvement in range sidelobe suppression. The performance of the above methods are verified by simulation results. Our future research avenues may concern practical experiments of adaptive thresholding (such as constant false alarm rate, CFAR) and target detection under this waveform scheme, as well as some robust waveform optimization methods for complex target detection [16,17].
Author Contributions: Conceptualization and supervision, X.H.; methodology, J.Z. and Y.S.; software validation, J.Z., N.J. and Z.X.; writing-original draft preparation J.Z.; writing-review and editing, N.J., Z.X. and C.F.; funding acquisition, J.Z. All authors have read and agree to the published version of the manuscript.