Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Study of Earthquake Landslide Hazard by Defining Potential Landslide Thickness Using Excess Topography: A Case Study of the 2014 Ludian Earthquake Area, China
Next Article in Special Issue
2D DOA and Polarization Estimation Using Parallel Synthetic Coprime Array of Non-Collocated EMVSs
Previous Article in Journal
Cascaded-Filter-Based Reverberation Suppression Method of Short-Pulse Continuous Wave for Active Sonar
Previous Article in Special Issue
Analysis of Electromagnetic Wave and Multipath Suppression from Overhead Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Perturbation Transmit Beamformer Based Fast Constant Modulus MIMO Radar Waveform Design

1
College of Electronic Information Engineering, Inner Mongolia University, Hohhot 010010, China
2
National Key Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
3
The 601 Institute, The 6th Academy of China Aerospace Science and Industry CORP, Hohhot 010010, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(16), 2950; https://doi.org/10.3390/rs16162950
Submission received: 14 June 2024 / Revised: 28 July 2024 / Accepted: 6 August 2024 / Published: 12 August 2024

Abstract

:
In this paper, a fast method to generate a constant-modulus (CM) waveform for a multiple-input, multiple-output, (MIMO) radar is proposed. To simplify the optimization process, the design of the transmit waveform is decoupled from the design of transmit beamformers (TBs) and subpulses. To further improve the computational efficiency, the TBs’ optimization is conducted in parallel, and a linear programming model is proposed to match the desired beampattern. Additionally, we incorporate the perturbation vectors into the TBs’ optimization so that the TBs can be adjusted to satisfy the CM constraint. To quickly generate the CM subpulses with the desired range-compression (RC) performance, the classical linear frequency modulation (LFM) signal and non-LFM (NLFM) are adopted as subpulses. Meanwhile, to guarantee the RC performance of the final angular waveform, the selection of LFM signal parameters is analyzed to achieve a low cross-correlation between subpulses. Numerical simulations verify the transmit beampattern performance, RC performance, and computational efficiency of the proposed method.

Graphical Abstract

1. Introduction

In recent years, due to the increasing complexity of electromagnetic environments and the growing demand for detection, how to utilize the waveform-diversity capability of a radar’s transmitting end has attracted much attention [1]. Each transmitting antenna of a multiple-input, multiple-output, (MIMO) radar can transmit different waveforms. Therefore, it is a typical waveform-diversity radar. Compared to a traditional phased-array radar, a MIMO radar has advantages, such as the flexible design of the transmit beampattern, a high spatial resolution, and a strong multiple-target tracking ability [2,3,4,5,6,7,8,9].
A MIMO radar is usually divided into two types, based on the distance between radar antennas, namely a distributed MIMO radar and a colocated MIMO radar. In this paper, we focus on the study of a colocated MIMO radar, which can improve spatial resolution, flexibility in the transmit beampattern design, and multiple-target tracking capability. The transmit waveform is critical for realizing its advantages over a conventional phased-array radar [10]. For a colocated MIMO radar, the waveform design is usually divided into a quasi-orthogonal waveform design [11,12,13,14,15] and a partially correlated waveform design [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36]. The quasi-orthogonal waveform is defined as waveforms approaching orthogonality at all mutual delays so that the target echoes caused by different waveforms can be separated at the receiving end. Many works have been dedicated to lowering the auto-correlation and cross-correlation levels of transmit waveforms to obtain ideal quasi-orthogonal waveforms [11,12,13,14,15]. The MIMO radar-transmitting quasi-orthogonal waveforms can achieve the highest spatial resolution [37]. However, the cost of this high spatial resolution is an omnidirectionally transmit beampattern and a low transmission gain [35].
In order to meet the detection demand for different scenarios, such as multiple-target tracking [23,24,25,26], interference mitigation [27,28,29], and region-of-interest searches [20,30,36], a compromise between the transmit beampattern and transmission gain should be achieved. Therefore, some work has begun to focus on a partially correlated waveform design. The partially correlated waveform design can be classified into two sub-classes. One approach aims at obtaining the maximization output of the signal-to-interference-plus-noise ratio [17,18,19]. Another approach is to focus on the design of the transmit beampattern and optimize the temporal characteristic of the transmit waveform simultaneously [16,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,36]. The first sub-class not only utilizes the optimization resources of the transmitting end but also uses the optimization resources of the receiving end. However, in these works, prior knowledge of interference is assumed. Additionally, when considering the constant-modulus (CM) constraint that guarantees the radio frequency amplifier’s operating at the maximum efficiency, joint optimization is an intractable problem, and it cannot be solved within the acceptable computation time for a radar. For the second sub-class, some work has only involved the transmit beampattern design and ignored the temporal characteristics of the transmit waveform [28,31,36]. In [31,38], the classical two-step scheme was adopted, which involves optimizing the covariance matrix first and then synthesizing the CM waveform to approximate the optimal covariance matrix. It should be noted that covariance matrix optimization is a semi-definite programming (SDP) problem with a heavy computational burden [39]. The methods proposed in [36] achieve high computational efficiency, but the generated waveform does not meet the CM constraint. In [28], a parallel design strategy was adopted to improve computational efficiency. However, the sidelobe control was ignored.
It is well known that radar detection performance is closely related to the temporal characteristics of waveforms. Therefore, some work has investigated the joint design of waveforms’ spatial and temporal characteristics for a MIMO radar. Because of their good Doppler tolerance and ease of generation, in [20,30], linear frequency modulation (LFM) waveforms were used directly. However, to match a desired beampattern, both the methods in [20,30] involved nonlinear programming to optimize the parameter of LFM, which was solved using a sequential quadratic programming algorithm with high computational complexity [30]. In addition to adopting the LFM signal directly, the similarity constraint was also adopted to optimize the ambiguity function (AF) of the transmit waveform [21,32]. All the above work has focused on the temporal characteristics of each waveform out of each transmitting antenna. Due to the different signal-processing procedures, some work has been dedicated to designing the temporal characteristics of angular waveforms (AWs), which refer to the synthesized waveforms of the transmit waveform in different directions [22,26,33,34]. The work [22,26,34] applied advanced, non-convex optimization algorithms to the joint design of the transmit beampattern and AWs. However, under the CM constraint, the computational burden of these algorithms is heavy, and the convergence speed is slow. Time-diversity technology was utilized in the work [16] so that the CM transmit waveform could be generated directly. However, in [16], the flexibility of a MIMO radar transmit beampattern was not fully explored.
Motivated by the above issues, in this paper, we focus on the fast method to generate CM waveforms, which makes sense in rapidly changing detection environments. Additionally, the joint design of a transmit beampattern and AWs is considered under the proposed signal model. Specifically, the main contributions of this paper can be summarized as follows:
(1) 
Direct CM waveform synthesis based on proposed transmission strategy:
In the proposed transmission strategy, the transmit waveform can be obtained via a product of a transmit beamformer (TB) matrix and a subpulse vector. Therefore, under the proposed transmission strategy, once the TBs and subpulses that satisfy the CM constraints are obtained, the final CM transmit waveform is achieved.
(2) 
Fast and flexible transmit beampattern design based on perturbation TB:
In order to quickly realize the desired transmit beampattern, we formulated the beampattern design as a linear programming (LP) problem and a set of parallel TBs’ optimization. Both the above problems can be solved efficiently. Additionally, the perturbation vectors are introduced in the TBs’ design so that the TBs can be adjusted to meet CM constraints. Combining the LP and TB design, the final beampattern is flexible, and most scenarios are included, such as multiple beams, a wide beam, sidelobe control, and notches’ generation.
(3) 
Fast and low cross-correlation CM subpulses’ generation based on analytical expression:
In order to quickly generate the CM subpulses, we directly adopted the classical radar waveforms as the subpulses, such as the LFM signal and the non-LFM (NLFM) signal. Both of them have analytical expressions, and they are easy to generate. Moreover, both of them have an excellent AF feature. In the proposed signal model, we show that the temporal characteristic of AWs is consistent with that of subpulses. In order to distinguish between the different subpulses, the relation between cross-correlation and LFM parameters was analyzed. Then, a proper parameter setting was proposed to produce low-cross-correlation subpulses.
The rest of the paper is organized as follows. Section 2 introduces the signal model and clarifies the transmission strategy and signal processing scheme. In Section 3, the problem model is formulated, and the proposed method is shown in Section 4. The simulation results are presented in Section 5. Conclusions are drawn in Section 6.
Notation: Standard case letters stand for scalars; bold, uppercase letters denote matrices; and bold, lowercase letters denote vectors. The notations * , T , and  H denote the conjugate, the transpose, and the conjugate transpose, respectively. C L × L denotes the L × L complex space. ⊙ denotes the Hadamard product. The subscript 2 denotes the Euclidean vector norm. denotes the modulus operation. For a vector, the operation is imposed in an element-wise way. a n g l e ( ) denotes the phase of each element in a vector. d i a g x represents a diagonal matrix with its main diagonal filled with x , and j denotes an imaginary unit.

