Next Article in Journal
Integrating Multiple Models Using Image-as-Documents Approach for Recognizing Fine-Grained Home Contexts
Next Article in Special Issue
A Multi-Scale U-Shaped Convolution Auto-Encoder Based on Pyramid Pooling Module for Object Recognition in Synthetic Aperture Radar Images
Previous Article in Journal
Optical Fiber FP Sensor for Simultaneous Measurement of Refractive Index and Temperature Based on the Empirical Mode Decomposition Algorithm
Previous Article in Special Issue
On the Slow-Time k-Space and its Augmentation in Doppler Radar Tomography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compressive Sensing-Based Bandwidth Stitching for Multichannel Microwave Radars

1
Intelligence, Surveillance and Space Division, Defence Science and Technology Group, Edinburgh, SA 5111, Australia
2
Maritime Division, Defence Science and Technology Group, Edinburgh, SA 5111, Australia
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(3), 665; https://doi.org/10.3390/s20030665
Submission received: 18 November 2019 / Revised: 17 January 2020 / Accepted: 19 January 2020 / Published: 24 January 2020
(This article belongs to the Special Issue Recent Advancements in Radar Imaging and Sensing Technology)

Abstract

:
The problem of obtaining high range resolution (HRR) profiles for non-cooperative target recognition by coherently combining data from narrowband radars was investigated using sparse reconstruction techniques. If the radars concerned operate within different frequency bands, then this process increases the overall effective bandwidth and consequently enhances resolution. The case of unknown range offsets occurring between the radars’ range profiles due to incorrect temporal and spatial synchronisation between the radars was considered, and the use of both pruned orthogonal matching pursuit and refined l 1 -norm regularisation solvers was explored to estimate the offsets between the radars’ channels so as to attain the necessary coherence for combining their data. The proposed techniques were demonstrated and compared using simulated radar data.

1. Introduction

The construction of high range resolution profiles (HRRP) of targets is a precursor to feature extraction for automatic target recognition (ATR), and normally requires the employment of a high-bandwidth waveform following detection by a lower resolution radar mode. Examples of recent papers in the non-cooperative target recognition (NCTR) literature focusing on feature extraction for ATR following HRRP construction are [1,2,3]. This paper considers the problem of HRRP construction, but using low resolution radars operating in different frequency bands for the purpose of combining their signals to achieve a higher resolution, and examines the problem of their data not being mutually coherent.
The ability to acquire high resolution range profiles of targets has improved over time as hardware capability has developed, with higher resolution being achieved by increasing the time-bandwidth product. In the early approaches, for narrowband radars with very limited instantaneous bandwidth, stepped-frequency waveforms were used with a single I,Q (that is, baseband quadrature) signal sample received after each frequency step. The set of samples is effectively used for the Fourier transform of the slant range profile, enabling the range response to be obtained simply by implementing an inverse Fourier transform (see, e.g., [4]). The greater the frequency range, the higher the range resolution, but the downside is that the burst of pulses can be so long that a scatterer may migrate between range cells, causing smearing of the range profile, and therefore requiring range compensation. An example of a recent paper involving the use of stepped-frequency waveforms is [5], and recent papers which have investigated the effects of target motion and aspect sensitivity are [6,7].
An advance on the stepped-frequency approach is to increase the time-bandwidth product using stretch processing, whereby a wideband LFM waveform is transmitted, and pulse-compression is achieved in hardware by mixing the received signal with an extended replica of the transmitted waveform. A point target at a particular range will manifest itself as a single frequency which is proportional to its range. The result, which is digitally sampled in time in order to facilitate the identification of the frequency components, is therefore a superposition of discrete frequencies, each corresponding to a point target at a different range. An application of the inverse Fourier transform again recovers the range profile (see, e.g., [4,8,9]).
A spectral analysis technique to improve range resolution was proposed in [10,11] based on autoregressive linear prediction. The main idea is to combine the mutually coherent signals received from multiple waveforms transmitted sequentially or concurrently, which have widely separate carrier frequencies or may even occupy entirely different frequency bands. Viewed in the spectral domain, the received signals from individual waveforms may be seen to occupy discrete wavebands which are separate or contiguous. If contiguous, then they can potentially be coherently combined to synthesise the signals that would have been received from a single wider bandwidth waveform in the manner presented in [12]. If separate, then presumably this coherent combination of signals would still be feasible, as would be the interpolation of the frequency response in the gaps between the bands under the a priori assumption that the signals are returned from discrete scatterers using spectral estimation techniques (see, e.g., [13]). Alternatively, the signals from different frequency bands can be jointly processed without explicitly filling the gaps between the bands. Since no new synthetic frequency band is actually constructed, this approach can be referred to as bandwidth stitching to distinguish it from bandwidth interpolation and extrapolation. The main challenge of this approach lies in the presence of phase errors in different frequency bands resulting from post-processed motion compensation which is often carried out separately for each frequency band. In this paper, we focus on the problem of bandwidth stitching for radar high-resolution range profiling and explore the use of sparse reconstruction to deal with the phase error problem.
It is convenient to formulate these ideas in the spectral domain, within which point scatterers appear as discrete sinusoids and which are amenable to analysis by spectral estimation techniques such as autoregression, as presented in [10,11]. Compressive sensing and sparse reconstruction, however, provide for the possibility of alternative signal representations, potentially allowing for greater flexibility and discriminating between signals of physical origin and receiver noise [13,14,15,16,17,18,19,20]. Instances of non-sinusoidal signals are waveforms in fast-time and signals returned from rotating objects when the angle of rotation is large. Compressive sensing and sparse reconstruction can also handle the situation corresponding to data being non-uniformly sampled in time or space, such as non-uniform PRF (pulse repetition frequency) waveforms and random sparse arrays.
Compressive sensing and sparse reconstruction were exploited in [17,19,21,22] to address the problem of gaps in the data both in slow-time and in frequency for inverse synthetic aperture radar (ISAR) imaging. However, these works assumed that the data were coherent across different sub-bands and that there were no model uncertainties. The work [23] took account of the possible lack of mutual coherence between the radars operating on the different sub-bands arising from incorrect timing synchronisation, or, equivalently, errors in antenna phase’s centre-relative locations. This is achieved by fitting an ultra-wideband all-pole signal model to the mutually-coherent sub-bands, which is then used for bandwidth interpolation and extrapolation prior to recovering the range profile by means of an inverse Fourier transform. This paper, however, proposes the use of compressive sensing and sparse reconstruction to deal with the non-coherence problem between different sub-bands.
To address the sub-band non-coherence problem, two different approaches were explored: (i) greedy pursuit and (ii) l 1 -norm regularisation. In the first approach, pruned orthogonal matching pursuit (POMP) [24], which was originally developed for micro-Doppler parameter estimation, is adopted to deal with the dictionary mismatch which is due to the phase errors in each sub-band resulting from the motion-compensation post processing. The main idea is to parameterise the dictionary as a function of the phase errors and to construct multiple realisations of the dictionary. A selective learning process is then used to discard the dictionaries which correspond to incorrect values of phase errors. Since a straight application of the POMP algorithm to the problem under consideration would have been computationally expensive, we first applied the POMP algorithm pairwise to sub-bands in order to estimate the phase errors, and then utilised the conventional OMP algorithm to determine the range profile based on the estimated phase error values. In the second approach, an l 1 -norm regularisation problem can be solved jointly both for the range profile vector and the phase errors [25]. The work [25] offers two solutions for joint synthetic aperture radar (SAR) imaging and phase error correction. The first solution is not applicable to the problem under consideration because a constant phase error in each sub-band is assumed. On the other hand, the second solution considers general arbitrary phase errors, and thus can be applied to our problem. Since the phase errors within each sub-band are a linear function of the range error coming from post-processing motion compensation, we also present refined variants of the second solution of [25] to take into account this underlying structure of the phase error.
The paper is organised as follows. Section 2 formulates the problem of bandwidth stitching for HRRP in the presence of phase errors. The POMP algorithm is applied in Section 3 to the bandwidth stitching problem under consideration. Section 4 presents l 1 -norm regularisation solvers. Numerical performance comparisons are provided in Section 5 and conclusions are drawn in Section 6.

