Next Article in Journal
A New Missing Data Imputation Algorithm Applied to Electrical Data Loggers
Previous Article in Journal
Theory and Application of Magnetic Flux Leakage Pipeline Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compressive Sensing Based Design of Sparse Tripole Arrays

1
Department of Automatic Control and Systems Engineering, University of Sheffield, Sheffield S1 3JD, UK
2
Department of Electronic and Electrical Engineering, University of Sheffield, Sheffield S1 3JD, UK
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(12), 31056-31068; https://doi.org/10.3390/s151229849
Submission received: 17 November 2015 / Revised: 3 December 2015 / Accepted: 7 December 2015 / Published: 10 December 2015
(This article belongs to the Section Physical Sensors)

Abstract

:
This paper considers the problem of designing sparse linear tripole arrays. In such arrays at each antenna location there are three orthogonal dipoles, allowing full measurement of both the horizontal and vertical components of the received waveform. We formulate this problem from the viewpoint of Compressive Sensing (CS). However, unlike for isotropic array elements (single antenna), we now have three complex valued weight coefficients associated with each potential location (due to the three dipoles), which have to be simultaneously minimised. If this is not done, we may only set the weight coefficients of individual dipoles to be zero valued, rather than complete tripoles, meaning some dipoles may remain at each location. Therefore, the contributions of this paper are to formulate the design of sparse tripole arrays as an optimisation problem, and then we obtain a solution based on the minimisation of a modified l 1 norm or a series of iteratively solved reweighted minimisations, which ensure a truly sparse solution. Design examples are provided to verify the effectiveness of the proposed methods and show that a good approximation of a reference pattern can be achieved using fewer tripoles than a Uniform Linear Array (ULA) of equivalent length.

1. Introduction

Sensor arrays are used in a wide range of application areas including Direction Of Arrival (DOA) estimation and beamforming [1], vehicle health maintenance [2] and many others. When considering fixed beamformer design, it is well known that an adjacent antenna separation of less than half a wavelength is required if using Uniform Linear Arrays (ULAs) in order to avoid the introduction of grating lobes [1]. This can become problematic when considering arrays with large apertures due to the cost associated with the number of antennas required. As a result, sparse arrays have become a desirable alternative as they allow a given aperture to be implemented using fewer antennas. The reason for this is that their nonuniform adjacent antenna separations allow mean separations greater than half of the wavelength, while still avoiding the introduction of grating lobes [3].
The tradeoff is that the unpredictable sidelobe behaviour of sparse arrays means that some optimisation of antenna locations is required in order to ensure an acceptable performance is achieved. For example, we could consider minimising the Peak Sidelobe Level (PSL) or matching the response to an acceptable reference response. Stochastic optimisation methods such as Genetic Algorithms (GAs) [4,5,6,7] and Simulated Annealing (SA) [8,9] have been shown to be able to solve such a problem. Improved performances can also be achieved by combining GAs with Difference Sets (DS) [10] and Almost Difference Sets (ADS) [11]. However, the disadvantages of GAs and similar methods are the potentially long computation times and the possibility of failing to converge to the optimal solution.
Instead, Compressive Sensing (CS) [12] (or alternatively the Bayesian formulation [13]) tells us that, when certain conditions are met, we can recover signals from fewer measurements than are traditionally required. CS-based methods have been applied to various array signal processing problems such as DOA estimation [14] and in the design of sparse antenna arrays [15,16,17,18,19]. For CS-based sparse array design methods, the aim is to use as few antennas as possible, while still getting an exact, or almost exact, match to a reference response. It has further been shown that the sparsity of the solution can be improved by converting the problem into a series of iteratively solved reweighted minimisations [20,21]. This is achieved by making the problem a closer approximation of an l 0 norm minimisation by adding a reweighting term that penalises small weight coefficients more heavily [22].
Previous work involving CS and sparse array design has tended to focus on the case of isotropic array elements. This means that the polarisation of the signal of interest is not accounted for, potentially degrading the performance of an array. Alternatively, we can consider arrays consisting of vector sensors allowing measurement of different components of the received electromagnetic waveform [23,24,25,26]. One such vector sensor is a tripole, which consists of three orthogonally orientated dipoles [27,28,29,30,31,32,33,34,35]. Such vector sensors have been used in various configurations for DOA and polarisation estimation [36,37,38,39,40,41,42,43]. However, in these works, fixed sensor locations are considered rather than addressing the location optimisation problem associated with sparse array design.
The novelty of this paper arises in the proposed solution to the problem of designing sparse tripole arrays using CS for the first time. Unlike for isotropic array elements (consisting of array elements of a single antenna), we now have three complex valued weight coefficients (one for each co-located dipole) that have to be simultaneously minimised for each potential tripole location. If this is not done, there is the potential for only the coefficients of single dipoles (rather than complete tripoles) to be made zero-valued. As a result, we can obtain a solution with at least one dipole at each location on the grid of potential sensor locations rather than obtaining a truly sparse solution. Instead, to ensure the desired sparsity is achieved, we propose two CS-based methods of solving this problem: (i) the optimisation can be formulated as a modified l 1 norm minimisation; (ii) the sparsity can be further improved by converting the problem into a series of iteratively solved reweighted minimisations. We provide broadside and off-broadside design examples to show that the proposed methods can successfully obtain a good approximation of a reference pattern using fewer tripoles than a ULA of equivalent length.
This work differs from previous work of the first and third authors as follows: (i) In [44,45,46], we considered the problem of designing sparse sensor arrays for both narrowband and wideband beamforming. Unlike the current work, the arrays consisted of isotropic array elements with a single sensor at each location. In [47], we dealt with the design of sparse linear arrays based on crossed dipoles. However, [47] deals with a quaternionic signal model, which is not utilised in this work; (ii) [48] deals with DOA estimation using a fixed spatially stretched tripole (no co-located tripole) rather than the sparse array design problem.
The remainder of this paper is structured as follows, Section 2 gives details of the proposed design methods including: the array model in Section 2.1 and the CS based formulations in Section 2.2. Design examples are presented in Section 3, with conclusions being drawn in Section 4.