2. Signal Model

2.1. Transmission Strategy

A colocated MIMO radar with N transmitting antennas in a uniform linear array (ULA) was considered. Let s n t denote the waveform emitted via the n-th transmitting antenna, and the transmit waveform can be written as
S t = s 1 t , s 2 t , , s N t T
where 0 t T p , T p is the pulse duration.
The AW at direction θ can be expressed as
y t , θ = a T ( θ ) S ( t )
where
a θ = e j 2 π N 1 2 d sin ( θ ) / λ , e j 2 π N 1 2 1 d sin ( θ ) / λ , , e j 2 π N 1 2 d sin ( θ ) / λ T / N
denotes the N-dimensional transmission-steering vector. d is the inter-transmitter spacing, and λ is the signal wavelength.
The average power received at direction θ is expressed as
P θ = 1 T p 0 T p a T θ S t S H t a * θ d t
where P ( θ ) is the transmit beampattern.
To simplify the optimization complexity, the transmission strategy in [28] is introduced here. The transmit waveform is resolved into a product of a TB matrix, X C N × K , and a subpulse vector, f ( t ) C K × 1 :
S t = Xf t = x 1 , x 2 , , x K f 1 t f 2 t T 1 f K t k = 1 K 1 T k
where x k = x k 1 x k 2 x k N T is an N-dimensional TB, and the subpulse vector is composed of the time shift of subpulses f k ( t ) with support 0 , T k . Obviously, there is k = 1 K T k = T p .
In this way, the transmit waveform can be constructed as follows:
S t = x 1 f 1 t + k = 2 K x k f k t k ˜ = 1 k 1 T k ˜
Similarly, the transmit beampattern can be rewritten as
P θ = 1 T p k = 1 K P k θ 0 T k | f k t | 2 d t
where
P k θ = a T θ x k x k H a * θ
represents the beampattern of the k-th TB.

2.2. Signal Processing Scheme

As we know, the signal processing scheme of the receiving end determines which properties of the transmit waveform we need to optimize. Here, the signal processing scheme in [33] is introduced. In this scheme, the spatial receiving channel (SRC) consists of a receive beamformer and a range compressor (RC). The details of each SRC and the introduced signal processing scheme are given in Figure 1.
In Figure 1, M is the number of receiving antennas, and w θ = w 1 θ , w 2 θ , , w M θ T denotes the receiving beamformer for the direction θ . Review (2) and (5); under the proposed transmission strategy, the AW at direction θ can be rewritten as
y ( t , θ ) = a T θ x 1 f 1 t + k = 2 K x k f k t k ˜ = 1 k 1 T k ˜ = α 1 θ f 1 t + k = 2 K α k θ f k t k ˜ = 1 k 1 T k ˜
where α k θ = a T θ x k . Then, at the receiving array, we have
y ˜ ( t , θ ) = a ˜ ( θ ) y ( t , θ )
where
a ˜ θ = e j 2 π M 1 2 d sin ( θ ) / λ , e j 2 π M 1 2 1 d sin ( θ ) / λ , , e j 2 π M 1 2 d sin ( θ ) / λ T / M
denotes the M-dimensional receiving steering vector.
The output of the SRC for the direction θ 0 can be given by
y ˜ r b t , θ 0 = w H θ 0 y ˜ t , θ = α ˜ θ 0 , θ y t , θ
where α ˜ θ 0 , θ = w H θ 0 a ˜ ( θ ) . Let θ N denote the natural first null beamwidth after the received beamforming. It can be noticed that, when θ θ 0 θ N / 2 , θ 0 + θ N / 2 , we have α ˜ ( θ 0 , θ ) 1 . In the mainlobe region, the signals have similar temporal characteristics [40]. Therefore, after the received beamforming, only the signals from direction θ 0 are considered [26], and the output of the SRC for the direction θ 0 is approximated as
y ˜ r b t , θ 0 = α ˜ θ 0 , θ 0 y t , θ 0
It can be seen that after, the received beamforming, the signal’s temporal characteristic entirely depends on the AW at the beam-pointing direction. If  α k ¯ ( θ 0 ) / α k ( θ 0 ) 1 , k = 1 , , K , k k ¯ , the AW in (8) can be approximated as
y t , θ 0 = α k ¯ θ 0 f k ¯ t k ˜ = 1 k ¯ 1 T k ˜
In this situation, the subpulses f k ¯ t , which were carried via beamformer x k , decide the characteristics of AW at direction θ 0 , and the corresponding RC can adopt f k ¯ t .
If α k ¯ ( θ 0 ) / α k ^ ( θ 0 ) 1 , k ¯ k ^ , the AW in (8) can be approximated as
y t , θ 0 = α k ¯ θ 0 f k ¯ t k ˜ = 1 k ¯ 1 T k ˜ + α k ^ θ 0 f k ^ t k ˜ = 1 k ^ 1 T k ˜
In this situation, the cross-correlation between f k ¯ t  and f k ^ ( t ) should be designed so that the subpulses carried via different TBs can be separated. In this way, the RC can use f k ¯ t or f k ^ ( t ) , and the characteristics of AW in direction θ 0 are decided via   f k ¯ t   or   f k ^ ( t ) .

3. Problem Formulation

According to (6), it can be noticed that the spatial properties of a transmit waveform are decided according to the beampattern of TBs and the corresponding subpulses’ energy. Moreover, from (12) and (13), we can find that the temporal properties of a transmit waveform are determined according to subpulses. Through the above analysis, the transmit waveform design is divided into two parts: the transmit beampattern design and the subpulses’ correlation design.

3.1. The Transmit Beampattern Design

Let P d ( θ ) denote the desired transmit beampattern, and let  β k = 0 T k f k ( t ) 2 d t / T p denote the proportion of the power of the k-th subpulse. When considering the CM constraint, without losing generality, it is supposed that f k t = 1 , t 0 , T k , k = 1 , , K . Then, it is natural that k = 1 K β k = 1 , and the transmit beampattern (6) can be rewritten as
P θ = k = 1 K P k θ β k
In summary, the transmit beampattern design can be formulated as the following optimization problem:
min β k , P k ( θ ) l = 1 L u l P d θ l k = 1 K β k P k θ l s . t . β k 0 , k = 1 , . . . , K k = 1 K β k = 1 P k θ = a T θ x k x k H a * ( θ )
where L is the spatial grid number, u l is the weight for the l-th grid point, and P d θ is the desired transmit beampattern. The objective function of problem (15) denotes the absolute error between the optimized transmit beampattern and the desired transmit beampattern. The constraints in problem (15) can be concluded according to the above derivation about TB beampatterns. In problem (15), the CM constraint is ignored, and it will be considered in the following solution process. Nonetheless, the problem (15) is still a quartic programming problem, and it is difficult to solve efficiently. In the following subsection, a simplified method is proposed to solve the problem (15).

3.2. The Subpulses Correlation Design