2. Problem Formulation

Consider a multistatic radar system consisting of M radar channels on different and distinct frequency sub-bands, approximately co-located, and illuminating a common target, such that their radar lines of sight (LoS) coincide but their range profiles are out of alignment. Each channel can individually produce a one dimensional range profile of the target, but with relatively coarse resolution. The bandwidth stitching problem can be briefly stated as follows: for the M generally non-coherent channels, the aim is to coherently combine, or “stitch,“ the channels together so that they can effectively produce a single range profile with resolution corresponding to the combined overall signal bandwidth.
Let f m , n ( n = 1 , , N m ) denote the n th frequency bin of the m th channel ( m = 1 , , M ). Here, N m is the number of frequency bins in the m th channel. We base the formulation on the point-scatterer model and assume that the target can be defined as consisting of K scattering centres at local line-of-sight coordinates x k (or local “down ranges”) and having complex-valued reflectivity coefficients α k , which are also assumed frequency-independent. The down-converted, pulse-compressed, motion-compensated signals received in each channel, in the frequency domain, can be written as
S m = [ , S m , n , ] n = 1 , , N m T ,
where superscript T denotes the transpose operation, and
S m , n = | A ( f m , n ) | 2 exp 4 π j f m , n c Δ R m k = 1 K α k exp 4 π j f m , n c x k .
Here, A ( f n , m ) represents the transmit radar waveform, the squared amplitude resulting from pulse compression processing; constant c denotes the speed of light; and Δ R m accounts for the range errors in the motion-compensation processing. Bandwidth stitching in this context amounts to estimating these phase errors as accurately as possible.
The signals S m in (1) can be rewritten in a more compact form as
S m = Λ m F m α ,
where
(4) Λ m = diag , exp 4 π j f m , n c Δ R m , n = 1 , , N m (5) F m = [ , F m , k , ] k = 1 , , K (6) F m , k = , exp 4 π j f m , n c x k , n = 1 , , N m T (7) α = [ , α k , ] k = 1 , , K T
Here, “diag” denotes a diagonal matrix; Λ m is referred to as the phase error matrix, of dimension N m × N m ; F m and α are respectively, dimensions N m × K and K × 1 ; S m is a column vector of dimension N m × 1 ; and the dagger symbol † refers to the K actual scatterers on the target.
To apply the sparse representation techniques of compressive sensing, we discretise the target’s local range coordinate x using a regularly-spaced range grid { x l } for l = 1 , , L x , with L x K , and construct the N m × L x dictionary matrices
F m = [ , F m , l , ] l = 1 , , L x ,
where
F m , l = , exp 4 π j f m , n c x l , n = 1 , , N m T
Are the “atoms” of dictionary F m in the frequency domain. The corresponding range profile vector
α = [ , α l , ] l = 1 , , L x T
Spans over the range grid { x l } . The received signal S m can also be written as
S m = Λ m F m α .
Since the target usually contains only a small number of dominant scattering centres relative to the total number of range resolution cells, the range profile α can be considered sparse (i.e., containing a small number of non-zero elements).
In the presence of unknown noise, (11) becomes
S ˜ m = Λ m F m α + n m .
where n m is the additive noise for channel m. Stacking up the individual channel signals S ˜ m , m = 1 , , M , gives
S ˜ = Λ F α + n ,
where
(14) S ˜ = [ , S ˜ m T , ] m = 1 , , M T (15) Λ = diag { , Λ m , } m = 1 , , M (16) F = [ , F m T , ] m = 1 , , M T (17) n = [ , n m T , ] m = 1 , , M T .
Note that S ˜ and n are column vectors of size ( m N m ) × 1 ; diagonal phase error matrix Λ is of size ( m N m ) × ( m N m ) ; dictionary matrix F is ( m N m ) × L x ; and the range profile α is again a column vector of size L x × 1 . Stacking the received signals amounts to a vertical stacking of the dictionary matrices from all channels and a diagonal concatenation of the corresponding phase error matrices. The stacking of multiple channels in this manner can improve the estimation accuracy for α , as will be demonstrated later in the paper.
A problem statement can thus be expressed as follows: given S ˜ as the measured signal,
find α and Λ , subject to S ˜ Λ F α , α is sparse .
The estimation of α over { x l } is the process of range profiling, giving the main desired output, whereas the estimation of the phase error matrix Λ is really only a necessary intermediate result; it is a function of Δ R 1 , Δ R 2 , , and Δ R M (recall that Δ R m is the range estimation error resulting from the motion-compensation process for channel m). Furthermore, since these errors arise from a lack of precise knowledge of the relative locations of the radar channels’ phase centres and are small relative to a range resolution cell, we may assume, without loss of generality, that Δ R 1 = 0 .

3. Greedy Pursuit Solutions

