Next Article in Journal
Tracking Ground Targets with a Road Constraint Using a GMPHD Filter
Next Article in Special Issue
EEG-Based Emotion Recognition Using Quadratic Time-Frequency Distribution
Previous Article in Journal
IMU Signal Generator Based on Dual Quaternion Interpolation for Integration Simulation
Previous Article in Special Issue
A Gyroscope Bias Estimation Algorithm Based on Map Specific Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Design of Space-Time Transmit and Receive Weights for Colocated MIMO Radar

1
School of Electronic and Information Engineering, Beihang University, Beijing 100083, China
2
Department of Electronic & Electrical Engineering, University of Sheffield, Sheffield S1 4ET, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 27 June 2018 / Revised: 12 August 2018 / Accepted: 13 August 2018 / Published: 18 August 2018
(This article belongs to the Special Issue Sensor Signal and Information Processing)

Abstract

:
Compared with single-input multiple-output (SIMO) radar, colocated multiple-input multiple-output (MIMO) radar can detect moving targets better by adopting waveform diversity. When the colocated MIMO radar transmits a set of orthogonal waveforms, the transmit weights are usually set equal to one, and the receive weights are adaptively adjusted to suppress clutter based on space-time adaptive processing technology. This paper proposes the joint design of space-time transmit and receive weights for colocated MIMO radar. The approach is based on the premise that all possible moving targets are detected by setting a lower threshold. In each direction where there may be moving targets, the space-time transmit and receive weights can be iteratively updated by using the proposed approach to improve the output signal-to-interference-plus-noise ratio (SINR), which is helpful to improve the precision of target detection. Simulation results demonstrate that the proposed method improves the output SINR by greater than 13 dB.

Graphical Abstract

1. Introduction

For colocated multiple-input multiple-output (MIMO) radar, the antennas are close to each other, and each antenna serves as a transceiver. Compared with single-input multiple-output (SIMO) radar, colocated MIMO radar can achieve better detection of moving targets by adopting waveform diversity [1,2,3].
To improve the performance of clutter mitigation, orthogonal waveforms are widely used in colocated MIMO radar [4,5]. Each transceiver transmits one unique and orthogonal waveform, and then receives and separates all the return echoes by using the appropriate matched filter bank. Compared with SIMO radar, more virtual phase centers are formed, and more degrees of freedom (DOFs) are achieved. Combined with space-time adaptive processing (STAP), MIMO radar with orthogonal waveforms can obtain better clutter mitigation performance and identify more target parameters. To realize good correlation properties of the transmitted signals, the common solution is to solve an optimization problem with constraints. The existing orthogonal waveform design algorithms can be divided into two categories. The first, which includes Cross Entropy [6], simulated annealing [7], and the Genetic Algorithm algorithms [8,9], optimizes the polyphase code sequences to achieve low aperiodic auto-correlation sidelobe peaks and cross-correlation peaks. The second, which includes the efficient cyclic algorithm [10] and the monotonic minimizer for integrated sidelobe level [11] and the coordinate-descent algorithm [12], optimizes the unimodular sequences to achieve a good integrated or peak sidelobe level of the auto-correlation function.
Although they provide more DOFs, orthogonal waveforms induce the omnidirectionality of the antenna, which reduces the signal-to-noise ratio (SNR). One solution is to adopt beamforming techniques [13,14,15,16,17], which design partial coherent transmitted signals by optimizing the signal covariance matrix with the constraint of the desired beampattern. In order to improve the signal-to-interference-plus-noise ratio (SINR), the joint design of the transmit signal and receive filter has been proposed [18,19,20,21,22]. In [18], constrained optimization procedures are devised based on phase-only waveforms to sequentially improve the SINR. In [19], the priori knowledge about geographical information is used in the optimization. In [20], it is assumed the target Doppler frequency is unknown and the worst-case SINR at the output of the filter bank is considered in the robust joint design. For non-cooperative radar networks, design of coded coherent waveforms has been presented to improve the SINR of each active radar by resorting to the theory of potential games [21]. For MIMO radar, the design of a space-time transmit code and space-time receive filter has been presented to improve the worst case SINR by optimizing the coded coherent signals [22]. For colocated MIMO radar, transmitted waveform covariance matrix has been optimized to focus the transmit beampattern into the target direction to improve the SINR [23]. The design of a space-time transmit code and receive filter has been presented to improve the SINR by optimizing the partial coherent signals under different constraints, such as a similarity constraint and a constant modulus requirement [24], integrated sidelobe level and peak sidelobe level constraints at the pulse compression output [25], or for generating BPSK waveforms [26].
The existing literature about improving the output SINR for colocated MIMO radar focuses on optimizing the covariance matrix of partially coherent signals, which reduces the dimension of the measurement space and results in the loss of abilities to suppress high-rank clutter. To preserve the full spatial degrees of freedom and improve the output SINR, a joint design of transmit and receive weights is proposed for colocated MIMO radar using orthogonal waveforms, and an iterative framework is determined to update the weights in real time. Initially, all possible moving targets are detected by setting a lower threshold, which may cause high false alarm rate. Then in each direction where there may be moving targets, the space-time transmit and receive weights are iteratively updated by exploiting the clutter covariance matrix and moving target status including locations and velocities. Through iteration, the output SINR is gradually increased, which is beneficial to improve detection performance and reduce the false alarm rate.
This paper is structured as follows. Section 2 introduces the signal model for colocated MIMO radar, and illustrates the optimization problem with constraints on the transmit and receive weights for improving the output SINR. Section 3 presents the joint design method of space-time transmit and receive weights. By exploiting the information about the dynamic environment contained in the return echoes, the transmit-receive weights can be updated in real time. The simulation results are demonstrated and analyzed in Section 4, and Section 5 concludes the paper.

