Range-Ambiguous Clutter Suppression via FDA MIMO Planar Array Radar with Compressed Sensing

Range-ambiguous clutter is an inevitable issue for airborne forward-looking array radars, especially with the high pulse repetition frequency (PRF). In this paper, a method to suppress the range-ambiguous clutter is proposed in an FDA-MIMO radar with a forward-looking planar array. Compressed sensing FDA technology is used to suppress the range-ambiguous clutter and the forward-looking non-uniformity short-range clutter of radar. Specifically, first, the range ambiguous clutter in different regions is separated by the characteristics of the planar array radar elevation dimension and FDA radar range coupling. Meanwhile, regarding the issue of the FDA radar main lobe moving between coherent pulses, a main lobe correction (MLC) algorithm proposes a solution for the issue, where the FDA radar cannot coherently accumulate signals in the case of non-full angle illumination. Finally, compressed sensing technology and elevation dimension filtering are utilized to suppress the range ambiguous clutter at the receiver, with the approach alleviating the range dependence of clutter in the observation region. A small number of clutter snapshots can obtain an approximately ideal clutter covariance matrix through compressed sensing sparse recovery. The method not only reduces the number of training samples, but also overcomes the problem of clutter non-uniformity in the forward-looking array. Therefore, the clutter suppression problems faced by the high repetition frequency airborne radar forward-looking array structure are solved. At the analysis stage, a comparison among the conventional MIMO and FDA methods is carried on by analyzing the improvement factor (IF) curves. Numerical results verify the effectiveness of the proposed method in range-ambiguous clutter suppression.


Introduction
When an airborne radar detects ground targets in the downward-looking mode, the strong ground clutter power will submerge the target and it is difficult to detect the moving target. Therefore, the clutter suppression of airborne radar is the core issue in radar target detection. Space-time adaptive processing (STAP) technology has advantages and shows considerable development in solving airborne radar clutter suppression; it has been used in a variety of early warning radars. It requires training samples that meet the conditions of independent and identical distribution (IID), the number of which is larger than or equal to twice the radar degrees of freedom [1,2]. However, according to actual requirements, the radar array is not necessarily consistent with the flight direction of the carrier platform. There is a particular angle between the flight direction and the radar array plane. When the included angle is 90 degrees, the radar changes from a side view array to a forwardlooking array configuration. For example, most fire control radars are forward-looking array structures. Under this array structure, the relationship between Doppler and angular dimensional power spectrum. Then, by selecting a small number of snapshots around the range gate to be measured for sparse processing, the clutter support set approximate to a single range gate can be obtained, followed by recovery ofs the clutter covariance matrix according to the dictionary and support set, so as to overcome the issue of power spectrum clutter spread due to non-uniform clutter snapshot data.
Finally, the problems of range-ambiguous clutter and clutter range dependence of forward-looking array radar are solved. Compared with the performances of conventional planar array MIMO radar and FDA radar, simulation results prove the effectiveness of the proposed method. The innovations of this method are summarized as follows: (1) In view of the range dependence of the forward-looking array radar clutter, this paper uses compressed sensing technology to divide the dictionary grid and construct a sparse dictionary according to the characteristics of the forward-looking radar clutter steering vector in the spatial Doppler two-dimensional power spectrum structure. Combined with the FDA technology to resist the range ambiguous clutter, this paper overcomes the problem that the conventional forward-looking array multiple range ambiguous clutters lack sparsity. The clutter is transformed into a sparsely recoverable target. A small number of echo snapshots are received near the detection range gate through compressed sensing to reconstruct the clutter covariance matrix, which is not approximately affected by the range dependence. It solves the contradiction that the conventional STAP requires both training samples to meet the IID condition and the selection of enough samples. Especially when the forward-looking clutter range dependence of training samples is serious, this method improves the performance of radar clutter spreading and clutter suppression and the moving target detection ability of forward-looking array radar. (2) Aiming to address the issue of main lobe moving in the transmitting pattern of FDA MIMO pulse radar, the main lobe rectification method of equivalent transmitting pattern is proposed. During the same coherent processing interval, the initial phase of an FDA radar pulse is compensated and restored at the transmitting and receiving ends, respectively. Although the coherent processing pulse irradiates different angle areas, they are all able to contain the desired target angle, so that the coherent processing pulses can irradiate the same target, improving the signal-to-noise ratio and completing the coherent accumulation. On the basis of this, the subsequent signal processing of pulse framework FDA radar is carried out. (3) Addressing the clutter ambiguous issue of airborne forward-looking array radar, clutters from different ranges have different power spectrum distributed characteristics because of the phase term of range coupling by combining FDA radar and the elevation dimension freedom of planar array. Different from the conventional ambiguous clutter, this method makes the target signal and clutter of the desired region maintain high gain, while the range-ambiguous clutter of the undesired region has a wide area low gain distribution, which is suppressed. Combined with clutter range dependent compensation technology, radar target parameter searching is realized.