In this section, we adopt the pruned OMP (POMP) technique, which was originally proposed for micro-Doppler parameter estimation, [24], for the problem of bandwidth stitching for range profiling. We start with the simplest case of two channels and then generalise it to the multiple channel case.

3.1. The Two-Channel Case

For this case, the signal model in (13) can be expressed as
S ˜ = Λ ( Δ R 2 ) F α + n ,
where Λ ( Δ R 2 ) is a function of the single unknown relative range error Δ R 2 ,
Λ ( Δ R 2 ) = diag { I N 1 , Λ 2 ( Δ R 2 ) } ,
with
Λ 2 ( Δ R 2 ) = diag , exp 4 π j f 2 , n c Δ R 2 , n = 1 , , N 2 .
In addition to the sparse range profile vector α , Δ R 2 is the only additional unknown parameter to be estimated. Let us rewrite (19) as
S ˜ = Φ ( Δ R 2 ) α + n
where
Φ ( Δ R 2 ) = Λ ( Δ R 2 ) F .
In this form, the problem can be viewed as a joint sparse reconstruction and parameter estimation problem with the parametric dictionary Φ ( Δ R 2 ) itself a function of the parameter Δ R 2 . This can be considered as a special dictionary learning problem where the objective is to solve simultaneously for both the sparse solution of α and the range error Δ R 2 .
To solve this problem, we adopt the POMP technique [24], which embeds a pruning operation into the iterative process of OMP. The main idea of POMP is to construct multiple realisations of the dictionary Φ based on a number L Δ of candidate values of Δ R 2 ; the OMP algorithm is applied to each dictionary realisation to find the atom which correlates most strongly with the current residual for that dictionary, and to recompute the residual with that atom’s contribution to the residual removed. To overcome possible excessive computations arising from outlier candidate values of Δ R 2 , a pruning operation is performed to exclude the half of the dictionaries which yield the largest residual errors, until a single dictionary realisation remains. The OMP iterations for the remaining dictionary are continued until a termination criterion is satisfied. The candidate value of Δ R 2 corresponding to this dictionary gives the final estimate of Δ R 2 . The basic POMP algorithm is summarised in Table 1 and its computational cost is shown in Appendix C to be of the order of L Δ N m L x .

3.2. The Multi-Channel Case

For this case, the noisy signal model (19) becomes
S ˜ = Φ ( Δ R 2 , , Δ R M ) α + n
where
Φ ( Δ R 2 , , Δ R M ) = diag { I N 1 , Λ 2 ( Δ R 2 ) , , Λ M ( Δ R M ) } F .
Here, the dictionary matrix is a function of the ( M 1 ) unknowns Δ R 2 , Δ R 3 , and Δ R M .
The POMP algorithm could be extended to multiple channels by computing candidate dictionaries based on a multi-dimensional grid of candidate values for Δ R 2 , Δ R 3 , , and Δ R M . The grid would consist of a total of ( M 1 ) dimensions, where the mth dimension corresponds to the unknown range error Δ R m + 1 of the ( m + 1 ) channel. Note that only a one-dimensional grid for Δ R 2 is required for the case of two channels. However, the cardinality of the dictionary set is exponentially dependent on the number of available channels; i.e., V ( M 1 ) , where V denotes the number of grid points in each parameter dimension. As a result, although this extension would be simple and straightforward, it is computationally expensive.
To alleviate this computational burden, we instead apply the POMP algorithm pairwise to channels in order to estimate the range errors Δ R 2 , , Δ R M relative to the first channel, and then utilise the conventional OMP algorithm to determine the range profile vector α based on these estimated values of Δ R 2 , , Δ R M . The procedure is summarised as below:
STEP 1: Estimation of range errors.
  • For each pair between the 1st and m th channel, m { 2 , , M } :
    -
    Calculate input signal:
    S ˜ 1 m = [ S ˜ 1 T , S ˜ m T ] T
    F 1 m = [ F 1 T , F m T ] T .
    -
    Construct candidate dictionaries based on a grid of L candidate values of Δ R m ( l ) ( l = 1 , , L ):
    Φ 1 m , l = diag { I N 1 , Λ m ( Δ R m ( l ) ) } F 1 m .
    -
    Perform POMP given S ˜ 1 m and Φ 1 m , 1 , , Φ 1 m , L to obtain an estimate of Δ R m (denoted as Δ R ^ m ).
 End for.
STEP 2: Estimation of range profile vector.
  • Compute signal and dictionary.
    (26) S ˜ = [ , S ˜ m T , ] m = 1 , , M T (27) F = [ , F m T , ] m = 1 , , M T (28) Λ = diag { , Λ m ( Δ R ^ m ) , } m = 1 , , M (29) Φ = Λ F .
  • Estimate α using OMP given S ˜ and Φ .
The computational cost of the general POMP algorithm is shown in Appendix C to be of the order of ( M 1 ) L Δ N m L x .

4. L 1 -Norm Regularisation Approach