2. Proposed Design Methods

2.1. Array Model

Figure 1 shows the tripole array structure being considered. There are M tripole locations spread over the y-axis with an adjacent separation distance d. For each tripole, there are three orthogonally orientated dipoles, one parallel to each axis. Also shown is a signal with its DOA defined by the angles θ and ϕ, which are limited as follows: 0 θ π / 2 and - π / 2 ϕ π / 2 . A plane-wave signal model is assumed, i.e., the signal impinges upon the array from the far field.
Figure 1. An array model with M potential tripole locations shown.
Figure 1. An array model with M potential tripole locations shown.
Sensors 15 29849 g001
The spatial steering vector of the array is given by
s s ( θ , ϕ ) = [ 1 , e - j 2 π d sin θ sin ϕ / λ , , e - j 2 π ( M - 1 ) d sin θ sin ϕ / λ ]
where λ is the wavelength of the signal of interest. For tripoles, the spatial-polarization coherent vector contains information about a signal’s polarisation and is given by [27]:
s p ( θ , ϕ , γ , η ) = sin γ cos θ cos ϕ e j η - cos γ sin ϕ sin γ cos θ sin ϕ e j η - cos γ cos ϕ - sin γ sin θ e j η = s p , x ( θ , ϕ , γ , η ) s p , y ( θ , ϕ , γ , η ) s p , z ( θ , ϕ , γ , η )
where 0 γ π / 2 is the auxiliary polarization angle and - π η < π is the polarization phase difference.
Now, the array can be split into three sub-arrays, one parallel to each axis. The steering vector of each of these sub-arrays is complex-valued and given by:
s x ( θ , ϕ , γ , η ) = s p , x ( θ , ϕ , γ , η ) s s ( θ , ϕ )
s y ( θ , ϕ , γ , η ) = s p , y ( θ , ϕ , γ , η ) s s ( θ , ϕ )
and
s z ( θ , ϕ , γ , η ) = s p , z ( θ , ϕ , γ , η ) s s ( θ , ϕ )
respectively.
The response of the array is given by
p ( θ , ϕ , γ , η ) = s ( θ , ϕ , γ , η ) w H
with
w = [ w x , 1 , w y , 1 , w z , 1 , , w x , M , w y , M , w z , M ]
where w i , m is the weight coefficient for the m t h dipole orientated parallel to the i-axis ( i = { x , y , z } ) and { . } H denotes the Hermitian transpose. Similarly, we also have
s ( θ , ϕ , γ , η ) = [ s x , 1 ( θ , ϕ , γ , η ) , s y , 1 ( θ , ϕ , γ , η ) , s z , 1 ( θ , ϕ , γ , η ) , , s x , M ( θ , ϕ , γ , η ) , s y , M ( θ , ϕ , γ , η ) , s z , M ( θ , ϕ , γ , η ) ]
where s i , m ( θ , ϕ , γ , η ) is the contribution of the m t h dipole to the overall steering vector parallel to the i-axis.
In what follows, we address the problem of designing a sparse linear tripole array, where there are two issues to be solved. Firstly, we have to find the optimal tripole locations, as there will no longer be a uniform adjacent tripole separation of d. Secondly, we have to find the weight coefficients that give an acceptable beam response. We propose to use design methods based on CS to solve these problems.

2.2. Compressive Sensing Based Design of Sparse Tripole Arrays

