Next Article in Journal
Multiscale Analysis of Reflected Radiation on Lunar Surface Region Based on MRRT Model
Previous Article in Journal
Exploiting the Sensitivity of Dual-Frequency Smartphones and GNSS Geodetic Receivers for Jammer Localization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Data-Independent Phase-Only Beamforming of FDA-MIMO Radar for Swarm Interference Suppression

1
School of Air Defense and Missile Defense, Air Force Engineering University, Xi’an 710051, China
2
College of Information and Communication, National University of Defense Technology, Wuhan 430010, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(4), 1159; https://doi.org/10.3390/rs15041159
Submission received: 3 January 2023 / Revised: 17 February 2023 / Accepted: 18 February 2023 / Published: 20 February 2023
(This article belongs to the Topic Radar Signal and Data Processing with Applications)

Abstract

:
This paper proposes two data-independent phase-only beamforming algorithms for frequency diverse array multiple-input multiple-output radar against swarm interference. The proposed strategy can form a deep null at the interference area to achieve swarm interference suppression by tuning the phase of the weight vector, which can effectively reduce the hardware cost of the receiver. Specifically, the first algorithm imposes constant modulus constraint and sidelobe level constraint, and the phase-only weight vector is solved. The second algorithm performs a constant modulus decomposition of the weight vector to obtain two phase-only weight vectors, and uses two parallel phase shifters to synthesize one beamforming weight. Both methods can obtain the phase-only weight to realize suppression for swarm interference. Simulation results demonstrate that our strategy shows superiority in beam shape, output signal-to-interference-noise ratio, and phase shifter quantization performance, and has the potential for use in many applications, such as radar countermeasures and electronic defense.

Graphical Abstract

1. Introduction

Beamforming is an essential technique in array signal processing and has been used to achieve interference suppression by forming a deep null in interference directions [1,2,3,4,5,6]. Currently, due to its powerful ability in signal processing, the technique has exhibited enormous potential in communication, radar, sonar and other fields. Unfortunately, the rapid development of unmanned aerial vehicles has made swarm interference a serious threat to airport aviation security [7,8]. It is important to suppress swarm interference to improve radar detection performance in the complex electromagnetic environment.
Over the past decades, we have witnessed explosive studies on numerous beamforming methods for phased array radar interference suppression, such as the classic linearly constrained minimum variance (LCMV) beamformer [9], minimum variance distortionless response (MVDR) beamformer, sampling matrix inversion (SMI) beamformer [10] and interference plus noise covariance matrix (INCM) reconstruction beamformer [11,12]. Recently, numerous studies have been conducted on data-independent beamforming methods [13,14,15,16,17,18,19]. Such methods do not require the acquisition of an accurate INCM matrix, and the performance of the beamformer is not affected by the sampling covariance matrix (SCM). Nevertheless, swarm interference is divided into mainlobe swarm interference and sidelobe swarm interference. The above phased array beamforming methods can only form a deep null at a specific angle and cannot suppress the mainlobe swarm interference by forming a null in a specific area.
Recently, due to the advantages of two-dimensional (2D) degrees of freedom (DOFs) in range and angle, frequency diverse array multiple-input multiple-output (FDA-MIMO) radar has been widely investigated [20,21,22,23,24,25,26]. Specifically, it can form a range-angle-dependent beampattern by adding the frequency offset between the elements of the transmitter. Thereby, the swarm interference can be suppressed by controlling the beampattern to form a two-dimensional deep null in the interference area [27].
Using FDA-MIMO radar, Ref. [28] proposed a preset broadened nulling beamformer (PBN-BF) by placing artificial interferences with appropriate powers around the nulls of the equivalent transmit beampattern. The enlarged widening of the nulls guarantees the effective suppression of mainlobe deceptive interference. However, a band-shaped null was generated based on the PBN-BF algorithm at the interference transmit spatial frequency, which made it impossible to form a 2D regional null. In Ref. [29], a data-independent beamforming method was proposed, which could form a regional null by assigning artificial interferences with prescribed powers in the joint transmit–receive spatial frequency domain. It could realize point-by-point successive null broadening control (SNBC) and multi-point concurrent null broadening control (CNBC). Based on the above strategies, these data-independent methods provide a way to suppress interference by adding artificial virtual interference to the interference area. Nevertheless, the swarm interference distributed in an area, excessively adding virtual interferences, brings beam distortion in the non-controlled area when forming a two-dimensional null. In addition, these methods ignore the consideration of output signal-to-interference-noise ratio (SINR).
Furthermore, the phase-only weight vector cannot be obtained by the aforementioned FDA-MIMO radar beamforming methods, and an amplitude adjustment unit with a high dynamic range is required behind each matched filter. Therefore, the hardware architecture is complicated, and high costs are involved when implementing these methods [30,31]. In the past, a few approaches to phase-only beamforming for phased array have been reported [32,33,34], such as using neural networks [32], numerical optimization techniques [33] and geometric approaches [34]. However, to date, the method to suppress swarm interference using a low-cost FDA-MIMO radar system has been rarely reported, which requires the weight vector in the receiver to be constant in amplitude but different in phase.
In order to fill the above gap, we propose two data-independent phase-only beamforming methods for FDA-MIMO radar against swarm interference using the powerful convex optimization theory. Specifically, we consider two different technical approaches to obtain the phase-only weight vector. The first method imposes constant modulus constraint and sidelobe level constraint to solve the phase-only weight vector with the maximum output SINR as the objective function. The second decomposes the weight vector into constant modulus complex numbers from a different perspective and uses the dual phase receiver to synthesize the one beamforming weight. Moreover, the constant modulus decomposition method can convert other controlling complex weight beamforming methods to phase-only beamforming. Simulation experiments show that both proposed methods can form a better beam shape in the interference area and obtain higher output SINR than the conventional beamforming methods. We summarize the main contributions of this study as follows:
(1)
We propose two data-independent phase-only beamforming methods for FDA-MIMO to suppress the swarm interference.
(2)
Our proposed algorithms can achieve interference suppression by only tuning the phase of the weight vector, which efficiently reduces the hardware cost of the FDA-MIMO radar.
(3)
Our proposed dual-phase shifter receiver for FDA-MIMO can convert other complex weight vector beamforming methods to phase-only beamforming.
(4)
The proposed method can obtain excellent output SINR with a small number of snapshots.
This paper is organized as follows. In Section 2, the signal model is established, and the related parameters are defined. Section 3 recalls the adaptive beamforming methods. In Section 4, two data-independent phase-only beamforming methods are introduced. In Section 5, numerical simulations are performed to demonstrate the performance of the proposed method, and Section 6 concludes the paper.
Notations: In this paper, · , · T and · H respectively denote the conjugate, the transpose and the conjugate transpose of the matrix or vector. · 2 refers to the Euclidian norm of the vector. The modulus of the complex number w is denoted by w . j = 1 denotes the imaginary unit. I m represents the m × m identity matrix, and m × n and m × n represent the sets of m × n complex matrix and real matrix, respectively. We use E · to represent the expected value operator. diag · indicates returning the elements of the main diagonal of the matrix as a vector. and represent the Hadamard product and Kronecker product, respectively. A B means A B is a positive semidefinite, i.e., the eigenvalues of A B are nonnegative. Finally, tr A represents the trace of a square matrix A .

