Transmit Beampattern Synthesis with Constant Beamwidth and Sidelobe Control for Wideband MIMO Radar

A beampattern synthesis approach is proposed to design the power spectral density matrix (PSDM), which is chosen to achieve a given transmit beampattern in wideband multiple-input multiple-output (MIMO) radar systems. The proposed approach focuses on transmit beampattern synthesis with constant beamwidth and sidelobe control. Moreover, the design problem is further converted to a convex optimization problem, which is solved efficiently via the modeling system CVX. In comparison to these recently developed wideband MIMO beampattern synthesis methods, the proposed approach maintains a constant beamwidth across the entire frequency band and provides a great improvement in sidelobe control. Numerical simulation results are obtained to validate the effectiveness of this approach.


Introduction
Transmit beampattern synthesis is a well-studied topic in standard phased-array radars [1][2][3].Recently, several transmit beampattern design approaches for narrowband MIMO radars have been reported [4][5][6][7][8][9].The main idea of these approaches is to design the proper probing signals according to the desired transmit beampattern.A relationship between the spatial beampattern and the covariance matrix of the transmit waveforms derived in [4,5] shows that it is possible to design a covariance matrix to synthesize the desired beampattern.This work is further extended in [6,7] to solve the beampattern design problem by a semidefinite quadratic programming (SQP) algorithm.
However, the discussion is limited to transmit beampattern synthesis for wideband MIMO radar.Antonio and Fuhrmann have designed the signal cross-spectral density matrixes (CSDM) in order to approximate the desired transmit beampattern in [10].He et al. [11] also have proposed a wideband beampattern formation via iterative techniques (WBFIT) to synthesize wideband MIMO waveforms, which satisfy some practical constraints such as constant-modulus or low peak-to-average power ratio (PAPR).In WBFIT, the beamwidth of the designed beampattern tends to shrink as the frequency increases for one narrow mainbeam.The method in [12] performs similarly as WBFIT, but with much lower computational costs.
Although these methods [10][11][12] have discussed the wideband beampattern synthesis problem, so little attention has been paid to the case of the constant beamwidth and the sidelobe attenuation.In this paper, we propose a beampattern synthesis approach to design power spectral density matrix (PSDM) in order to match a transmit energy distribution in both space and frequency for wideband MIMO radar.In our proposed approach, the mainlobe of the designed beampattern has an almost constant beamwidth for different frequencies and the sidelobe level of the designed beampattern is strictly below a desired sidelobe value.In addition, the proposed approach is converted to second-order cone programming (SOCP) problem, which is easily solved by the convex optimization toolbox CVX.Simulation results show that the proposed approach significantly outperforms the existing WBFIT approach in terms of both frequency invariant property and sidelobe attenuation.
The rest of this paper is organized as follows.The problem is formulated in Section 2. In Section 3, the proposed Notation 1.We denote vectors and matrices by lowercased and uppercased letters, respectively.(⋅)  denotes the vector or matrix conjugate transpose operation, (⋅)  denotes the transpose operation, (⋅) * denotes the conjugate operation, diag(⋅) denotes diagonals of a matrix, | ⋅ | denotes the absolute value operation, ‖ ⋅ ‖ 2 denotes the Euclidean norm, and rank(⋅) refers to the rank operation.

Problem Formulation
When a MIMO radar system employs  transmit antennas with the interelement spacing , the signal transmitted by the th antenna is given by where   denotes the carrier frequency of the transmitted signal and   () is the baseband signal, whose time support is [0, ].
Then ( 8) can be rewritten as In practical applications, the baseband signal and the spatial angle need to be discretized, so that ( 6) and ( 7) can be rewritten, respectively, as follows: where where  denotes the number of samples and is assumed to be even.So, the beampattern in (8) can be expressed on the discrete angle-frequency grid as where In this paper, when the transmit beampattern and a desired beampattern are defined as (13) and   , respectively, the goal is to choose the appropriate PSDM () ( ∈ [−/2, . . ., 0, . . ., /2 − 1]), under some constraints, so that the beampattern in (13) matches the desired one   .

The Proposed Transmit Beampattern Design
In this section, how to synthesize wideband transmit beampattern with constant beamwidth across the entire frequency band and sidelobe control is described.Here, the beampattern synthesis goal is to minimize the error between the designed beampattern and the desired one over the mainlobe region and to maintain the sidelobe region to be lower than a given threshold value.Mathematically, such beampattern synthesis problem can be formulated as follows: where    ∈ Θ    (  = 1, . . .,   ) and   ∈ Θ  ( = 1, . . ., ) that denote that a fine grid (uniform or nonuniform) of points covers mainlobe region and sidelobe region of interest, respectively.    and  denote the mainlobe region and the desired sidelobe value of a desired beampattern, respectively.And  > 0 is a variable to control the magnitude of the desired beampattern.
In (15), the PSDM cannot be chosen freely and should satisfy the following constraints: where / is the transmitted power from each antenna, and  is a total transmit power for  antennas.The first constraint  1 means that transmit power from all the antenna elements should be the same, and the second constraint  2 ensures that () is a positive semidefinite matrix.
Using the Parseval equality, the first constraint  1 on {  ()} imposes the following constraint on {  ()}: Therefore,  1 in (16a) can be rewritten as Combining (15) with the above PSDM constraints in (16a) and (16b), the beampattern design problem can be formulated as Next, it will be shown that the optimal design problem (19) can be converted to a second-order cone programming (SOCP) [15] and can be efficiently solved by the wellestablished interior point method [16], such as public domain software CVX toolbox [17].
Introducing  It can be observed that, by introducing new variables, the optimal design problem in (19) can be converted to the SOCP problem (22), which is a convex optimization problem [15].The second-order cone,   = {(, ) ∈   ×  | ‖‖ 2 ≤ }, can be solved using CVX in MATLAB with the code given in "{, } ⟨In⟩ complex lorentz (m)." We therefore emphasize more on convex problem formulation (22) and the computational complexity analysis for the solution is not provided.Now, if () is of rank one, an optimal solution to ( 14) can be obtained.However, if the corresponding rank is greater than one, we need to resort to randomization techniques to extract a feasible solution [18].
Hence the proposed approach works as follows.Firstly, the frequency band occupied by the wideband signal is divided into a certain number of narrowband bins.Secondly, PSDM of the signal vectors at each frequency bin is designed via proposed approach in (22).Then, a feasible solution to (14) from a solution () is extracted.

Numerical Results
This section focuses on demonstrating the performances of the proposed approach using an  = 10 MIMO radar system with uniform linear array (ULA), where the interelement spacing  is given by half-wavelength.The parameters for simulations are set as follows:   = 1 GHz,  = 200 MHz,  = 64,  = 180, and the total transmit power is  = 1.Although the design of multiple beams is not considered in these simulations, it can be achieved via the proposed approach.So, it is sufficient only to discuss the single beam design.Firstly, the desired beampattern which is synthesized using the proposed approach denotes the following: Here the mainlobe region is set to Θ    = [80 ∘ , 100 ∘ ]; that is, the beamwidth of the desired beampattern is 20 ∘ , and the sidelobe region is with 1 ∘ increment between adjacent grid points.Figure 1(a) shows the obtained beampatterns at all the 64 frequencies using the WBFIT approach of [11] and the proposed approach in (19).In this case, the desired sidelobe value is set to  = −35 dB.From the plots, it can been seen that the beampattern obtained by proposed approach shows a clean mainlobe at 90 ∘ across the entire frequency range, and the main beamwidth synthesized by proposed approach keeps fixed for all in-band frequency.Figure 1(b) compares the beampattern designed by the proposed approach with the WBFIT approach at center frequency.From Figure 1(b), it is noted that the designed beampattern using the proposed approach has a sidelobe level of −35 dB.
Secondly, the following desired beampattern with a wider mainlobe is considered: Similarly, Figure 2(a) shows the obtained beampatterns at all the 64 frequencies using the WBFIT approach and the proposed approach.It can be clearly seen that the designed beampattern, using the proposed approach, maintains a constant beamwidth across the entire frequency band from Figure 2(a), and more importantly, the designed sidelobe level is strictly below −40 dB from Figure 2(b).
From these simulation examples, compared with the existing WBFIT approach, the proposed approach has achieved significant improvement in both frequency invariant property and sidelobe attenuation.Therefore, the merit of the proposed approach is validated.

Conclusions
In this paper, a modified approach has been proposed to successfully synthesize the PSDM according to the desired transmit beampattern in wideband MIMO radar systems.With the proposed approach, the mainlobe of the designed beampattern has an almost constant beamwidth for different frequencies and the sidelobe level of the designed beampattern is strictly below a desired sidelobe value (e.g., −35 dB and −40 dB).Numerical results were provided to illustrate the effectiveness and validity of the proposed approach.

2 International
Journal of Antennas and Propagation approach is developed.Design examples and simulation results are shown in Section 4 and concluded remarks are given in Section 5.

Figure 1 :
Figure 1: Comparison of the synthesized beampattern.The desired beampattern is given in (23).(a) Plot at all the 64 frequencies using the WBFIT approach and the proposed approach.(b) Plot at center frequency.

Figure 2 :
Figure 2: Comparison of the synthesized beampattern.The desired beampattern is given in (24).(a) Plot using the WBFIT approach and the proposed approach.(b) Plot at center frequency.