Suppose P r ( θ , ϕ , γ , η ) is a reference response which we wish to achieve. First, consider Figure 1 as being a grid of potential tripole locations. In this instance, ( M - 1 ) d is the potential aperture of the array and M is a large number. Sparseness is then introduced by selecting the weight coefficients to give as few active tripoles as possible, or, in other words, as few non-zero valued weight coefficients as possible. This has to be done while still giving a designed response that is close to the desired one.
This problem is formulated as
min | | w | | 1 subject to | | p r - S w H | | 2 α
where | | w | | 1 is the l 1 norm of the weight coefficients and used as an estimate of the number of nonzero weight coefficients in w , i.e., an approximation to the l 0 norm, p r is the vector holding the desired beam response at the sampled angular and polarisation points of interest, S is the matrix composed of the corresponding steering vectors and α places a limit on the allowed difference between the desired and the designed responses. Therefore, by minimising | | w | | 1 in Equation (9), we are minimising the number of tripoles we have to use, while the constraint ensures that we still get a good approximation of the reference response.
In detail, p r and S are respectively given by
p r = [ P r ( θ 1 , ϕ 1 , γ 1 , η 1 ) , , P r ( θ L , ϕ L , γ L , η L ) ] T
S = [ s ( θ 1 , ϕ 1 , γ 1 , η 1 ) , , s ( θ L , ϕ L , γ L , η L ) ] T
where L is the number of points sampled at each dimension of the desired beam response and { . } T denotes the transpose operation. In this work, p r is the ideal response, i.e., a value of one for the mainlobe and zeros for the other entries. It would also be possible to use other reference responses if desired, for example, a response which defines a maximum sidelobe level or one with nulls at a given direction.
Note, we have to choose L to be large enough to ensure that all angular and polarisation points of interest are being considered. However, the larger L is, the longer the computation time will be. An acceptable compromise seems to be to sample the sidelobe regions every 1 ° . We can also help improve the sparsity of the solution by increasing the value of M. A larger value of M gives a denser grid which is more likely to include the optimal locations. As a result, the same amount of error between designed and reference responses can be obtained using fewer tripoles. However, we will get to the point where the optimal locations will be already included on the grid and further increases will only serve to lengthen the computation time. This means there is again a compromise to be achieved when selecting the value of M, and this will be explored in the design examples we present in this paper.
One problem with the above formulation is that the three complex weight coefficients associated with each tripole can not be guaranteed to be minimised simultaneously. As a result, we may not be able to obtain a sparse solution even if we have a sparse weight vector w after the minimisation. In other words, it would be possible to end up only removing a single dipole from each potential location rather than truly introducing sparsity by removing complete tripoles. Therefore, we further modify the formulation in Equation (9) by converting it into a modified l 1 norm minimisation.
First, we rewrite Equation (9) as
min t ϵ R + subject to | | p r - S w H | | 2 α , | w | 1 t
where
| w | 1 = m = 1 M | | w m | | 2
and w m = [ w x , m , w y , m , w z , m ] . Here w m contains the weight coefficients for all the dipoles that make up the tripole at the m t h grid location. Therefore, by minimising the value of t the second constraint in Equation (12) ensures the combined weight coefficients for the tripoles are also minimised.
Now, we decompose t to t = m = 1 M t m , t m ϵ R + . This allows us to minimise the combined weight coefficients for each tripole separately. In vector form, we have
t = [ 1 , , 1 ] t 1 t M = 1 T t .
Then, Equation (12) can be rewritten as
min t 1 T t subject to | | p r - S w H | | 2 α | | w m | | 2 t m , m = 1 , , M .
Now, define
w ^ = [ t 1 , R ( w x , 1 ) , R ( w y , 1 ) , R ( w z , 1 ) , - I ( w x , 1 ) , - I ( w y , 1 ) , - I ( w z , 1 ) , t 2 , , - I ( w z , M ) ] ,
c ^ = [ 1 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , , 1 , 0 , 0 , 0 , 0 , 0 , 0 ] T ,
p ^ r = [ R ( p r ) , I ( p r ) ] T
and
S ^ = 0 0 R ( s x , 1 ) I ( s x , 1 ) R ( s y , 1 ) I ( s y , 1 ) R ( s z , 1 ) I ( s z , 1 ) - I ( s x , 1 ) R ( s x , 1 ) - I ( s y , 1 ) R ( s y , 1 ) - I ( s z , 1 ) R ( s z , 1 ) 0 0 - I ( s z , M ) R ( s z , M ) T
where s i , m is the contribution of the m t h dipole parallel to the i - axis to the steering vector of the array, R ( . ) gives the real component of a scalar and I ( . ) gives the imaginary component of a scalar and R ( . ) gives the real components of a vector/matrix and I ( . ) gives the imaginary components of a vector/matrix. Here, we have included the values of t m in w ^ as we want to automatically find them with the weight coefficients instead of them being predetermined. We then use c ^ to select the values of t m for minimisation, and the zeros are added in S ^ to ensure they do not contribute to the error between the reference and designed responses.
Then, the final formulation is as follows
min w ^ w ^ c ^ subject to | | p ^ r - S ^ w ^ T | | 2 α | | w m | | 2 t m , m = 1 , , M
which can be solved using cvx, a package for specifying and solving convex programs [49,50].
As with the traditional CS based design methods, the solution can be improved by considering the problem as a series of iteratively solved reweighted minimisations to bring the solution closer to that of the l 0 norm minimisation [20,21,22]. This is achieved by the addition of a reweighting term which is found from the previous iteration and penalises small non-zero valued weight coefficients more heavily, meaning all non-zero valued coefficients are treated in a more uniform manner (as for the l 0 norm minimisation). As a result, these small non-zero valued weight coefficients are less likely to be repeated in the next iteration and the sparsity of the solution is improved.
In this form, the problem is formulated as follows
min w ^ w ^ c ^ subject to | | p ^ r - S ^ w ^ T | | 2 α δ m k | | w m | | 2 t m , m = 1 , , M
where we now have
c ^ = [ δ 1 k , 0 , 0 , 0 , 0 , 0 , 0 , δ 2 k , 0 , , δ M k , 0 , 0 , 0 , 0 , 0 , 0 ] T
and
δ m k = ( | w m k - 1 | + ϵ ) - 1 .
This is then solved iteratively until the number of non-zero valued weight coefficients has remained constant for a few iterations.
Note, the term ϵ is added for numerical stability and should be set slightly smaller than the minimum implemented weight coefficient for a single tripole. However, its addition means that a zero-valued weight coefficient in one iteration will not be guaranteed to be repeated in the next. This has the result of meaning a solution is not always guaranteed, but when a solution is possible, it will usually be achieved in less than 10 iterations. For the first iteration, Equation (20) can be solved and the reweighting terms introduced from the second iteration onwards. Finally, w m gives the current estimate of the weight coefficients for the m t h tripole, whereas w m k - 1 gives the estimate from the previous iteration.