Transmitting of MLC-FDA Radar
In this paper, we will discuss the range-ambiguous clutter suppression method of airborne forward-looking planar array radar. The geometry of FDA radar is shown in Figure 1a, which is in the right-hand coordinate system. The radar pulse-repetition interval is T = 1/f PRF , where f PRF is the pulse repetition frequency. The number of coherent processing pulses is K. The speed of the radar platform is V, which is parallel to the y-axis. The planar array is perpendicular to the speed direction, and the height of the platform is H. The azimuth, elevation and spatial cone angle corresponding to the clutter block are expressed as ϕ, θ, ψ, respectively. Consequently, we have cos ψ = cos θ· cos ϕ. The planar array is M × N-dimensional array element distribution, in which the carrier frequencies transmitted by each row are the same, and the transmission waveforms of different rows Remote Sens. 2022, 14, 1926 5 of 29 keep orthogonality to each other, as shown in Figure 1b. The carrier frequencies emitted by the first to M-th array elements in the same column are shown in the following equation: f m = f 0 + (m−1)∆f m = 1,2, . . . , M, (1) where ∆f is the frequency increment between adjacent rows. It must keep much smaller than the carrier frequency f 0 . It should be noted that the large value size of ∆f will have the effect of decoherence of the target reflection coefficient [31,32].
azimuth, elevation and spatial cone angle corresponding to the clutter block are expressed as φ, θ, ψ, respectively. Consequently, we have cos cos cos The planar array is M × N-dimensional array element distribution, in which the carrier frequencies transmitted by each row are the same, and the transmission waveforms of different rows keep orthogonality to each other, as shown in Figure 1b. The carrier frequencies emitted by the first to M-th array elements in the same column are shown in the following equation: where ∆f is the frequency increment between adjacent rows. It must keep much smaller than the carrier frequency f0. It should be noted that the large value size of ∆f will have the effect of decoherence of the target reflection coefficient [31,32]. The corresponding normalized Doppler frequency, normalized elevation vertical spatial frequency and azimuth horizontal spatial frequency are, respectively, expressed as: The signal transmitted by the array element in row m can be expressed as: where Φ m (t) represents the orthogonal waveform corresponding to row m-th. Here, discussion is based on an ideal orthogonal waveform premise. Therefore, the conditions satisfied are: where τ represents time delay. ∀q, p = 1, 2, . . . , M. The conventional pulse system FDA radar pattern is characterized by S-shaped moving within the pulse width. When the pulse width is not large enough, the high gain irradiated zone of its main lobe is limited and the full angle coverage of [−90, 90] cannot be realized. Due to the frequency increment between array elements, the initial phases accumulated by each array element are also different after a pulse period, T. As a result, because each pulse has different initial phases, the high gain irradiated angular region of the main lobe within the pulse width is also different, i.e., the main lobe moving problem of FDA radar. Therefore, it is impossible to irradiate the target at the same angle with high gain within the coherent processing interval, let alone carry out signal processing such as coherent accumulation or clutter cancellation for the same target. The radar cannot complete its normal operation. Figure 2 is the schematic diagram of FDA radar range (time domain)-angle two-dimensional transmitting pattern. The red 0 degree is the desired target position. In the figure, three pulses transmitted by the radar are drawn. The gray box represents the pulse width, the main lobe distribution of FDA radar pattern in the box is represented by a blue dotted line, and the initial phase of the array element of the three pulses is also marked. In the figure, ∆f = 1/(2T). We note that the range of the main lobe irradiation angle of the first pulse pattern is [−30, 0], while each array element of the second pulse accumulates the initial phase of 2π∆fT(m−1), compared with the first pulse. M represents the array element serial number. The zone of the main lobe irradiation angle of the second pulse becomes [−60, −90] and [70,90]. The accumulated phase of the third pulse accumulation is 2π∆f 2T(m−1). As a result of ∆f = 1/(2T) in the figure, the initial phase of the third pulse accumulation is an integral multiple of 2π. Therefore, the third pulse and the first pulse have the same main lobe irradiation zone. In this paper, the phase weighting between pulses of transmitting array elements is used to adjust the beampattern main lobe position, so as to have a solution of main lobe moving of FDA radar. When transmitting the signal, each array element adds the main lobe moving correction weight, corresponding to pulses. Thus, the issue that the main lobe of coherent pulses cannot irradiate the target due to the accumulation of ∆f phase is adjusted, as shown in Figure 3, and the angle areas of the main lobe of the corrected FDA pulse pattern can irradiate the same target. In the coherent accumulation time, the radar can achieve high gain illumination of the same target. The red correction phase part in the figure shows the effect of the main lobe correction weight after the array element is weighted. In the weighted intra pulse pattern, the first pulse main lobe is irradiated in the area of [−30, 0], the second pulse is irradiated in the area of [−20, 20], and the third pulse main lobe is irradiated in the area of [−5, 38]. Three pulse main lobes all cover the desired target position, 0 degree.  In this paper, the phase weighting between pulses of transmitting array elements is used to adjust the beampattern main lobe position, so as to have a solution of main lobe moving of FDA radar. When transmitting the signal, each array element adds the main lobe moving correction weight, corresponding to pulses. Thus, the issue that the main lobe of coherent pulses cannot irradiate the target due to the accumulation of ∆f phase is adjusted, as shown in Figure 3, and the angle areas of the main lobe of the corrected FDA pulse pattern can irradiate the same target. In the coherent accumulation time, the radar can achieve high gain illumination of the same target. The red correction phase part in the figure shows the effect of the main lobe correction weight after the array element is weighted. In the weighted intra pulse pattern, the first pulse main lobe is irradiated in the area of [−30, 0], the second pulse is irradiated in the area of [−20, 20], and the third pulse main lobe is irradiated in the area of [−5, 38]. Three pulse main lobes all cover the desired target position, 0 degree.  In this paper, the phase weighting between pulses of transmitting array elements is used to adjust the beampattern main lobe position, so as to have a solution of main lobe moving of FDA radar. When transmitting the signal, each array element adds the main lobe moving correction weight, corresponding to pulses. Thus, the issue that the main lobe of coherent pulses cannot irradiate the target due to the accumulation of ∆f phase is adjusted, as shown in Figure 3, and the angle areas of the main lobe of the corrected FDA pulse pattern can irradiate the same target. In the coherent accumulation time, the radar can achieve high gain illumination of the same target. The red correction phase part in the figure shows the effect of the main lobe correction weight after the array element is weighted. In the weighted intra pulse pattern, the first pulse main lobe is irradiated in the area of [−30, 0], the second pulse is irradiated in the area of [−20, 20], and the third pulse main lobe is irradiated in the area of [−5, 38]. Three pulse main lobes all cover the desired target position, 0 degree.  The main lobe rectification weighting vector added to each array element is constructed as follows when transmitting signals: The main lobe rectification weighting vector added to each array element is constructed as follows when transmitting signals: where ⊗ represents the Kronecker product, g C ∈ C MK×1 , g C (mk) = e −j2π(m−1)∆ f (k−1+z k−1 )/ f PRF and z k−1 are the compensation coefficient corresponding to the k-th pulse, which determine the main lobe movement of the k-th pulse relative to that of the first pulse. Here, we take z k−1 = (k−1)x/180, with x the offset coefficient to reflect the angle difference of the main lobe illumination area of the adjacent pulse beampattern (since the FDA pattern is S-shaped, there is a case where individual pulse positions z k−1 need to be adjusted nonlinearly). Substitute the above equation into Equation (3), and the signal of the k-th pulse transmitted by the m-th array element reaching the clutter block can be written as: where ε is the product of target scattering coefficient and signal propagation gain. τ r = 2R/c, being a signal reaching the clutter block, the range delay term only produces one-way delay τ r /2; τ sm = (m − 1)d sin ϕ/c = (m − 1) f hm / f m is the phase delay term generated by the spatial position of the array element. c represents the speed of light and n c is the Gaussian white noise in the propagation process.cos ψ = cos θ· cos ϕ, f dm = 2v cos ψ f m /c represents the Doppler frequency, corresponding to f m generated by platform motion. In Equation (6), it is only affected by the velocity of the transmitting array, so we have f dmt = v cos ψ f m /c. Since the radar transmission is a narrow-band signal, there The last term −j2π(m − 1)∆ f (k − 1 + z k−1 )/ f PRF in the above equation is the compensation term substituted into the transmitted signal Equation (3) by main lobe correction weighting, which corresponds to the red correction phase part in Figure 3.

