Next Article in Journal
Recent Progress on Modeling Land Emission and Retrieving Soil Moisture on the Tibetan Plateau Based on L-Band Passive Microwave Remote Sensing
Next Article in Special Issue
Robust MIMO Waveform Design in the Presence of Unknown Mutipath Return
Previous Article in Journal
Asian Rice Calendar Dynamics Detected by Remote Sensing and Their Climate Drivers
Previous Article in Special Issue
A Robust Sparse Imaging Algorithm Using Joint MIMO Array Manifold and Array Channel Outliers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Radar Phase-Coded Waveform Design with Local Low Range Sidelobes Based on Particle Swarm-Assisted Projection Optimization

School of Information Science and Engineering, Harbin Institute of Technology, Weihai 264209, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(17), 4186; https://doi.org/10.3390/rs14174186
Submission received: 7 July 2022 / Revised: 12 August 2022 / Accepted: 23 August 2022 / Published: 25 August 2022
(This article belongs to the Special Issue Theory and Applications of MIMO Radar)

Abstract

:
In modern electronic warfare, cognitive radar with knowledge-aided waveforms would show significant flexibility in anti-interference. In this paper, a novel method, named particle swarm-assisted projection optimization (PSAP), is introduced to design phase-coded waveforms with multi-level low range sidelobes, which mainly considers the stability for randomized initialization under the unimodular constraint. Firstly, the mathematical problem corresponding to avoid the range sidelobe masking from multiple non-cooperative targets or interference is formulated by giving different threat levels. Then, based on the alternating direction decomposition idea, the original problem is divided into triple-variable ones where these non-linear approximations can be solved via alternating projections along with FFT. Furthermore, the PSAP method with swarm intelligence, learning factor, and particle-assisted projection could ensure the optimization convergence in a parallel way, which could relax the non-convex constraint and enhance the global exploiting performance. Finally, simulations for several typical scenarios and numerical results are all provided to assess the waveforms generated by PSAP and other prevalent ones.

Graphical Abstract

1. Introduction

Waveform diversity aided by high-performance radar hardware has received considerable attention and even made a great step forward to cognitive radar (CR) [1,2,3]. Most radars transmit a modulated waveform and then use some type of matched filter (MF) to enhance the signal-to-noise ratio (SNR) of the return echo. In mathematical sense, the MF output is usually the convolution between the received signal and the time-reversed replica of the transmitted signal, which is also regarded as the aperiodic auto-correlation [4,5,6]. Generally, to suppress range sidelobes which might obscure small targets of interest (especially, some dot targets or low RCS target) and further improve the anti-interference ability, the transmitting waveforms with desirable auto-correlation property are imperative by using some prior information [7,8,9,10]. Moreover, in engineering, the constant modulus (CM) of waveform (in most cases, i.e., unimodular property) could maximize the transmitter’s efficiency, but also makes the mathematical problem of generating waveform be more non-convex [11,12,13].
In the past 10 years, to achieve waveforms with ideal range sidelobes, minimizing integrated sidelobe level (ISL) and weighted integrated sidelobe level (WISL) have been developed as the classical metrics [14,15,16]. Therein, typical methods, such as cyclic algorithm new (CAN) [4], iterative spectral approximation algorithm (ISAA) [12], coordinate descent method (CD) [15], majorization minimization (MM) [16], weighted CAN (WeCAN) [17], alternating direction method of multipliers (ADMM) [18,19], etc., have also been presented and have received much attention. Note that, for mathematical problems under convex constraints, CANs could give some asymptotic convergence and make a big difference, while for the non-convex case CANs might stagnate into a suboptimum or local area [20,21]. Authors in [22] discussed the successive application of MM and phase gradient algorithm to synthesize the waveforms with low sidelobes. Additionally, authors adopted the limited memory Broyden–Fletcher–Goldfarb–Shanno algorithm (L-BFGS) to solve a fourth order polynomial formula and design unimodular sequences, but their ideas might incur invalidation for waveforms with large size [23]. Meanwhile, waveform optimization based on simulated annealing and stochastic neighborhood searching mechanism have also been presented, but these heuristic algorithms may be restricted to the modulus intricacy and large size of the waveforms [24,25].
Note that, when discussing the non-convex optimization with random initialization, these algorithms above would lead to a different terminus and only guarantee a local convergence [26]. Especially, to design the random phase-coded waveforms, how to tackle the random initialization for non-convex case has always been the key issue in engineering. Parallel optimization based on swarm intelligence maybe a good choice to improve the robustness. Unlike the optimization methods regarding each phase-coded unit, the waveform sample-based projection optimization using FFT will make a difference. In this paper, we use the idea of swarm particle intelligence [27,28,29] and combine alternating projection and particle swarm intelligence together to improve the global exploiting. To this end, our particle swarm-assisted projection optimization (PSAP) method is introduced to design the waveforms with desirable range sidelobes. Firstly, the mathematical problem is formulated to tackle the multiple non-cooperative targets or interference. Furthermore, using the alternating direction idea, the original problem is divided into some triple-variable ones considering different constraints. Then, the spectrum approximation in the sense of F-norm can be transformed into multi-variable alternating optimization cases. Finally, with the help of particle swarm intelligence, phase retrieval, learning factor and accelerated projections, PSAP method and its accelerated version have been formulated.
The remainder of paper is organized as follows. In Section 2, the system model is shown, and the formulated optimization problem for minimizing sidelobes is derived. In Section 3, PSAP as a novel alternative optimization mechanism based on swarm intelligence and FFT is presented. In Section 4, the performance of the proposed algorithms is evaluated and a series of numerical examples are also provided. Finally, in Section 5, the concluding remarks are provided.