3. Design Examples

We will now consider broadside and off-broadside design examples to verify the effectiveness of the two design methods. All examples are implemented on a computer with an Intel Xeon CPU E3-1271 (3.60 GHz), 16GB of RAM and Windows 7 operating system, using Matlab R2014a and cvx [49,50]. In all of the figures that follow positive values of θ indicate the value range θ [ 0 ° , 90 ° ] with ϕ = 90 ° , while negative values of θ [ - 90 ° , 0 ° ] indicate an equivalent range of θ [ 0 ° , 90 ° ] with ϕ = - 90 ° .

3.1. Broadside Design Examples

First, we will consider broadside ( θ M L = 0 ° , ϕ M L = 90 ° ) design examples for both methods and show how the value of M effects their performances. The length of the array is set to be 10 λ and the values α = 0 . 5 , γ = 55 ° and η = 100 ° selected. Finally, the sidelobe regions are defined as θ S L = [ 10 ° , 90 ° ] for ϕ S L = ± 90 ° .
Table 1 compares the performances of the proposed methods for values of M = 101 , 301 and 501. For each case, the results for the non-reweighted example are found using Equation (20) and for the reweighted example Equation (21) is used. Here, we are comparing the aperture of the designed array, number of tripoles required, mean adjacent tripole separation ( Δ d ¯ ), the minimum adjacent tripole separation, the computation time, designed mainlobe location ( θ ^ M L ) and PSL. The first three give an indication of the level of sparsity introduced by the design methods, the computation time is used to indicate the implementation complexity and the final two measures give an indication of how desirable the final response is. Note, we are not using | | p r - S w H | | 2 as a performance measure as the same value of α is used in each case, meaning the value of | | p r - S w H | | 2 will be the same for both methods.
Table 1. Broadside performance comparison.
Table 1. Broadside performance comparison.
ExampleAperture / λ Number of Tripoles Δ d ¯ / λ Min( Δ d )Computation Time (s) θ ^ M L PSL (dB)
Non-reweighted
M= 10110160.670.155.61 90 ° −32.28
Reweighted
M= 1015.2080.740.6025.92 90 ° −23.25
Non-reweighted
M= 30110140.770.5723.86 90 ° −32.53
Reweighted
M= 3015.1380.730.6088.88 90 ° −22.94
Non-reweighted
M= 50110140.770.5849.90 90 ° −32.69
Reweighted
M= 5015.1280.730.60168.52 90 ° −22.89
The first thing to note is that, for all values of M, both methods have successfully obtained a good approximation of the reference pattern and achieved a desirable beam response. This is evident from the fact that the mainlobes are always in the correct location with sufficient sidelobe attenuation simultaneously achieved.
When looking at the results for the non-reweighted examples, we can see that it is possible, but not guaranteed, that increasing the value of M can help increase the sparsity. This can be explained by considering that the denser the grid of potential tripole locations, the more likely it is that the grid will include the optimal tripole locations. As a result, it would be possible to achieve the same error between reference and designed responses using fewer tripoles. However, we will get to the point where the grid is dense enough to include the optimal locations and further increasing the value of M has little or no effect. This is evident from the observation that increasing M from 101 to 301 has increased the mean and minimum adjacent separations as well as decreasing the number of tripoles required. In other words, the sparsity has been increased. However, when increasing M further to 501 only the minimum adjacent sensor separation has increased, showing less of an improvement has been achieved. It is also worth noting that increasing M has been shown to always increase the computation time. If we take the computation time for M = 101 as a reference, we can see that there is a 425 % increase when M = 301 and a 978 % increase in computation time when M = 501 . Similarly, for the reweighted examples, we have increases of 343 % for M = 301 and 650 % for M = 501 . As a result, a balance has to be achieved rather than continually increasing M. In this case, M = 301 seems to be a good choice.
From the results for M = 101 , we can see that using the reweighted method has improved all measures of sparsity of the result (i.e., number of tripoles, mean and minimum adjacent tripole separation). For M = 301 and M = 501 we can see that the reweighted method still gives a significant improvement in the number of tripoles that are required. As a result, a greater cost saving is made. However, the iterative nature of this method means that in each case there has been an increase in the computation time, as would be expected. For M = 101 we have a 462 % increase, 373 % for M = 301 and 338 % for M = 501 .
As an example, the tripole locations for each method when M = 301 are shown in Table 2 and Table 3, respectively, with Figure 2 showing the designed responses. For completeness a ULA with adjacent tripole separation of λ / 2 and aperture of 10 λ is also shown for a comparison. The weight coefficients for this comparison ULA can be found using the following:
min | | p ^ r - S ^ U L A w ^ U L A T | | 2 subject to S ^ U L A , θ M L w ^ U L A T = 1
where S U L S is given by
S ^ U L A = R ( s U L A , x , 1 ) I ( s U L A , x , 1 ) R ( s U L A , y , 1 ) I ( s U L A , y , 1 ) R ( s U L A , z , 1 ) I ( s U L A , z , 1 ) - I ( s U L A , x , 1 ) R ( s U L A , x , 1 ) - I ( s U L A , y , 1 ) R ( s U L A , y , 1 ) - I ( s U L A , z , 1 ) R ( s U L A , z , 1 ) - I ( s U L A , z , M ) R ( s U L A , z , M ) T
s U L A , i , m is the contribution of the m t h dipole parallel to the i-axis, S ^ U L A , θ M L is S ^ U L A just for the mainlobe,
w ^ U L A = [ R ( w U L A , x , 1 ) , R ( w U L A , y , 1 ) , R ( w U L A , z , 1 ) , - I ( w U L A , x , 1 ) , - I ( w U L A , y , 1 ) , - I ( w U L A , z , 1 ) , , - I ( w U L A , z , M ) ]
and w U L A , i , m is the weight coefficient for the m t h dipole parallel to the i-axis. We can view this as finding the weight coefficients w ^ U L A that minimise the difference between the reference and designed responses, subject to keeping a unity response for the mainlobe direction. In practice, this has the effect of ensuring a distortionless response at the mainlobe while minimising the sidelobes.
Table 2. Locations for the non-reweighted broadside design example, M = 301 .
Table 2. Locations for the non-reweighted broadside design example, M = 301 .
n1234567891011121314
d n / λ 00.801.622.423.254.024.725.285.986.757.588.389.2010
Table 3. Locations for the reweighted broadside design example, M = 301 .
Table 3. Locations for the reweighted broadside design example, M = 301 .
n12345678
d n / λ 2.433.234.034.705.305.976.777.57
Figure 2. Responses for broadside design example, M = 301 .
Figure 2. Responses for broadside design example, M = 301 .
Sensors 15 29849 g002
Note, Figure 2 shows that the responses for the arrays designed using the two proposed methods are slightly less desirable than the comparison ULA. This is to be expected, as significantly less tripoles are being used (21 for the ULA compared to 14 or 8 for the proposed methods). However, they still give acceptable approximations of the reference response. Here, the important thing is the reduction in the number of tripoles as was the aim of the proposed methods.