Receiving of MLC-FDA Radar
After being reflected by the clutter block, the k-th pulse echo signal received by the n-th array element can be written as: where τ sn = sin θ cos ϕ(n − 1)d/c = (n − 1) f v / f m . It can be seen from the above equation it is necessary to separate the signals of each channel, and then carry out subsequent signal processing. The processing process is shown in Figure 4. After the received signal reaches each receiving array element, it is presented as s n . Each channel is separated into M channels through matched filtering and down conversion. Therefore, the received signal is changed from N channels to MN channels. After adding a total of K coherent pulses, the received data is MNK dimension. nel transmitting array elements. After the superimposed signals enter the receiver, it is necessary to separate the signals of each channel, and then carry out subsequent signal processing. The processing process is shown in Figure 4. After the received signal reaches each receiving array element, it is presented as sn. Each channel is separated into M channels through matched filtering and down conversion. Therefore, the received signal is changed from N channels to MN channels. After adding a total of K coherent pulses, the received data is MNK dimension. The signal of the k-th pulse, emitted by the m-th one and received by the n-th element, is presented as: where The signal of the k-th pulse, emitted by the m-th one and received by the n-th element, is presented as: where Since Equation (1) ∆f << f 0 , we have: Remote Sens. 2022, 14,1926 10 of 29 The phase term containing ∆f in the FDA expression can be simplified to: The received and processed data is MNK*1-dimensional data, which can be written in the form of a steering vector, as: where Their specific forms are given below: where The main lobe rectification compensation weight vector is written as: where According to Equation (12), when there are U times of range ambiguous areas and N C clutter blocks in the same illumination range, the total echoes received by the receiver can be written as:

MLC-FDA Range-Ambiguous Clutter Suppression Method Based on Compressed Sensing
According to Equation (13), the spatial angular frequency of composite vertical transmitting is expressed as For the desired target located at R 0 , according to [25,33], and especially Equation (15) of Article [25], we construct the range compensation steering vector: where are all one-column vectors. After substituting the above equation into Equation (15) for compensation, the composite vertical transmitting spatial angular frequencies of the unambiguous region and the 1st ambiguous region are expressed as: where R 0 and R 1 represent the oblique ranges of the unambiguous range region and the 1st range-ambiguous region, respectively. R 1 = R 0 + R u . The maximum unambiguous range can be expressed as R u = c/(2f PRF ). Figure 5 shows the situation where the signals of different range-ambiguous clutter are superimposed into the radar receiver. Because the reflected clutter blocks differ by a maximum unambiguous range R u , and the difference between pulses is a pulse period T = 1/f PRF , taking the range unambiguous region and the 1st range-ambiguous region as an example, the superimposed pulses differ by one pulse. For the k-th transmitting pulse, the sequence number of the received pulse corresponding to the range unambiguous range area is k unchanged. Here we take coherent pulse number K = 5 as an example. The sequence number corresponding to the superimposed pulse in the 1st range ambiguous area is k 1 = mod 5 |(k − 1) + 4|+1 . When k = 1, 2, 3, 4, 5, the echo sequence numbers received from the 1st range ambiguous area are k 1 = 5, 1, 2, 3, 4, respectively.
OR PEER REVIEW 12 of 29 where R0 and R1 represent the oblique ranges of the unambiguous range region and the 1st range-ambiguous region, respectively. R1 = R0 + Ru. The maximum unambiguous range can be expressed as Ru = c/(2fPRF). Figure 5 shows the situation where the signals of different range-ambiguous clutter are superimposed into the radar receiver. Because the reflected clutter blocks differ by a maximum unambiguous range Ru, and the difference between pulses is a pulse period T = 1/fprf, taking the range unambiguous region and the 1st range-ambiguous region as an example, the superimposed pulses differ by one pulse. For the k-th transmitting pulse, the sequence number of the received pulse corresponding to the range unambiguous range area is k unchanged. Here we take coherent pulse number K = 5 as an example. The sequence number corresponding to the superimposed pulse in the 1st range ambiguous area is k1 = Since we compensate the main lobe movement of the FDA pattern during transmission, we need to eliminate the effect of main lobe rectification compensation on the signal when receiving the signal. Therefore, we carry out corresponding main lobe correction reduction compensation for the receiving signal, and the reduction compensation vector is constructed as follows: Since we compensate the main lobe movement of the FDA pattern during transmission, we need to eliminate the effect of main lobe rectification compensation on the signal when receiving the signal. Therefore, we carry out corresponding main lobe correction reduction compensation for the receiving signal, and the reduction compensation vector is constructed as follows: This is substituted into Equation (15) (taking the signal of the k-th pulse, emitted by the m-th one and received by the n-th element as an example). The process of twice compensation is shown as follows. The echo signal in the range unambiguous region is written as follows: The echo from the 1st range-ambiguous region can be expressed as: where ∆z = z mod 5 |k−1+4| − z k−1 . After two compensations, because of R u = c/(2 f PRF ), the 1st ambiguous region echo of Equation (20) can be further written as: β = mod 5 |(k − 1) + 4|−(k − 1) + 1 + ∆z = mod 5 |k + 3|−k + 2 + ∆z . According to the above equation, through twice compensation, the signal from the unambiguous region is the same as MIMO radar signal form, and its FDA range coupling phase term and main lobe correction phase term are compensated and eliminated. It can be denoted as: The signal form of the 1st ambiguous clutter region can be expressed as: From Equations (22) and (23), we can see that the spatial angular frequency of composite vertical transmitting of the range unambiguous area is f cv0 = d/c sin ϕ f 0 . The spatial angular frequency of composite vertical transmitting of the 1st range ambiguous area is Compared with the range unambiguous region, the composite vertical transmission spatial angular frequency of the range ambiguous region has one more term, −∆ f β/ f PRF , and the β corresponding to different range ambiguous regions is different. Meanwhile, for the same range ambiguous region, when the value of k is different, the corresponding value of β is also different. The coherent pulse number is denoted as K. Taking K = 5 with the offset coefficient x = 3 as an example, when k equals 1, 2, 3, 4, and 5, respectively, the β(k) values are 76/15, −1/60, −1/60, −1/60, and −1/60. Therefore, during power spectrum scanning, for the same ambiguous area, the spatial angular frequencies of the composite vertical transmitting are different. The clutter power spectrum is low gain with dispersed distribution, and a high gain clutter ridge can be difficult to form.
During the spectrum scanning, the clutter of the range unambiguous region is the same as the conventional MIMO radar signal form. Under the geometric framework of the forward-looking array, it presents a positive elliptical distribution in the elevation angle, horizontal angle and Doppler three-dimensional power spectrum. Since the radar illumination angle zone is [0, 180], the clutter presents a semi-elliptical shape. As shown in Figure 6, there is a three-dimensional clutter spectrum distribution of clutter in different areas under a single range cell. Figure 6a represents MIMO radar. The blue represents the clutter power spectrum distribution of the unambiguous area. It corresponds to a large elevation angle, with the clutter distribution characteristics leading to a small ellipse radius. The corresponding elevation angles of echoes from range ambiguous areas with different ranges are different. And the change in elevation angle decreases with increase in range. Therefore, the difference in corresponding elevation angles of clutter from range ambiguous areas is small. Their clutter distribution overlaps each other. Different orange semi-ellipses in the figure represent the clutter distribution of the 1st and the 2nd range-ambiguous regions. Figure 6b shows the distribution characteristics of range ambiguous clutter of conventional FDA radar. Because of the range coupled phase terms of elevation spatial frequency, clutter corresponding to different ranges are separated from one another. Figure 6c shows the range-ambiguous clutter of MLC-FDA radar. Through twice compensation, the phase term of the clutter unambiguous region is the same as that of conventional MIMO and presents a semi-elliptical distribution. Due to the existence of the −∆ f β/ f PRF term in the spatial angular frequency of composite vertical transmission of the ambiguous region, the clutter energy of the range ambiguous region is challening to concentrate to form a high gain clutter ridge during power spectrum scanning. Instead, it presents a low gain distribution in the whole three-dimensional space.   Substituting the two compensation vectors g RC , e R into Equation (15), it can be denoted as: Subsequently, after elevation dimension filtering, the clutter distribution in the horizontal receive angle Doppler two-dimensional domain of the signal after two compensations is shown in Figure 7. Substituting the two compensation vectors RC g , R e into Equation (15), it can be denoted as: Subsequently, after elevation dimension filtering, the clutter distribution in the horizontal receive angle Doppler two-dimensional domain of the signal after two compensations is shown in Figure 7.  Due to the forward-looking array configuration, the clutter after elevation dimension filtering still has serious range dependence. We use the compressed sensing method and a small number of snapshots data around the desired range gate to recover the clutter covariance matrix of the desired position. Firstly, we construct the horizontal reception angle and Doppler two-dimensional domain dictionary Ω. N s is the horizontal reception spatial dispersion coefficient and N d is the Doppler dispersion coefficient. The horizontal reception angle and Doppler frequencies are f h ∈ [−0.5, 0.5], f d ∈ [−0.5, 0.5], respectively. A complete clutter spectrum dictionary is constructed as: and the dictionary vector can be written as: where The two-dimensional power spectrum is meshed by sparse dictionary Ω. For L range gate echo data: According to the reconstruction process of sparse recovery, we have: where µ is the sparsity constraint parameter. Here, we take µ as one third of the freedom degree NK of the planar array. We use the orthogonal matching pursuit (OMP) method for sparse recovery. The calculation process will not be repeated; CVX toolbox or an iterative algorithm can be used. The algorithm can be referred to in the relevant literature. It is worth noting that we set two iteration termination conditions in the calculation process; one is that y − Ωγ 2 2 is less than or equal to the noise power, and the other is that the number of iterations is greater than or equal to the sparsity constraint parameter µ. The OMP algorithm will terminate the iteration if any condition is met. After obtaining the sparse coefficient matrix γ, the clutter covariance matrixR = Ωγ * (Ωγ) H of compressed sensing estimation is recovered, and the clutter adaptive weight is calculated according to the clutter covariance matrix. The method performance is discussed in the next section. Figure 8 shows the entire compressed sensing main lobe rectification FDA radar signal processing flow. The main lobe rectification vector g C ∈ MK × 1 is weighted at the transmitting end to correct the transmitting main lobe. The K coherent pulses received by the N receiving array elements are converted into MNK-dimensional signal data through orthogonal filtering and down-conversion. The secondary range compensation and main lobe rectification restoration compensations are carried out. The echo of the range unambiguous region after compensation is consistent with the conventional MIMO signal form. The transmission spatial angular frequency of the signal from the rangeambiguous region contains the range coupling phase term −∆ f β/ f PRF , and can therefore be different from the signal of the unambiguous region in the power spectrum. At the same time, the separation of the ambiguous clutter is completed. Then signal is converted into NK dimension data through elevation dimension filtering. Then, the data is processed by compressed sensing covariance matrix sparse recovery, adaptive filtering and other subsequent signal processing. mension data through elevation dimension filtering. Then, the data is processed by compressed sensing covariance matrix sparse recovery, adaptive filtering and other subsequent signal processing.

Simulation Results
In this section, the algorithms above will be simulated and discussed. According to the signal processing sequence, there are three parts. The simulation mainly discusses the correction of the pattern main lobe after transmission, the three-dimensional power spectrum distribution of clutter after twice compensation, and the final clutter suppression by compressed sensing sparse recovery after elevation dimension filtering. The radar system parameters are listed in Table 1, unless otherwise specified.

MLC FDA Radar Main Lobe Correction Emission Pattern
The main lobe irradiation angle areas of three pulse FDA emission patterns are shown in Figure 9. It is assumed that the target is at 0 degree, the pulse width set as 4.76 × 10 5 s, and the pulse period is T= 6.67 × 10 4 s, ∆f = 1/(2T). Because FDA radar contains the phase term of range coupling, the irradiation angle zone of the main lobe within the pulse

Simulation Results
In this section, the algorithms above will be simulated and discussed. According to the signal processing sequence, there are three parts. The simulation mainly discusses the correction of the pattern main lobe after transmission, the three-dimensional power spectrum distribution of clutter after twice compensation, and the final clutter suppression by compressed sensing sparse recovery after elevation dimension filtering. The radar system parameters are listed in Table 1, unless otherwise specified.

MLC FDA Radar Main Lobe Correction Emission Pattern
The main lobe irradiation angle areas of three pulse FDA emission patterns are shown in Figure 9. It is assumed that the target is at 0 degree, the pulse width set as 4.76 × 10 5 s, and the pulse period is T= 6.67 × 10 4 s, ∆f = 1/(2T). Because FDA radar contains the phase term of range coupling, the irradiation angle zone of the main lobe within the pulse width changes with the range, showing an S-shaped moving. It can be observed from pulses in the figure that the irradiation position of the main lobe tilts in the direction of the angle decreasing with the change in range (time). Three continuous pulses are shown in Figure 9a. In the conventional FDA pattern, the main lobe irradiation region moves around between pulses. From the observation of Figure 9a, we can see that the 1st pulse main lobe irradiation region is [−30 • , 0 • ]. Due to the existence of ∆f, each array element accumulates different initial phases when transmitting the 2nd pulse after one pulse period. Therefore, the main lobe irradiation angle is affected by the initial phase, and the irradiation zone of the main lobe changes to [−68 • , −90 • ] and [26 • , 90 • ]. At the 3rd pulse, because of ∆f = 1/(2T), the initial phase of each array element ∆f undergoing 2T accumulation is an integer multiple of 2π, so the 3rd pulse main lobe irradiation zone is the same as that of the 1st pulse. From the Figure 9a, we can see that the angle of the main lobe irradiated from different pulses is not the same. The conventional FDA cannot maintain the position of the main lobe beampattern. Therefore, during pulse coherent accumulation processing, the desired target echo signal received by the radar cannot maintain high gain, so it is impossible to carry on subsequent signal processing. To overcome these problems, we use the main lobe correction weight. The elements are weighted between coherent pulses in Figure 9b, so as to suppress the influence of the initial phase accumulation generated by ∆f on the main lobe irradiation angle. Therefore, the main lobe with different beam directions within pulse width can irradiate the same angle of the target, so that the target echo signal has high gain of the main lobe. As shown in Figure 9b, the 0 • position of the desired target is within the irradiation areas of the main lobe of the different pulses. According to Equation (5), the 1st pulse period initial phase compensation is 0, so the irradiation position of its main lobe is the same as that in Figure 9a. The m-th array element compensation amount for the 2nd pulse period is 2π(m − 1)∆ f (1 + z 1 )/ f PRF . After main lobe correction compensation, the 2nd pulse main lobe area of the beampattern changes from [−34 • , 4 • ] of the 1st pulse, to [−16.5 • , 11.5 • ]. The 3rd pulse compensation is 2π(m − 1)∆ f (2 + z 2 )/ f PRF . After compensation, the main lobe irradiation area is [−1 • , 38 • ]. The target position 0 • has been covered in the irradiation areas of coherent pulse main lobes, so as to ensure that the radar can process the signals at the same angle. Figure 9c shows the situation after compensation with coherent pulse number K = 5. Similar to the above, it will not be discussed in detail here. main lobe irradiation region is [−30°, 0°]. Due to the existence of ∆f, each array element accumulates different initial phases when transmitting the 2nd pulse after one pulse period. Therefore, the main lobe irradiation angle is affected by the initial phase, and the irradiation zone of the main lobe changes to [−68°, −90°] and [26°, 90°]. At the 3rd pulse, because of ∆f = 1/(2T), the initial phase of each array element Δf undergoing 2T accumulation is an integer multiple of 2π, so the 3rd pulse main lobe irradiation zone is the same as that of the 1st pulse. From the Figure 9a, we can see that the angle of the main lobe irradiated from different pulses is not the same. The conventional FDA cannot maintain the position of the main lobe beampattern. Therefore, during pulse coherent accumulation processing, the desired target echo signal received by the radar cannot maintain high gain, so it is impossible to carry on subsequent signal processing. To overcome these problems, we use the main lobe correction weight. The elements are weighted between coherent pulses in Figure 9b, so as to suppress the influence of the initial phase accumulation generated by ∆f on the main lobe irradiation angle. Therefore, the main lobe with different beam directions within pulse width can irradiate the same angle of the target, so that the target echo signal has high gain of the main lobe. As shown in Figure 9b, the 0° position of the desired target is within the irradiation areas of the main lobe of the different pulses. According to Equation (5), the 1st pulse period initial phase compensation is 0, so the irradiation position of its main lobe is the same as that in Figure 9a. The m-th array element compensation amount for the 2nd pulse period is . After main lobe correction compensation, the 2nd pulse main lobe area of the beampattern changes from [−34°, 4°] of the 1st pulse, to [−16.5°, 11.5°]. The 3rd pulse compensation is . After compensation, the main lobe irradiation area is [−1°, 38°]. The target position 0° has been covered in the irradiation areas of coherent pulse main lobes, so as to ensure that the radar can process the signals at the same angle. Figure 9c shows the situation after compensation with coherent pulse number K = 5. Similar to the above, it will not be discussed in detail here.
(a)  Figures 10 and 11 show the power spectrum distribution of the conventional MIMO, FDA radars and the proposed approach in the horizontal, vertical spatial angular frequencies and Doppler three-dimensional space before the elevation dimension filtering. In order to illustrate the clutter ridge structure, the simulation in the figures shows the clutter power spectrum distribution of a single range cell of three range ambiguous regions. The main slant ranges of the range ambiguous areas are 9 km, 84 km, 159 km, respectively,  Figures 10 and 11 show the power spectrum distribution of the conventional MIMO, FDA radars and the proposed approach in the horizontal, vertical spatial angular frequencies and Doppler three-dimensional space before the elevation dimension filtering. In order to illustrate the clutter ridge structure, the simulation in the figures shows the clutter power spectrum distribution of a single range cell of three range ambiguous regions. The main slant ranges of the range ambiguous areas are 9 km, 84 km, 159 km, respectively, corresponding to the range unambiguous region, the 1st range ambiguous region, and the 2nd ambiguous region. The spatial power spectra with power greater than or equal to −10 dB and −30 dB are selected for explanation. Observing the three groups of data in Figure 10, we can see that the vertical spatial angular frequency corresponding to the unambiguous region in Figure 10a is f cv0 = 0.3. The corresponding frequencies of the clutters corresponding to the 1st and 2nd range ambiguous regions are f cv1 = 0.1222 and f cv2 = 0.0778, respectively. Since the difference is small, it can be seen from Figure 10a that they overlap each other in spectrum. Figure 10b shows the power spectrum of conventional FDA radar after full angle irradiation or solving the issue of the main lobe moving. Clutter of different regions are separated from each other, and the separation degree is affected by radar parameter ∆f. The elevation dimension filter and the positions of the separated clutter influence the suppression result of range ambiguous clutter. The proposed method, clutter power spectrum, is shown in Figure 10c. As described in the above section, after twice compensation, the spatial angular frequency of unambiguous signal echo composite vertical transmission f cv = d/c sin ϕ f 0 is the same as that of conventional MIMO. Therefore, the clutter ridge distribution is consistent with that of conventional MIMO. The corresponding vertical spatial angular frequency is f cv0 = 0.3. For range ambiguous areas, the spatial angular frequency of the composite vertical transmission is f cv = −∆ f β/ f PRF + d/c sin ϕ f 0 . Due to the phase term of the range region and pulse number coupling, the clutter echo in the range ambiguous region difficultly forms a clutter ridge with high gain during power spectrum scanning, but is distributed in a wide area with low gain in three-dimensional space. Therefore, the three-dimensional clutter power spectrum is not displayed. Therefore, we distinguish the clutter in the range unambiguous region from the clutter in the range ambiguous region. By observing the −30 dB three-dimensional power spectrum, it can be seen that the display zones of the conventional MIMO radar and FDA radar clutter power spectrum in Figure 11a,b become larger. While in Figure 11c, the clutter in the ambiguous region of the proposed method is irregularly distributed in three-dimensional space. Compared with the first two methods, the clutter energy is dispersed in the whole three-dimensional space and the gain is low. uous region in Figure 10a is fcv0 = 0.3. The corresponding frequencies of the clutters corresponding to the 1st and 2nd range ambiguous regions are fcv1 = 0.1222 and fcv2 = 0.0778, respectively. Since the difference is small, it can be seen from Figure 10a that they overlap each other in spectrum. Figure 10b shows the power spectrum of conventional FDA radar after full angle irradiation or solving the issue of the main lobe moving. Clutter of different regions are separated from each other, and the separation degree is affected by radar parameter ∆f. The elevation dimension filter and the positions of the separated clutter influence the suppression result of range ambiguous clutter. The proposed method, clutter power spectrum, is shown in Figure 10c. As described in the above section, after twice compensation, the spatial angular frequency of unambiguous signal echo composite vertical transmission