2. Signal Model of FDA-MIMO Radar

Consider a colocated FDA-MIMO radar, as shown in Figure 1, in which the transmitter and receiver are uniform linear arrays (ULA) with N and M elements, respectively. The element spacing of the transmitter and the receiver is d = λ 0 / 2 , λ 0 = c / f 0 is the wavelength, c represents the speed of light, and f 0 is the reference carrier frequency. The frequency offset across the transmitter elements is Δ f .
Thus, the carrier frequency of the n-th transmitter element can be written as
f n = f 0 + n Δ f , n = 0 , 1 , N 1 .
The transmit signal of the n-th element can be expressed as
s n t = E N ψ n t e j 2 π f 0 + n Δ f t , 0 t T .
where E denotes the transmitted energy, T is the radar pulse duration and ψ n t is the baseband envelope of the n-th transmit element, i.e., 0 T ψ n 1 t ψ n 1 t d t = 1 , which is satisfied with the orthogonality condition
0 T ψ n 1 t ψ n 2 t τ e j 2 π Δ f n 2 n 1 t d t = 0 , n 1 n 2 , τ ,
where, τ is the delay time, n 1 , n 2 = 0 , 1 , N 1 .
Assuming that a far-field target is located at r 0 , θ 0 , τ n , m = 2 r 0 c n d sin θ 0 c m d sin θ 0 c is the echo delay time from the n-th transmit element to the m-th receive element, and we set τ 0 = 2 r 0 c . Then, the received signal of the m-th element can be expressed as
x m ( t ) = ξ n = 0 N 1 ψ n t τ n , m e j 2 π f 0 + n Δ f t τ n , m = ξ n = 0 N 1 ψ n t τ n , m e j 2 π f 0 + n Δ f t τ 0 + n d sin θ 0 c + m d sin θ 0 c = ξ n = 0 N 1 ψ n t τ n , m e j 2 π f 0 t τ 0 e j 2 π f 0 n d sin θ 0 c × e j 2 π f 0 m d sin θ 0 c e j 2 π n Δ f t τ 0 e j 2 π n Δ f n d sin θ 0 c + m d sin θ 0 c ξ e j 2 π f 0 t τ 0 e j 2 π f 0 m d sin θ 0 c n = 0 N 1 ψ n t τ 0 e j 2 π f 0 n d sin θ 0 c e j 2 π n Δ f t τ 0
where the approximation is due to n Δ f c n d sin θ 0 + m d sin θ 0 1 [35,36]. ξ represents the target echo complex coefficient (determined by the target reflection coefficient, space propagation coefficient, amplitude, phase of the transmit signal, etc.). It should be pointed out that this paper mainly studies the beamforming method, which is independent of the target doppler. The doppler factor is ignored for the sake of simplicity.
The signal processing chain of the receiver is shown in Figure 2. The received signal x m t by the m-th receiving element is first mixed with e j 2 π f 0 t in the analog device, and the output signal after mixed is
x ˜ m t = ξ e j 2 π f 0 τ 0 e j 2 π f 0 m d sin θ 0 c n = 0 N 1 ψ n t τ 0 e j 2 π f 0 n d sin θ 0 c e j 2 π n Δ f t τ 0 .
After the mixed output signal x ˜ m t is matched with ϕ n t in the matched filter, ϕ n t is defined as
ϕ n t = ψ n t e j 2 π n Δ f t .
According to Equation (3), after the matched filtering, the n-th output signal of the m-th receive element can be expressed as
x ^ n , m t = ζ t e j 2 π f 0 τ 0 e j 2 π f 0 m d sin θ 0 c e j 2 π f 0 n d sin θ 0 c e j 2 π n Δ f τ 0 ,
where ζ represents the complex-valued coefficient after matched filtering.
After matched filtering is performed on the M elements, the FDA-MIMO radar received target echo signal at time t can be expressed as
x s t = x ^ 0 , 0 t , , x ^ 0 , N 1 t , x ^ 1 , 1 t , , x ^ M 1 , N 1 t T = γ 0 Γ 0 b r θ 0 a t r 0 , θ 0
where γ 0 = e j 2 π f 0 τ 0 and Γ 0 N M × 1 is the complex vector after matched filtering. a t r 0 , θ 0 N × 1 and b r θ 0 M × 1 indicate the transmit and receive steering vectors, respectively, which have the forms of
a t r 0 , θ 0 = 1 , e j 2 π f t 0 , , e j 2 π N 1 f t 0 T ,
  b r θ 0 = 1 , e j 2 π f r 0 , , e j 2 π M 1 f r 0 T ,