3.2. Off-Broadside Design Examples

We will now present off-broadside design examples for both of the proposed methods. Here, we will use the parameters from the broadside design example but with the value of M = 301 fixed. The performance will be considered over a range of desired mainlobe locations ( θ M L varies for a fixed ϕ M L = 90 ° ) and, in each case, there will be a transition region of 10 ° on either side of the mainlobe with the remaining regions being the sidelobes.
Table 4 summarises the performance of the proposed methods over the range of mainlobes being considered. We can see that an acceptable approximation of the reference response is achieved in all cases (i.e., mainlobe is within 1 ° of what is designed and sufficient sidelobe attenuation is achieved). In addition, for each example, using the reweighted method has reduced the number of tripoles required, but at the expense of an increased computation time ( 321 % for θ = 10 ° , 295 % for θ = 30 ° , 582 % for θ = 50 ° and 544 % for θ = 70 ° ). In the case of the θ M L = 50 ° and θ M L = 70 ° examples, this reduction in tripoles has even made the difference between whether there are any savings in number of tripoles as compared to an equivalent length ULA or not. However, for the non-reweighted method, we can also get a result with improved sparsity by increasing the value of α, although this will mean that the designed response will be a worse approximation of the reference response as a result. For example, see Table 5, where a value of α = 0 . 95 is shown to give a reduction in the number tripoles required when using the non-reweighted method.
Table 4. Off-broadside performance comparison, for α = 0 . 5 .
Table 4. Off-broadside performance comparison, for α = 0 . 5 .
θ M L Aperture / λ Number of Tripoles Δ d ¯ / λ Min( Δ d )Computation Time (s) θ ^ M L PSL (dB)
Non-reweighted 10 ° 10150.710.6822.17 10 ° −35.18
Reweighted 10 ° 5.7390.720.7071.17 10 ° −21.27
Non-reweighted 30 ° 10180.590.5523.40 30 ° −30.92
Reweighted 30 ° 6.53120.590.5769.05 29 ° −23.27
Non-reweighted 50 ° 10230.450.3524.88 50 ° −24.69
Reweighted 50 ° 8.2061.641.57144.89 49 ° −18.12
Non-reweighted 70 ° 10240.430.1721.18 69 ° −19.61
Reweighted 70 ° 1081.430.17115.22 70 ° −17.51
Table 5. Off-broadside performance comparison, for α = 0 . 95 .
Table 5. Off-broadside performance comparison, for α = 0 . 95 .
θ M L Aperture / λ Number of Tripoles Δ d ¯ / λ Min( Δ d )Computation Time (s) θ ^ M L PSL (dB)
Non-reweighted 50 ° 10180.590.4825.25 50 ° −15.23
Non-reweighted 70 ° 10200.530.3325.44 70 ° −7.59
Here, we show the responses for θ M L = 10 ° as an example in Figure 3. Again, the response for a ULA of aperture 10 λ and adjacent tripole separation of λ / 2 is shown for a comparison. As for the broadside example with M = 301 , we can see that the designed responses are again slightly less desirable (wider mainlobe and increased PSL). However, this is an acceptable cost for achieving the desired reduction in the number of tripoles used. For completeness, the resulting tripole locations are shown in Table 6 and Table 7, respectively.
Figure 3. Responses for off-broadside design example, θ M L = 10 ° .
Figure 3. Responses for off-broadside design example, θ M L = 10 ° .
Sensors 15 29849 g003
Table 6. Locations for the non-reweighted off-broadside design example, θ M = 10 ° .
Table 6. Locations for the non-reweighted off-broadside design example, θ M = 10 ° .
n123456789101112131415
d n / λ 00.681.402.122.853.554.2855.726.457.157.888.609.3210
Table 7. Locations for the reweighted off-broadside design example, θ M = 10 ° .
Table 7. Locations for the reweighted off-broadside design example, θ M = 10 ° .
n123456789
d n / λ 2.132.833.574.2755.736.437.177.87