2. The Signal Model and Problem Description

In this section, we discuss the range sidelobes masking effect from some strong RCS scatters which might obscure small targets of interest (such as dot target or low RCS target), and mainly focus on the waveform design for the static target detection in masking scenario and ignore the Doppler effect. For general description, the pre-modulated transmitting waveform in time domain has the discrete form of base band sequence with N code elements, i.e., s = [ s 1 s 2 s N ] T . The received signal is down converted to base band and undergoes the MF at the receiver [30]. The vector format of autocorrelation function which can be regarded as the MF output at the zero Doppler shift, has been listed as follows
α ( s ) = [ α N + 1 ( s ) α 1 ( s ) α 0 ( s ) α 1 ( s ) α N 1 ( s ) ] T
α n ( s ) = n ˜ = 1 N n s n ˜ s n ˜ + n * = α n * ( s )
As shown in [12], suppose that a strong point scatterer (or the interference) with echo power σ t 2 ( q ) exists in the q-th range cell, and a weak target of interest with echo power σ t 2 ( r ) exists in the r-th range cell, the noise plus range sidelobe interference for the r-th range cell can be represented as
σ I 2 ( r ) = Ω ( N 1 | q r | ) | σ t ( q ) α q r ( s ) / N | 2 + σ n 2 ,
where Ω ( n ) = { 1 n 0 0 n < 0 , σ n 2 is the power of the thermal noise, and α q r ( s ) denotes the (q-r)-th lag of the aperiodic autocorrelation of (1). Then the target Signal power to Interference plus Noise Ratio (i.e., SINR) can be expressed as
SINR = σ t 2 ( r ) σ I 2 ( r ) = = σ t 2 ( r ) Ω ( N 1 | q r | ) | σ t ( q ) α q r ( s ) / N | 2 + σ n 2
here, of all the variables in (4), α q r ( s ) is the only one under the control of radar transmitter.
Despite the phase-coded waveform s with arbitrary amplitude, the CM waveform could maximize the transmitter’s efficiency [13,16,31]. Let x = [ x 1 x 2 x N ] T N × 1 denote CM phase-coded one with N discrete phase-coded units, i.e., x n = e j ψ n where ψ n denotes the n-th phase-unit extracted from [ 0 , 2 π ] . ( ) T , ( ) * means the operation of transpose, complex conjugate, respectively. Similarly, we use α n ( x ) = n ˜ = 1 N n x n ˜ x n ˜ + n * = α n * ( x ) to denote the range sidelobes of MF output in lieu of α n ( s ) in (1) and (2). In modern electronic countermeasures scenarios, range sidelobes units occupied by some powerful interference or extended-scatters which mask the interesting target need to be suppressed [8,10,17]. Namely, local low range sidelobes is more convenient for weak targets detection and anti-masking effect as discussed in (3) and (4). Suppose that a powerful dot-scatter exists in the q-th range unit and inevitably impacts the target detection of the r-th one. With the help of some prior information R r , i.e., R r = { ± | n q | : n Z ^ r } \ { N 1 , N + 1 } where Z ^ r denotes the area with a foreseeable target, we could further describe the range interval R r to be suppressed. We use an indicating vector 𝔃 = [ z 1 , z 2 z 2 N ] T to formulate the area-mapping of these pre-suppressed range sidelobes interval, i.e., z n = 1 when n R r , and z n = 0 when n R r . The classical weighted ISL given some prior information has been discussed in [4,7,17], i.e.,
WISL ( x ) = min x 2 n = 1 N 1 w n | α n ( x ) | 2 , w n 0
Borrowed the idea of (5), we denote { ϖ n } n = 1 2 N as the weight corresponding to each unit in 𝔃 = [ z 1 , z 2 z 2 N ] T , i.e., ϖ n = δ 1 1 when z n = 1 , and ϖ n = 1 when z n = 0 . We further define x ˜ as the desirable waveform with ideal local low sidelobes, i.e., lim x x ˜ α ( x ˜ ) = ϖ α ( x ) , denotes the Hadamard element-wise product. Let v denote the approximation vector with v ϖ α ( x ) . To design waveform with desirable property, namely, we should make α ( x ) and α ( x ˜ ) be more approximate. Here we use the norm-metric α ( x ) v 2 to denote the approximation level of them. Finally, the objective function can be formulated as
min x α ( x ) α ( x ˜ ) 2 = min x α ( x ) v 2
According to the “Parseval-type equality” in [17], i.e., FC T x 2 = x 2 , the objective function of (6) is equivalent to
min x FC T α ( x ) FC T v 2
FC T v = f ˜ f ˜ * 2 N × 1 ,   FC T α ( x ) = f f * 2 N × 1
where C denotes the extend or cutoff matrix with C = [ I N × N 0 N × N ] , f ˜ and f = FC T x refers to the frequency spectrum of x ˜ and the designing one, respectively. The DFT matrix F 2 N × 2 N is constituted by unity exponential factor, i.e., f ^ k H = [ e j w k e j 2 N w k ] with w k = 2 π k / 2 N . Next, the optimization problem in (6) can be transformed into the spectrum approximation as following
min x ( FC T x ) ( FC T x ) * f ˜ f ˜ * 2 min x k = 1 2 N | ( | ( f ^ k H C T x ) | 2 ( ( f ^ k H C T v ) ) 2 ) | 2 s . t . | x n | = 1 ,   n = 1 , 2 , , N
Note that the problem of (9) is a quartic function of { x n } n = 1 N , using the auxiliary phase vector ϕ = [ ϕ 1 ϕ 2 ϕ 2 N ] T , (9) can be justified ‘almost equivalent’ to the following quadratic function of { x n } n = 1 N [17], i.e.,
min x , ϕ k = 1 2 N | ( f ^ k H C T v ) e j ϕ k ( f ^ k H C T x ) | 2
and the ideal frequency spectrum vector f ˜ can be further expressed as
f ˜ = ( FC T v ) 1 / 2 exp ( j ϕ )
As v 2 N × 1 ϖ α ( x ) , (10) implies that once given x , the ideal v satisfies
min v v ϖ α ( x ) 2 s . t . | x n | = 1 , n = 1 , 2 , , N
namely, v ϖ α ( x ) , so that ( f ^ k H C T v ) will be constant once given x, then (10) in the vector format has
min ϕ ( FC T v ) 1 / 2 exp ( j ϕ ) ( FC T x ) 2 s . t . { v ϖ α ( x ) ϕ n [ 0 , 2 π ] , n = 1 , 2 , , 2 N
For brevity, defining a novel operator diag ( ) which rearranges the column vector to be a square diagonal matrix, i.e., ( FC T v ) 1 / 2 exp ( j ϕ ) = diag ( exp ( j ϕ ) ) ( FC T v ) 1 / 2 , and the objective function (13) has
min ϕ diag ( j ϕ ) ( FC T v ) 1 / 2 FC T x 2 s . t . { v ϖ α ( x ) ϕ n [ 0 , π ] , n = 1 , 2 , , N
Similarly, given ϕ and v , the quadratic optimization problem of (10) has
min x ( ( FC T v ) 1 / 2 exp ( j ϕ ) ) FC T x 2 s . t . | x n | = 1 , n = 1 , 2 , , N
Let J = ( FC T v ) 1 / 2 2 N × 1 , then diag ( exp ( j ϕ ) ) J = diag ( J ) exp ( j ϕ ) , (14) can be rewritten as
min ϕ diag ( J ) exp ( j ϕ ) FC T x 2
Meanwhile, diag ( J ) is also an invertible matrix, the estimated ϕ in (16) could be given by
ϕ = a n g ( ( ( diag ( J ) ) H diag ( J ) ) 1 ( diag ( J ) ) H FC T x )
where a n g ( ) = i m ( ln ( diag ( ) ) ) represents the phase extractor from its vector argument, and im ( ) denotes the imaginary component extraction operator. Similarly, recall “Parseval-type equality” again, the objective function in (15) could be expressed as:
min x CF H ( ( FC T v ) 1 / 2 exp ( j ϕ ) ) x 2 = CF H ( diag ( exp ( j ϕ ) ) ( FC T v ) 1 / 2 ) x 2 s . t . | x n | = 1 , n = 1 , 2 , , N
On account of CM property, we only retain the phase section of the estimated x of (18). Then the designed waveform by the phase-retrieval operation of [20] has
x = exp ( j a n g ( CF H ( diag ( exp ( j ϕ ) ) ( FC T v ) 1 / 2 ) ) )

3. The Proposed Particles Swarm Assisted Projection Framework

In recent years, CAN [4], RISAAP [5], PONLP [12], CD [15], WeCAN [17], and ADMM [18] are all presented to deal with the Non-deterministic Polynomial-time hard (NP-hard) problem. They have the common trait, i.e., iterative optimization mechanism. Therein, each iteration of them requires to solve a non-convex problem under unimodular constraint, no matter by virtue of handling the bisection gradient optimization or FFT-based one. However, the initialized random phases where each phase unit is distributed in [ 0 , 2 π ] , would incur a different terminus when using different random initialization. Namely, each Monte-Carlo trial could obtain different solution and remain some local convergence. To tackle these, we borrow the idea of parallel optimization to combine the particles swarm intelligence and projection optimization together, where the novel particles projection mechanism will avoid the local area in the statistical sense. Next, we present the PSAP framework in lieu of the traditional evolution mechanism of PSO or DE [28,29], and thus could enhance the global exploiting for non-convex phase-coded problem. The detailed descriptions of PSAP have been listed as follows
Step 1. Formulate the waveform set rather than one single sequence, i.e., X = [ x 1 x m x M ] M × N with x m = [ x m ( 1 ) x m ( n ) x m ( N ) ] T N × 1 . Note that each sequence has x m ( n ) = e j ψ n , where ψ n denotes the independent phase-coded variables extracted from [ 0 , 2 π ] . Additionally, the waveform set could also be initialized by Frank or Barker sequence.
Step 2. Define a novel metric to assess the sidelobe performance of waveform in the range interval R l as
f i t n e s s ( x m ) = 1 n u m _ R r k R r α k ( x m )
where n u m _ R r denotes the number of pre-suppressed sidelobes units.
Step 3. Using (20) as the fitness function to evaluate each waveform x m of set X , and select the best fitness function and its corresponding waveform p t .
Step 4. For the t-th iteration, select M ^ < M waveforms from waveform set X to formulate the novel subset. Here we should update the subset by some rules that if the best fitness function and its corresponding waveform have not been incorporated, then use it to replace the worst one in current subset.
Step 5. For each waveform x m at the t-th iteration, utilize the oversampling FFT to get x m f m , and then formulate the relaxing factor δ ( x m ) , relational factor d m , and also projection vector v m by using (22)~(24) respectively, as follows
f m = ( F ϖ α ( x m ) ) ( 1 / 2 ) ( exp ( j a n g ( FC T x m ) ) )
δ ( x m ) = ( F ϖ α ( x m ) ) ( 1 / 2 ) ( exp ( j a n g ( FC T x m ) ) ) x m ( exp ( j a n g ( C T F f m ) ) ) x m
d m = δ ( x m ) ( exp ( j a n g ( C T F f m ) ) x m )
v m = w d m + c 1 rand ( y m x m ) + c 2 rand ( p x m )
where d m denotes the relational factor of the m-th waveform at the t-th iteration to the t + 1-th iteration, rand denotes one random value of [0, 1], w [ 1 , 1.5 ] represent the inertia factor, c 1 [ 0.5 , 1 ] is the learning factor related to the best waveform y m which is selected from initial iteration to the current iteration { x m 1 x m t } . c 2 [ 0 , 0.5 ] represents the learning factor related to the best waveform p t at the current t-th iteration.
Step 6. Use the multi-particles projection of (25)~(26) to achieve unimodular waveforms at the t + 1-th iteration, as follows
x m = x m + v m
x ^ m = ( exp ( j a n g ( x m ) ) )
where x ^ m with CM property has been obtained by the phase retrieval operation in (26).
Step 7. Consider the remaining subset, update them by
ϕ m = a n g ( diag ( ( ( FC T ( ϖ α ( x m ) ) ) 1 / 2 ) 1 ) FC T x m )
x ^ m = exp ( j a n g ( CF H ( diag ( exp ( j ϕ m ) ) ( FC T ( ϖ α ( x m ) ) ) 1 / 2 ) ) )
Step 8. Merge the subset and remaining part, use the following selection rules to get y m , i.e., y m = x ^ m when f i t n e s s ( x ^ m ) < f i t n e s s ( x m ) , or not y m = x m , and select the best fitness function and its corresponding waveform p t .
Step 9. Repeat step 3~8 until | f i t n e s s ( p t ) f i t n e s s ( p t 1 ) | < ε or num > iter_num, then output p t . Otherwise, update t = t + 1 and continue iterating.
Furthermore, we could incorporate the gradient steepest idea into PSAP. The objective function of (15) can be expressed as:
min x J ( x ) = FC T x | FC T diag ( 𝔃 ) α ( x ) | a n g ( f ˜ ) 2
Let ψ = a n g ( f ˜ ) , then the two-variable optimization problem has
min ψ , x J ( x , ψ ) = FC T x | FC T diag ( 𝔃 ) α ( x ) | ψ 2
Given the latest iterative x , ψ can be achieved by the follow gradient optimization, and its gradient matrix has
ψ J ( x ) = abs ( F [ x 0 N × 1 ] ) exp ( j a n g ( FC T x ) ) | FC T diag ( 𝔃 ) α ( x ) | ψ 2 ψ = j diag ( abs ( FC T x ) ) diag ( | FC T diag ( 𝔃 ) α ( x ) | ) ψ ) ( exp ( j ( a n g ( FC T x ) ψ ) ) exp ( j ( ψ a n g ( FC T x ) ) ) )
2 abs ( FC T x ) exp ( j a n g ( FC T x ) ) | FC T diag ( 𝔃 ) α ( x ) | ψ 2 ψ ψ T = diag ( aba ( FC T x ) ) diag ( v ) diag ( exp ( j ( a n g ( FC T x ) ψ ) ) + exp ( j ( ψ a n g ( FC T x ) ) ) )
Let (32) equal to 0, then ψ has:
ψ = a n g ( FC T x ) + π ξ
where ξ , as known, only if ξ is the even value, the hessian matrix of (33) can be positive, and (30) will get the minimum. Moreover, given ψ , the phase vector of x has:
a n g ( x ) = a n g ( CF H ( | FC T diag ( 𝔃 ) α ( x ) | ψ ) ) + π ξ
Finally, the detailed descriptions of particle swarm-assisted projection with optimizing mechanism (named as PSAPOM) have been listed as follow:
Step 1~Step 3 are similar to PSAP.
Step 4. Define x ¯ ( 1 ) = p t , then calculate the gradient direction d 0 = g 0 = ψ J ( x ¯ ( 1 ) ) , then search the best length ϑ which satisfies
J ( diag ( exp ( j ϑ d k ) ) x ¯ ( k ) ) J ( diag ( exp ( j η d k ) ) x ¯ ( k ) ) , η 0
define x ¯ ( k + 1 ) = diag ( exp ( j ϑ d k ) ) x ¯ ( k ) , then
g k + 1 = ψ J ( x ¯ ( k + 1 ) )
d k + 1 = ( g k + 1 + ( g k + 1 g k ) T g k + 1 d k / g k 2 )
if | f i t n e s s ( x ¯ ( k + 1 ) ) f i t n e s s ( x ¯ ( k ) ) | < 10 5 or k > 100 , output the initiation p t = x ¯ ( k + 1 ) ;
Step 5. For the current iteration, select M ^ < M waveforms from set X to formulate the novel subset. Here we should update the subset by some rules that if the best fitness function and its corresponding waveform have not been incorporated, then use it to replace the worst one in current subset.
Step 6. For each waveform x m of subset, utilize the oversampling FFT to get x m f m , and then formulate the relaxing factor δ ( x m ) , relational factor d m , and also projection vector v m by using (39)~(41) respectively, as follows
f m = ( F ϖ α ( x m ) ) ( 1 / 2 ) ( exp ( j a n g ( FC T x m ) ) )
δ ( x m ) = ( F ϖ α ( x m ) ) ( 1 / 2 ) ( exp ( j a n g ( F C T x m ) ) ) x m ( exp ( j a n g ( C F H f m ) ) ) x m
d m = δ ( x m ) ( ( exp ( j a n g ( C F H f m ) ) ) x m )
v m = w d m + c 1 rand ( y m x m ) + c 2 rand ( p t x m )
where d m denotes the relational factor of the m-th waveform at the t-th iteration to the t + 1-th iteration, rand denotes one random value in [0, 1], w [ 1 , 1.5 ] represent the inertia factor, c 1 [ 0.5 , 1 ] is the learning factor related to the best waveform y m which is chosen from initial iteration to the current iteration { x m 1 x m t } . c 2 [ 0 , 0.5 ] represents the learning factor related to the best waveform p t at the current t-th iteration.
Step 7. Use the multi-particles projection to achieve unimodular waveforms at t + 1-th iteration,
x m = x m + v m
x ^ m = ( exp ( j a n g ( x m ) ) )
where x ^ m with CM property has been obtained by the phase retrieval operation of (43).
Step 8. For the remaining subset, update them by
ϕ m = a n g ( diag ( ( ( F C T ( ϖ α ( x m ) ) ) 1 / 2 ) 1 ) FC T x m )
x ^ m = exp ( j a n g ( C F H ( diag ( exp ( j ϕ m ) ) ( F C T ( ϖ α ( x m ) ) ) 1 / 2 ) ) )
Step 9. Merge the subset and the remaining part, use the following rules to get y m , i.e.,
i f fitness ( x ^ m ) < fitness ( x m ) y m = x ^ m e l s e y m = x m
then select the best fitness function and its corresponding waveform p t .
Step 10. Repeat above-mentioned steps 5~9 until num > iter_num or | f i t n e s s ( p t ) f i t n e s s ( p t 1 ) | < ε , then output p t . Otherwise, update t = t + 1 and continue iterating.

4. Simulations and Performance Analysis

Note that the non-convex optimization problem under different initialization is usually providing a different terminus, and hard to obtain the global solution within polynomial time [5,11,26]. For the phase-coded CM waveform design, selecting an effective technique has always been the focus [12,13,14]. In this section, to further assess PSAP’s performance, we firstly assume code length of waveform N = 150 , then the initialized waveform set has
s e t 0 = [ x 1 0 x m 0 x M 0 ] N × M
x m 0 = [ x m 0 ( 1 ) x m 0 ( n ) x m 0 ( N ) ] T N × 1
where M = 20 , x m ( n ) = e j ψ n . In addition, the inertia weight, individual learning factor as well as group learning factor are set as w = 1.5 , c 1 = 0.5 , and c 2 = 0.5 , respectively. The iterations number iter_num is set as 2000, threshold value of | f i t n e s s ( p t ) f i t n e s s ( p t 1 ) | is set as ε = 10 10 .
Next, we take the scenario of suppressing one single area for comparison, i.e., single interval R r = [ 2 : 30 ] . PSAP will be compared with WeCAN, ISAA, RISAAP, and PONLP, by 20 Monte Carlo (MC) trials. Here, we define the averaging autocorrelation sidelobe level (Aver-ACL) and local PSL (LPSL) in suppressed regions as the metric, as follows
Aver - ACL = 1 n u m _ R r k R r 20 log 10 | α k ( x ) | | α 0 ( x ) | ( d B )
LPSL = 20 log 10 ( max ( | α k ( x ) | N ) ) ( d B ) , k R r
for the sake of comparison, all methods or algorithms would be initialized by random phase-coded sequence. Performance comparisons have been shown in Table 1 and Figure 1. Simulations are performed on a PC with 3.40 GHz i7 CPU.
In Table 1 and Figure 1, our proposed PSAP with parallel optimization mechanism has achieved the best performance, PONLP with stochastic gradient optimization has achieved only −139 dB of Aver-ACL and −130 dB of LPSL, while WeCAN achieves −48 dB of Aver-ACL and −35 dB of LPSL. Moreover, ISAA and RISAAP using the mechanism of alternating projection has obtained −228 dB and −158 dB of Aver-ACL, −216 dB and −144 dB of LPSL, respectively. We need to declare that, with 20 log 10 ( ) as the mathematical metric referring to [4,5,12,17], these methods or algorithms conduct their iterating or optimizing with the same stop criteria/thresholds (iter_num is set as 2000, threshold value is set as ε = 10 10 ). Here, the phenomenon with –228 dB or −333 dB may be unnecessary to have so low sidelobe levels for practical engineering, but these low values in mathematical sense would demonstrate some quickly converging or optimizing level of our PSAP framework even for the future electronic countermeasure scenario.
To further discuss the computation complexity of them, here we mainly demonstrate the number of iterations in convergence graph (as shown in Figure 2), and the CPU time consumption in Table 1. In Figure 2, we use log | f i t n e s s ( x t ) f i t n e s s ( x t 1 ) | as the convergence metric of the y-coordinate. Figure 2 has shown obvious converging difference. PSAP uses the idea of particles swarm intelligence, i.e., M = 20, to establish thus cooperative optimization (in CPU model), also occupies much more time than RISAAP, ISAA and PONLP. WeCAN has consumed the longest time 5.11 s. By 20 MC trails, we can see that PSAP has achieved the best robustness performance, which is own to the parallel cooperative mechanism. In addition, these trails and simulations are all based on CPU; when given the GPU condition, PSAP will occupy much less time than others. As seen in Table 1 and Figure 2, WeCAN has slow convergence and PONLP with the steepest descent gradient might stagnate into the local area. ISAA and RISAAP have oscillations in Figure 2 which might attribute to the alternating projection between multi-local areas for the non-convex case.
In addition, we also discuss the Aver-ACL and LPSL comparisons with different code lengths N (i.e., 150, 200, 250, 300) which are shown in Figure 3. No matter N = 150, 200, 250, and 300, our proposed PSAP could obtain the best Aver-ACL and LPSL in different code lengths. Meanwhile, the results of these methods in Figure 3 have shown some similarity that WeCAN might lose the performance for engineering application.
Meanwhile, we also consider the sidelobe suppression in multiple-area case, i.e., R l = [ 2 : 10 ] [ 30 : 40 ] . Here, we assign different suppressing levels with δ 1 = 0 , δ 2 = 10 4 where the former corresponds to the range sidelobe area next to the mainlobe, and the second refers to the farther one. Note that, the interferences near the mainlobe would also produce more effect than the farther one, and affect the detecting performance. Namely, we arrange two suppressing areas with different weights to demonstrate the suppressing levels for different non-cooperative targets. When discussing δ 1 = 0 , δ 2 = 10 4 , the former δ 1 = 0 with relative-low weights means more emphasis on the first area. In Table 2 and Figure 4, we could see that the first area has achieved more excellent performance than the second one, which is due to the weights of the different non-cooperative targets.
In Table 2 with recorded time-consumption data, obviously, our proposed PSAP with 0.84 s has achieved the best Aver-ACL and LPSL using 20 MC trials, but with more time than RISAAP (with 0.16 s), ISAA (with 0.34 s) and PONLP (with 0.45 s). Moreover, WeCAN (with 5.94 s) has lost the performance and cost more time than others.
Moreover, we also consider the multiple-case, i.e., R l = [ 2 : 10 ] [ 30 : 40 ] and arrange two suppressing areas with same weights δ 1 = 0 ,   δ 2 = 0 . In Table 3 and Figure 5, we could see that these two areas have achieved same excellent performance. Table 3 has shown same characteristics as shown in Table 1 and Table 2. Our proposed PSAP has achieved the best result by 20 MC trials. WeCAN has lost its performance.
Note that the initialization of algorithms is indeed significant no matter for cyclic algorithms or alternating projection. As the Gradient Descent (GD) mechanism could accelerate the exploiting for local optimization, and we combine the GD and PSAP together to formulate the PSAPOM and enhance the global robustness. To further assess the performance of PSAPOM and PSAP, despite the random phase-coded sequence, here we assume that these algorithms have been initialized by the Frank-coded or Barker-coded sequence, respectively. For the length N = N ˜ 2 (N = 196), the Frank sequence is given by:
x ( n N ˜ + n ˜ + 1 ) = e j 2 π n n ˜ / N ˜ ; n , n ˜ = 0 , 1 , 2 , N ˜ 1
we also assume another waveform sequence (N = 169) formulated by the 13 Barker sequence. As shown in Figure 6 and Table 4, for both the Frank sequence and Barker sequence, PSAPOM also achieved the same performance as PSAP, but its time consumption, 0.4210 s and 0.2502 s, was less than PSAP.
In these comparisons, we could draw a basic conclusion that our improved PSAP algorithms have a remarkable performance compared to WeCAN, PONLP, ISAA, and RISAAP, which will have a large influence on future practical applications. Moreover, given the sophisticated scenarios, i.e., the unimodular constraint and multiple suppressing areas, our PSAP and PSAPOM have shown a more powerful convergence than these methods. Meanwhile, our simulations have demonstrated that PSAP and PSAPOM have excellent robustness and stability. We may attribute these traits to FFT leverage and swarm particle intelligence due to alternating projection mechanisms.

5. Conclusions

In this paper, a PSAP framework (such as PSAP and PSAPOM) is introduced to design a unimodular CM phase-coded waveform with local low range sidelobes. Therein, PSAP with learning factor and particle-assisted projection could improve the convergence in the non-convex case. Numerical trails and simulations have also provided plenty of analysis to assess the waveforms generated by PSAP, WeCAN, ISAA, RISAAP, and PONLP. Regarding statistical performance, PSAP and PSAPOM via swarm intelligence and parallel optimization idea have achieved outstanding results. Additionally, in this paper, as we only discussed the masking scenarios of detecting static targets despite the Doppler effect of moving targets, in our further research, we will continue designing other phase-coded waveforms considering the non-zero Doppler effect. Moreover, we will also use GPU to accelerate the distributed parallel optimization.

Author Contributions

Conceptualization, X.F. and F.L.; Data curation, X.F.; Formal analysis, F.L., W.C. and Z.Z.; Funding acquisition, Z.Z.; Investigation, F.L. and Y.Z.; Methodology, X.F., F.L., W.C., Z.Z. and Y.Z.; Project administration, X.F. and Y.Z.; Resources, Z.Z.; Supervision, Z.Z.; Validation, F.L. and W.C.; Writing—original draft, X.F. and W.C.; Writing—review and editing, X.F. and Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 42127804.

Data Availability Statement

Not applicable.

Acknowledgments

We also appreciate the anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

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. Guolong, C.; Xianxiang, Y.; Jing, Y.; Yue, F.; Lingjiang, K. An Overview of Waveform Optimization Methods for Cognitive Radar. J. Radars 2019, 8, 537–557. [Google Scholar]
  3. Li, Y.; Vorobyov, S.A. Fast algorithms for designing unimodular waveform(s) with good correlation properties. IEEE Trans. Signal Process. 2017, 66, 1197–1212. [Google Scholar] [CrossRef]
  4. 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]
  5. Feng, X.; Zhao, Y.; Zhou, Z.; Zhao, Z.-F. Waveform design with low range sidelobe and high Doppler tolerance for cognitive radar. Signal Process. 2017, 139, 143–155. [Google Scholar] [CrossRef]
  6. Kajenski, P.J. Design of low-sidelobe phase-coded waveforms. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 2891–2898. [Google Scholar] [CrossRef]
  7. Bu, Y.; Yu, X.; Yang, J.; Fan, T.; Cui, G. A new approach for design of constant modulus discrete phase radar waveform with low WISL. Signal Process. 2021, 187, 108145. [Google Scholar] [CrossRef]
  8. Fan, W.; Liang, J.; Yu, G.; So, H.C.; Lu, G. Minimum local peak sidelobe level waveform design with correlation and/or spectral constraints. Signal Process. 2020, 171, 107450. [Google Scholar] [CrossRef]
  9. Bolhasani, M.; Mehrshahi, E.; Ghorashi, S.A.; Alijani, M.S. Constant envelope waveform design to increase range resolution and SINR in correlated MIMO radar. Signal Process. 2019, 163, 59–65. [Google Scholar] [CrossRef]
  10. Thakur, A.; Saini, D.S. MIMO radar sequence design with constant envelope and low correlation side-lobe levels. AEU-Int. J. Electron. Commun. 2021, 136, 153769. [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] [Green Version]
  12. Li, F.; Zhao, Y.; Qiao, X. A waveform design method for suppressing range sidelobes in desired intervals. Signal Process. 2014, 96, 203–211. [Google Scholar] [CrossRef]
  13. Ge, P.; Cui, G.; Karbasi, S.M.; Kong, L.; Yang, J. A template fitting approach for cognitive unimodular sequence design. Signal Process. 2016, 128, 360–368. [Google Scholar] [CrossRef]
  14. Aubry, A.; De Maio, A.; Jiang, B.; Zhang, S. Ambiguity function shaping for cognitive radar via complex quartic optimization. IEEE Trans. Signal Process. 2013, 61, 5603–5619. [Google Scholar] [CrossRef]
  15. Kerahroodi, M.A.; Aubry, A.; De Maio, 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] [CrossRef]
  16. Esmaeili-Najafabadi, H.; Leung, H.; Moo, P.W. Unimodular waveform design with desired ambiguity function for cognitive radar. IEEE Trans. Aerosp. Electron. Syst. 2019, 56, 2489–2496. [Google Scholar] [CrossRef]
  17. He, H.; Li, J.; Stoica, P. Waveform Design for Active Sensing Systems: A Computational Approach; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  18. Cheng, Z.; Liao, B.; He, Z.; Li, J.; Han, C. A nonlinear-ADMM method for designing MIMO radar constant modulus waveform with low correlation sidelobes. Signal Process. 2019, 159, 93–103. [Google Scholar] [CrossRef]
  19. Liu, Y.; Jiu, B.; Liu, H. ADMM-based transmit beampattern synthesis for antenna arrays under a constant modulus constraint. Signal Process. 2020, 171, 107529. [Google Scholar] [CrossRef]
  20. Patton, L.K.; Rigling, B.D. Phase retrieval for radar waveform optimization. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 3287–3302. [Google Scholar] [CrossRef]
  21. Zhang, X.; Wang, X. Waveform design with controllable modulus dynamic range under spectral constraints. Signal Process. 2021, 189, 108285. [Google Scholar] [CrossRef]
  22. Wu, Z.J.; Zhou, Z.Q.; Wang, C.X.; Li, Y.-C.; Zhao, Z.-F. Doppler resilient complementary waveform design for active sensing. IEEE Sens. J. 2020, 20, 9963–9976. [Google Scholar] [CrossRef]
  23. Wang, Y.C.; Dong, L.; Xue, X.; Yi, K.-C. On the design of constant modulus sequences with low correlation sidelobes levels. IEEE Commun. Lett. 2012, 16, 462–465. [Google Scholar] [CrossRef]
  24. Xue-Bin, S.; Zhan-Min, L.; Cheng-Lin, Z.; Zheng, Z. Cognitive UWB Pulse Waveform Design Based on Particle Swarm Optimization. Adhoc Sens. Wirel. Netw. 2012, 16, 215–228. [Google Scholar]
  25. Wang, S. Efficient heuristic method of search for binary sequences with good aperiodic autocorrelations. Electron. Lett. 2008, 44, 731–732. [Google Scholar] [CrossRef]
  26. Zhao, D.; Wei, Y.; Liu, Y. Design of unimodular sequence train with low central and recurrent autocorrelations. IET Radar Sonar Navig. 2018, 13, 45–49. [Google Scholar] [CrossRef]
  27. Tanweer, M.R.; Suresh, S.; Sundararajan, N. Self regulating particle swarm optimization algorithm. Inf. Sci. 2015, 294, 182–202. [Google Scholar] [CrossRef]
  28. Zhang, H.; Xie, J.; Zong, B. Bi-objective particle swarm optimization algorithm for the search and track tasks in the distributed multiple-input and multiple-output radar. Appl. Soft Comput. 2021, 101, 107000. [Google Scholar] [CrossRef]
  29. Lin, R.; Soltanalian, M.; Tang, B.; Li, J. Efficient design of binary sequences with low autocorrelation sidelobes. IEEE Trans. Signal Process. 2019, 67, 6397–6410. [Google Scholar] [CrossRef]
  30. Levanon, N.; Mozeson, E. Radar Signals; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  31. Klauder, J.R.; Price, A.C.; Darlington, S.; Albersheim, W.J. The theory and design of chirp radars. Bell Syst. Tech. J. 1960, 39, 745–808. [Google Scholar] [CrossRef]