where f t 0 = d sin θ 0 λ 2 r 0 Δ f c and f r 0 = d sin θ 0 λ denote the transmit and receive spatial frequencies, respectively.
Suppose the swarm interference is located in area Θ J , which consists of numerous mutually independent deceptive interferences (as shown the Figure 3). Similarly, for the interference located at ( r j , θ j ) ( ( r j , θ j ) Θ J ,   j = 1 , , J ) , the received signal after matched filtering can be expressed as
x j t = γ j Γ j b r θ j a t r j , θ j ,
where γ j = e j 2 π f 0 τ j and Γ j N M × 1 denotes the j-th complex vector after matched filtering.
In this circumstance, the received signal of the FDA-MIMO radar at time t can be expressed as
x t = x s t + j = 1 J x j t + n t ,
where n t N M × 1 denotes the white Gaussian noise signal with zero mean and covariance matrix σ n 2 I .

3. Analysis of Adaptive Beamforming in Interference Suppression

In order to present the data-independent phase-only beamforming method, we briefly recall the adaptive beamforming method in this section.
According to Equation (12), the echo covariance matrix R , the target signal covariance matrix R s and the interference plus noise covariance matrix R j + n can be expressed respectively as
R = E x t x H t ,
R s = E x s t x s t H = σ s 2 a θ 0 , r 0 a H θ 0 , r 0 ,
R j + n = E j = 1 J x j t + n t j = 1 J x j t + n t H = j = 1 J σ j 2 a θ j , r j a H θ j , r j + σ n 2 I ,
where σ s 2 and σ n 2 represent the energy of the target and noise, respectively, σ j 2 indicates the energy of the j-th interference, and a θ 0 , r 0 = b r θ 0     a t r 0 , θ 0 and a θ j , r j = b r θ j     a t r j , θ j denote the steering vector of the target and the j-th interference, respectively. Then, R can be written as
R = σ s 2 a θ 0 , r 0 a H θ 0 , r 0 + j = 1 J σ j 2 a θ j , r j a H θ j , r j + σ n 2 I .
For a given complex weight vector w = w 1 , , w N M T N M × 1 , the output signal of the radar can be expressed as Y t = w H x t . Hence, the energy of the output signal is given by
E Y t 2 = w H R s w + w H R j + n w ,
where w H R s w and w H R j + n w denote the energy of the output target signal and the interference plus the noise signal after beamforming, respectively. Therefore, the output SINR after beamforming can be defined as
SIN R out = w H R s w w H R j + n w .
In order to maximize the SIN R out , the optimal weight vector w opt is given by
w opt = μ R j + n 1 a θ 0 , r 0 ,
where, μ is a normalization factor that does not affect the SIN R out . In practical applications, R j + n cannot be accurately obtained. It is usually replaced by the sample covariance matrix R ^ , which is given by
R ^ = 1 L l = 1 L x l x H l ,
where L is the number of sampling snapshots. The sample covariance matrix R ^ , however, contains the target signal, which causes the SMI beamformer to be distorted in the target position and leads to a significant decrease in SIN R out . To obtain a better weight vector when including the target signal, the INCMR method was proposed in [11]. Nevertheless, this method still relies on the calculation of the echo covariance matrix. When the number of sampling snapshots is small, it is difficult to obtain the precise INCM information. Furthermore, the phase-only weight vector cannot be obtained, the hardware architecture is complicated and high cost is required when implementing this method.
These discussions motivate us to design a beamforming method in the swarm interference scenario, which needs to satisfy the following conditions: (i) does not depend on INCM acquisition; (ii) achieves interference suppression when only the interference area is known; (iii) maximizes the output SINR; (iv) the weight vector is phase-only.

4. Data-Independent Phase-Only Beamforming

In the previous section, we briefly reviewed the conventional adaptive beamforming techniques. To design the optimal weight vector, this problem can be equivalent to
max w w H R s w w H R j + n w ,
s .   t .   F θ s , r s θ 0 , r 0 ρ , θ s , r s Θ s ,  
F θ j , r j θ 0 , r 0 γ , θ j , r j Θ J ,
where Θ s denotes the area around the target, ρ indicates the lower bound level constraint to guarantee the mainlobe gain, and γ represents the upper bound level constraint to suppress the swarm interference. F θ , r θ 0 , r 0 represents the normalized array response at θ , r , which is defined as
F θ , r θ 0 , r 0 = w H a θ , r 2 w H a θ 0 , r 0 2 .
Therefore, the problem (21) can be considered equivalent to
max w w H R s w w H R j + n w ,
s .   t .   w H a θ s , r s 2 w H a θ 0 , r 0 2 ρ , θ s , r s Θ s ,
w H a θ j , r j 2 w H a θ 0 , r 0 2 γ , θ j , r j Θ J .

4.1. Constant Modulus Constraint