Three-Dimensional Distribution Structure of Received Clutter Signal Power Spectrum
is the same as that of conventional MIMO. Therefore, the clutter ridge distribution is consistent with that of conventional MIMO. The corresponding vertical spatial angular frequency is fcv0 = 0.3. For range ambiguous areas, the spatial angular frequency of the composite vertical transmission is . Due to the phase term of the range region and pulse number coupling, the clutter echo in the range ambiguous region difficultly forms a clutter ridge with high gain during power spectrum scanning, but is distributed in a wide area with low gain in three-dimensional space. Therefore, the three-dimensional clutter power spectrum is not displayed. Therefore, we distinguish the clutter in the range unambiguous region from the clutter in the range ambiguous region. By observing the −30 dB threedimensional power spectrum, it can be seen that the display zones of the conventional MIMO radar and FDA radar clutter power spectrum in Figure 11a,b become larger. While in Figure 11c, the clutter in the ambiguous region of the proposed method is irregularly distributed in three-dimensional space. Compared with the first two methods, the clutter energy is dispersed in the whole three-dimensional space and the gain is low.
The above mainly discusses the proposed method's spatial distribution characteristics of clutter in range unambiguous and ambiguous areas after twice compensation, and compares it with the clutter structure of conventional MIMO radar and FDA radar.

Power Spectrum Clutter Suppression Effect and Analysis
As shown in Figure 8, after twice compensation and elevation dimension filtering, the MKN*1-dimensional echo data becomes KN*1-dimensional. Through different elevation dimension filter design, the radar has different suppression effects on clutter at dif-

Power Spectrum Clutter Suppression Effect and Analysis
As shown in Figure 8, after twice compensation and elevation dimension filtering, the MKN*1-dimensional echo data becomes KN*1-dimensional. Through different elevation dimension filter design, the radar has different suppression effects on clutter at dif- The above mainly discusses the proposed method's spatial distribution characteristics of clutter in range unambiguous and ambiguous areas after twice compensation, and compares it with the clutter structure of conventional MIMO radar and FDA radar.