4. Conclusions

This paper has addressed the problem of designing sparse linear tripole arrays for the first time. Unlike for isotropic array elements, each antenna now has three complex valued weight coefficients (one for each co-located dipole) which have to be simultaneously minimised in order to guarantee a truly sparse solution. Therefore, we have proposed two CS-based solutions to this problem: (i) a modified l 1 norm minimisation; (ii) a series of iteratively solved reweighted minimisations, which give an improved level of sparsity. Both broadside and off-broadside design examples have been provided to verify the effectiveness of the proposed methods. In both cases, the results show that a good approximation of the reference pattern can be achieved using fewer tripoles than a ULA of equivalent length. Note, in this paper, we have focused on the case of a single polarisation of interest, as we were more interested in showing it was possible to design sparse tripole arrays using CS rather than dealing with multiple polarisations. In order to handle multiple polarisations, it may be necessary to extend the proposed methods to designing planar arrays, which is seen as an area for future research.

Acknowledgments

We appreciate the support of the UK Engineering and Physical Sciences Research Council (EPSRC) via the project Bayesian Tracking and Reasoning over Time (BTaRoT) grant EP/K021516/1.

Author Contributions

Matthew Hawes developed the proposed design methods and wrote the manuscript. Wei Liu formulated the original problem being considered. Wei Liu and Lyudmila Mihaylova were involved in the theoretical development of the solutions, in editing and structuring the manuscript. All authors were involved in discussing potential solutions, the results and determining the best test scenarios to be used.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van Trees, H.L. Optimum Array Processing, Part IV of Detection, Estimation, and Modulation Theory; John Wiley & Sons, Inc.: New York, NY, USA, 2002. [Google Scholar]
  2. Rodger, J.A. Toward reducing failure risk in an integrated vehicle health maintenance system: A fuzzy multi-sensor data fusion kalman filter approach for IVHMS. Expert Syst. Appl. 2012, 39, 9821–9836. [Google Scholar] [CrossRef]
  3. Jarske, P.; Saramaki, T.; Mitra, S.K.; Neuvo, Y. On properties and design of nonuniformly spaced linear arrays. IEEE Trans. Acoust. Speech Signal Process. 1988, 36, 372–380. [Google Scholar] [CrossRef]
  4. Haupt, R.L. Thinned arrays using genetic algorithms. IEEE Trans. Antennas Propag. 1994, 42, 993–999. [Google Scholar] [CrossRef]
  5. Yan, K.K.; Lu, Y. Sidelobe reduction in array-pattern synthesis using genetic algorithm. IEEE Trans. Antennas Propag. 1997, 45, 1117–1122. [Google Scholar]
  6. Chen, K.; He, Z.; Han, C. Design of 2-dimensional sparse arrays using an improved genetic algorithm. In Proceedings of the IEEE Workshop on Sensor Array and Multichannel Processing, Waltham, MA, USA, 12–14 July 2006; pp. 209–213.
  7. Cen, L.; Ser, W.; Yu, Z.L.; Rahardja, S. An improved genetic algorithm for aperiodic array synthesis. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, Las Vagas, NV, USA, 31 March–4 April 2008; pp. 2465–2468.
  8. Trucco, A.; Murino, V. Stochastic optimization of linear sparse arrays. IEEE J. Ocean. Eng. 1999, 24, 291–299. [Google Scholar] [CrossRef]
  9. Repetto, S.; Trucco, A. Designing superdirective microphone arrays with a frequency-invariant beam pattern. IEEE Sens. J. 2006, 6, 737–747. [Google Scholar] [CrossRef]
  10. Caorsi, S.; Lommi, A.; Massa, A.; Pastorino, M. Peak sidelobe level reduction with a hybrid approach based on gas and difference sets. IEEE Trans. Antennas Propag. 2004, 52, 1116–1121. [Google Scholar] [CrossRef]
  11. Oliveri, G.; Massa, A. Genetic algorithm (GA)-enhanced almost difference set (ADS)-based approach for array thinning. IET Microwav. Antennas Propag. 2011, 5, 305–315. [Google Scholar] [CrossRef]
  12. Candes, E.; Romberg, J.; Tao, T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory 2006, 52, 489–509. [Google Scholar] [CrossRef]
  13. Ji, S.; Xue, Y.; Carin, L. Bayesian compressive sensing. IEEE Trans. Signal Process. 2008, 56, 2346–2356. [Google Scholar] [CrossRef]
  14. Wang, W.; Wu, R. High resolution direction of arrival (DOA) estimation based on improved orthogonal matching pursuit (OMP) algorithm by iterative local searching. Sensors 2013, 13, 11167–11183. [Google Scholar] [CrossRef] [PubMed]
  15. Li, L.; Zhang, W.; Li, F. The design of sparse antenna array. CoRR 2008, 11, 110–117. [Google Scholar]
  16. Carin, L. On the relationship between compressive sensing and random sensor arrays. IEEE Antennas Propag. Mag. 2009, 51, 72–81. [Google Scholar] [CrossRef]
  17. Cen, L.; Ser, W.; Cen, W.; Yu, Z.L. Linear sparse array synthesis via convex optimization. In Proceedings of the IEEE International Symposium on Circuits and Systems, Paris, France, 30 May–2 June 2010; pp. 4233–4236.
  18. Prisco, G.; D’Urso, M. Exploiting compressive sensing theory in the design of sparse arrays. In Proceedings of the IEEE Radar Conference, Kansas City, MO, USA, 23–27 May 2011; pp. 865–867.
  19. Oliveri, G.; Carlin, M.; Massa, A. Complex-weight sparse linear array synthesis by bayesian compressive sampling. IEEE Trans. Antennas Propag. 2012, 60, 2309–2326. [Google Scholar] [CrossRef]
  20. Fuchs, B. Synthesis of sparse arrays with focused or shaped beampattern via sequential convex optimizations. IEEE Trans. Antennas Propag. 2012, 60, 3499–3503. [Google Scholar] [CrossRef]
  21. Prisco, G.; D’Urso, M. Maximally sparse arrays via sequential convex optimizations. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 192–195. [Google Scholar] [CrossRef]
  22. Candes, E.J.; Wakin, M.B.; Boyd, S.P. Enhancing sparsity by reweighted l1 minimization. J. Fourier Anal. Appl. 2008, 14, 877–905. [Google Scholar] [CrossRef]
  23. Williams, D.B.; Madisetti, V. (Eds.) Digital Signal Processing Handbook, 1st ed.; CRC Press, Inc.: Boca Raton, FL, USA, 1997.
  24. Zhong, X.; Premkumar, A.; Wang, H. Multiple wideband acoustic source tracking in 3-D space using a distributed acoustic vector sensor array. IEEE Sens. J. 2014, 14, 2502–2513. [Google Scholar] [CrossRef]
  25. Zheng, G. A novel spatially spread electromagnetic vector sensor for high-accuracy 2-D DOA estimation. Multidimens. Syst. Signal Process. 2015, 16, 1–26. [Google Scholar] [CrossRef]
  26. Wang, Y.; Xia, W.; He, Z.; Sun, Y. Polarimetric detection for vector-sensor processing in quaternion proper gaussian noises. Multidimens. Syst. Signal Process. 2015, 1–22. [Google Scholar] [CrossRef]
  27. Compton, R.T. The tripole antenna: An adaptive array with full polarization flexibility. IEEE Trans. Antennas Propag. 1981, 29, 944–952. [Google Scholar] [CrossRef]
  28. Leprettre, B.; Martin, N.; Glangeaud, F.; Navarre, J.P. Three-component signal recognition using time, time-frequency, and polarization information-application to seismic detection of avalanches. IEEE Trans. Signal Process. 1998, 46, 83–102. [Google Scholar] [CrossRef]
  29. Bansal, R. Tripole to the rescue [antennas]. IEEE Antennas Propag. Mag. 2001, 43, 106–107. [Google Scholar] [CrossRef]
  30. Wong, K. Direction finding/polarization estimation-dipole and/or loop triad(s). IEEE Trans. Aerosp. Electron. Syst. 2001, 37, 679–684. [Google Scholar] [CrossRef]
  31. Wong, K. Blind beamforming/geolocation for wideband-FFHs with unknown hop-sequences. IEEE Trans. Aerospace Electron. Syst. 2001, 37, 65–76. [Google Scholar] [CrossRef]
  32. Lundback, J.; Nordebo, S. Analysis of a tripole array for polarization and direction of arrival estimation. In Proceedings of the IEEE Workshop on Sensor Array and Multichannel Signal Processing, Barcelona, Spain, 18–21 July 2004; pp. 284–288.
  33. Au-Yeung, C.K.; Wong, K. CRB: sinusoid-sources’ estimation using collocated dipoles/loops. IEEE Trans. Aerospace Electron. Syst. 2009, 45, 94–109. [Google Scholar] [CrossRef]
  34. Yuan, X.; Wong, K.; Xu, Z.; Agrawal, K. Various compositions to form a triad of collocated dipoles/loops, for direction finding and polarization estimation. IEEE Sens. J. 2012, 12, 1763–1771. [Google Scholar] [CrossRef]
  35. Getu, B. Theoretical capacity of tripole array MIMO systems within a fixed area. In Proceedings of the Wireless and Mobile Networking Conference (WMNC), Dubai, UAE, 23–25 April 2013; pp. 1–4.
  36. Zoltowski, M.; Wong, K. Closed-form eigenstructure-based direction finding using arbitrary but identical subarrays on a sparse uniform cartesian array grid. IEEE Trans. Signal Process. 2000, 48, 2205–2210. [Google Scholar] [CrossRef]
  37. Zoltowski, M.; Wong, K. ESPRIT-based 2-D direction finding with a sparse uniform array of electromagnetic vector sensors. IEEE Trans. Signal Process. 2000, 48, 2195–2204. [Google Scholar] [CrossRef]
  38. Wong, K.; Zoltowski, M. Closed-form direction finding and polarization estimation with arbitrarily spaced electromagnetic vector-sensors at unknown locations. IEEE Trans. Antennas Propag. 2000, 48, 671–681. [Google Scholar] [CrossRef]
  39. Wong, K.; Li, L.; Zoltowski, M. Root-MUSIC-based direction-finding and polarization estimation using diversely polarized possibly collocated antennas. IEEE Antennas Wirel. Propag. Lett. 2004, 3, 129–132. [Google Scholar] [CrossRef]
  40. Wong, K.; Yuan, X. Vector cross-product direction-finding with an electromagnetic vector-sensor of six orthogonally oriented but spatially noncollocating dipoles/loops. IEEE Trans. Signal Process. 2011, 59, 160–171. [Google Scholar] [CrossRef]
  41. He, J.; Ahmad, M.O.; Swamy, M. Extended-aperture angle-range estimation of multiplefresnel-region sources with a linear tripole array using cumulants. Signal Process. 2012, 92, 939–953. [Google Scholar] [CrossRef]
  42. Zheng, G.; Wu, B.; Ma, Y.; Chen, B. Direction of arrival estimation with a sparse uniform array of orthogonally oriented and spatially separated dipole-triads. IET Radar Sonar Navig. 2014, 8, 885–894. [Google Scholar] [CrossRef]
  43. Song, Y.; Yuan, X.; Wong, K. Corrections to vector cross-product direction-finding with an electromagnetic vector-sensor of six orthogonally oriented but spatially noncollocating dipoles/loops [jan 11 160-171]. IEEE Trans. Signal Process. 2014, 62, 1028–1030. [Google Scholar] [CrossRef]
  44. Hawes, M.; Liu, W. Location optimization of robust sparse antenna arrays with physical size constraint. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 1303–1306. [Google Scholar] [CrossRef]
  45. Hawes, M.; Liu, W. Sparse array design for wideband beamforming with reduced complexity in tapped delay-lines. IEEE Trans. Acoust. Speech Lang. Process. 2014, 22, 1236–1247. [Google Scholar] [CrossRef]
  46. Hawes, M.; Liu, W. Compressive sensing-based approach to the design of linear robust sparse antenna arrays with physical size constraint. IET Microwav. Antennas Propag. 2014, 8, 736–746. [Google Scholar] [CrossRef] [Green Version]
  47. Hawes, M.; Liu, W. Design of fixed beamformers based on vector sensor arrays. Int. J. Antennas Propag. 2015, 2015, 1–9. [Google Scholar] [CrossRef]
  48. Zhang, X.; Liu, Z.; Liu, W.; Xu, Y. Quasi-vector-cross-product based direction finding algorithm with a spatially stretched tripole. In Proceedings of the 2013 IEEE International Conference of IEEE Region 10 (TENCON13), Xi’an, China, 22–25 October 2013; pp. 1–4.
  49. Research, C. CVX: Matlab Software for Disciplined Convex Programming, Version 2.0 Beta. Available online: http://cvxr.com/cvx (accessed on 20 September 2015).
  50. Grant, M.; Boyd, S. Graph implementations for nonsmooth convex programs. In Recent Advances in Learning and Control (a tribute to M. Vidyasagar; Blondel, V., Boyd, S., Kimura, H., Eds.; Springer-Verlag Limited: Berlin Heidelberg, Germany, 2008; pp. 95–110. [Google Scholar]

Share and Cite

MDPI and ACS Style

Hawes, M.; Liu, W.; Mihaylova, L. Compressive Sensing Based Design of Sparse Tripole Arrays. Sensors 2015, 15, 31056-31068. https://doi.org/10.3390/s151229849

AMA Style

Hawes M, Liu W, Mihaylova L. Compressive Sensing Based Design of Sparse Tripole Arrays. Sensors. 2015; 15(12):31056-31068. https://doi.org/10.3390/s151229849

Chicago/Turabian Style

Hawes, Matthew, Wei Liu, and Lyudmila Mihaylova. 2015. "Compressive Sensing Based Design of Sparse Tripole Arrays" Sensors 15, no. 12: 31056-31068. https://doi.org/10.3390/s151229849

APA Style

Hawes, M., Liu, W., & Mihaylova, L. (2015). Compressive Sensing Based Design of Sparse Tripole Arrays. Sensors, 15(12), 31056-31068. https://doi.org/10.3390/s151229849

Article Metrics

Back to TopTop