The sparse reconstruction problem (18) can be solved via the following l 1 regularised optimisation:
min α , Λ S Λ F α 2 2 + μ α 1 ,
where μ is a regularisation parameter. It should be emphasized that this is not a conventional l 1 regularisation formulation because of the unknown phase error matrix Λ resulting from the estimation error of the motion-compensation phase. Therefore, Λ must be jointly estimated with α :
{ α ^ , Λ ^ } = arg min α , Λ S Λ F α 2 2 + μ α 1 .
Two solutions for this joint estimation problem were presented in [25]. The first solution assumes that the phase error matrix for the m th sub-band is modelled as
Λ m = exp { j ϕ m } I N m × N m .
In other words, the phase errors for different frequency bins of a particular sub-band are identical. However, this assumption is invalid in the problem under consideration because the phase error is a function of frequency and thus has different values for different frequency bins. Therefore, that solution is not applicable in this case. The second solution considers a general phase error matrix
Λ m = diag { , exp { j ϕ m , n } , } n = 1 , , N m
where the phase errors ϕ m , n can be arbitrary. Although this solution can be used, it does not exploit the underlying structure of the phase errors; i.e., ϕ m , n = 4 π f m , n c Δ R m . In what follows, we will also present other refined versions, building on the second solution of [25], while exploiting prior knowledge of the structure of the phase error.
The l 1 norm can be approximated as [26,27,28,29]:
α 1 l = 1 L | α l | 2 + δ 1 / 2
In order to overcome the nondifferentiably of the l 1 norm at the origin. Here, δ is a small non-negative parameter. Using this approximation, the minimisation problem in (31) becomes
{ α ^ , Λ ^ } = arg min α , Λ S Λ F α 2 2 + λ l = 1 L | α l | 2 + δ 1 / 2 .
The solution of (35) tends to the solution of (31) as δ approaches zero. Therefore, a small value of δ should be used to ensure the validity of this approximation. The quasi-Newton approach can be adopted to solve the modified l 1 regularised optimisation (35), as below.
The gradient of the objective function of (35) is given by
( α ) = H ( α ) α 2 F H Λ H S ,
where the superscript H denotes the Hermitian transpose operation. Here, H is the Hessian matrix given by
H ( α ) = 2 F H Λ H Λ F + λ W ( α ) = 2 F H F + λ W ( α ) ,
where
W ( α ) = diag , | α l | 2 + δ 1 / 2 , .
Since the Hessian matrix is a function of the unknown α , the minimisation (35) is solved iteratively. Given the estimates α ^ ( i ) and Λ ^ ( i ) from a previous iteration i, the new solutions at iteration i + 1 are obtained in the following two steps.
  • Calculate α ^ ( i + 1 ) by setting ( α ) = 0 given H ( α ^ ( i ) ) and Λ ^ ( i ) :
    α ^ ( i + 1 ) = 2 ( H ( α ^ ( i ) ) ) 1 F H ( Λ ^ ( i ) ) H S = F H F + 1 2 λ W ( α ^ ( i ) ) 1 F H ( Λ ^ ( i ) ) H S .
  • Calculate Λ ^ ( i + 1 ) given α ^ ( i + 1 ) . The phase error matrix Λ ^ ( i + 1 ) is obtained by solving:
    Λ ^ ( i + 1 ) = arg min Λ S Λ F α ^ ( i + 1 ) 2 2
    or equivalently
    Λ ^ m ( i + 1 ) = arg min Λ m S m Λ m F m α ^ ( i + 1 ) 2 2
    for m = 2 , , M . Note that Λ 1 = I N 1 (since Δ R 1 = 0 ); thus, no estimation is required for Λ 1 .
The algorithm may be halted when the objective function falls below a threshold, or when a maximum number of iterations is reached, or when the relative change in the objective function falls below a threshold.
Various methods for calculating the phase error matrix Λ ^ ( i + 1 ) in Step 2 are given in the following sections.

4.1. Unstructured Approach

Ignoring the underlying structure of the phase errors, i.e., ϕ m , n = 4 π f m , n c Δ R m , Λ ^ m can be considered as a diagonal matrix with arbitrary elements ϕ m , n :
Λ m = diag { , exp { j ϕ m , n } , } n = 1 , , N m .
Therefore, ϕ m , n can be estimated as [25]
ϕ ^ m , n ( i + 1 ) = tan 1 { S m , n Y ^ m , n * ( i + 1 ) } { S m , n Y ^ m , n * ( i + 1 ) }
where { · } and { · } denote operations to extract the imaginary and real parts of a complex number, and tan 1 stands for a four-quadrant arctangent operation. Here, Y ^ m , n ( i + 1 ) is the n th element of Y ^ m ( i + 1 ) which is defined as Y ^ m ( i + 1 ) = F m α ^ ( i + 1 ) . As a result, we obtain:
Λ ^ m ( i + 1 ) = diag { , exp { j ϕ ^ m , n ( i + 1 ) } , } n = 1 , , N m .

4.2. Gauss–Newton Approach

Taking into account the underlying structure of the phase errors, Λ m ( Δ R m ) is a function of Δ R m , and the minimisation (41) can be re-expressed as
Δ R ^ m ( i + 1 ) = arg min Δ R m S m Λ m ( Δ R m ) F m α ^ ( i + 1 ) 2 2 .
By letting e m = S m Λ m ( Δ R m ) F m α ^ ( i + 1 ) , we have
e m = [ , e m , n , ] n = 1 , , N m T
where
e m , n = S m , n U m , n exp 4 π j f m , n c Δ R m l = 1 L α ^ l ( i + 1 ) exp 4 π j f m , n c x l
and
U m , n = exp { j ϕ m , n } .
Here α ^ l ( i + 1 ) is the l th element of α ^ ( i + 1 ) . As we are estimating real quantities, it is more convenient to reformulate the problem as the minimisation of a real function in order to apply the Gauss–Newton. The details of the Gauss–Newton algorithm for updating Δ R ^ m ( i + 1 ) are given in Appendix A.
Using Δ R ^ m ( i + 1 ) , we obtain the phase error matrix as
Λ ^ m ( i + 1 ) = diag , exp 4 π j Δ R ^ m ( i + 1 ) c f m , n , n = 1 , , N m .

4.3. Linear Regression-Based Approach

By noting that ϕ m , n = 4 π Δ R m c f m , n , the gradient 4 π Δ R m c can be calculated via a linear least squares estimator using ϕ ^ m , n obtained from (43) [30]. Specifically, we have
4 π Δ R ^ m ( i + 1 ) c , ϕ ^ m ( i + 1 ) T = ( A m T A m ) 1 A m T b m ( i + 1 )
where
(51a) A m = [ , A m , n , ] n = 1 , , N m T , with A m , n = [ f m , n , 1 ] (51b) b m ( i + 1 ) = [ , ϕ ^ m , n unwrapped ( i + 1 ) , ] n = 1 , , N m T .
Note that ϕ ^ m , n unwrapped ( i + 1 ) is the unwrapped version of ϕ ^ m , n ( i + 1 ) and ϕ ^ m ( i + 1 ) is an estimate for the initial phase ϕ m ( i + 1 ) which results from the unwrapping process. From (50), Δ R ^ m ( i + 1 ) is obtained and can then be used for computing Λ ^ m ( i + 1 ) as in (49).

4.4. Differenced-Phase-Based Approach

Subtracting the estimated phase errors of two successive frequency bins, we obtain:
4 π ( f m , n + 1 f m , n ) c Δ R m = ϕ m , n + 1 ϕ m , n .
Therefore, Δ R m can be estimated as [30]
Δ R ^ m ( i + 1 ) = c 4 Δ f m ( N m 1 ) n = 1 N m 1 Δ ϕ m , n ( i + 1 )
where
Δ ϕ m , n ( i + 1 ) = tan 1 sin ( ϕ ^ m , n + 1 ( i + 1 ) ϕ ^ m , n ( i + 1 ) ) cos ( ϕ ^ m , n + 1 ( i + 1 ) ϕ ^ m , n ( i + 1 ) ) .
Note that the four-quadrant arctangent has been used here to handle the phase wrapping. An estimate of the phase error matrix Λ ^ m ( i + 1 ) is now obtained as in (49) using Δ R ^ m ( i + 1 ) in (53).