Power Spectrum Clutter Suppression Effect and Analysis
As shown in Figure 8, after twice compensation and elevation dimension filtering, the MKN*1-dimensional echo data becomes KN*1-dimensional. Through different elevation dimension filter design, the radar has different suppression effects on clutter at different angles, which will not be discussed in depth here. We use the ordinary target beamforming results as the elevation dimension filter. The structure of the clutter spectrum under a single range cell is shown above. Here, we use all range cells to illustrate the final suppression effect. Figure 12 illustrates two-dimensional clutter spectrum scanning in the horizontal space domain and Doppler domain of conventional MIMO, FDA and the proposed approach. The range dependence of forward-looking array clutter is serious. Different clutter range cell echoes spreading in the power spectrum pollute a large target detection area, and the performance of radar parameter estimation decreases sharply. As shown in Figure 12a, the clutter of MIMO radar is seriously spread. The range ambiguous clutter makes the radar observation zone almost entirely occupied by the range dependent clutter. In Figure 12b, FDA radar improves the pollution of the clutter power spectrum, compared with MIMO, by separating the clutter of different range-ambiguous regions. Figure 12c shows the proposed approach. Due to the ambiguous clutter phase term coupled with ambiguous multiplicity and pulse number, the clutter has a low gain distribution in a wide area of three-dimensional space. In addition to the elevation dimension filtering, the clutter energy of the range ambiguous regions is further filtered. Therefore, in the two-dimensional power spectrum in the horizontal spatial domain and Doppler domain, the clutter power is the least affected by range-ambiguous clutter.
Sens. 2022, 14, x FOR PEER REVIEW 23 of 2 ferent angles, which will not be discussed in depth here. We use the ordinary target beam forming results as the elevation dimension filter. The structure of the clutter spectrum under a single range cell is shown above. Here, we use all range cells to illustrate the fina suppression effect. Figure 12 illustrates two-dimensional clutter spectrum scanning in the horizonta space domain and Doppler domain of conventional MIMO, FDA and the proposed ap proach. The range dependence of forward-looking array clutter is serious. Different clut ter range cell echoes spreading in the power spectrum pollute a large target detection area and the performance of radar parameter estimation decreases sharply. As shown in Figure  12a, the clutter of MIMO radar is seriously spread. The range ambiguous clutter makes the radar observation zone almost entirely occupied by the range dependent clutter. In Figure 12b, FDA radar improves the pollution of the clutter power spectrum, compared with MIMO, by separating the clutter of different range-ambiguous regions. Figure 12c shows the proposed approach. Due to the ambiguous clutter phase term coupled with ambiguous multiplicity and pulse number, the clutter has a low gain distribution in a wide area of three-dimensional space. In addition to the elevation dimension filtering, the clutter energy of the range ambiguous regions is further filtered. Therefore, in the two dimensional power spectrum in the horizontal spatial domain and Doppler domain, the clutter power is the least affected by range-ambiguous clutter.  Figure 13 gives us the improvement factor (IF) curves of three approaches at 0° posi tion. The black, blue and red represent the MIMO, FDA and the proposed approach, re spectively. Compared with MIMO and FDA radars, it can be seen from Figure 13 that the most serious clutter spreading phenomenon of the power spectrum is MIMO radar, fol lowed by FDA radar. The clutter ridge distribution of the two methods is relatively wide  Figure 13 gives us the improvement factor (IF) curves of three approaches at 0 • position. The black, blue and red represent the MIMO, FDA and the proposed approach, respectively. Compared with MIMO and FDA radars, it can be seen from Figure 13 that the most serious clutter spreading phenomenon of the power spectrum is MIMO radar, followed by FDA radar. The clutter ridge distribution of the two methods is relatively wide on the Doppler dimension power spectrum. The clutter spreading image of this method is the lightest and the clutter ridge is the narrowest. Meanwhile, the main lobe region of the coherent pulse beampattern can be maintained. The performance of radar target detection is superior.
PEER REVIEW 25 of 29 Figure 13. IF curves.
Then we use the compressed sensing method to restore the unambiguous range dependent clutter covariance matrix. Figure 14a,b show the reception space and Doppler domains' two-dimensional power spectrum scanning figures, corresponding to the compressed sensing restored clutter covariance matrix of 10 snapshot data and the ideal clutter covariance matrix. Compared with Figure 12c, it can be seen that, after compressed sensing restoration, the clutter spreading phenomenon is greatly improved. The clutter covariance matrix after compressed sensing is very close to the ideal covariance matrix of the desired range gate in Figure 14b, and this method only needs 10 training samples. It not only solves the problem of few training samples meeting the IID condition, but also obtains the range independent covariance matrix. It reduces the serious range dependence of radar clutter and discards the influence of clutter in range ambiguous areas of target Then we use the compressed sensing method to restore the unambiguous range dependent clutter covariance matrix. Figure 14a,b show the reception space and Doppler domains' two-dimensional power spectrum scanning figures, corresponding to the compressed sensing restored clutter covariance matrix of 10 snapshot data and the ideal clutter covariance matrix. Compared with Figure 12c, it can be seen that, after compressed sensing restoration, the clutter spreading phenomenon is greatly improved. The clutter covariance matrix after compressed sensing is very close to the ideal covariance matrix of the desired range gate in Figure 14b, and this method only needs 10 training samples. It not only solves the problem of few training samples meeting the IID condition, but also obtains the range independent covariance matrix. It reduces the serious range dependence of radar clutter and discards the influence of clutter in range ambiguous areas of target detection. The radar surveillance area will no longer face the pollution of a large number of clutter signals. This method improves the range ambiguous clutter overlap and clutter range dependence of high-speed moving radar in the forward-looking planar array system, and improves the radar target detection performance. affected by range ambiguous clutter. And the improvement effect after compensation is limited. Comparing with Figure 12, we notice that the clutter power spectrum spreading caused by the clutter range-dependent problem of the three radar systems has been improved. From the improved clutter power spectra, we can more clearly observe the influence of clutter in the ambiguous area on the clutter power spectra. The parameter estimation performance of the radar power spectrum for the MIMO and FDA methods is reduced because of range ambiguous clutter. The range compensation method is affected by the range-ambiguous clutter. At the same time, the clutter range dependence leads to widening of the clutter ridge, and the improvement effect after compensation is limited. Therefore, the proposed method has the best performance.  Figure 15 is the performance evaluation curve of the improvement factor (IF) of radar clutter in the Doppler domain with a horizontal receiving angle of 0 degree. The black is the IF curve of MIMO radar. It can be seen from the curve that there are two notches at −0.5 and −0.25. We can see that they are caused by range unambiguous region clutter and the 1st range ambiguous region clutter, respectively. In these two areas and their middle Subsequently, through comparison with the DW (Doppler warping) method, we illustrate the advantages of this method on forward-looking array range dependence suppression. The DW method compensates the sample data of different range cells by estimating the Doppler dimensional frequency differences between the fixed angle expected range cell and other training samples. After elevation dimension filtering, the conventional MIMO and FDA methods use the DW method to suppress clutter spreading, and a comparison is made with the method in this paper. The compensated power spectra are shown in Figure 14c,d. Figure 14c shows that of conventional MIMO radar. Compared with Figure 12a, the influence of clutter range dependence in the range unambiguous area on power spectrum spreading is significantly reduced. The clutter of the range ambiguous area still has an impact on the radar power spectrum after elevation dimension filtering. In Figure 14c, the area with f d /f r = −0.2 can be observed. The clutter from the 1st range-ambiguous region generates −10 to −12 dB clutter power here, polluting the radar detectable area. Figure 14d shows the clutter power spectrum corresponding to FDA radar after compensation. In addition, in the f d /f r = −0.2 zone, the distribution of clutter ridges in the 1st ambiguous region can be seen, and the power is between [−12, −16] dB.
Compared with conventional MIMO radar, the influence of clutter power gain is reduced. Due to the existence of range-ambiguous clutter in the two methods, the available zone of parameter estimation of the radar power spectrum decreases. The range compensation method is affected by range ambiguous clutter. And the improvement effect after compensation is limited. Comparing with Figure 12, we notice that the clutter power spectrum spreading caused by the clutter range-dependent problem of the three radar systems has been improved. From the improved clutter power spectra, we can more clearly observe the influence of clutter in the ambiguous area on the clutter power spectra. The parameter estimation performance of the radar power spectrum for the MIMO and FDA methods is reduced because of range ambiguous clutter. The range compensation method is affected by the range-ambiguous clutter. At the same time, the clutter range dependence leads to widening of the clutter ridge, and the improvement effect after compensation is limited. Therefore, the proposed method has the best performance. Figure 15 is the performance evaluation curve of the improvement factor (IF) of radar clutter in the Doppler domain with a horizontal receiving angle of 0 degree. The black is the IF curve of MIMO radar. It can be seen from the curve that there are two notches at −0.5 and −0.25. We can see that they are caused by range unambiguous region clutter and the 1st range ambiguous region clutter, respectively. In these two areas and their middle area of Doppler domain, the radar performance has decreased significantly. The blue is the IF curve of the conventional FDA radar. There are two notches at −0.5 and −0.25, the same as that of the MIMO radar, remaining the clutter echo of the range unambiguous and the 1st range ambiguous areas. As a result, similar to the MIMO radar curve, the radar performance degrades significantly at and around the ambiguous clutter notches. In the case of range unambiguous clutter, there is only one notch in the corresponding IF curve. DW technology can approximately align the spectrum peak in the Doppler frequency domain, so as to further improve clutter suppression performance. Conventional radar can effectively suppress the clutter ridge spreading and maintain good target detection performance. However, in the presence of range ambiguous clutter, the performance of these kinds of methods will decline significantly. Range ambiguous clutter and unambiguous clutter have different clutter spreading degrees in the Doppler domain. The compensation methods cannot compensate multiple range ambiguous clutter at the same time. Therefore, clutter pollution still exists in the radar monitoring area. The method in this paper, shown by the red line in Figure 15, eliminates the influence of ambiguous clutter on the radar power spectrum because of the different clutter distribution characteristics of different regions. Therefore, there is only one clutter at −0.5, and the radar performance is no longer affected by range ambiguous clutter. Because it is not affected by clutter in other areas, the power spectrum has a sparse characteristic. The clutter in the unambiguous range area is heavily dependent on range. The compressed sensing method restores the clutter covariance matrix with a small amount of snapshot data, so that the unambiguous clutter spreading phenomenon is the lightest, and there is no conflict between clutter compensation in different range regions, so as to improve the performance of target detection. longer affected by range ambiguous clutter. Because it is not affected by clutter in other areas, the power spectrum has a sparse characteristic. The clutter in the unambiguous range area is heavily dependent on range. The compressed sensing method restores the clutter covariance matrix with a small amount of snapshot data, so that the unambiguous clutter spreading phenomenon is the lightest, and there is no conflict between clutter compensation in different range regions, so as to improve the performance of target detection.