2. Optimization Model for Colocated MIMO Radar

Figure 1a shows the observation geometry. A right-handed Cartesian coordinate system is established, where the origin is located at the center of the observed area and the x-axis, y-axis, and z-axis represent the broadside direction, flight direction, and height direction, respectively. The radar moves along the y-axis. A uniform linear array with N identical antennas is mounted on the radar. The spacing between the two adjacent antennas is half of a wavelength. A target moves at a speed of v t , and v r a d denotes the radial component of v t . v s denotes the platform speed. φ and θ are the azimuth and depression angles, respectively. In this paper, “speed” is loosely used to mean “magnitude of the velocity vector”. According to [27], angles and speeds are defined as follows:
  • Platform speed: This is the speed of the platform along the flight path.
  • Target speed: This is the speed of one moving target along the moving direction, which can be any direction in the free space.
  • Target radial speed: This is the speed of one moving target along the line-of-sight direction.
  • Nadir: The nadir is the point on the area directly below the radar phase center, so that the “normal” to the area at the nadir passes through the radar phase center.
  • Radar track: As the nadir point moves along the area, it traces out the radar track.
  • Depression angle: This is the remaining angle of the off-line angle, which is the angle from the target to the nadir relative to the radar phase center.
  • Azimuth angle: This is the clockwise angle from the direction of the radar track to the direction of interest relative to the nadir point viewed from above (i.e., projected to the ground plane).
Usually in radar systems, the echoes corresponding to a group of multiple pulses are coherently processed, and a declaration is made whether targets exist or not. The duration of these multiple pulses is called the coherent processing interval (CPI). Figure 1b describes the transmit and receive processes. In each CPI, each antenna transmits the orthogonal pulse train consisting of M rectangular pulses, which are p n ( t m T ) , m = 1 , 2 , , M , where T denotes the pulse repetition interval, n = 1 , 2 , , N represents the n-th antenna, and
T p / 2 T p / 2 p i ( t ) p j * ( t ) d t = { 1 i = j 0 e l s e    i = 1 , , N ;   j = 1 , , N
where T p is the pulse width and ( ) denotes the conjugate operator. At the m-th transmission moment, all the pulses are weighted by the space time transmit weights ω t ( m , n ) to form the signal s t , m , n ( t ) which is radiated in the free space. s t , m , n ( t ) equals
s t , m , n ( t ) = ω t ( m , n ) p n ( t m T ) exp [ j 2 π f c ( t m T ) ]
where f c is the carrier frequency, all the space time transmit weights ω t ( m , n ) will form the space time transmit weight vector
ω t = [ ω t ( 1 , 1 ) , ω t ( 2 , 1 ) , , ω t ( M , 1 ) , ω t ( 1 , 2 ) , ω t ( 2 , 2 ) , , ω t ( M , N ) ] T
The return echoes reflected from the clutters and moving target are received by each antenna and pass through the matched filter bank. The filter bank is composed of N different matched filters, and each filter corresponds to one transmit waveform. Thus, the echoes are separated. By applying the space-time receive weights ω r , the clutter is adaptively suppressed and moving targets are detected.
To focus on the optimization of transmit-receive weights, the simplified assumptions are adopted in the paper, which are not ambiguous returns as well as not returns from adjacent range cells. After demodulation and matched filtering, the target’s echo and the clutter’s echoes received by the n-th antenna transmitted by the i-th antenna at the m-th pulse repetition interval can be expressed as
x t , m , i , n ( ω t ) = k = 1 K A k ( t ) ω t ( m , i ) exp [ j 2 π λ ( m 1 ) T ( 2 v s cos φ k ( t ) cos θ + 2 v r a d , k ) ] × exp [ j 2 π λ ( ( i 1 ) λ 2 cos θ cos φ k ( t ) + ( n 1 ) λ 2 cos θ cos φ k ( t ) ) ] , m = 1 , , M i , n = 1 , , N
and
x c , m , i , n ( ω t ) = φ = 0 2 π A ( c ) ( φ ) ω t ( m , i ) exp [ j 2 π λ ( m 1 ) T ( 2 v s cos φ cos θ ) ] × exp [ j 2 π λ ( ( i 1 ) λ 2 cos θ cos φ + ( n 1 ) λ 2 cos θ cos φ ) ] d φ , m = 1 , , M i , n = 1 , , N
As there is a one-to-one functional mapping relationship between the depression angle θ and some range gate [12], and each range gate is a minimum range sampling unit, the independence of return echoes on θ has been omitted to simplify the expressions. The vector format of the target’s echo and clutter’s echoes can be expressed as
x t ( ω t ) = [ x t , 1 , 1 , 1 , x t , 2 , 1 , 1 , , x t , M , 1 , 1 , x t , 1 , 2 , 1 , x t , 2 , 2 , 1 , , x t , M , N , N ] T x c ( ω t ) = [ x c , 1 , 1 , 1 , x c , 2 , 1 , 1 , , x c , M , 1 , 1 , x c , 1 , 2 , 1 , x c , 2 , 2 , 1 , , x c , M , N , N ] T
In Equations (3) and (4), ω t = [ ω t ( 1 , 1 ) , ω t ( 1 , 2 ) , , ω t ( M , N ) ] T , K is the number of targets and equals to 1 in this paper, λ is the wavelength, v s is the platform speed, θ is the depression angle, i denotes the i-th transmit antenna, and n denotes the n-th receive antenna. A k ( t ) , v r a d , k and φ k ( t ) are the complex amplitude, radial velocity, and relative azimuth angle of the k-th target, respectively. A ( c ) ( φ ) is the clutter complex amplitude and is assumed to be characterized by a circularly symmetric complex Gaussian distribution [28,29].
In the temporal domain, the angular frequency is a reflection of the rate of phase change. Its counterpart in spatial domain can be defined similarly [30]. The normalized spatial frequency is a reflection of the rate of phase change with the spatial antennas in Hertz, and can be defined as f s = d cos θ cos φ λ , where d is the spacing between two adjacent antennas which is half of a wavelength in this paper, therefore the normalized spatial frequency is f s = cos θ cos φ 2 . Let s s f and s t d denote the space steering vector and time steering vector, which represent the spatial phase and temporal phase in the return echoes, respectively. They are defined as
s s f = [ 1 exp ( j 2 π f s ) exp ( j 2 π f s ( N 1 ) ) ] T
and
s t d = [ 1 exp ( j 2 π f d ) exp ( j 2 π f d ( M 1 ) ) ] T
where f s and f d denote the normalized spatial frequency and the normalized Doppler frequency, respectively, and ( ) Τ denotes the transpose operator.
Define s t = s s f s t d and s r ( ω t ) = s s f ( ω t s t ) as the transmit and receive space-time steering vectors, respectively, where is the Hadamard product and is the Kronecker product. s t represents the unweighted echo phase term vector for the transmit. s r ( ω t ) represents the weighted echo reflected by one target whose reflection amplitude is one. Therefore, the echo from one range gate consisting of moving targets is expressed as
x ( ω t ) = k = 1 K A k ( t ) s r ( t , k ) ( ω t ) + n + 0 2 π A ( c ) s r ( c ) ( φ , ω t ) d φ  
where n is the vector of circularly symmetric white Gaussian noise samples with zero-mean and covariance matrix σ n 2 [25], and s r ( t , k ) ( ω t ) and s r ( c ) ( φ , ω t ) are the receive space-time steering vectors for the k-th target and clutter, respectively.
To improve the SINR, the interference-plus-noise power after receive weighting should be minimized under the condition that the target power is preserved. Therefore, the following optimization model should be satisfied:
min ω t , ω r ω r H Q q ( ω t ) ω r , s . t . ω r H s r ( ω t ) = 1 ω t H s t = g ,
where
Q q ( ω t ) = E [ ( n + 0 2 π A ( c ) s r ( c ) ( φ , ω t ) d φ ) ( n + 0 2 π A ( c ) s r ( c ) ( φ , ω t ) d φ ) H ] ,
E [ ] and ( ) H are the expectation and the conjugate transpose operators, respectively, and ω r denotes the space-time receive weights. g = 1 M N H s t , where 1 M N is a column vector of all-ones with the dimension of M N . When radar starts to illuminate some area, there is no target information as well as clutter information, then transmit weights will be set to 1 to achieve the original clutter suppression ability. g is to keep the radiation power reaching the target constant, and will hold constant during the illumination time as g = 1 M N H s t ,where s t indicates the target’s direction where the received radiation power will keep constant.
The objective function in Equation (8) is meant to minimize the interference-plus-noise output power. The first constraint indicates that the target power after receive weighting remains constant and the second constraint guarantees that the radiation power reaching the target is also constant. The space-time transmit and receive weights will be jointly and adaptively designed to satisfy Equation (8).

3. Joint Design of Transmit and Receive Weights

Joint design of transmit and receive weights is implemented iteratively, which is illustrated in Figure 2. Before starting to illuminate the observation area, the radar system works in the passive mode to records noise data [31,32]. And the noise covariance matrix can be obtained as
Q n = 1 L l = 1 L n l n l H .
where n l represents the space-time data of the l -th range gate, and L is the number of range gates [12]. The noise covariance matrix is assumed to be unchanged during illumination.
In the first CPI, the transmitter sends out orthogonal waveforms weighted by the initial space-time transmit weights (STTWs) which are set to all ones. The clutter analyzer acquires the clutter pulse noise covariance matrix by processing the return echoes. If the echo in the l 0 -th range gate is processed, the clutter-plus-noise covariance matrix can be obtained by analyzing echoes in other range gates
Q q = 1 L 1 l = 1 , l l 0 L x l x l H .
where x l represents the echo in the l -th range gate. Since the clutter and noise are statistically independent, the clutter covariance matrix can be achieved Q c = Q q Q n .
The space-time-receive weights generator optimally designs the space-time receive weights (STRWs) based on the space time adaptive technology [33,34]. In the moving targets analyzer, the echoes are weighted by the STRWs to suppress the clutter, and a lower threshold is set to detect moving targets. Although the lower threshold results in a higher false alarm rate, it is guaranteed that no targets are missed. Then the locations and velocities of each possible target are extracted.
In the subsequent CPIs, the system works in a closed-loop mode. For one direction where there may be a target of interest, the STTWs and STRWs are adaptively updated by applying the proposed approach. From one CPI to the next CPI, the interference-plus-noise output power is suppressed more and more, which leads to an improvement of the SINR. The moving targets analyzer raises the threshold with the SINR to better extract the information of the target. When the difference of SINRs between adjacent CPIs is little, another direction is selected to perform the closed-loop.
If the information about clutter and target is known to the system designer, the optimization problem about the transmit-receive weights in Equation (8) will have optimal solutions, since the optimization problem is a quadratic objective function with two linear constraints whose hessian matrix is positive definite [35]. If the statistic characteristics of clutter is stationary and the target keeps the same motion state, the proposed method will converge, which is guaranteed by the fact that the optimization problem listed in Equation (8) has an optimal solution. The transmit-receive weights are designed when the proposed iterative method converges and will be the optimal solution mentioned above.
In the joint design, the STRWs are calculated based on the deterministic STTWs. Therefore, Equation (8) can be simplified as
min ω r ω r H Q q ω r , s . t . ω r H s r = 1 ,
where only STRWs remain in the constraint. The STRWs can be obtained by applying the linear weighting as [34]
ω r = Q q 1 s r s r H Q q 1 s r .
By solving the optimization problem (8) based on the Lagrange multiplier method (see Appendix A), the optimal STRWs can be obtained
ω t = ( [ ( D 22 g D 12 ) Q c , r 1 ( 1 M N ) s 1 + ( g D 11 D 21 ) Q c , r 1 ( 1 M N ) s 2 ] D 11 D 22 D 12 D 21 ) ,
where
s 1 = [ ( i = 1 N s s f , i ω r , i ) s t ]
s 2 = s t
and
D 11 = s 1 H Q c , r 1 ( 1 M N ) s 1 D 12 = s 1 H Q c , r 1 ( 1 M N ) s 2 D 21 = s 2 H Q c , r 1 ( 1 M N ) s 1 D 22 = s 2 H Q c , r 1 ( 1 M N ) s 2 .
In Equations (14) and (15), s s f , i = exp ( j 2 π f s ( i 1 ) ) represents the i-th element of the space steering vector; ω r , i represents the i-th sub vector of the STRWs vector ω r = [ ω r , 1 Τ , ω r , 2 Τ , , ω r , N Τ ] Τ ; Q c , r ( 1 M N ) = F Q c ( 1 M N ) F H , where F = [ D i a g { ω r , 1 } , D i a g { ω r , 2 } , , D i a g { ω r , N } ] , and D i a g { a } denotes the diagonal matrix whose diagonal elements are the elements of vector a ; and the clutter covariance matrix Q c ( 1 M N ) equals
Q c ( 1 M N ) = Q q ( 1 M N ) Q n ,
where Q q ( 1 M N ) and Q n denote the clutter pulse noise covariance matrix and the noise covariance matrix, respectively. Q q ( 1 M N ) can be estimated from the secondary data which contains only clutter and noise [36]. Q n can be estimated from the collected noise data when the radar system works in the passive mode.

4. Simulation and Analysis

4.1. Simulation Parameters

Simulation parameters are listed in Table 1, which refer to the standard parameters set in [34]. The radar platform flies parallel to the ground at a uniform speed, while the target moves perpendicular to the flight direction, as shown in Figure 1. Additionally, four orthogonal transmit waveforms generated by the genetic algorithm [8] and corresponding receive filters are used in the simulation. The true interference and target parameters are assumed perfectly known in the simulation results. To valid the proposed method, the conventional method serves as a standard performance where the STTWs are all set equal to one and STTWs are optimized by the space-time adaptive technology.

4.2. Simulation Results and Analysis

This paper adopts orthogonal transmit waveforms to preserve high DOFs and joint design of transmit-receive weights to improve the output SINR. Here, the improvement factor (IF) is used to evaluate the performance of the colocated MIMO radar and is defined as the ratio of the output SINR to the input SINR [34], i.e.,
I F = ω r H ( f S T , f S R , f D ) s r ( f S T , f S R , f D ) s r H ( f S T , f S R , f D ) ω r ( f S T , f S R , f D ) ω r H ( f S T , f S R , f D ) Q q ω r ( f S T , f S R , f D ) × t r ( Q q ) s r H ( f S T , f S R , f D ) s r ( f S T , f S R , f D ) .
According to Equation (17), IF is a function of the space transmit frequency f S T , space receive frequency f S R , and Doppler frequency f D . Usually, normalized frequencies are more convenient and are thus commonly used in analysis. The normalized space transmit and receive frequencies are defined as the space transmit and receive frequencies normalized to the reciprocal of the corresponding sensor spacing. The normalized Doppler frequency is defined as the Doppler frequency normalized to the pulse repetition frequency.
Figure 3 demonstrates IFs achieved by the proposed method and the conventional method. Figure 3a is the result in the second CPI. In the first CPI, the STTWs are all set equal to one, and the STRWs are determined by (13). In the second CPI, the STTWs are calculated according to (14) and STRWs are calculated according to (13), which indicates that STTWs and STRWs are updated in the iteration. Then IF is calculated by (17) and Figure 3a is generated according to IF. In Figure 3b, the IF is achieved by setting STTWs to one and calculating STRWs according to (13). It indicates that Figure 3b is the same as IF achieved by the proposed method in the first CPI.
The distributions of IFs in Figure 3a,b are similar. In each figure, the narrow trench is formed along the clutter trajectory, which indicates that the clutter can be suppressed. Everywhere outside the clutter trench, the IF plateau exists, where detection of moving targets is optimum. The difference between the two figures is that the IF plateau of the proposed method is higher than that of the conventional method, which means that the proposed method has better performance. Figure 4 illustrates this issue more clearly by displaying three profiles of the stereograms in Figure 3. The variations of IF with normalized frequencies show that the minimums of the notches achieved by the two methods are nearly the same. Outside the notches, the IF achieved by the proposed method is more than 13 dB better than that of the conventional method.
Figure 5 demonstrates the IF curves with respect to the normalized Doppler frequency as the number of CPIs increases. Initially, the STTWs are all set to one in the first CPI. Therefore, the IF curve depending on the first CPI also represents the performance of the conventional method. Once the radar system has received echoes, the STTWs are adaptively updated by exploiting the information contained in the previous CPI. From one CPI to the next CPI, the interference-plus-noise output power is suppressed more and more, which leads to improvement of the SINR. The moving targets analyzer raises the threshold with the SINR to better extract the information of the target. Therefore, the IF curve improves as the CPI increases. Thus, a few CPIs are needed to allow the system to reach the optimal state. As shown in Figure 5, the proposed method approximately converges after the fifth CPI. Compared with the conventional method, the improvement of IF is greater than 13 dB, which is equivalent to the improvement of the output SINR under the condition that the input SINR is constant.

5. Conclusions

This paper has proposed a joint design of space-time transmit and receive weights. By exploiting the information of the clutter and target in the return echoes, the output SINR of colocated MIMO radar can be improved. Moreover, because of rapid convergence, the proposed approach is suitable for real-time applications. Although the proposed approach is designed for the single moving target, it can be extended to more complicated situations. Future research will focus on designing robust strategies following and accounting for target steering mismatches that can arise due to multipath and phase noise [37,38,39]. And the upper bound to the IF will also be studied.

Author Contributions

All authors contributed equally to this paper in ideas, design of methodology and writing. Additionally, Z.Y. was responsible for project administration, data analysis and experiments validation. S.W. accomplished software programming, data analysis and presentation.

Funding

This research was funded by China Scholarship Council grant number 201706025038.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Equation (8) is rewritten as follows:
min ω t ω r H Q q ( ω t ) ω r , s . t . ω r H s r ( ω t ) = 1 ω t H s t = g .
However, the objective function and the first constraint are implicit functions of ω t . Therefore, to obtain the solution to ω t , they must be changed.
The objective function equals
min ω t ω r H Q q ( ω t ) ω r = min ω t ω r H Q c ( ω t ) ω r = min ω t ω r H E [ ( φ = 0 2 π A ( c ) ( φ ) s r ( c ) ( φ , ω t ) d φ ) × ( φ = 0 2 π A ( c ) ( φ ) s r ( c ) ( φ , ω t ) d φ ) H ] ω r = min ω t ω r H E [ ( φ = 0 2 π A ( c ) ( φ ) s r ( c ) ( φ , 1 M N ) ( 1 N ω t ) d φ ) × ( φ = 0 2 π A ( c ) ( φ ) s r ( c ) ( φ , 1 M N ) ( 1 N ω t ) d φ ) H ] ω r = min ω t ω r H D i a g { 1 N ω t } E [ ( φ = 0 2 π A ( c ) ( φ ) s r ( c ) ( φ , 1 M N ) d φ ) × ( φ = 0 2 π A ( c ) ( φ ) s r ( c ) ( φ , 1 M N ) d φ ) H ] D i a g { 1 N ω t } ω r = min ω t ω r H D i a g { 1 N ω t } E [ x c ( 1 M N ) x c H ( 1 M N ) ] D i a g { 1 N ω t } ω r = min ω t ( 1 N ω t ) Τ D i a g { ω r } Q c ( 1 M N ) D i a g { ω r } ( 1 N ω t ) = min ω t ω t Τ Q c , r ( 1 M N ) ω t ,
where
Q c , r ( 1 M N ) = F Q c ( 1 M N ) F H F = [ D i a g { ω r , 1 } , D i a g { ω r , 2 } , , D i a g { ω r , N } ] Q c ( ω t ) = E [ ( φ = 0 2 π A ( c ) s r ( c ) ( φ , ω t ) d φ ) ( φ = 0 2 π A ( c ) s r ( c ) ( φ , ω t ) d φ ) H ] .
The superscripts ( ) Τ and ( ) denote the transpose and conjugate operators, respectively; 1 N and 1 M N are column vectors of all-ones with dimensions of N and M N , respectively; and D i a g { a } denotes the diagonal matrix whose diagonal elements are the elements of the vector a . The STRWs vector can be expressed as ω r = [ ω r , 1 Τ , ω r , 2 Τ , , ω r , N Τ ] Τ , and ω r , i represents the i-th sub vector with the dimension of M N . Additionally, s s f , i = exp ( j 2 π f s ( i 1 ) ) represents the i-th element of the space steering vector.
The first constraint can be expressed as
1 = ω r H s r ( ω t ) = ω r H [ s s f ( s t ω t ) ] = i = 1 N s s f , i ω r , i H ( s t ω t ) = [ ( i = 1 N s s f , i ω r , i H ) s t Τ ] ω t .
Let s 1 = [ ( i = 1 N s s f , i ω r , i ) s t ] and s 2 = s t . Then the optimization problem becomes
min ω t ω t T Q c , r ( 1 M N ) ω t , s . t . s 1 H ω t = 1 s 2 H ω t = g .
By using the Lagrange multiplier method, the solution is
ω t = ( [ ( D 22 g D 12 ) Q c , r 1 ( 1 M N ) s 1 + ( g D 11 D 21 ) Q c , r 1 ( 1 M N ) s 2 ] D 11 D 22 D 12 D 21 ) ,
where
D 11 = s 1 H Q c , r 1 ( 1 M N ) s 1 D 12 = s 1 H Q c , r 1 ( 1 M N ) s 2 D 21 = s 2 H Q c , r 1 ( 1 M N ) s 1 D 22 = s 2 H Q c , r 1 ( 1 M N ) s 2 .

References

  1. Forsythe, K.W.; Bliss, D.W.; Fawcett, G.S. Multiple-input multiple-output (MIMO) radar: Performance issues. In Proceedings of the 38th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 7–10 November 2004; Volume 1, pp. 310–315. [Google Scholar]
  2. Li, J.; Stoica, P. MIMO radar with colocated antennas. IEEE Signal Process. Mag. 2007, 24, 106–114. [Google Scholar] [CrossRef]
  3. Stoica, P.; Li, J.; Xie, Y. On probing signal design for MIMO radar. IEEE Trans. Signal Process. 2007, 55, 4151–4161. [Google Scholar] [CrossRef]
  4. Bliss, D.W.; Forsythe, K.W. Multiple-input multiple-output (MIMO) radar and imaging: Degrees of freedom and resolution. In Proceedings of the 37th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 9–12 November 2003; Volume 1, pp. 54–59. [Google Scholar]
  5. Robey, F.C.; Coutts, S.; Weikle, D.; McHarg, J.C.; Cuomo, K. MIMO radar theory and experimental results. In Proceedings of the 38th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 7–10 November 2004; Volume 1, pp. 300–304. [Google Scholar]
  6. Khan, H.A.; Zhang, Y.; Ji, C.; Stevens, C.J.; Edwards, D.J.; O’Brien, D. Optimizing Polyphase Sequences for Orthogonal Netted Radar. IEEE Signal Process. Lett. 2006, 13, 589–592. [Google Scholar] [CrossRef]
  7. Deng, H. Polyphase code design for Orthogonal Netted Radar systems. IEEE Trans. Signal Process. 2004, 52, 3126–3135. [Google Scholar] [CrossRef]
  8. Liu, B.; He, Z.; Zeng, J.; Liu, B. Polyphase orthogonal code design for MIMO radar systems. In Proceedings of the 2006 CIE International Conference on Radar, Shanghai, China, 16–19 October 2006; pp. 1–4. [Google Scholar]
  9. Liu, B.; He, Z.; He, Q. Optimization of Orthogonal Discrete Frequency-Coding Waveform Based on Modified Genetic Algorithm for MIMO Radar. In Proceedings of the 2007 International Conference on Communications, Circuits and Systems, Kokura, Japan, 11–13 July 2007; pp. 966–970. [Google Scholar]
  10. Stoica, P.; He, H.; Li, J. New Algorithms for Designing Unimodular Sequences with Good Correlation Properties. IEEE Trans. Signal Process. 2009, 57, 1415–1425. [Google Scholar] [CrossRef]
  11. Song, J.; Babu, P.; Palomar, D.P. Optimization Methods for Designing Sequences with Low Autocorrelation Sidelobes. IEEE Trans. Signal Process. 2015, 63, 3998–4009. [Google Scholar] [CrossRef]
  12. Aubry, A.; Maio, A.D.; Kerahroodi, M.A.; Naghsh, M.M.; Modarres-Hashemi, M. A Coordinate-Descent Framework to Design Low PSL/ISL Sequences. IEEE Trans. Signal Process. 2017, 65, 5942–5956. [Google Scholar] [Green Version]
  13. Fuhrmann, D.R.; Antonio, G.S. Transmit beamforming for MIMO radar systems using partial signal correlation. In Proceedings of the 38th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 7–10 November 2004; Volume 1, pp. 295–299. [Google Scholar]
  14. Stoica, P.; Li, J.; Zhu, X. Waveform synthesis for diversity-based transmit beampattern design. IEEE Trans. Signal Process. 2008, 56, 2593–2598. [Google Scholar] [CrossRef]
  15. Wang, Y.C.; Wang, X.; Liu, H.; Luo, Z.Q. On the design of constant modulus probing signals for MIMO radar. IEEE Trans. Signal Process. 2012, 60, 4432–4438. [Google Scholar] [CrossRef]
  16. Jiu, B.; Liu, H.; Wang, X.; Zhang, L.; Wang, Y.; Chen, B. Knowledge-Based Spatial-Temporal Hierarchical MIMO Radar Waveform Design Method for Target Detection in Heterogeneous Clutter Zone. IEEE Trans. Signal Process. 2014, 63, 543–554. [Google Scholar] [CrossRef]
  17. Aubry, A.; Maio, A.D.; Huang, Y. MIMO Radar Beampattern Design via PSL/ISL Optimization. IEEE Trans. Signal Process. 2016, 64, 3955–3967. [Google Scholar] [CrossRef]
  18. Aubry, A.; Maio, A.D.; Piezzo, M.; Farina, A.; Wicks, M. Cognitive design of the receive filter and transmitted phase code in reverberating environment. IET Radar Sonar Navig. 2012, 6, 822–833. [Google Scholar] [CrossRef]
  19. Aubry, A.; Maio, A.D.; Farina, A.; Wicks, M. Knowledge-Aided (Potentially Cognitive) Transmit Signal and Receive Filter Design in Signal-Dependent Clutter. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 93–117. [Google Scholar] [CrossRef]
  20. Aubry, A.; Maio, A.D.; Naghsh, M.M. Optimizing Radar Waveform and Doppler Filter Bank via Generalized Fractional Programming. IEEE J. Sel. Top. Signal Process. 2015, 9, 1387–1399. [Google Scholar] [CrossRef]
  21. Piezzo, M.; Aubry, A.; Buzzi, S.; Maio, A.D.; Farina, A. Non-cooperative code design in radar networks: A game-theoretic approach. EURASIP J. Adv. Signal Process. 2013, 2013, 63. [Google Scholar] [CrossRef]
  22. Karbasi, S.M.; Aubry, A.; Carotenuto, V.; Naghsh, M.M.; Bastani, M.H. Knowledge-based design of space–time transmit code and receive filter for a multiple-input–multiple-output radar in signal-dependent interference. IET Radar Sonar Navig. 2015, 9, 1124–1135. [Google Scholar] [CrossRef]
  23. Haghnegahdar, M.; Imani, S.; Ghorashi, S.A.; Mehrshahi, E. SINR Enhancement in Colocated MIMO Radar Using Transmit Covariance Matrix Optimization. IEEE Signal Process. Lett. 2017, 24, 339–343. [Google Scholar] [CrossRef]
  24. Cui, G.; Yu, X.; Carotenuto, V.; Kong, L. Space-Time Transmit Code and Receive Filter Design for Colocated MIMO Radar. IEEE Trans. Signal Process. 2017, 65, 1116–1129. [Google Scholar] [CrossRef]
  25. Imani, S.; Nayebi, M.M.; Ghorashi, S.A. Colocated MIMO Radar SINR Maximization under ISL and PSL Constraints. IEEE Signal Process. Lett. 2018, 25, 422–426. [Google Scholar] [CrossRef]
  26. Bolhasani, M.; Mehrshahi, E.; Ghorashi, S.A. Waveform Covariance Matrix Design for Robust Signal-Dependent Interference Suppression in Colocated MIMO Radars. Signal Process. 2018, 152, 311–319. [Google Scholar] [CrossRef]
  27. Cumming, I.G.; Wong, F.H. Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation; Artech House: Norwood, MA, USA, 2005. [Google Scholar]
  28. Caves, R.G.; Harley, P.J.; Quegan, S. Clutter statistics along edge features in synthetic aperture radar imagery. J. Electromagn. Waves Appl. 1995, 9, 327–353. [Google Scholar]
  29. Grossi, E.; Lops, M.; Venturino, L. Track-before-detect for multiframe detection with censored observations. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 2032–2046. [Google Scholar] [CrossRef]
  30. Wang, B.; Wang, Y.; Hui, C. Spatial wavelet transform preprocessing for direction of arrival estimation. In Proceedings of the IEEE Antennas and Propagation Society International Symposium, San Antonio, TX, USA, 16–21 June 2002; pp. 672–675. [Google Scholar]
  31. Brenner, T.; Weiss, G.; Klein, M.; Kuschel, H. Signals and data fusion in a deployable multiband passive-active radar (DMPAR). In Proceedings of the IET International Conference on Radar Systems, Glasgow, UK, 22–25 October 2012; pp. 1–6. [Google Scholar]
  32. O’Hagan, D.W.; Ummenhofer, M.; Kuschel, H.; Heckenbach, J. A passive/active dual mode radar concept. In Proceedings of the 2013 14th International Radar Symposium, Dresden, Germany, 19–21 June 2013; Volume 1, pp. 136–142. [Google Scholar]
  33. Ward, J. Space-time adaptive processing for airborne radar. In Proceedings of the IEE Colloquium on Space-Time Adaptive Processing, London, UK, 6 April 1998; pp. 2/1–2/6. [Google Scholar]
  34. Klemm, R. Principles of Space-Time Adaptive Processing, 3rd ed.; The Institution of Engineering and Technology: London, UK, 2006. [Google Scholar]
  35. Chong, E.K.P.; Żak, S.H. An Introduction to Optimization, 4th ed.; John Wiley & Sons: New York, NY, USA, 2013. [Google Scholar]
  36. Kelly, E.J. An Adaptive Detection Algorithm. IEEE Trans. Aerosp. Electron. Syst. 1986, AES-22, 115–127. [Google Scholar] [CrossRef]
  37. Karbasi, S.M.; Aubry, A.; Maio, A.D.; Bastani, M.H. Robust Transmit Code and Receive Filter Design for Extended Targets in Clutter. IEEE Trans. Signal Process. 2015, 63, 1965–1976. [Google Scholar] [CrossRef]
  38. Aubry, A.; Maio, A.D.; Foglia, G.; Orlando, D. Diffuse Multipath Exploitation for Adaptive Radar Detection. IEEE Trans. Signal Process. 2015, 63, 1268–1281. [Google Scholar] [CrossRef]
  39. Cheng, X.; Aubry, A.; Ciuonzo, D.; Maio, A.D.; Wang, X. Robust Waveform and Filter Bank Design of Polarimetric Radar. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 370–384. [Google Scholar] [CrossRef]
Figure 1. Illustration of the multiple-input multiple-output (MIMO) radar with uniform linear arrays. (a) shows the observation geometry, and (b) describes the transmit and receive processes.
Figure 1. Illustration of the multiple-input multiple-output (MIMO) radar with uniform linear arrays. (a) shows the observation geometry, and (b) describes the transmit and receive processes.
Sensors 18 02722 g001aSensors 18 02722 g001b
Figure 2. Joint design of space-time transmit and receive weights.
Figure 2. Joint design of space-time transmit and receive weights.
Sensors 18 02722 g002
Figure 3. Variations of improvement factor (IF) with normalized frequencies. Results achieved by the proposed and conventional methods are shown in (a) and (b), respectively.
Figure 3. Variations of improvement factor (IF) with normalized frequencies. Results achieved by the proposed and conventional methods are shown in (a) and (b), respectively.
Sensors 18 02722 g003aSensors 18 02722 g003b
Figure 4. IF curves of the proposed and conventional methods. IF curves with respect to the normalized Doppler frequency, space transmit frequency, and space receive frequency are shown in (a), (b), and (c), respectively. In each figure, two other normalized frequencies are zero.
Figure 4. IF curves of the proposed and conventional methods. IF curves with respect to the normalized Doppler frequency, space transmit frequency, and space receive frequency are shown in (a), (b), and (c), respectively. In each figure, two other normalized frequencies are zero.
Sensors 18 02722 g004aSensors 18 02722 g004b
Figure 5. IF curves with respect to the normalized Doppler frequency.
Figure 5. IF curves with respect to the normalized Doppler frequency.
Sensors 18 02722 g005
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParameterValue
carrier frequency10 GHz
pulse repetition frequency12 KHz
inter-CPI interval0.125 s
platform speed90 m/s
altitude3 km
number of transmit elements4
number of receive elements4
transmit sensor spacing0.015 m
receive sensor spacing0.015 m
target speed56.61 m/s
target initial range10 km
target initial direction90°
signal/noise ratio−10 dB
clutter/noise ratio20 dB
number of temporal samples8

Share and Cite

MDPI and ACS Style

Yu, Z.; Wang, S.; Liu, W.; Li, C. Joint Design of Space-Time Transmit and Receive Weights for Colocated MIMO Radar. Sensors 2018, 18, 2722. https://0-doi-org.brum.beds.ac.uk/10.3390/s18082722

AMA Style

Yu Z, Wang S, Liu W, Li C. Joint Design of Space-Time Transmit and Receive Weights for Colocated MIMO Radar. Sensors. 2018; 18(8):2722. https://0-doi-org.brum.beds.ac.uk/10.3390/s18082722

Chicago/Turabian Style

Yu, Ze, Shusen Wang, Wei Liu, and Chunsheng Li. 2018. "Joint Design of Space-Time Transmit and Receive Weights for Colocated MIMO Radar" Sensors 18, no. 8: 2722. https://0-doi-org.brum.beds.ac.uk/10.3390/s18082722

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