In this subsection, we present a data-independent phase-only beamforming method based on constant modulus constraint (CMC). When the modulus of the weight vector is constrained to h , the problem (23) can be converted into
max w w H R s w w H R j + n w ,
s .   t .   w H a θ s , r s 2 w H a θ 0 , r 0 2 ρ , θ s , r s Θ s ,
w H a θ j , r j 2 w H a θ 0 , r 0 2 γ , θ j , r j Θ J ,
w i = h , i = 1 , , N M .
According to the definition of R s and R j + n , the SIN R out can be expressed as
SINR out = w H R s w w H R j + n w = σ s 2 w H a θ 0 , r 0 2 σ n 2 w 2 2 + j = 1 J σ j 2 w H a θ j , r j 2 = 1 σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 + j = 1 J σ j 2 σ s 2 w H a θ j , r j 2 w H a θ 0 , r 0 2
To maximize the SIN R out , we need to set a very small γ to suppress the interference. Then, the problem of maximizing SIN R out can be equivalent to
max w SIN R out max w w H a θ 0 , r 0 2 .
Proof. 
See Appendix A. □
For (26), we expect max w w H a θ 0 , r 0 2 = 1 to guarantee the mainlobe gain. Therefore, the maximization problem (26) is converted into a minimization problem min w 1 w H a θ 0 , r 0 2 . The problem (24) can be considered equivalent to
min w 1 w H a θ 0 , r 0 2 ,
s .   t .   w H a θ s , r s 2 ρ , θ s , r s Θ s ,
w H a θ j , r j 2 γ , θ j , r j Θ J ,
w i = h , i = 1 , , N M .
We define A θ , r = a θ , r a H θ , r N M × N M and W = w w H N M × N M , and using the semidefinite relaxation (SDR) technique [37], we obtain
w H a θ , r 2 = w H a θ , r a H θ , r w = tr a θ , r a H θ , r w w H = tr A θ , r W
In such a way, (27b) and (27c) can be converted to tr A θ s , r s W ρ and tr A θ j , r j W γ , respectively. Furthermore, Equation (27d) is constrained by W = w w H and w i = p , and it can be written as diag W i = h 2 . Thus, problem (27) can be expressed as
min w 1 tr A θ 0 , r 0 W ,
s .   t .   tr A θ s , r s W ρ , θ s , r s Θ s ,
tr A θ j , r j W γ , θ j , r j Θ J ,
diag W i = h 2 , i = 1 , , N M ,
tr W α ,
W 0 .
where α is a small positive value to prevent the solution from being zero. The optimal W can be obtained by using the CVX toolbox to solve problem (29). Then, let the eigendecomposition of W be
W = i = 1 N M λ i q i q i T .
where λ 1 λ 2 λ N M > 0 denotes the eigenvalue after eigendecomposition, and q i N M × 1 is the eigenvector corresponding to λ i . Then, the optimal weight vector w is given by [38]
w = λ 1 q 1 .
Therefore, the main steps of the proposed CMC algorithm are summarized in Algorithm 1. The computational complexity of the CMC algorithm is mainly derived from the calculation of the convex optimization, and the computational complexity of solving the convex optimization of the N M × 1 dimensional vector is O ( N 3 M 3 ) . Thus, the computational complexity of the CMC algorithm is O ( N 3 M 3 ) .
Algorithm 1: Proposed CMC Algorithm.
1:
Input:  Θ s , Θ J , ρ , γ and h .
2:
Use the convex optimization toolbox to solve problem (29) and obtain W .
3:
Using (30) and (31) yields the optimal weight vector w .
4:
Output: the optimal weight vector w .

4.2. Constant Modulus Decomposition

In the previous subsection, we introduced the CMC algorithm. In this subsection, from a different perspective, we decompose the weight vector into constant modulus complex number and use dual phase receive array to synthesize the phase-only beampattern.
Based on (A8) in Appendix A, we constrain w H a θ 0 , r 0 2 = 1 to guarantee the mainlobe gain, and the problem max w SIN R out can be equivalent to min w w 2 2 . Thus, the problem (23) can be written as
min w w 2 2 ,
s .   t .   w H a θ 0 , r 0 2 = 1 ,
w H a θ s , r s 2 ρ , θ s , r s Θ s ,
w H a θ s , r s 2 γ , θ j , r j Θ J .
Similarly, using the SDR technique, (32) can be transformed into
min w tr W ,
s .   t .   tr A θ 0 , r 0 W = 1 ,  
tr A θ s , r s W ρ , θ s , r s Θ s ,  
tr A θ j , r j W γ , θ j , r j Θ j ,  
tr W α ,  
W 0 .  
Likewise, the problem (33) can be solved using the CVX toolbox to acquire the optimal W . The optimal weight vector w can be obtained through Equations (30) and (31). Notice that the optimal weight vector w is not phase-only, and we intend to convert it into phase-only beamforming by designing a dual-phase shifter receiver. Before that, we present the following lemma.
Lemma 1.
For any complex number p = g e j ω 0 g 2 b , ω and g denote the phase and modulus of the complex number p , respectively; p can be decomposed into
p = b e j α 1 + b e j α 2 ,
where α 1 = ω cos 1 g 2 b and α 2 = ω + cos 1 g 2 b .
Proof 
. See Appendix B. □
The above analysis shows that any complex number can be decomposed into two complexes with equal modulus and unequal phase. This means that we can use a dual-phase shifter at the receiver to synthesize the initial weight vector equivalently. We use Lemma 1 (set b = w m ) to decompose w i , which can be expressed as
w i = b e j α p 1 + b e j α p 2 ,  
where w m = max w 1 , w 2 , , w N M ,
α p 1 = ω i cos 1 w i 2 b ,  
α p 2 = ω i cos 1 w i 2 b ,  
ω i denotes the phase of the complex number w i . After all elements of w have been decomposed, we obtain
w = w p 1 + w p 2 ,  
where w p 1 and w p 2 are phase-only weight vectors. Therefore, we can design the signal processing flow shown in Figure 4 at the receiver.
It should be pointed out that after applying Equation (38), other beamforming methods such as CNBC, SNBC and PBN-BF can also be applied to this array structure. These beamforming methods can be transformed into phase-only beamforming. We implement this in the simulation. The procedure of the constant modulus decomposition (CMD) algorithm is given in Algorithm 2. Similarly, the computational complexity of the CMD algorithm is mainly derived from the calculation of the convex optimization. In addition, the dimensionality of the solved weight vector is N M × 1 . Thus, the computational complexity of the CMD algorithm is O ( N 3 M 3 ) .
Algorithm 2: Proposed CMD Algorithm.
1:
Input :   Θ s ,   Θ J ,   ρ   and   γ .
2:
Using the convex optimization toolbox to solve problem (33) yields W .
3:
Using (30) and (31) yields the optimal weight vector w .
4:
Use (35) and (38) to perform constant mode decomposition of w .
5:
Output: the phase-only weight vectors w p 1   and   w p 2 .
Comparing the CMC and CMD algorithms, we can observe that the CMC algorithm solves the weight vector using the convex optimization with the constant modulus constraint, so that the phase-only weight vector can be acquired directly. It can use the receiver as shown in Figure 2 to synthesize the beampattern. Comparatively, the CMD algorithm achieves the phase-only weight vector by constant modulus decomposition. It requires the dual-phase shifter receiver as shown in Figure 4 to synthesize the beampattern. Furthermore, since the CMC algorithm contains constant-mode constraints, its feasible domain is smaller than the CMD algorithm. Thus, the performance of the CMD algorithm is better than the CMC algorithm, as will be demonstrated in the simulation.