There are two parts that are considered in the subpulses’ correlation design, which are auto-correlation and cross-correlation. Auto-correlation is related to RC performance, and cross-correlation determines the isolation degree between different subpulses. Therefore, subpulse correlation design can be specifically considered according to the radar’s working mode. When the radar is in tracking mode, the different TBs are designed for targets at different directions, and we have α k ¯ ( θ ) / α k ( θ ) 1 , k = 1 , , K , k k ¯ . In this situation, according to (11), only the auto-correlation of subpulses can be considered. The efficiency of waveform generation was considered, and the classical radar signals with analytical expressions, such as LFM and NLFM, could be adopted here.
When a radar is in search mode, the scanning beams are close to each other, and the correlation between TBs cannot be ignored. In this situation, in addition to auto-correlation, the cross-correlation properties between subpulses should be considered so that the subpulses carried via different TBs can be separated. Accordingly, the following problem is proposed to optimize the cross-correlation between subpulses:
min f k ¯ ( t ) , f k ^ ( t ) k ¯ = 1 , K k ^ = 1 , K k ¯ k ^ f k ¯ t f k ^ * t + τ d t s . t . f k ¯ t = 1 , t 0 , T k ¯ , k ¯ = 1 , K f k ^ t = 1 , t 0 , T k ^ , k ^ = 1 , K , k ¯ k ^
Some existing work has proposed methods to solve the above problem [11,12,13,14,15]. However, because the above problem is NP-hard, most of these methods cannot obtain a solution within the accepted computation time for a radar. For example, the proposed method in [15] involved singular value decomposition, and the computation complexity was around O I N 3 . In [14], though the majorization–minimization algorithm was used, the computation complexity was around O I N 2 . In the above, I represents the number of iterations. In the following section, we analyze the cross-correlation between different-parameter LFM signals, and we propose a method to obtain LFM signals with a low cross-correlation level and without complicated computations.

4. Problem Analysis and Proposed Method

In Section 3, the transmit waveform design is divided into the transmit beampattern design and the subpulses’ correlation design. In this section, problem (15) and problem (16) are analyzed successively, and the corresponding waveform design methods are proposed.

4.1. Analysis of Transmit Beampattern Design Problem (15)

Problem (15) is a quartic programming problem about the power proportion of subpulses and the TB beampattern. It is hard to get a satisfactory result directly with an accepted computation cost. However, it can be noticed that, when P k ( θ ) , k = 1 , , K are fixed, problem (15) can be simplified as an LP problem concerning β k , k = 1 , , K :
min β k , η l l = 1 L u l η l s . t . β k 0 , k = 1 , , K k = 1 K β k = 1 P d θ l k = 1 K β k P k θ l η l k = 1 K β k P k θ l P d θ l η l
where η l , l = 1 , , L are auxiliary variables. Problem (17) can be solved efficiently via some classical algorithms, e.g., the simplex method and dual simplex method. Based on the above analysis, an alternating optimization algorithm was employed for problem (15); in particular, we optimized β k , k = 1 , , K for fixed TB beampatterns first, and then the TBs were optimized for the fixed β k , k = 1 , , K . For clarity, the proposed strategy is summarized in Algorithm 1.
Algorithm 1 Proposed strategy for transmit beampattern design.
Input: Initial TBs x k ( 0 ) , k = 1 , , K , = 0
Output: The power proportion of subpulses β k ( ) , k = 1 , K and TBs x k ( ) , k = 1 , , K
1:
For fixed TBs x k ( ) , k = 1 , , K , solve problem (17) using dual simplex method, and obtain the optimized β k ( ) , k = 1 , K ;
2:
According to (14), calculate the transmit beampattern P ( ) ( θ ) , and obtain θ max = max θ l P ( ) θ l P d θ l , l = 1 , , L ;
3:
Let = + 1 ;
4:
Optimize x k ( ) , k = 1 , , K to satisfy P k ( ) θ max = P d θ max , k = 1 , , K ;
5:
Repeat steps 1–4 until some termination conditions are satisfied.

4.2. Perturbation TB-Based Transmit Beampattern Design

It can be seen that the remaining problem of Algorithm 1 is how to optimize x k , k = 1 , , K to satisfy P k θ max = P d θ max , k = 1 , , K . Additionally, another key issue is that the TBs should satisfy the following CM constraint:
x k = ξ , k = 1 , , K
where ξ denotes a constant vector, and all adjusted TB modulus values will be equal to it. Combining it with f k t = 1 , t 0 , T k , k = 1 , , K , we actually require the transmit waveform of each transmitting antenna to meet the CM constraint.
To satisfy Equation (18), the idea of perturbation TBs is proposed. In array signal processing, the array error cannot be ignored. Therefore, many works have proposed a robust beamforming algorithm under an array error case [41,42,43,44,45,46,47]. Almost all of the above works described the array error as a norm-bounded uncertainty set, and they proposed methods such as the worst-case optimization technique, convex optimization algorithm, and closed-form solution iteration to solve this problem. In this paper, we do not consider the array error, but we regard this error as an artificial perturbation to force the TB satisfying the CM constraint. Suppose that there is a set of perturbation vectors, λ k , k = 1 , , K and λ k 2 = ε x k 2 , k = 1 , , K , ε 1 . The key of the proposed idea is adding different perturbation vectors for each TB to satisfy the following equation:
x k + λ k = ξ , k = 1 , , K
Let V k r θ and V k u θ stand for the normalized transmit magnitude response and its upper boundary, respectively, and they can be expressed as
V k r θ = x k + λ k H a * θ / x k + λ k H a * θ k
V k u θ = x k H a * θ / x k H a * θ k + ε x k 2 / x k H a * θ k 1 ε x k 2 / x k H a * θ k
where θ k denotes the mainlobe direction of x k and the derivation of the upper boundary; see Appendix A.
Review Algorithm 1; after perturbation, when θ max is in the sidelobe region of P d θ , P k θ max = P d θ max , k = 1 , , K is equivalent to
V k r 2 θ max = P d θ max , k = 1 , , K
where it is assumed that P d θ and P k θ have been normalized. For the sidelobe region, (22) can be relaxed as
V k u 2 θ max = P d θ max , k = 1 , , K
It is assumed that the TBs are designed to generate a “fan” of beams jointly covering the target region [28], as shown in Figure 2. Based on this assumption, for the mainlobe region of the desired beampattern, the beampattern-matching problem is converted to maximize the mainlobe gain of each TB and allocate the power proportion for them to approach the desired mainlobe shape. This simplification cannot guarantee the optimal solution for problem (17), but it can decouple the optimization between TBs and improve the computational efficiency. After perturbation, the gain of the TBs can be written as
γ k = x k + λ k H a * θ k 2 x k + λ k 2 2 , k = 1 , , K
Using the triangle inequality, we have the lower boundary of the perturbation TB gain:
γ k 1 1 + ε 2 x k H a * θ k x k 2 ε 2 , k = 1 , , K
When θ max is in the mainlobe region of P d θ , we only focus on the transmission gain of each perturbation TB. According to inequality (25), the problem for the k-th TB can be relaxed as
max x k x k H a * ( θ k ) x k 2
In summary, for the k-th TB, the following problem is proposed:
max x k x k H a * θ k x k 2 s . t . V k u θ max = P d θ max
Problem (27) can be solved by using C 2 W O R D efficiently [47]. Compared with the traditional adaptive beamforming algorithms, which need covariance matrix inversion, the computational complexity of C 2 W O R D to obtain a closed-form solution is only O N 2 .
After all TBs are obtained, the CM constraint is considered. Suppose that x ˜ k , k = 1 , , K is the solution to (27), and we have
x ¯ c = k = 1 K x ˜ k K
Then, the final TBs are set as x k = x ¯ c a n g l e ( x ˜ k ) , k = 1 , , K . It can be proven that the modulus difference in solutions obtained from C 2 W O R D is small. The details are seen in Appendix B. Finally, the algorithm concerning perturbation TBs’ design is summarized in Algorithm 2.
After the above analysis, the proposed strategy in Algorithm 1 can be refined into the above Algorithm 3, which is called the perturbation TB-based transmit beampattern design (PTB-TBD) here.
Algorithm 2 Perturbation vector-based CM TBs’ design.
Input: Initial TBs x k ( ) , k = 1 , , K , the mainlobe directions of TBs θ k , k = 1 , , K , the maximum error direction θ max , the desired beampattern P d θ , and the perturbation parameter ε ;
Output: Optimized TBs x k ( + 1 ) , k = 1 , , K ;
1:
In parallel, solve problem (27), and obtain the solution, x ˜ k , k = 1 , , K ;
2:
Obtain x ¯ c via (28);
3:
Calculate x k ( + 1 ) = a n g l e x ˜ k x ¯ c ;
4:
Calculate the distance d k = | | x ¯ c | x ˜ k | | | 2 , k = 1 , , K ;
5:
If d k ε x k 2 , k = 1 , , K , output d k = x ¯ c x ˜ k 2 , k = 1 , , K ; otherwise, increase ε , and repeat above steps.
Algorithm 3 PTB-TBD.
Input: Initial TBs x k ( 0 ) , k = 1 , , K , = 0 , and  θ κ , k = 1 , , K ;
Output: The power proportion of subpulses β k ( ) , k = 1 , K and TBs x k ( ) , k = 1 , , K ;
1:
For fixed TBs x k ( ) , k = 1 , , K , solve problem (17) using the dual simplex method, and obtain the optimized β k ( ) , k = 1 , K ;
2:
According to (14), calculate the transmit beampattern P ( ) ( θ ) , and obtain θ max = max θ l P ( ) θ l P d θ l , θ l Ω s , where Ω s denotes the sidelobe region;
3:
Let = + 1 ;
4:
Optimize x k ( ) , k = 1 , , K using Algorithm 2;
5:
Repeat steps 1–4 until some termination conditions are satisfied.