Figure 1. ACL comparison for suppressing sidelobes in a single area.
Figure 1. ACL comparison for suppressing sidelobes in a single area.
Remotesensing 14 04186 g001
Figure 2. Convergence comparison of different methods.
Figure 2. Convergence comparison of different methods.
Remotesensing 14 04186 g002
Figure 3. Aver-ACL and LPSL comparison with different code lengths N. (a) Aver-ACL comparison of RISAAP, ISAA, WeCAN, PONLP, and PSAP. (b) LPSL comparison of RISAAP, ISAA, WeCAN, PONLP, and PSAP.
Figure 3. Aver-ACL and LPSL comparison with different code lengths N. (a) Aver-ACL comparison of RISAAP, ISAA, WeCAN, PONLP, and PSAP. (b) LPSL comparison of RISAAP, ISAA, WeCAN, PONLP, and PSAP.
Remotesensing 14 04186 g003aRemotesensing 14 04186 g003b
Figure 4. ACL comparison for suppressing sidelobes for multiple areas.
Figure 4. ACL comparison for suppressing sidelobes for multiple areas.
Remotesensing 14 04186 g004
Figure 5. ACL comparison of suppressing sidelobes for multiple areas.
Figure 5. ACL comparison of suppressing sidelobes for multiple areas.
Remotesensing 14 04186 g005
Figure 6. ACL performance comparison for suppressing multiple sidelobe areas. (a) Description of PSAPOM and PSAP initialized by Frank-coded sequence; (b) Description of PSAPOM and PSAP initialized by the 13 Barker-coded sequence.
Figure 6. ACL performance comparison for suppressing multiple sidelobe areas. (a) Description of PSAPOM and PSAP initialized by Frank-coded sequence; (b) Description of PSAPOM and PSAP initialized by the 13 Barker-coded sequence.
Remotesensing 14 04186 g006aRemotesensing 14 04186 g006b
Table 1. Performance comparison for suppressing single area sidelobes of different algorithms.
Table 1. Performance comparison for suppressing single area sidelobes of different algorithms.
AlgorithmsSingle Area
Aver-ACL (dB)LPSL (dB)Time Consumption (s)
RISAAP−228.4560−216.54900.14
ISAA−158.0292−144.88210.35
WeCAN−48.5622−35.65735.11
PONLP−139.6091−130.27620.48
PSAP−333.4510−321.84620.76
Table 2. Performance comparison for suppressing multiple area sidelobes with different weights.
Table 2. Performance comparison for suppressing multiple area sidelobes with different weights.
AlgorithmsMultiple Area
Aver-ACL
in 1# Area (dB)
LPSL
in 1# Area (dB)
Aver-ACL
in 2# Area (dB)
LPSL
in 2# Area (dB)
Time Consumption (s)
RISAAP−233.2183−220.6892−163.5196−163.51100.16
ISAA−119.9062−117.8060−100.1328−93.86270.34
WeCAN−61.2209−47.9854−41.1983−34.94655.94
PONLP−147.6981−141.6530−116.1224−112.52600.45
PSAP−336.3241−329.4050−169.5424−169.54200.84
Table 3. Performance comparison for suppressing multiple area sidelobes with same weights.
Table 3. Performance comparison for suppressing multiple area sidelobes with same weights.
AlgorithmsMultiple Area
Aver-ACL
in 1# Area (dB)
LPSL
in 1# Area (dB)
Aver-ACL
in 2# Area (dB)
LPSL
in 2# Area (dB)
RISAAP−238.1425−232.5312−236.9185−235.2602
ISAA−137.7802−134.5042−134.4622−131.5260
WeCAN−59.4995−50.2199−53.3124−48.0641
PONLP−151.2146−142.0220−149.6574−142.9440
PSAP−340.8439−326.2900−336.8021−323.6540
Table 4. Algorithm (initialized by the Frank or Barker sequence) comparison for suppressing range sidelobes.
Table 4. Algorithm (initialized by the Frank or Barker sequence) comparison for suppressing range sidelobes.
AlgorithmsFrank SequenceBarker Sequence
Aver-ACL in 1# Area (dB)Aver-ACL in 2# Area (dB)Aver-ACL in 1# Area (dB)Aver-ACL in 2# Area (dB)
PSAPOM−321.1162−320.7214−326.3520−323.2342
PSAP−321.4633−321.2361−325.7421−323.1652
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Feng, X.; Zhao, Z.; Li, F.; Cui, W.; Zhao, Y. Radar Phase-Coded Waveform Design with Local Low Range Sidelobes Based on Particle Swarm-Assisted Projection Optimization. Remote Sens. 2022, 14, 4186. https://doi.org/10.3390/rs14174186

AMA Style

Feng X, Zhao Z, Li F, Cui W, Zhao Y. Radar Phase-Coded Waveform Design with Local Low Range Sidelobes Based on Particle Swarm-Assisted Projection Optimization. Remote Sensing. 2022; 14(17):4186. https://doi.org/10.3390/rs14174186

Chicago/Turabian Style

Feng, Xiang, Zhanfeng Zhao, Fengcong Li, Wenqing Cui, and Yinan Zhao. 2022. "Radar Phase-Coded Waveform Design with Local Low Range Sidelobes Based on Particle Swarm-Assisted Projection Optimization" Remote Sensing 14, no. 17: 4186. https://doi.org/10.3390/rs14174186

APA Style

Feng, X., Zhao, Z., Li, F., Cui, W., & Zhao, Y. (2022). Radar Phase-Coded Waveform Design with Local Low Range Sidelobes Based on Particle Swarm-Assisted Projection Optimization. Remote Sensing, 14(17), 4186. https://doi.org/10.3390/rs14174186

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