5. Simulation Results

In this section, simulation experiments are performed to substantiate the performance of the proposed methods. The main simulation parameters are listed in Table 1. The convex optimization problem is solved by the CVX toolbox in MATLAB [39].

5.1. Beampattern Comparison for Different Algorithms

In this subsection, we compare the beampattern of different methods to verify the performance of the proposed algorithms. For comparison purposes, the INCMR-linear constraint sector suppressed (LCSS) algorithm [9], SMI-MVDR algorithm [10], SNBC algorithm [29], and CNBC algorithm [29] are simulated. The swarm interference is distributed in an area Θ J , Θ J = θ j , r j 30 θ j 40 , 10   km r j 20   km . The swarm interference power is 30 dB and the target power is 10 dB . The desired beampattern level is expected to be lower than −70 dB at the swarm interference area.
Figure 5 plots the 2D beampattern of different methods. Figure 6a,b shows cross-sectional views of the 2D beampattern of Figure 5 at θ = 35 and r = 45   km , respectively. It can be seen from Figure 5 that both proposed methods can form the regional deep null in the swarm interference area and effectively suppress the swarm interference. Furthermore, the beam shape of the proposed methods is better than the other methods. It can be seen from Figure 6 that the CMD method has the lowest beampattern level in the interference area and less distortion in the uncontrolled area. The mainlobe width of the two proposed methods is the same as the SNBC method but better than the CNBC method.
Moreover, we compare the beampattern obtained by different methods when the swarm interference is distributed in two areas, where Θ J = Θ J 1 Θ J 2 , Θ J 1 and Θ J 2 are given by Θ J 1 = θ j , r j 30 θ j 40 , 10   km r j 20   km , Θ J 2 = 5 θ j 5 , 10   km r j 20   km . Figure 7 shows the 2D beampattern of different methods. It can be seen that under the two swarm interference areas, two proposed methods can form the regional deep null in the interference area. The sidelobe level of the CMD method is lower than the CMC method.

5.2. Comparison of the Output SINR

To verify the performance of the proposed method, in this subsection, we compare the output SINR of different methods at different SNR and the number of snapshots. The swarm interference is distributed in an area Θ J = θ j , r j 30 θ j 40 , 10   km r j 20   km .
Figure 8a shows the output SINR of different beamformers in the SNR from −10 dB to 30 dB when the number of snapshots is 200. Furthermore, under the SNR is 20 dB. Figure 8b shows the output SINR of different beamformers in the number of snapshots from 200 to 1000. In Figure 8a, it can be found that the output SINR increases with the input SNR. From Figure 8b, it can be seen that the proposed methods can obtain excellent output SINR when the number of snapshots is small. In addition, since the proposed methods are data-independent, it can be found that the output SINR does not change with the number of snapshots and maintains good performance. Moreover, as the objective function is to maximize the output SINR, it can be observed that the two proposed methods outperform other methods. The performance of the CMD method is better than the CMC method because its feasible domain is bigger than the CMC algorithm.

5.3. Beampattern on the Different Quantization Bits

In this subsection, we compare the beampattern performance of phase-only methods on different quantization bits. Figure 9 and Figure 10 show the beampattern of the CMC and CMD method when the phase shifter has different quantization bits, respectively. In addition, Figure 11 and Figure 12 plot the beampattern of CNBC and SNBC after constant modulus decomposition on the dual-phase shifter receiver with different quantization bits, respectively. Furthermore, it can be found from Figure 9 and Figure 10 that the synthesis beampattern is significantly different from the original beampattern when the phase shifters have different quantization bits. Figure 11 and Figure 12 demonstrate that the CMD method can be applied to the conventional beamforming method and converted into phase-only beamforming. It can be observed that the beampattern can reach the performance of the original beampattern as the number of quantization bits of the phase shifter increases. When the quantization bit is 3, the beampattern of the interference area is different from the original beampattern, but the mainlobe beampattern remains suitable. With eleven quantization bits, the quantized beampattern is approximate to the original beampattern.

5.4. Output SINR on the Different Quantization Bits

In this subsection, we compare the output SINR of phase-only methods on different quantization bits. Figure 13 and Figure 14 plot the output SINR of the CMC, CMD, Phase-only (PO) CNBC and PO SNBC methods on different quantization bits, respectively. It can be seen that as the number of quantization bits of the phase shifter increases, the output SINR can approximate the original performance. When the quantization bit is 9, the quantized output SINR is very different from the original output SINR. The output SINR can approach the original performance when the quantization bits reach 13.

6. Conclusions

In this paper, we proposed two data-independent phase-only beamforming methods for FDA-MIMO radar to suppress the swarm interference. The two proposed algorithms can suppress swarm interference effectively. Firstly, we convert the output SINR into the data-independent form by imposing the constant modulus constraint and interference area level constraint. The phase-only weight vector is solved by maximizing the output SINR as the objective function. Secondly, from a different perspective, we proposed the CMD beamforming method. This method decomposed the weight vector into constant modulus complex number and used the dual-phase shifter receiver to synthesize the phase-only beampattern. It is worth emphasizing that the CMD method can also be applied to other complex weight beamforming methods to achieve phase-only beampattern synthesis. Simulation experiments demonstrate the superiority of the proposed method. As a future work, we shall consider how to reduce the computational complexity of the proposed methods.

Author Contributions