5. Simulation and Discussion

Numerical simulations are presented in this section to evaluate the performance of the methods described in previous sections.

5.1. Scenario 1: Two Sub-Bands

We consider a synthetic scenario with two sub-bands at carrier frequencies of f 1 = 6 GHz and f 2 = 8 GHz, each having a bandwidth of B = 300 MHz and 64 frequency steps (i.e., N = N 1 = N 2 = 64 ). The range profile is discretised over a grid with a length of ( N 1 ) c / ( 2 B ) = 31.5 m and a grid step of Δ Grid = c / ( 10 B ) = 0.1 m. We consider a far-field target consisting of six point scatterers which are aligned with the grid. Figure 1 plots the true range profile of the target. We set Δ R 2 = 2.78 Δ Grid for the case of existing phase errors. The signal-to-noise ratio is set to 10 dB.
Figure 2 compares the reconstructed range profiles obtained by the conventional OMP algorithm without and with the presence of phase errors. The OMP is terminated when the signal residual reaches the noise level or after 15 iterations have been carried out. We observe that OMP successfully reconstructs the range profile of the target by correctly identifying the scatterers of the target with accurate range and coefficient estimates when no phase errors exist. However, OMP provides unsatisfactory results in the presence of phase errors, where the reconstructed image is observed as exhibiting many spurious scatterers. Similar observations are obtained for the results obtained by the conventional l 1 -norm regularised optimisation solver (without phase error correction), as shown in Figure 3. Here, we set δ = 10 5 and λ = 0.001 max | c l | where c l is the l th element of c = F H S ˜ . The l 1 -norm regularised optimisation solver is stopped if the relative change in the l 2 -norm of the range profile vector α falls below 10 5 or after it reaches 500 iterations. The performance degradation of these conventional sparse reconstruction techniques is not unexpected, since they were not originally developed to cope with dictionary mismatch arising from the presence of the phase errors.
Figure 4 shows the reconstructed range profile obtained by the POMP. POMP constructs candidate dictionaries based on a grid of Δ R 2 with a grid step size of Δ Grid / 100 . The same stopping criteria of OMP is used for POMP. We observe that POMP produces a range profile which is almost identical to the ground truth, thereby demonstrating the effectiveness of POMP in terms of dealing with the phase errors between different sub-bands.
Figure 5 shows the results obtained by different l 1 -norm regularised optimisation solvers with phase error correction, as presented in Section 4. The same parameters and stopping criteria of the conventional l 1 -norm regularised optimisation solver as described above are used in the simulations. Although these algorithms exhibit some improvements over the conventional l 1 -norm regularisation (i.e., without phase error correction), they provide poorer results compared to that of POMP. Specifically, the peaks of the reconstructed range profiles obtained by these algorithms only appear close to but not exactly at the true scatterer positions. In addition, the magnitudes of the peaks are much smaller than the ground truth values.
The inferior performances of these algorithms can be explained by noting that the l 1 -norm regularised optimisation in (31) is a nonconvex problem due to the phase error matrix Λ . Figure 6 plots the objective function of (31) as a function of Δ R 2 assuming that α is perfectly known. We observe that this objective function has many local maxima and minima, confirming the non-convexity of the l 1 -norm regularised optimisation in (31). The reason for this non-convexity is explained in Appendix B. Due to this nonconvexity, the iterative solvers presented in Section 4 are prone to converge to local minima; thus, limiting the effectiveness of this approach.
The performance of the OMP, POMP, and l 1 -norm regularised optimisation methods are now compared using the earth mover’s distance (EMD) between the true and reconstructed range profiles. EMD [31] is a metric estimating the distance between two distributions or equivalently the minimal amount of work required to transform one distribution to the other. Figure 7 and Figure 8 show the EMD performance of the OMP, POMP, and l 1 -norm regularised optimisation methods, averaged over 100 Monte Carlo runs, versus different levels of SNR (signal-to-noise ratio) and phase error, respectively. It is observed that the POMP method yields the smallest EMD values amongst all algorithms considered. Since a smaller value of EMD corresponds to a higher level of similarity between the true and reconstructed range profiles, this observation indicates that the reconstructed range profile obtained by POMP is closer to the ground-truth range profiles than those obtained from the OMP and l 1 -norm regularised optimisation methods. This verifies the performance advantage of the POMP method from a statistical point of view.

5.2. Scenario 2: Four Sub-Bands

We now consider another scenario with four sub-bands at carrier frequencies of f 1 = 6 GHz, f 2 = 8 GHz, f 3 = 10 GHz, and f 4 = 12 GHz, each having a bandwidth of B = 300 MHz and 64 frequency steps (i.e., N = N 1 = N 2 = N 3 = N 4 = 64 ). The range errors are set to Δ R 2 = 2.78 Δ Grid , Δ R 3 = 1.33 Δ Grid , and Δ R 4 = 3.69 Δ Grid . Other simulation parameters and the true range profile of the target remain unchanged as in the previous simulation example.
Figure 9 compares the reconstructed range profiles obtained by the conventional OMP algorithm and the POMP algorithm presented in Section 3.2. OMP results in an unsatisfactorily reconstructed range profile with many spurious peaks, as expected, because it ignores the phase errors between different sub-bands. In contrast, the POMP is capable of reconstructing the true range profile with a high accuracy thanks to the use of dictionary learning with a pruning process. Note that, given the inferior performance of the l 1 -norm regularised optimisation approach compared to POMP, as demonstrated in the previous simulation scenario, this approach is excluded from the comparison here.

6. Conclusions

This paper explores the use of the POMP algorithm and l 1 -norm regularisation solvers for the problem of sparsity-driven HRRP with bandwidth stitching in the presence of phase errors. We observe that the l 1 -norm regularisation solvers do not provide significant performance improvement over the conventional sparse reconstruction algorithms due to the nonconvexity of l 1 -norm regularised optimisation when phase errors exists. In contrast, POMP is observed to be capable of effectively dealing with the phase errors and thus be able to reconstruct the range profile of the target with high accuracy. Simulation results show a significant performance improvement by POMP over OMP and the conventional and refined l 1 -norm regularisation. In future work, we propose using experimental data for a more general scenario where the true scatterers constituting the target are located in off-grid positions with respect to the dictionary grid, and the true range errors have off-grid values. We shall also consider the more general case of frequency-dependence of scatterer RCS (radar cross-section). A potential approach for this is to exploit the framework of spectral compressive sensing [32,33].