4.3. Subpulses’ Cross-Correlation Design and CM Waveform Synthesis

As described in Section 3.2, the cross-correlation between subpulses should be optimized when the radar is in search mode. However, problem (16) is difficult to obtain a solution for within the accepted computation time. Given that the strong Doppler tolerance signal is more proper in search mode, two LFM signals with different parameters are analyzed in this section.
The LFM signals with different carrier frequencies are written as
s i t = exp j 2 π f i t + j π μ t 2 , t T k 2 , i = 1 , 2
where μ denotes the chirp rate. Without a loss of generality, it is assumed that the LFM signals have an equal pulse width.
When 0 τ T k , the cross-correlation of the above signals is
R 12 τ = s 1 t s 2 * t + τ d t = T k 2 T k 2 τ exp j 2 π f 1 t + j π μ t 2 exp j 2 π f 2 ( t τ ) j π μ ( t τ ) 2 d t = exp 3 j π f 2 τ 2 j π μ τ 2 j π f 1 τ × sin c ( f 1 f 2 + μ τ ) ( T k τ ) × ( T k τ )
When T k τ 0 , the cross-correlation of the above signals is
R 12 τ = s 1 t s 2 * t + τ d t = T k 2 τ T k 2 exp j 2 π f 1 t + j π μ t 2 exp j 2 π f 2 ( t τ ) j π μ ( t τ ) 2 d t = exp 3 j π f 2 τ 2 j π μ τ 2 j π f 1 τ × sin c ( f 1 f 2 + μ τ ) ( T k + τ ) × ( T k + τ )
In summary, the cross-correlation between the above LFM signals can be expressed as
R 12 ( τ ) = sin c f 1 f 2 + μ τ T k | τ | × T k | τ | , T k τ T k
From (32), it can be seen that R 12 ( τ ) decreases as f 1 f 2 increases. Based on this conclusion, we can generate low cross-correlation LFM signals by allocating different carrier frequencies. In this way, though the optimal solution of (16) is not obtained, the expected cross-correlation level can still be obtained by expanding the difference between carrier frequencies.
Finally, the overall fast CM waveform synthesis process is summarized in the flowchart Figure 3. After a review of the entire CM waveform synthesis process, it can be seen that the key to achieving fast performance through the proposed method lies in the following four points:
  • Matching the transmit beampattern design via LP solving;
  • Obtaining TBs via a closed-form solution without matrix inversion;
  • The parallel optimization of TBs and the CM operation based on perturbation vectors;
  • Generating low-cross-correlation CM subpulses based on analytical expressions.

5. Simulation Results

This section consists of three parts: transmit beampattern performance, RC performance, and a comparison with existing methods. In the transmit beampattern part, different desired transmit beampatterns are adopted, and the effectiveness of the proposed LP is shown. In the RC performance part, the auto- and cross-correlation of the designed subpulses are shown, and the final RC performance is presented. Finally, in the comparison part, we compare the waveform performance and computation time of the proposed method with those of other existing methods, and the numerical results verify the fast generation performance of the proposed method.
In all the simulations, a colocated MIMO ULA with N = 16 transmitting antennas, separated by a half-wavelength, was considered. For Algorithm 3, the termination condition was set such that the iteration number reached 500 or θ Ω s , P θ P d θ 0 , and the weights for the spatial grid point were set to 1. Unless otherwise specified, the simulation parameters remained unchanged. All the simulations were performed in the MATLAB 2017b/Windows 10 environment on a computer with a 2.1 GHz CPU and 32 GB of RAM.

5.1. Transmit Beampattern Design

In this subsection, we present the results of numerical simulations that were conducted to evaluate the proposed algorithm for the transmit beampattern design. The search mode was considered first, and the desired beampattern is expressed as
P d θ = 1 , θ 30 , 30 0.01 , θ 30 , 30
where a wide beam with a −20 dB sidelobe level was adopted. In order to avoid missing the effective TBs, the spacing between θ k was set to half of the 3-dB beamwidth. The specific value for each θ k is listed in Table 1. In Figure 4, the initial TBs’ beampattern and the desired beampattern are shown. It can be seen that the mainlobe region is covered with the initial transmit beams. Let the perturbation parameter ε = 0.03 , and let us conduct the proposed PTB-TBD. In Table 1, we represent the final power proportion of each subpulse. It can be found that some TBs were not used. Figure 5 shows the optimized TBs beampattern and the final transmit beampattern, respectively. It can be seen that the beams of TBs were concentrated in the mainlobe region, and the sidelobe level was lower than −20 dB. After power proportion allocation, we can see that the final transmit beampattern could approach the desired mainlobe shape, and the sidelobe satisfied the expected level. To verify the influence of the perturbation parameter ε , the final transmit beampatterns under different ε are shown in Figure 6. It can be seen that the impact of the perturbation parameter on the transmit beampattern performance was minimal. This result may imply that the optimal CM solution is included in the uncertain set ε 0.03 . Finally, Figure 7 depicts the amplitude of each TB obtained using the proposed method. It can be found that, for each transmitting antenna, the TBs satisfied the CM constraint.
Next, the multiple-target tracking scene was considered. The desired beampattern was obtained by adding three traditional single beams whose beam directions were 35 , 0 , and 35 , respectively, and the desired sidelobe level was set to −16 dB. Because the target directions are known in the tracking mode, the θ k , k = 1 , 2 , 3 were set to 35 , 0 and 35 , respectively. The perturbation parameter remained unchanged. Figure 8 shows the transmit beampattern performance. We can see that the performance was similar to that in Figure 5. In Figure 9, we show the influence of the perturbation parameter on the transmit beampattern with multiple beams. It can be seen that the result is similar to that in search mode. In Figure 10, the iterative process is shown. It can be found that the direction with the maximum sidelobe level was optimized in each iteration. Additionally, the sidelobe level gradually decreased as the iteration progressed.
Finally, the transmit beampattern design in the interference scenario was considered. It makes sense to form a transmit beampattern with notches at the interference directions [48]. Therefore, an interference region was added to the desired beampattern mentioned above, respectively. For the search mode, the interference region was 35 , 40 , and the notches’ level was −30 dB. For the tracking mode, the interference region was 15 , 20 , and the notches’ level was −26 dB. Figure 11 shows the transmit beampattern performance. It can be observed that the transmit beampattern produced using the proposed method could simultaneously form sidelobes and notches lower than the expected levels.

5.2. Subpulses’ Design and RC Performance