G.C.: methodology, software, writing—original draft. C.W.: conceptualization, investigation—review & editing. J.G.: formal analysis, funding acquisition, writing—review & editing. M.T.: methodology, supervision, validation, writing—review & editing. Y.L.: validation, writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part by the National Nature Science Foundation of China under Grant 62201580 and by the Natural Science Foundation of Shaanxi Province under Grants 2021JM-222 and 2023-JC-YB-553.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof 
( max SIN R out max w H a θ 0 , r 0 2 ).
Firstly, the maximization problem (25) is converted into a minimization problem:
min w σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 + j = 1 J σ j 2 σ s 2 w H a θ j , r j 2 w H a θ 0 , r 0 2 .
According to constraint (24c), we have
j = 1 J σ j 2 σ s 2 w H a θ j , r j 2 w H a θ 0 , r 0 2 γ j = 1 J σ j 2 σ s 2 .  
Thus, we can obtain
σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 + j = 1 J σ j 2 σ s 2 w H a θ j , r j 2 w H a θ 0 , r 0 2 σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 + γ j = 1 J σ j 2 σ s 2 .
Let us set g w = σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 + j = 1 J σ j 2 σ s 2 w H a θ j , r j 2 w H a θ 0 , r 0 2 , f w = σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 . When g w and f w are respectively taken as objective functions, they have the same feasible domain due to the same constraints. Suppose w is the optimal solution to f w , then we have the inequation
g w f w + γ j = 1 J σ j 2 σ s 2 .
As f w g w , we have
g w f w .  
Subsequently, we can obtain
g w g w f w f w + γ j = 1 J σ j 2 σ s 2 .
Owing to f w f w 0 , Equation (A6) is rewritten as
g w g w γ j = 1 J σ j 2 σ s 2 .  
Therefore, w is the optimal solution that satisfies the error γ j = 1 J σ j 2 σ s 2 in the feasible region of g w . Therefore, f w can be used as the objective function when γ is very small. In beamforming, we usually set a very small γ to suppress the interference. Hence, the objective function (A1) can be equivalent to
min w σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 + j = 1 J σ j 2 σ s 2 w H a θ j , r j 2 w H a θ 0 , r 0 2 min w σ n 2 σ s 2 w 2 2 w H a θ 0 , r 0 2 .  
where σ s 2 σ n 2 denotes the input SNR, which is a constant. Furthermore, w is constrained by (24d), w 2 2 is also a constant. Hence, we can obtain
max SIN R out max w H a θ 0 , r 0 2 .

Appendix B

Proof of Lemma 1.
For different values of g, we can divide it into the following three cases to prove.
(1)
When g = 0 , there is α 1 α 2 = π , and the above formula is obviously established.
(2)
When g = 2 b , there is α 1 = ω , α 2 = ω , and then b e j α 1 + b e j α 2 = g e j ω .
(3)
Consider the common case that 0 < g < 2 b , we plot p = g e j ω in the complex plane, as shown in Figure A1, where O P denotes the complex p and O P = g . After that, find the midpoint Q of the O P such that Q P = O Q = g 2 . We make a vertical line of O P through point Q and intersect the circle of radius b at points A and B, i.e., O P A B . Therefore, we obviously find O P = O A + O B .
According to the geometric relationship in Figure A1, the angle δ between O A and O Q can be given by
δ = cos 1 O Q O A = cos 1 g 2 b .
Therefore, the phase angle of O A is α 1 = ω cos 1 g 2 b , and the phase angle of O B is α 2 = ω + cos 1 g 2 b . Using the fact that O P = O A + O B , we obtain
p = b e j ω cos 1 g 2 b + b e j ω + cos 1 g 2 b ,  
where b can take any value within b g / 2 . In practical implementation, it is more suitable to set b = g / 2 .
Combining the above three cases, Lemma 1 is proved. □
Figure A1. The diagram of complex number decomposition. (a) 0 < g < b ; (b) b g < 2 b .
Figure A1. The diagram of complex number decomposition. (a) 0 < g < b ; (b) b g < 2 b .
Remotesensing 15 01159 g0a1