Author Contributions

Conceptualisation, P.B., N.H.N., and H.-T.T.; methodology, P.B., N.H.N., and H.-T.T.; simulation, N.H.N.; validation, P.B., N.H.N., and H.-T.T.; writing—original draft preparation, N.H.N.; writing—review and editing, P.B., N.H.N., and H.-T.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Defence Science and Technology Group, Australia.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of Gauss–Newton Algorithm for Estimation of Δ R ^ m ( i + 1 )

We rewrite e m , n as
e m , n = S m , n exp 4 π j f m , n c Δ R m Z m , n
with
Z m , n = U m , n l = 1 L α ^ l ( i + 1 ) exp 4 π j f m , n c x l .
By noting that
exp 4 π j f m , n c Δ R m Z m , n = Z m , n R cos 4 π f m , n c Δ R m + Z m , n I sin 4 π f m , n c Δ R m + j Z m , n I cos 4 π f m , n c Δ R m Z m , n R sin 4 π f m , n c Δ R m
where Z m , n R and Z m , n I are the real and imaginary components of Z m , n , we can decouple and stack the real and imaginary components of e m to form a real-valued vector as
ϵ m = e m R T , e m I T T
where
e m R = , e m , n R , n = 1 , , N m T
e m I = , e m , n I , n = 1 , , N m T
and
e m , n R = Z m , n R cos 4 π f m , n c Δ R m + Z m , n I sin 4 π f m , n c Δ R m
e m , n I = Z m , n I cos 4 π f m , n c Δ R m Z m , n R sin 4 π f m , n c Δ R m .
Using (46)–(A8), the minimization (45) becomes
Δ R ^ m ( i + 1 ) = arg min Δ R m ϵ ( Δ R m ) 2 2 .
By adopting the Gauss–Newton algorithm, Δ R ^ m ( i + 1 ) can be computed from Δ R ^ m ( i ) as
Δ R ^ m ( i + 1 ) = Δ R ^ m ( i ) J m T ( i ) J m ( i ) 1 J m T ( i ) ϵ ( Δ R ^ m ( i ) )
where ϵ ( Δ R ^ m ( i ) ) is an estimated version of ϵ computed from Δ R ^ m ( i ) and J m ( i ) is the Jacobian of ϵ ( Δ R m ) with respect to Δ R m evaluated at Δ R m = Δ R ^ m ( i ) .
The expression for the Jacobian J m of ϵ ( Δ R m ) with respect to Δ R m is given by
J m = J m R T , J m I T T
where
J m R = , J m , n R , n = 1 , , N m T
J m I = , J m , n I , n = 1 , , N m T
and
J m , n R = 4 π f m , n c Z m , n R sin 4 π f m , n c Δ R m + Z m , n I cos 4 π f m , n c Δ R m
J m , n I = 4 π f m , n c Z m , n I sin 4 π f m , n c Δ R m + Z m , n R cos 4 π f m , n c Δ R m .

Appendix B. Analysis of the Nonconvexity of the l 1 -Norm Regularised Optimisation Problem (31)

From Equation (31)
{ α ^ , Λ ^ } = arg min α , Λ S Λ F α 2 2 + μ α 1 ,
or
{ α ^ , Λ ^ } = arg min α arg min Λ S Λ F α 2 2 + μ α 1 .
S Λ F α 2 2 = S Λ F α H S Λ F α
S Λ F α 2 2 = S H S 2 e [ ( F H Λ H S ) H α ] + α H F H F α .
Only the term ( F H Λ H S ) H α is dependant on Δ R 2 , so the objective function is minimised with respect to Δ R 2 when the correlation between F H Λ H S and α is maximised. Now
F H Λ H S = F 1 H S 1 + F 2 H Λ 2 H S 2
And the terms F 1 H S 1 and F 2 H Λ 2 H S 2 represent the conventional pulse compression (i.e., transformation from frequency domain to range domain) for each of the radars with a range offset Δ R 2 . These are non-sparse and low resolution due to the oversampling in range. So F H Λ H ( Δ R 2 ) S represents the linear superposition of the two conventional range profiles. When these range profiles from the two radars are correctly aligned in range, they will better correlate with the true range profile α . Also if the range offset Δ R 2 is such that a scatterer for one radar is superimposed upon a different scatterer for the other radar, a local minimum will occur. That explains Figure 6 and the reason for its non-convexity.

Appendix C. Analysis of Operation Count for POMP

Consider first the case of M = 2 . Let the size of the grid for the parameter Δ R 2 be L Δ = 2 N Δ , so that a dictionary is constructed for each of these 2 N Δ values of Δ R 2 . OMP is implemented with the number of dictionaries halved at each stage; hence, the name “pruned” OMP. At stage k of OMP, the number of complex multiplications and divisions for a single dictionary is denoted C o m p ( k ) . Here k = 1 , , N Δ + 1 with 2 N Δ k + 1 dictionaries considered at stage k. The purpose of this section is to estimate the dependence of the computational cost of POMP on the size of the dictionary L x and the grid size L Δ for the parameter Δ R 2 .
With reference to Table 1, the significant costs for POMP are associated with the Identify and Update steps. At each stage of OMP for a given dictionary, the Identify step performs atom/residual correlations which require O ( N m L x ) complex multiplications. The Update step performs a linear least squares estimation requiring Gaussian elimination which, at stage k of OMP, has an operation count O ( k 3 ) .
Due to the halving of the number of dictionaries at each stage, the total operational count required until only one dictionary is left (although more OMP steps may be required for that dictionary until the residual is sufficiently small) is therefore of order
k = 1 N Δ + 1 ( N m L x + k 3 ) 2 N Δ k + 1
or
( 2 N Δ + 1 1 ) N m L x + ( N Δ + 1 ) 3 + 2 N Δ + 1 k = 1 N Δ k 3 z k
with z = 1 2 .
The finite sum k = 1 n k 3 z k is referred to as a low-order polylogarithm, for which a formula may be derived [34]. This formula can be shown to have a leading term of order n 3 z n + 3 so that the overall operation count for POMP is O ( ( 2 N Δ + 1 1 ) N m L x + N Δ 3 ) . As the size of the grid for Δ R 2 is L Δ = 2 N Δ , N Δ = log L Δ log 2 , and the operation count in terms of L Δ is O L Δ N m L x + log L Δ log 2 3 . We see that to leading order the computational cost is proportional to the size of the grid for Δ R 2 .
For the case of general M this cost is multiplied by ( M 1 ) .

References

  1. Ai, X.; Li, Y.; Wang, X.; Xiao, S. Some results on characteristics of bistatic high-range resolution profiles for target classification. IET Radar Sonar Navig. 2012, 6, 379–388. [Google Scholar] [CrossRef]
  2. Ai, X.; Zou, X.; Liu, J.; Li, Y.; Xiao, S. Bistatic high range resolution profiles of precessing cone-shaped targets. IET Radar Sonar Navig. 2013, 7, 615–622. [Google Scholar] [CrossRef]
  3. Liu, H.; Chen, F.; Du, L.; Bao, Z. Robust radar automatic target recognition algorithm based on HRRP signature. Front. Electr. Electron. Eng. 2012, 7, 49–55. [Google Scholar] [CrossRef]
  4. Tait, P. Introduction to Radar Target Recognition; IET: London, UK, 2005. [Google Scholar]
  5. Mao, Z.; Wei, Y. Interpulse-frequency-agile and intrapulse-phase-coded waveform optimisation for extend-range correlation sidelobe suppression. IET Radar Sonar Navig. 2017, 11, 1530–1539. [Google Scholar] [CrossRef]
  6. Bai, X.; Zhou, F.; Bao, Z. High-resolution radar imaging of space targets based on HRRP series. IEEE Trans. Geosci. Remote. Sens. 2013, 52, 2369–2381. [Google Scholar] [CrossRef]
  7. Ajorloo, A.; Hadavi, M.; Bastani, M.H.; Nayebi, M.M. Radar HRRP modeling using dynamic system for radar target recognition. Radioengineering 2014, 23, 121–127. [Google Scholar]
  8. Rihaczek, A.W. Principles of High-Resolution Radar; Artech House: Norwood, MA, USA, 1996. [Google Scholar]
  9. Wehner, D.R. High Resolution Radar; Artech House, Inc.: Norwood, MA, USA, 1987; 484p. [Google Scholar]
  10. Cuomo, K. A Bandwidth Extrapolation Technique for Improved Range Resolution Of Coherent Radar Data; Project Report CJP-60, Rev. 1; Massachusetts Inst of Tech: Lexington, MA, USA, 1992. [Google Scholar]
  11. Borison, S.; Bowling, S.B.; Cuomo, K.M. Super-resolution methods for wideband radar. Linc. Lab. J. 1992, 5, 441–461. [Google Scholar]
  12. French, A. Improved high range resolution profiling of aircraft using stepped-frequency waveforms with an S-band phased array radar. In Proceedings of the 2006 IEEE Conference on Radar, Verona, NY, USA, 24–27 April 2006; pp. 69–75. [Google Scholar]
  13. Eldar, Y.C.; Kutyniok, G. (Eds.) Compressed Sensing: Theory and Applications; Cambridge University Press: New York, NY, USA, 2012. [Google Scholar]
  14. Amin, M. (Ed.) Compressive Sensing for Urban Radar; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  15. Tropp, J.A.; Wright, S.J. Computational methods for sparse solution of linear inverse problems. Proc. IEEE 2010, 98, 948–958. [Google Scholar] [CrossRef] [Green Version]
  16. Potter, L.C.; Ertin, E.; Parker, J.T.; Cetin, M. Sparsity and compressed sensing in radar imaging. Proc. IEEE 2010, 98, 1006–1020. [Google Scholar] [CrossRef]
  17. Tomei, S.; Bacci, A.; Giusti, E.; Martorella, M.; Berizzi, F. Compressive sensing-based inverse synthetic radar imaging imaging from incomplete data. IET Radar Sonar Navig. 2016, 10, 386–397. [Google Scholar] [CrossRef]
  18. Nguyen, N.H.; Dogancay, K.; Tran, H.T.; Berry, P. An image focusing method for sparsity-driven radar imaging of rotating targets. Sensors 2018, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Qiu, W.; Giusti, E.; Bacci, A.; Martorella, M.; Berizzi, F.; Zhao, H.; Fu, Q. Compressive sensing-based algorithm for passive bistatic ISAR with DVB-T signals. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 2166–2180. [Google Scholar] [CrossRef]
  20. Nguyen, N.H.; Dogancay, K.; Tran, H.; Berry, P. Nonlinear least-squares post-processing for compressive radar imaging of a rotating target. In Proceedings of the International Conference on Radar (RADAR), Brisbane, Australia, 27–30 August 2018; pp. 1–6. [Google Scholar]
  21. Winkler, V.; Edrich, M. Multiband radar signal processing. In Proceedings of the 2015 16th International Radar Symposium (IRS), Dresden, Germany, 24–26 June 2015; pp. 108–113. [Google Scholar]
  22. Giusti, E.; Tomei, S.; Bacci, A.; Martorella, M.; Berizzi, F. Autofocus for CS based ISAR imaging in the presence of gapped data. In Proceedings of the 2nd International Workshop on Compressed Sensing Applied to Radar (CoSeRa), Bonn, Germany, 17–19 September 2013. [Google Scholar]
  23. Cuomo, K.M.; Piou, J.E.; Mayhan, J.T. Ultra-wideband coherent processing. Linc. Lab. J. 1997, 10, 203–222. [Google Scholar]
  24. Li, G.; Varshney, P.K. Micro-Doppler parameter estimation via parametric sparse representation and pruned orthogonal matching pursuit. IEEE J. Sel. Top. Appl. Earth Obs. Remote. Sens. 2014, 7, 4937–4948. [Google Scholar] [CrossRef]
  25. Onhon, N.O.; Cetin, M. A sparsity-driven approach for joint SAR imaging and phase error correction. IEEE Trans. Image Process. 2012, 21, 2075–2088. [Google Scholar] [CrossRef]
  26. Cetin, M.; Karl, W.C. Feature-enhanced synthetic aperture radar image formation based on nonquadratic regularization. IEEE Trans. Image Process. 2001, 10, 623–631. [Google Scholar] [CrossRef] [Green Version]
  27. Zhang, L.; Qiao, Z.J.; Xing, M.; Li, Y.; Bao, Z. High-resolution ISAR imaging with sparse stepped-frequency waveforms. IEEE Trans. Geosci. Remote. Sens. 2011, 49, 4630–4651. [Google Scholar] [CrossRef]
  28. Vogel, C.R.; Oman, M.E. Fast, robust total variation-based reconstruction of noisy, blurred images. IEEE Trans. Image Process. 1998, 7, 813–824. [Google Scholar] [CrossRef] [Green Version]
  29. Acar, R.; Vogel, C.R. Analysis of bounded variation penalty methods for ill-posed problems. Inverse Probl. 1994, 10, 1217–1229. [Google Scholar] [CrossRef]
  30. Liao, Y. Phase and Frequency Estimation: High-Accuracy and Low-Complexity Techniques. Master’s Thesis, Worcester Polytechnic Institute, Worcester, MA, USA, 2011. [Google Scholar]
  31. Rubner, Y.; Tomasi, C.; Guibas, L.J. The Earth Mover’s Distance as a Metric for Image Retrieval. Int. J. Comput. Vis. 2000, 40, 99–121. [Google Scholar] [CrossRef]
  32. Duarte, M.F.; Baraniuk, R.G. Spectral compressive sensing. Appl. Comput. Harmon. Anal. 2013, 35, 111–129. [Google Scholar] [CrossRef]
  33. Liu, S.; Zhang, Y.D.; Shan, T.; Tao, R. Structure-aware Bayesian compressive sensing for frequency-hopping spectrum estimation with missing observations. IEEE Trans. Signal Process. 2018, 66, 2153–2166. [Google Scholar] [CrossRef]
  34. Olver, F.W.J.; Olde Daalhuis, A.B.; Lozier, D.W.; Schneider, B.I.; Boisvert, R.F.; Clark, C.W.; Miller, B.R.; Saunders, B.V.; Cohl, H.S.; McClain, M.A. (Eds.) NIST Digital Library of Mathematical Functions; Release 1.0.25 of 2019-12-15; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2019. Available online: http://dlmf.nist.gov/ (accessed on 23 January 2020).
Figure 1. True range profile of synthetic target under consideration.
Figure 1. True range profile of synthetic target under consideration.
Sensors 20 00665 g001
Figure 2. Performance of conventional OMP.
Figure 2. Performance of conventional OMP.
Sensors 20 00665 g002
Figure 3. Performance of conventional l 1 -norm regularised optimisation solver.
Figure 3. Performance of conventional l 1 -norm regularised optimisation solver.
Sensors 20 00665 g003
Figure 4. Performance of POMP in the presence of phase errors.
Figure 4. Performance of POMP in the presence of phase errors.
Sensors 20 00665 g004
Figure 5. Performance of l 1 -norm regularised optimisation solver with phase error correction.
Figure 5. Performance of l 1 -norm regularised optimisation solver with phase error correction.
Sensors 20 00665 g005
Figure 6. Illustration of the nonconvexity of the l 1 -norm regularised optimisation problem (31). The objective function of (31) is plotted against Δ R 2 assuming that α is perfectly known.
Figure 6. Illustration of the nonconvexity of the l 1 -norm regularised optimisation problem (31). The objective function of (31) is plotted against Δ R 2 assuming that α is perfectly known.
Sensors 20 00665 g006
Figure 7. Earth mover’s distance (EMD) performance of the OMP, POMP, and l 1 -norm regularised optimisation methods versus various of SNRs ( Δ R 2 = 2 . 78 Δ Grid ).
Figure 7. Earth mover’s distance (EMD) performance of the OMP, POMP, and l 1 -norm regularised optimisation methods versus various of SNRs ( Δ R 2 = 2 . 78 Δ Grid ).
Sensors 20 00665 g007
Figure 8. EMD performance of the OMP, POMP, and l 1 -norm regularised optimisation methods versus various levels of phase error ( SNR = 10 dB).
Figure 8. EMD performance of the OMP, POMP, and l 1 -norm regularised optimisation methods versus various levels of phase error ( SNR = 10 dB).
Sensors 20 00665 g008
Figure 9. Performance comparison between OMP and POMP for Simulation Scenario 2 (with four sub-bands).
Figure 9. Performance comparison between OMP and POMP for Simulation Scenario 2 (with four sub-bands).
Sensors 20 00665 g009
Table 1. The pruned orthogonal matching pursuit (POMP) algorithm (M = 2).
Table 1. The pruned orthogonal matching pursuit (POMP) algorithm (M = 2).
INPUT:
  • Noisy signal data vector S ˜
  • Candidate dictionaries Φ 1 , Φ 2 , , Φ L Δ , corresponding to L Δ candidate values of Δ R 2
PROCEDURE:
  • Initialization:
    -
    set the initial indexes of active dictionaries to Θ 1 = { 1 , , L Δ } ;
    -
    set the corresponding residual vectors to r 1 = = r L Δ = S ˜ ;
    -
    set the initial support Λ to ∅, the null set;
  • for i = 1 ; i : = i + 1 until | Θ i | = = 1 (the cardinality of Θ i ) and | r l | < ϵ for l Θ i ,
    • for every l Θ i , perform OMP as follows
      -
      Identify:
      • c l = Φ l H r l
      • j l = arg max j | c j |
      -
      Merge supports:
      • Λ l = Λ l j l
      -
      Update * :
      • α ^ l , Λ l = Φ l , Λ l H Φ l , Λ l 1 Φ l , Λ l H S ˜
      • r l = S ˜ Φ l , Λ l α ^ l , Λ l
    • end for
    • if | Θ i | > 1
      • Remove indices of | Θ i | / 2 candidate dictionaries that correspond to | Θ i | / 2 largest residual errors from | Θ i |
    • end if
  end for
OUTPUT:
  • The range profile estimate α ^ l , Λ l where l is the last remaining element of Θ i
  • The estimate of Δ R 2 is the value of Δ R 2 corresponding to Φ l
* Φ l , Λ l consists of the columns of Φ l with indices belonging to Λ l and α ^ Λ l consists of the elements of α ^ l with indices belonging to Λ l .

Share and Cite

MDPI and ACS Style

Berry, P.; Nguyen, N.H.; Tran, H.-T. Compressive Sensing-Based Bandwidth Stitching for Multichannel Microwave Radars. Sensors 2020, 20, 665. https://doi.org/10.3390/s20030665

AMA Style

Berry P, Nguyen NH, Tran H-T. Compressive Sensing-Based Bandwidth Stitching for Multichannel Microwave Radars. Sensors. 2020; 20(3):665. https://doi.org/10.3390/s20030665

Chicago/Turabian Style

Berry, Paul, Ngoc Hung Nguyen, and Hai-Tan Tran. 2020. "Compressive Sensing-Based Bandwidth Stitching for Multichannel Microwave Radars" Sensors 20, no. 3: 665. https://doi.org/10.3390/s20030665

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