In this part, some simulations are presented that were conducted to verify the effectiveness of the proposed subpulse design method and the proposed signal processing scheme. Search mode was considered first, and the transmit beampattern result in Figure 5 was adopted here. From Table 1, we can see the useful TBs and the corresponding power proportion. It was assumed that the pulse width was 100 μ s, and the width of each subpulse could be obtained via the power proportion. For search mode, the LFM signal was adopted, and the signal expression for each subpulse was
s k ˜ t = exp j 2 π f k ˜ t + j π μ k ˜ t 2 , t T k ˜ 2 , k ˜ = 1 , , K ˜
where f k ˜ , μ k ˜ , and T k ˜ denote the start frequency, chirp rate, and subpulse width of the k ˜ th subpulse, respectively. K ˜ is the number of useful TBs whose power proportion is non-zero. From Table 1, we can see that K ˜ = 11 . To ensure the same range resolution for each subpulse, we set μ k ˜ = B / T k ˜ and B = 5 MHz . Finally, in order to obtain a low cross-correlation, according to (32), the start frequency was set to
f k ˜ + 1 f k ˜ = B , k ˜ = 1 , , K ˜ 1
Figure 12 shows the frequency spectrum of each subpulse. In Figure 12, subpulses 1 to 11 are carried via the useful TBs x k ˜ , k ˜ = 1 , , 11 , respectively. It can be observed that the interval between each start frequency is B, and the frequency spectrum is non-overlapping between subpulses. Figure 13 shows the correlation performance of subpulses. It can be seen that the cross-correlation level between subpulses is lower than −25 dB, which can satisfy the detection requirement of a radar.
In Figure 14, we verified the RC performance of the proposed signal processing scheme. The RC process was implemented using
RC τ = y ˜ r b t , θ k f k * t + τ d t
where θ k denotes the mainlobe direction of the k-th TB, and f k ( t ) is the subpulse carried via the k-th TB. According to Table 1, it can be understood that the corresponding RC for the AW at 26 . 3 is subpulse 1. Figure 14a shows the RC result of the AW at 26 . 3 and the auto-correlation of subpulse 1. It can be seen that the RC performance was almost the same as that of the expected LFM signal. It can be noticed that some TBs were not used, for example, the TB pointing to 30 . Given that the signals within the mainlobe region were highly correlated, subpulse 1 was also adopted as the RC for the AW at 30 . Figure 14b shows the RC result for the AW at 30 and the auto-correlation of subpulse 1. It can be found that results similar to those of 26 . 3 could be obtained.
In the tracking mode, only the auto-correlation of subpulses was considered, and the transmit beampattern result in Figure 8 was adopted. For the result in Figure 8, the power proportion for each TB was 0.3330, 0.3334, and 0.3336, respectively. It was assumed that the pulse width was 10 μ s , and the width of each subpulse could be obtained using the power proportion. For the tracking mode, the NLFM signal was adopted, and the phase expression of each subpulse was
φ k ˜ t = π ξ k ˜ t 2 T k ˜ π ζ k ˜ T k ˜ 1 4 t 2 / T k ˜ 2 2 , k ˜ = 1 , 2 , 3
where
ξ 1 = 30 / T 1 , ξ 2 = 30 / T 2 , ξ 3 = 60 / T 3 ζ 1 = 8 / T 1 , ζ 2 = 13 / T 2 , ζ 3 = 18 / T 3
and, for the TB pointing to 35 , 0 and 35 , the width of the subpulse was T 1 = 3.33 μ s, T 2 = 3.334 μ s, and T 3 = 3.336 μ s, respectively. Figure 15 shows the auto-correlation performance of subpulses. It can be seen that, because of the NLFM form, the subpulses had a sidelobe level lower than −20 dB. With the AW at 30 taken as an example, Figure 16 shows the RC result and the auto-correlation of subpulse 1. It can be seen that the RC performance was almost the same as that of the expected NLFM signal.

5.3. Comparison with Existing Methods

For comparison, the alternating-direction method of multipliers (ADMM) in [26] is considered here. For the transmit beampattern, the work [26] aimed at concentrating the transmission power in the mainlobe regions and maximizing the transmission gain, but it ignored the sidelobe level control. For temporal characteristics, the work [26] focused on the correlation of the beam directions’ AW, and the AW was adopted as the RC. In the proposed method, the subpulse was adopted as the RC, and the correlation of subpulses was considered. It is difficult to solve problem (15) directly. In order to verify the improvement in the solving efficiency using the proposed method, we converted problem (15) into the following SDP problem (37), and we solved it using CVX. According to the optimal solution to the SDP problem, the CM waveform was synthesized via the cyclic algorithm (CA) [38].
min R l = 1 L u l P d θ l a T θ l R a * ( θ l ) s . t . R 0
where R = k = 1 K β k x k x k H .
In the following simulations, the multiple-target tracking scene was considered. For the proposed method and SDP, the desired beampattern was obtained by adding two traditional single beams whose beam directions were 20 and 30 , respectively, and the desired sidelobe level was set to −16 dB. For the ADMM, the approaching points were calculated according to the description in [26]. For the temporal design, the proposed method adopted the NLFM signal as subpulses, and the parameter in (36) was changed as follows:
ξ 1 = 45 / T 1 , ξ 2 = 45 / T 2 ζ 1 = 13 / T 1 , ζ 2 = 13 / T 2
and the pulse width T p = T 1 + T 2 = 3.2 μ s. For the ADMM, the same NLFM signal was adopted as the desired AW. It should be noted that the pulse width in ADMM was 1.6 μ s in order to maintain the same range resolution as that in the proposed method. After temporal sampling, the pulse code lengths of all methods were 320.
In Figure 17, we set N = 10 and show the transmit beampattern result of the proposed method, SDP, CA, and ADMM. Additionally, we also show the result of the proposed method without considering the CM constraint. It can be seen that the proposed method and ADMM achieved similar mainlobe performance, but the mainlobe generated via SDP and CA was slightly wider. This is because SDP and CA use the criterion of minimizing the matching error, while the proposed method and ADMM maximize the mainlobe gain during waveform generation. It can be found that the proposed method achieved a lower sidelobe level than the other methods because only the proposed method focuses on controlling sidelobe levels. Comparing the proposed method with and without a CM constraint, we can find that the presence or absence of a CM constraint has almost no effect on the transmit beampattern performance of the proposed method. Figure 18 shows the RC performance of different methods. In Figure 18, the RC results of AW in beam directions are shown. It can be seen that the proposed method has a lower temporal sidelobe level than the other methods. This is because the proposed method directly uses NLFM signals, while ADMM optimizes phase-coded signals to approximate the desired NLFM signals, and CA does not consider AW properties.
In Table 2, we show the execution time of per iteration (ETPI) under the scenario in Figure 17 to compare the computational efficiency of all methods. In Table 2, except for the transmitting antenna number N, all the other settings are the same as those in Figure 17. The results of each method are the average of 100 trials with different random initial points. From Table 2, we can find that the ETPI of the proposed method is much less than that of ADMM and SDP, and the gap in the ETPI became larger as the transmitting antenna number increased. Moreover, the simulation experience revealed that the typical iteration number of the proposed method is about 10∼20, while the typical iteration number of ADMM is about 1000∼2000. Although CA is more efficient with a small number of elements, it is a subsequent step of SDP, and the combination of the two is necessary to obtain a CM waveform. In Table 3, the temporal sampling frequency was changed to evaluate the influence of the code length on the computational efficiency. In Table 3, except for the code length, all the other settings are the same as those in Figure 17. Because the subpulses are produced via analytical expressions, it can be seen that the ETPI of the proposed method remained almost unchanged with a code-length increase. Additionally, due to the fact that SDP optimization only involves the number of array elements, its ETPI is almost unaffected. However, the ETPI of ADMM increases significantly with the code length. According to the analysis in [26], we can understand that the computational complexity of ADMM is about O ( N L 2 ) , where L denotes the code length.

6. Conclusions

In this paper, we have proposed a fast method to produce a CM MIMO radar waveform. In order to improve the computational efficiency, a series of reasonable, simplified operations were proposed. The final simulations showed that the simplified operations can guarantee the beampattern and RC performance of the transmit waveform while improving the computational efficiency. In the proposed method, the key idea is the perturbation vector, which guarantees the CM property without a heavy computational burden. However, the boundary of perturbation parameter ε was not clarified, which may need further research. Moreover, for the transmit beampattern design, the improper desired beampattern may lead to no solution. Therefore, the transmit beampattern boundary of the proposed method should be investigated in the future.

Author Contributions

Conceptualization, H.Z. and J.Y.; methodology, H.Z. and H.W.; software, H.Z.; validation, H.Z. and H.W.; investigation, J.X. and Y.S.; writing—original draft preparation, H.Z.; writing—review and editing, Y.Z. and J.Y.; supervision, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Inner Mongolia Autonomous Region of China under Grants 2023QN06003 and 2021JQ07, the National Natural Science Foundation of China under Grants 62071345, 62361046, and 62371264, the Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region under Grant NJYT22109, the Training Plan for Young Innovative of Grassland Talents Project in Inner Mongolia Autonomous Region under Grant Q2022003, and the Innovation Capability Support Program of Shaanxi under Grant 2023KJXX-015.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

According to triangle inequality, we have
x k + λ k H a * θ = x k H a * θ + λ k H a * θ x k H a * θ + λ k H a * θ x k H a * θ + ε x k 2
and
x k + λ k H a * θ = x k H a * θ + λ k H a * θ x k H a * θ λ k H a * θ x k H a * θ ε x k 2
where we utilize the properties a ( θ ) 2 = 1 .
Then, according to the expression of V k r θ , we can obtain the upper boundary as
V k u θ = x k H a * θ + ε x k 2 x k H a * θ k ε x k 2 = x k H a * θ / x k H a * θ k + ε x k 2 / x k H a * θ k 1 ε x k 2 / x k H a * θ k

Appendix B

In the -th iteration, it is assumed that there are two TBs, w and w ˜ , produced via C 2 W O R D , and they can be expressed as
w = I P + β P w ( 1 )
w ˜ = I P + β ˜ P w ˜ ( 1 )
where I denotes the identity matrix, and P = a θ max a H θ max , w ( 1 ) and w ˜ ( 1 ) denote the TB result of the previous iteration. The modulus vector of w and w ˜ can be written as
w a b s = d i a g e j φ 1 w
w ˜ a b s = d i a g e j φ 2 w ˜
where φ 1 = a n g l e ( w ) and φ 2 = a n g l e ( w ˜ ) . Then, we have
w a b s w ˜ a b s 2 = d i a g e j φ 1 I P + β P w ( 1 ) d i a g e j φ 2 I P + β ˜ P w ˜ ( 1 ) 2 d i a g e j φ 1 w ( 1 ) d i a g e j φ 2 w ˜ ( 1 ) 2 + 1 β ˜ d i a g e j φ 2 P w ˜ ( 1 ) 1 β d i a g e j φ 1 P w ( 1 ) 2
Use the property
e j φ 1 = I P + β P w 1 I P + β P w 1
e j φ 2 = I P + β ˜ P w ˜ 1 I P + β ˜ P w ˜ 1
where the fraction operation is imposed in an element-wise way. Then, we have
d i a g e j φ 1 w ( 1 ) d i a g e j φ 2 w ˜ ( 1 ) 2 = I P * + β * P * w ( 1 ) * I P + β P w ( 1 ) w ( 1 ) I P * + β ˜ * P * w ˜ ( 1 ) * I P + β ˜ P w ˜ ( 1 ) w ˜ ( 1 ) 2 = w ( 1 ) 2 1 β * P * w ( 1 ) * w ( 1 ) w ( 1 ) 1 β P w ( 1 ) w ˜ ( 1 ) 2 1 β ˜ * P * w ˜ ( 1 ) * w ˜ ( 1 ) w ˜ ( 1 ) 1 β ˜ P w ˜ ( 1 ) 2
According to the procedure of Algorithm 3, θ max is located in the sidelobe region, and we have a H θ max w 1 1 and a H θ max w ˜ 1 1 . Additionally, according to the analysis in [49], we have 0 < β < 1 and 0 < β ˜ < 1 . Therefore, the influence of items ( 1 β ) Pw ( 1 ) and ( 1 β ˜ ) P w ˜ ( 1 ) is ignored, and Equation (A11) is rewritten as
d i a g e j φ 1 w ( 1 ) d i a g e j φ 2 w ˜ ( 1 ) 2 = 1 β ˜ * d i a g e j φ 2 ( 1 ) P * w ˜ ( 1 ) * 1 β * d i a g e j φ 1 ( 1 ) P * w ( 1 ) * 2
where the property w ( 1 ) = w ˜ ( 1 ) is used, φ 1 ( 1 ) = a n g l e w ( 1 ) , and φ 2 ( 1 ) = a n g l e w ˜ ( 1 ) .
Combining (A12) and inequality (A8), we have
w a b s w ˜ a b s 2 1 β ˜ * d i a g e j φ 2 ( 1 ) P * w ˜ ( 1 ) * 1 β * d i a g e j φ 1 ( 1 ) P * w ( 1 ) * 2 + 1 β ˜ d i a g e j φ 2 P w ˜ ( 1 ) 1 β d i a g e j φ 1 P w ( 1 ) 2
Let a H θ max w 1 = κ e j ω and a H θ max w ˜ 1 = κ ˜ e j ω ˜ . Then, inequality (A13) can be rewritten as
w a b s w ˜ a b s 2 κ ˜ 1 β ˜ * e j ω ˜ e j φ 2 ( 1 ) a * θ max κ 1 β * e j ω e j φ 1 ( 1 ) a * θ max 2 + κ ˜ 1 β ˜ e j ω ˜ e j φ 2 a θ max κ 1 β e j ω e j φ 1 a θ max 2 κ ˜ 1 β ˜ * + κ 1 β * + κ ˜ 1 β ˜ + κ 1 β
where the property a ( θ ) 2 = 1 is used. From the above analysis, we know that
κ ˜ 1 β ˜ * 1 κ 1 β * 1 κ ˜ 1 β ˜ 1 κ 1 β 1
Therefore, the difference w a b s w ˜ a b s 2 is small.

References

  1. Blunt, S.D.; Mokole, E.L. Overview of radar waveform diversity. IEEE Aerosp. Electron. Syst. Mag. 2016, 31, 2–42. [Google Scholar] [CrossRef]
  2. Fuhrmann, D.; San Antonio, G. Transmit beamforming for MIMO radar systems using partial signal correlation. In Proceedings of the Conference Record of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 7–10 November 2004; Volume 1, pp. 295–299. [Google Scholar] [CrossRef]
  3. Fishler, E.; Haimovich, A.; Blum, R.; Cimini, R.; Chizhik, D.; Valenzuela, R. Performance of MIMO radar systems: Advantages of angular diversity. In Proceedings of the Conference Record of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 7–10 November 2004; Volume 1, pp. 305–309. [Google Scholar] [CrossRef]
  4. Robey, F.; Coutts, S.; Weikle, D.; McHarg, J.; Cuomo, K. MIMO radar theory and experimental results. In Proceedings of the Conference Record of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 7–10 November 2004; Volume 1, pp. 300–304. [Google Scholar] [CrossRef]
  5. Yan, J.; Dai, J.; Pu, W.; Liu, H.; Greco, M. Target Capacity Based Resource Optimization for Multiple Target Tracking in Radar Network. IEEE Trans. Signal Process. 2021, 69, 2410–2421. [Google Scholar] [CrossRef]
  6. Yan, J.; Pu, W.; Zhou, S.; Liu, H.; Greco, M.S. Optimal Resource Allocation for Asynchronous Multiple Targets Tracking in Heterogeneous Radar Networks. IEEE Trans. Signal Process. 2020, 68, 4055–4068. [Google Scholar] [CrossRef]
  7. Yan, J.; Pu, W.; Zhou, S.; Liu, H.; Bao, Z. Collaborative detection and power allocation framework for target tracking in multiple radar system. Inf. Fusion 2020, 55, 173–183. [Google Scholar] [CrossRef]
  8. Zhang, R.; Cheng, L.; Wang, S.; Lou, Y.; Gao, Y.; Wu, W.; Ng, D.W.K. Integrated Sensing and Communication with Massive MIMO: A Unified Tensor Approach for Channel and Target Parameter Estimation. IEEE Trans. Wirel. Commun. 2024; early access. [Google Scholar] [CrossRef]
  9. Ma, Y.; Miao, C.; Long, W.; Zhang, R.; Chen, Q.; Zhang, J.; Wu, W. Time-Modulated Arrays in Scanning Mode Using Wideband Signals for Range-Doppler Estimation With Time-Frequency Filtering and Fusion. IEEE Trans. Aerosp. Electron. Syst. 2024, 60, 980–990. [Google Scholar] [CrossRef]
  10. Liu, H.; Zhou, S.; Zang, H.; Cao, Y. Two waveform design criteria for colocated MIMO radar. In Proceedings of the 2014 International Radar Conference, Lille, France, 13–17 October 2014; pp. 1–5. [Google Scholar] [CrossRef]
  11. Zhong, K.; Hu, J.; Pan, C. Constant Modulus MIMO Radar Waveform Design via Iterative Optimization Network Method. IEEE Trans. Instrum. Meas. 2023, 72, 8503911. [Google Scholar] [CrossRef]
  12. Qiu, Z.; Duan, K.; Wang, Y.; Liao, Z.; He, J. Designing Constant Modulus Approximate Binary Phase Waveforms for Multitarget Detection in MIMO Radar Using LSTM Networks. IEEE Trans. Geosci. Remote Sens. 2024, 62, 5103920. [Google Scholar] [CrossRef]
  13. Raei, E.; Alaee-Kerahroodi, M.; Babu, P.; Bhavani Shankar, M. Generalized waveform design for sidelobe reduction in MIMO radar systems. Signal Process. 2023, 206, 108914. [Google Scholar] [CrossRef]
  14. Song, J.; Babu, P.; Palomar, D.P. Sequence Set Design With Good Correlation Properties via Majorization-Minimization. IEEE Trans. Signal Process. 2016, 64, 2866–2879. [Google Scholar] [CrossRef]
  15. 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]
  16. Zou, X.; Jin, G.; He, F.; Zhang, Y. A New Waveform Design Method for Multi-Target Inverse Synthetic Aperture Radar Imaging Based on Orthogonal Frequency Division Multiplexing Chirp. Remote Sens. 2024, 16, 308. [Google Scholar] [CrossRef]
  17. Hu, J.; Zhang, W.; Zhu, H.; Zhong, K.; Xiong, W.; Wei, Z.; Li, Y. Constant modulus waveform design for MIMO radar via manifold optimization. Signal Process. 2022, 190, 108322. [Google Scholar] [CrossRef]
  18. Zhong, K.; Hu, J.; Zhao, Z.; Yu, X.; Cui, G.; Liao, B.; Hu, H. MIMO Radar Unimodular Waveform Design With Learned Complex Circle Manifold Network. IEEE Trans. Aerosp. Electron. Syst. 2024, 60, 1798–1807. [Google Scholar] [CrossRef]
  19. Chen, X.; Wu, H.; Cheng, Y.; Feng, W.; Guo, Y. Joint Design of Transmit Sequence and Receive Filter Based on Riemannian Manifold of Gaussian Mixture Distribution for MIMO Radar. IEEE Trans. Geosci. Remote Sens. 2023, 61, 5102213. [Google Scholar] [CrossRef]
  20. Li, H.; Zhao, Y.; Cheng, Z.; Feng, D. Correlated LFM Waveform Set Design for MIMO Radar Transmit Beampattern. IEEE Geosci. Remote Sens. Lett. 2017, 14, 329–333. [Google Scholar] [CrossRef]
  21. Yu, X.; Cui, G.; Zhang, T.; Kong, L. Constrained transmit beampattern design for colocated MIMO radar. Signal Process. 2018, 144, 145–154. [Google Scholar] [CrossRef]
  22. Wang, J.; Wang, Y. On the Design of Constant Modulus Probing Waveforms With Good Correlation Properties for MIMO Radar via Consensus-ADMM Approach. IEEE Trans. Signal Process. 2019, 67, 4317–4332. [Google Scholar] [CrossRef]
  23. Friedlander, B. On Transmit Beamforming for MIMO Radar. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 3376–3388. [Google Scholar] [CrossRef]
  24. Yan, J.; Liu, H.; Jiu, B.; Chen, B.; Liu, Z.; Bao, Z. Simultaneous Multibeam Resource Allocation Scheme for Multiple Target Tracking. IEEE Trans. Signal Process. 2015, 63, 3110–3122. [Google Scholar] [CrossRef]
  25. Yan, J.; Jiu, B.; Liu, H.; Chen, B.; Bao, Z. Prior Knowledge-Based Simultaneous Multibeam Power Allocation Algorithm for Cognitive Multiple Targets Tracking in Clutter. IEEE Trans. Signal Process. 2015, 63, 512–527. [Google Scholar] [CrossRef]
  26. Zheng, H.; Jiu, B.; Li, K.; Liu, H. Joint Design of the Transmit Beampattern and Angular Waveform for Colocated MIMO Radar under a Constant Modulus Constraint. Remote Sens. 2021, 13, 3392. [Google Scholar] [CrossRef]
  27. Zheng, H.; Jiu, B.; Liu, H.; Liu, Y.; Wang, X. Prior Knowledge Based Transmit Nulling Method for MIMO Radar. In Proceedings of the 2018 International Conference on Radar (RADAR), Brisbane, Australia, 27–31 August 2018; pp. 1–5. [Google Scholar] [CrossRef]
  28. Zheng, H.; Jiu, B.; Liu, H. Fast Constant Modulus MIMO Radar Waveform Design for Interference Mitigation. In Proceedings of the 2020 IEEE Radar Conference (RadarConf20), Florence, Italy, 21–25 September 2020; pp. 1–6. [Google Scholar] [CrossRef]
  29. Fan, W.; Liang, J.; Li, J. Constant Modulus MIMO Radar Waveform Design With Minimum Peak Sidelobe Transmit Beampattern. IEEE Trans. Signal Process. 2018, 66, 4207–4222. [Google Scholar] [CrossRef]
  30. Niu, B.; Zhao, Y.; Ge, B.; Zhang, T.; Zhang, M.; Tang, D. MIMO radar partially correlated waveform design based on chirp rate diversity. Signal Process. 2024, 220, 109450. [Google Scholar] [CrossRef]
  31. Stoica, P.; Li, J.; Xie, Y. On Probing Signal Design For MIMO Radar. IEEE Trans. Signal Process. 2007, 55, 4151–4161. [Google Scholar] [CrossRef]
  32. Cheng, Z.; He, Z.; Fang, M.; Li, J.; Xie, J. Spectrally Compatible Waveform Design for MIMO Radar Transmit Beampattern with Par and Similarity Constraints. In Proceedings of the 2018 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Calgary, AB, Canada, 15–20 April 2018; pp. 3286–3290. [Google Scholar] [CrossRef]
  33. Zhou, S.; Liu, H.; Hongtao Su, H.Z. Doppler sensitivity of MIMO radar waveforms. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 2091–2110. [Google Scholar] [CrossRef]
  34. Wang, Y.C.; Xu, W.; 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]
  35. Ahmed, S.; Thompson, J.S.; Petillot, Y.R.; Mulgrew, B. Finite Alphabet Constant-Envelope Waveform Design for MIMO Radar. IEEE Trans. Signal Process. 2011, 59, 5326–5337. [Google Scholar] [CrossRef]
  36. Bouchoucha, T.; Ahmed, S.; Al-Naffouri, T.; Alouini, M.S. DFT-Based Closed-Form Covariance Matrix and Direct Waveforms Design for MIMO Radar to Achieve Desired Beampatterns. IEEE Trans. Signal Process. 2017, 65, 2104–2113. [Google Scholar] [CrossRef]
  37. Li, J.; Stoica, P. MIMO Radar Diversity Means Superiority. In MIMO Radar Signal Processing; Wiley: Hoboken, NJ, USA, 2009; pp. 1–64. [Google Scholar] [CrossRef]
  38. 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]
  39. Luo, Z.Q.; Ma, W.K.; So, A.M.C.; Ye, Y.; Zhang, S. Semidefinite Relaxation of Quadratic Optimization Problems. IEEE Signal Process. Mag. 2010, 27, 20–34. [Google Scholar] [CrossRef]
  40. Bo, J.; Liu, H.; Xu, W.; Lei, Z.; Wang, Y.; Bo, C. 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]
  41. Tian, J.; Lei, S.; Lin, Z.; Gao, Y.; Hu, H.; Chen, B. Robust Pencil Beam Pattern Synthesis With Array Position Uncertainty. IEEE Antennas Wirel. Propag. Lett. 2021, 20, 1483–1487. [Google Scholar] [CrossRef]
  42. Dai, Y.; Sun, C.; Liu, X. On the 2D Beampattern Optimization of Sparse Group-Constrained Robust Capon Beamforming with Conformal Arrays. Remote Sens. 2024, 16, 421. [Google Scholar] [CrossRef]
  43. Shi, S.; Xu, Y.; Zhao, K.; Liu, Z. Polarised beampattern synthesis against array manifold mismatch. IET Microwaves Antennas Propag. 2021, 15, 910–924. [Google Scholar] [CrossRef]
  44. Shariati, N.; Zachariah, D.; Bengtsson, M. Minimum sidelobe beampattern design for MIMO radar systems: A robust approach. In Proceedings of the 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, 4–9 May 2014; pp. 5312–5316. [Google Scholar] [CrossRef]
  45. Cheng, Z.; He, Z.; Li, R.; Wang, Z. Robust transmit beampattern matching synthesis for MIMO radar. Electron. Lett. 2017, 53, 620–622. [Google Scholar] [CrossRef]
  46. Yu, X.; Cui, G.; Piezzo, M.; Iommelli, S.; Kong, L. Robust constrained waveform design for MIMO radar with uncertain steering vectors. EURASIP J. Adv. Signal Process. 2017, 2017, 2. [Google Scholar] [CrossRef]
  47. Zhang, X.; He, Z.; Zhang, X.; Xie, J. Robust Sidelobe Control via Complex-Coefficient Weight Vector Orthogonal Decomposition. IEEE Trans. Antennas Propag. 2019, 67, 5411–5425. [Google Scholar] [CrossRef]
  48. Neng-Jing, L.; Yi-Ting, Z. A survey of radar ECM and ECCM. IEEE Trans. Aerosp. Electron. Syst. 1995, 31, 1110–1120. [Google Scholar] [CrossRef]
  49. Zhang, X.; He, Z.; Liao, B.; Zhang, X.; Peng, W. Pattern Synthesis for Arbitrary Arrays via Weight Vector Orthogonal Decomposition. IEEE Trans. Signal Process. 2018, 66, 1286–1299. [Google Scholar] [CrossRef]
Figure 1. The introduced signal processing scheme.
Figure 1. The introduced signal processing scheme.
Remotesensing 16 02950 g001
Figure 2. A fan of beams jointly covering the target region.
Figure 2. A fan of beams jointly covering the target region.
Remotesensing 16 02950 g002
Figure 3. A flowchart of the proposed CM waveform synthesis method.
Figure 3. A flowchart of the proposed CM waveform synthesis method.
Remotesensing 16 02950 g003
Figure 4. The initial TBs beampattern (search mode).
Figure 4. The initial TBs beampattern (search mode).
Remotesensing 16 02950 g004
Figure 5. The transmit beampattern performance (search mode).
Figure 5. The transmit beampattern performance (search mode).
Remotesensing 16 02950 g005
Figure 6. The final transmit beampatterns under different ε values (search mode).
Figure 6. The final transmit beampatterns under different ε values (search mode).
Remotesensing 16 02950 g006
Figure 7. The amplitude of each TB.
Figure 7. The amplitude of each TB.
Remotesensing 16 02950 g007
Figure 8. The transmit beampattern performance (tracking mode).
Figure 8. The transmit beampattern performance (tracking mode).
Remotesensing 16 02950 g008
Figure 9. The final transmit beampatterns under different ε values (tracking mode).
Figure 9. The final transmit beampatterns under different ε values (tracking mode).
Remotesensing 16 02950 g009
Figure 10. The change in the transmit beampattern with the iterative process (tracking mode).
Figure 10. The change in the transmit beampattern with the iterative process (tracking mode).
Remotesensing 16 02950 g010
Figure 11. The transmit beampattern with notches.
Figure 11. The transmit beampattern with notches.
Remotesensing 16 02950 g011
Figure 12. The frequency spectrum of each subpulse.
Figure 12. The frequency spectrum of each subpulse.
Remotesensing 16 02950 g012
Figure 13. The auto- and cross-correlation of subpulses. (a) The auto-correlation of each subpulse. (b) The cross-correlation between subpulse 1 and other subpulses.
Figure 13. The auto- and cross-correlation of subpulses. (a) The auto-correlation of each subpulse. (b) The cross-correlation between subpulse 1 and other subpulses.
Remotesensing 16 02950 g013
Figure 14. The RC performance (search mode).
Figure 14. The RC performance (search mode).
Remotesensing 16 02950 g014
Figure 15. The auto-correlation of subpulses.
Figure 15. The auto-correlation of subpulses.
Remotesensing 16 02950 g015
Figure 16. The RC performance (tracking mode).
Figure 16. The RC performance (tracking mode).
Remotesensing 16 02950 g016
Figure 17. The transmit beampattern comparison.
Figure 17. The transmit beampattern comparison.
Remotesensing 16 02950 g017
Figure 18. The RC results’ comparison.
Figure 18. The RC results’ comparison.
Remotesensing 16 02950 g018
Table 1. The mainlobe directions of TBs and the power proportion of the corresponding subpulse.
Table 1. The mainlobe directions of TBs and the power proportion of the corresponding subpulse.
x 1 x 2 x 3 x 4 x 5 x 6 x 7 x 8 x 9
θ k 30 ° 26.3 ° 22.5 ° 18.8 ° 15 ° 11.3 ° 7.5 ° 3.8 °
β k 00.1488000.11940.03500.10380.04550.0981
x 10 x 11 x 12 x 13 x 14 x 15 x 16 x 17
θ k 3.8°7.5°11.3°15°18.8°22.5°26.3°30°
β k 0.04430.10380.03420.1206000.14640
Table 2. The ETPI (s) versus the transmitting antenna number.
Table 2. The ETPI (s) versus the transmitting antenna number.
N 1030507090
Proposed0.05330.05870.06040.06210.0743
ADMM1.31642.82887.345512.464718.6786
SDP0.12814.266210.561718.445533.3160
CA0.00120.01430.07170.49781.3980
Table 3. The ETPI (s) versus the code length.
Table 3. The ETPI (s) versus the code length.
Code Length 3264128256512
Proposed0.05160.05240.05270.05300.0539
ADMM0.03600.08650.21260.86684.4349
SDP0.12350.12470.12930.12870.1286
CA0.00020.00030.00060.00110.0015
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

Zheng, H.; Wu, H.; Zhang, Y.; Yan, J.; Xu, J.; Sun, Y. Perturbation Transmit Beamformer Based Fast Constant Modulus MIMO Radar Waveform Design. Remote Sens. 2024, 16, 2950. https://doi.org/10.3390/rs16162950

AMA Style

Zheng H, Wu H, Zhang Y, Yan J, Xu J, Sun Y. Perturbation Transmit Beamformer Based Fast Constant Modulus MIMO Radar Waveform Design. Remote Sensing. 2024; 16(16):2950. https://doi.org/10.3390/rs16162950

Chicago/Turabian Style

Zheng, Hao, Hao Wu, Yinghui Zhang, Junkun Yan, Jian Xu, and Yantao Sun. 2024. "Perturbation Transmit Beamformer Based Fast Constant Modulus MIMO Radar Waveform Design" Remote Sensing 16, no. 16: 2950. https://doi.org/10.3390/rs16162950

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