References

  1. Zheng, Z.; Yang, T.; Wang, H.C. Robust Adaptive Beamforming via Simplified Interference Power Estimation. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 3139–3152. [Google Scholar] [CrossRef]
  2. Zhang, X.; He, Z.; Liao, B.; Zhang, X.; Yang, Y. Pattern Synthesis via Oblique Projection-Based Multipoint Array Response Control. IEEE Trans. Antennas Propag. 2019, 67, 4602–4616. [Google Scholar] [CrossRef]
  3. Ge, Q.; Zhang, Y.; Wang, Y.; Zhang, D. Multi-Constraint Adaptive Beamforming in the Presence of the Desired Signal. IEEE Commun. Lett. 2020, 24, 2594–2598. [Google Scholar] [CrossRef]
  4. Li, Z.; Zhang, Y.; Ge, Q.; Guo, Y. Middle Subarray Interference Covariance Matrix Reconstruction Approach for Robust Adaptive Beamforming with Mutual Coupling. IEEE Commun. Lett. 2019, 23, 664–667. [Google Scholar] [CrossRef]
  5. Zhang, X.; Wang, X.; So, H.C. Linear Arbitrary Array Pattern Synthesis with Shape Constraints and Excitation Range Control. IEEE Antennas Wirel. Propag. Lett. 2021, 20, 1018–1022. [Google Scholar] [CrossRef]
  6. Wang, X.; Amin, M.; Cao, X. Analysis and Design of Optimum Sparse Array Configurations for Adaptive Beamforming. IEEE Trans. Signal Process. 2018, 66, 340–351. [Google Scholar] [CrossRef]
  7. Zhou, X.; Wen, X.; Wang, Z. Swarm of micro flying robots in the wild. Sci. Robot 2022, 7, eabm5954. [Google Scholar] [CrossRef]
  8. Soria, E. Swarms of flying robots in unknown environments. Sci. Robot 2022, 7, eabq2215. [Google Scholar] [CrossRef]
  9. Amar, A.; Doron, M. A linearly constrained minimum variance beamformer with a pre-specified suppression level over a pre-defined broad null sector. Signal Process. 2015, 109, 165–171. [Google Scholar] [CrossRef]
  10. Liu, J.; Liu, W.; Liu, H.; Chen, B.; Xia, X.; Dai, F. Average SINR Calculation of a Persymmetric Sample Matrix Inversion Beamformer. IEEE Trans. Signal Process. 2016, 64, 2135–2145. [Google Scholar] [CrossRef]
  11. Gu, Y.; Leshem, A. Robust Adaptive Beamforming Based on Interference Covariance Matrix Reconstruction and Steering Vector Estimation. IEEE Trans. Signal Process. 2012, 60, 3881–3885. [Google Scholar]
  12. Gu, Y.; Goodman, N.; Hong, S.; Li, Y. Robust adaptive beamforming based on interference covariance matrix sparse reconstruction. Signal Process. 2014, 96, 375–381. [Google Scholar] [CrossRef]
  13. Zhang, X.; He, Z.; Liao, B.; Zhang, X.; Cheng, Z.; Lu, Y. A2RC: An Accurate Array Response Control Algorithm for Pattern Synthesis. IEEE Trans. Signal Process. 2017, 65, 1810–1824. [Google Scholar] [CrossRef]
  14. Zhang, X.; He, Z.; Liao, B.; Zhang, X.; Peng, W. Pattern Synthesis with Multipoint Accurate Array Response Control. IEEE Trans. Antennas Propag. 2017, 65, 4075–4088. [Google Scholar] [CrossRef]
  15. Ai, X.; Gan, L. MPARC: A fast beampattern synthesis algorithm based on adaptive array theory. Signal Process. 2021, 189, 108259. [Google Scholar] [CrossRef]
  16. Zhuang, Y.; Zhang, X.; He, Z. Beampattern Synthesis for Phased Array with Dual-Phase-Shifter Structure. IEEE Trans. Antennas Propag. 2021, 69, 7053–7058. [Google Scholar] [CrossRef]
  17. Zhang, X.; He, Z.; Xia, X.G.; Liao, B.; Zhang, X.; Yang, Y. OPARC: Optimal and Precise Array Response Control Algorithm—Part I: Fundamentals. IEEE Trans. Signal Process. 2019, 67, 652–667. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, X.; He, Z.; Xia, X.G.; Liao, B.; Zhang, X.; Yang, Y. OPARC: Optimal and Precise Array Response Control Algorithm—Part II: Multi-Points and Applications. IEEE Trans. Signal Process. 2019, 67, 668–683. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, X.; He, Z.; Liao, B.; Yang, Y.; Zhang, J.; Zhang, X. Flexible Array Response Control via Oblique Projection. IEEE Trans. Signal Process. 2019, 67, 3126–3139. [Google Scholar] [CrossRef]
  20. Tan, M.; Wang, C.; Li, Z. Correction Analysis of Frequency Diverse Array Radar about Time. IEEE Trans. Antennas Propag. 2021, 69, 834–847. [Google Scholar] [CrossRef]
  21. Wen, C.; Tao, M.; Peng, J.; Wu, J.; Wang, T. Clutter Suppression for Airborne FDA-MIMO Radar using Multi-Waveform Adaptive Processing and Auxiliary Channel STAP. Signal Process. 2019, 154, 280–293. [Google Scholar] [CrossRef]
  22. Lan, L.; Rosamilia, M.; Aubry, A.; Maio, A.D.; Liao, G. Single-Snapshot Angle and Incremental Range Estimation for FDA-MIMO Radar. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 3705–3718. [Google Scholar] [CrossRef]
  23. Wang, W.; Dai, M.; Zheng, Z. FDA Radar Ambiguity Function Characteristics Analysis and Optimization. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 1368–1380. [Google Scholar] [CrossRef]
  24. Xu, J.; Liao, G.; Huang, L. Robust adaptive beamforming for fast-moving target detection with FDA-STAP radar. IEEE Trans. Signal Process. 2017, 65, 973–984. [Google Scholar] [CrossRef]
  25. Wen, C.; Huang, Y.; Peng, J.; Wu, J.; Zheng, G.; Zhang, Y. Slow-Time FDA-MIMO Technique with Application to STAP Radar. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 74–95. [Google Scholar] [CrossRef]
  26. Wen, C.; Ma, C.; Peng, J.; Wu, J. Bistatic FDA-MIMO radar space-time adaptive processing. Signal Process. 2019, 163, 201–212. [Google Scholar] [CrossRef]
  27. Tan, M.; Wang, C.; Xue, B. A Novel Deceptive Jamming Approach against Frequency Diverse Array Radar. IEEE Sen. J. 2021, 21, 8323–8332. [Google Scholar] [CrossRef]
  28. Lan, L.; Xu, J.; Liao, G.; Zhang, Y.; Fioranelli, F.; So, H.C. Suppression of Mainbeam Deceptive Jammer with FDA-MIMO Radar. IEEE Trans. Veh. Technol. 2020, 69, 11584–11598. [Google Scholar] [CrossRef]
  29. Lan, L.; Liao, G.; Xu, J.; Zhang, Y.; Liao, B. Transceive Beamforming with Accurate Nulling in FDA-MIMO Radar for Imaging. IEEE Trans. Geosci. Remote Sens. 2020, 58, 4145–4159. [Google Scholar] [CrossRef]
  30. Aslan, Y.; Puskely, J.; Roederer, A.; Yarovoy, A. Phase-Only Control of Peak Sidelobe Level and Pattern Nulls Using Iterative Phase Perturbations. IEEE Antennas Wirel. Propag. Lett. 2019, 18, 2081–2085. [Google Scholar] [CrossRef]
  31. Wen, C.; Huang, Y.; Peng, J.; Zheng, G.; Liu, W.; Zhang, J. Reconfigurable Sparse Array Synthesis with Phase-Only Control via Consensus-ADMM-Based Sparse Optimization. IEEE Trans. Veh. Technol. 2021, 70, 6647–6661. [Google Scholar] [CrossRef]
  32. Ghayoula, R.; Fadlallah, N.; Gharsallah, A.; Rammal, M. Phase-only adaptive nulling with neural networks for antenna array synthesis. IET Microw. Antennas Propag. 2009, 3, 154–163. [Google Scholar] [CrossRef]
  33. Luyen, T.V.; Giang, T.V.B. Interference suppression of ULA antennas by phase-only control using bat algorithm. IEEE Antennas Wireless Propag. Lett. 2017, 16, 3038–3042. [Google Scholar] [CrossRef]
  34. Zhang, X.; He, Z.; Liao, B.; Zhang, X. Fast Array Response Adjustment with Phase-Only Constraint: A Geometric Approach. IEEE Trans. Antennas Propag. 2019, 67, 6439–6451. [Google Scholar] [CrossRef]
  35. Wen, C.; Peng, J.; Zhou, Y.; Wu, J. Enhanced Three-Dimensional Joint Domain Localized STAP for Airborne FDA-MIMO Radar under Dense False-Target Jamming Scenario. IEEE Sen. J. 2018, 18, 4154–4166. [Google Scholar] [CrossRef]
  36. Chen, G.; Wang, C.; Gong, J.; Tan, M.; Liu, Y. DCFDA-MIMO: A dual coprime FDA-MIMO radar framework for target localization. Digital Signal Process. 2023, 133, 103811. [Google Scholar] [CrossRef]
  37. Luo, Z.; Ma, W.; So, A.; Ye, Y.; Zhang, S. Semidefinite Relaxation of Quadratic Optimization Problems. IEEE Signal Process. Magazine 2010, 27, 20–34. [Google Scholar] [CrossRef]
  38. Fuchs, B. Application of Convex Relaxation to Array Synthesis Problems. IEEE Trans. Antennas Propag. 2014, 62, 634–640. [Google Scholar] [CrossRef] [Green Version]
  39. MATLAB Software for Disciplined Convex Programming; CVX Res., Inc.: San Ramon, CA, USA, 2012.