Conclusions
This paper introduces a range ambiguous clutter suppression method for forwardlooking planar array FDA-MIMO radar. The FDA-MIMO radar through weighting the

Conclusions
This paper introduces a range ambiguous clutter suppression method for forwardlooking planar array FDA-MIMO radar. The FDA-MIMO radar through weighting the main lobe rectification phase at the transmitting end makes the main lobe illumination zones of coherent pulses cover the same target angle, so that the radar can coherently accumulate for a target at the same range and angle. It improves the signal-to-noise ratio. Subsequently, signal compensations are carried out at the receiving end. Because of the range coupling characteristics of FDA echoes, clutters of different range-ambiguous regions have different composite transmission spatial frequencies after compensation. The signal of the observed desired region is the same as the conventional MIMO signal, while the range-ambiguous echoes of non-observed regions present a low gain discrete distribution in the power spectrum. Therefore, the ambiguous clutter can be eliminated. Because it is not affected by clutter in other areas, the power spectrum has a sparse characteristic. The approximate ideal covariance matrix can be recovered by using only a few clutter snapshots through compressed sensing. At the same time, the problems of insufficient training samples and non-uniformity of clutter data in the conventional STAP method are solved. The detection performance of the conventional STAP method is improved. The simulation results show the effectiveness of this method compared with conventional methods. It is important to emphasize that this approach is applicable to general array geometries.