Figure 1. The FDA-MIMO array structure.
Figure 1. The FDA-MIMO array structure.
Remotesensing 15 01159 g001
Figure 2. The signal processing chain at the receiver array.
Figure 2. The signal processing chain at the receiver array.
Remotesensing 15 01159 g002
Figure 3. Swarm interference model.
Figure 3. Swarm interference model.
Remotesensing 15 01159 g003
Figure 4. Receiver processing flow of the CMD method.
Figure 4. Receiver processing flow of the CMD method.
Remotesensing 15 01159 g004
Figure 5. Two–dimensional beampattern of different methods when the swarm interference is distributed in an area. (a) Initial beampattern, (b) CMC, (c) CMD, (d) SNBC, (e) CNBC, (f) SMI-MVDR, and (g) INCMR-LCSS.
Figure 5. Two–dimensional beampattern of different methods when the swarm interference is distributed in an area. (a) Initial beampattern, (b) CMC, (c) CMD, (d) SNBC, (e) CNBC, (f) SMI-MVDR, and (g) INCMR-LCSS.
Remotesensing 15 01159 g005
Figure 6. The beampattern cross–sectional of different methods at target range and angle. (a) The beampattern on θ = 35 . (b) The beampattern on r = 45   km .
Figure 6. The beampattern cross–sectional of different methods at target range and angle. (a) The beampattern on θ = 35 . (b) The beampattern on r = 45   km .
Remotesensing 15 01159 g006
Figure 7. Two–dimensional beampattern of different methods when the swarm interference is distributed in two areas. (a) CMC, (b) CMD, (c) SNBC, (d) CNBC, (e) SMI-MVDR, and (f) INCMR–LCSS.
Figure 7. Two–dimensional beampattern of different methods when the swarm interference is distributed in two areas. (a) CMC, (b) CMD, (c) SNBC, (d) CNBC, (e) SMI-MVDR, and (f) INCMR–LCSS.
Remotesensing 15 01159 g007
Figure 8. The output SINR of different methods. (a) The output SINR versus SNR. (b) The output SINR versus the number of snapshots.
Figure 8. The output SINR of different methods. (a) The output SINR versus SNR. (b) The output SINR versus the number of snapshots.
Remotesensing 15 01159 g008
Figure 9. The CMC beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Figure 9. The CMC beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Remotesensing 15 01159 g009
Figure 10. The CMD beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Figure 10. The CMD beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Remotesensing 15 01159 g010
Figure 11. The constant modulus decomposition CNBC beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Figure 11. The constant modulus decomposition CNBC beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Remotesensing 15 01159 g011
Figure 12. The constant modulus decomposition SNBC beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Figure 12. The constant modulus decomposition SNBC beampattern on different quantization bits: (a) 3 bits, (b) 7 bits, (c) 11 bits, (d) cross–section of different quantization at target angle.
Remotesensing 15 01159 g012
Figure 13. The output SINR versus SNR on different quantization bits.
Figure 13. The output SINR versus SNR on different quantization bits.
Remotesensing 15 01159 g013
Figure 14. The output SINR versus the number of snapshots on different quantization bits.
Figure 14. The output SINR versus the number of snapshots on different quantization bits.
Remotesensing 15 01159 g014
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersSymbolsValueParametersSymbolsValue
Transmit elements N 15Target range r 0 45 km
Receive elements M 15Target angle θ 0 35°
Carrier frequency f 0 10 GHzFrequency offset Δ f 1500 Hz
Element spacing d 0.015 mCoefficient α 0.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, G.; Wang, C.; Gong, J.; Tan, M.; Liu, Y. Data-Independent Phase-Only Beamforming of FDA-MIMO Radar for Swarm Interference Suppression. Remote Sens. 2023, 15, 1159. https://doi.org/10.3390/rs15041159

AMA Style

Chen G, Wang C, Gong J, Tan M, Liu Y. Data-Independent Phase-Only Beamforming of FDA-MIMO Radar for Swarm Interference Suppression. Remote Sensing. 2023; 15(4):1159. https://doi.org/10.3390/rs15041159

Chicago/Turabian Style

Chen, Geng, Chunyang Wang, Jian Gong, Ming Tan, and Yibin Liu. 2023. "Data-Independent Phase-Only Beamforming of FDA-MIMO Radar for Swarm Interference Suppression" Remote Sensing 15, no. 4: 1159. https://doi.org/10.3390/rs15041159

APA Style

Chen, G., Wang, C., Gong, J., Tan, M., & Liu, Y. (2023). Data-Independent Phase-Only Beamforming of FDA-MIMO Radar for Swarm Interference Suppression. Remote Sensing, 15(4), 1159. https://doi.org/10.3390/rs15041159

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop