Next Article in Journal
Reducing Nervousness in Master Production Planning: A Systematic Approach Incorporating Product-Driven Strategies
Next Article in Special Issue
TransPCGC: Point Cloud Geometry Compression Based on Transformers
Previous Article in Journal
Reinforcement Learning Derived High-Alpha Aerobatic Manoeuvres for Fixed Wing Operation in Confined Spaces
Previous Article in Special Issue
Hardware Suitability of Complex Natural Resonances Extraction Algorithms in Backscattered Radar Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constant-Beamwidth LCMV Beamformer with Rectangular Arrays

by
Vitor Probst Curtarelli
1,* and
Israel Cohen
2
1
Graduate Program in Electrical Engineering, Universidade Federal de Santa Catarina, Florianópolis CEP 88040-900, SC, Brazil
2
Andrew and Erna Viterbi Faculty of Electrical and Computer Engineering, Technion–Israel Institute of Technology, Technion City, Haifa 3200003, Israel
*
Author to whom correspondence should be addressed.
Algorithms 2023, 16(8), 385; https://doi.org/10.3390/a16080385
Submission received: 12 July 2023 / Revised: 27 July 2023 / Accepted: 28 July 2023 / Published: 10 August 2023
(This article belongs to the Special Issue Digital Signal Processing Algorithms and Applications)

Abstract

:
This paper presents a novel approach utilizing uniform rectangular arrays to design a constant-beamwidth (CB) linearly constrained minimum variance (LCMV) beamformer, which also improves white noise gain and directivity. By employing a generalization of the convolutional Kronecker product beamforming technique, we decompose a physical array into virtual subarrays, each tailored to achieve a specific desired feature, and we subsequently synthesize the original array’s beamformer. Through simulations, we demonstrate that the proposed approach successfully achieves the desired beamforming characteristics while maintaining favorable levels of white noise gain and directivity. A comparative analysis against existing methods from the literature reveals that the proposed method performs better than the existing methods.

1. Introduction

Beamformers play a crucial role in diverse fields, such as telecommunications [1], acoustics [2,3], hearing aids [4], and others [5,6,7,8]. Among the array configurations used for beamforming, rectangular arrays are an interesting option to be explored [9,10,11], as they offer distinct advantages over linear arrays, providing enhanced spatial information regarding impinging sources [12,13] and reduced redundancy due to their asymmetry [14].
The development of robust adaptive beamformers with frequency-invariant characteristics has been a significant point of interest as, in this case, frequency does not affect the behavior of the beamformer. Some desired features are a constant beamwidth [15] and null steering [16]. One approach for null steering is using a linearly constrained minimum variance (LCMV) beamformer [17,18,19], which cancels interfering signals from given directions and steers the main beam toward the desired signal. However, it lacks a robust mechanism for maintaining a constant beamwidth. Constant-beamwidth (CB) beamforming [15,20,21] can be accomplished by using window-based beamforming techniques [22], but these methods cannot incorporate directional restrictions. CB-LCMV beamformers have been recently explored [23], however only in the context of linear sensor arrays, leaving space for their study in the context of different array configurations.
Other metrics that are desired to be enhanced are white noise gain and directivity factor, which can respectively be maximized by the delay-and-sum (DS) and superdirective (SD) beamformers [24]. Another quality that is often required for designing a beamformer is a distortionless response to the desired source or to the desired source direction. This ensures that the desired signal is unaltered by filtering processes.
While constructing a beamformer with multiple beamforming features is non-trivial, efforts have been made to combine multiple beamforming techniques for a single array of sensors. Two notable approaches are the Kronecker product (KP) method [25,26] and the linear convolutional Kronecker product (LCKP) method [23]. The LCKP method is valid only for linear arrays, however it allows for the virtual utilization of more sensors than are physically available [23]. Meanwhile, the KP method can be used in linear or rectangular arrays, but it does not increase the number of sensors available for beamforming. For both combination methods (the KP and LCKP), beampattern features and distortionless constraints are preserved in the combined beamformers. Therefore, by using beamformers with desired beampattern characteristics that respect the distortionless constraint, these are maintained in the end result.
In this paper, we propose a novel approach to constructing a CB-LCMV beamformer for rectangular arrays. To this end, we generalize the LCKP beamforming technique to the case of rectangular arrays. We synthesize beamformers for virtual subarrays and, with our proposed generalized technique, apply them to a full array to achieve the desired beamwidth and null placement. This is achieved without sacrificing white noise gain or the directivity factor. The performance of the proposed method is compared against that of beamformers obtained through the KP and LCKP methods. Our results demonstrate superior performance in terms of the beamwidth, white noise gain, and directivity of the beamformers obtained using the proposed method when compared to the performance of the beamformers obtained using the methods in the literature.
This paper is organized as follows: Section 2 presents the array and signal model considered for the problem; Section 3 shows the traditional beamforming techniques and methods that are used further. In Section 4, the newly proposed method for an array analysis for rectangular arrays is introduced and detailed, as well as its proposed usage for solving the problem at hand. In Section 5, we present the conducted simulations and discuss the results, comparing them to those in the literature. Finally, Section 6 concludes this paper, presenting an overview of the main contributions.

2. Signal and Array Model

Let S be a uniform rectangular array (URA) of sensors over the x y plane in an anechoic environment with desired and undesired sources. The URA comprises M x sensors spaced δ x apart along the x -axis and M y sensors spaced δ y apart along the y -axis, resulting in a total of M = M x M y sensors. Assume a source in the far field on the same plane as the sensor array (that is, with elevation ϕ = 0 o ), impinging on it from an azimuth angle θ . As it is unusual, in speech enhancement, for the desired and undesired sources to have the same azimuth, differing only by elevation, we assume the elevation to be 0 o . This constraint can be easily removed without affecting the developed mathematical framework.
Let D ( ω , θ ) denote the steering matrix of size M x × M y with elements { m x , m y } given by
[ D ( ω , θ ) ] m x , m y = exp { j w c r m x , m y cos ( θ ψ m x , m y ) } ,
where c = 340 m/s is the speed of sound, ( r m x , m y , ψ m x , m y ) are the polar coordinates of the sensor at ( m x δ x , m y δ y ) , ω = 2 π f is the angular frequency, f is the temporal frequency, and j = 1 is the imaginary unit. We let d ( ω , θ ) = V ( D ( ω , θ ) ) denote the M × 1 steering vector, with M = m x m y and V ( · ) being the vectorization operation, and we let D ( ω , θ ) = V 1 ( d ( ω , θ ) ; M y ) denote the inverse vectorization of d ( ω , θ ) .
The M × 1 observed signal vector y ( ω ) for all sensors in the frequency domain can be written as
y ( ω ) = d ( ω , θ d ) X ( ω ) + v ( ω ) ,
where X ( ω ) is the desired signal at the reference sensor, d ( ω , θ d ) is the steering vector of the desired source from direction θ d , and v ( ω ) is the additive noise signal vector. All signals are assumed to be zero-mean and uncorrelated. We can estimate the desired signal X ( ω ) as Z ( ω ) using the beamformer h ( ω ) (assumed to be a 2-D beamformer) through linear filtering,
Z ( ω ) = h ( ω ) H y ( ω ) ,
where the superscript H denotes the conjugate-transpose operator. In that case,
Z ( ω ) = X ( ω ) + v rn ( ω ) ,
where v rn ( ω ) = h ( ω ) H v ( ω ) is the residual noise at the beamformer’s output. The beamformer is called distortionless if it satisfies h ( ω ) H d ( ω , θ d ) = 1 for all ω . This constraint guarantees that the beamformer does not affect the desired signal, only altering the undesired noise signal.
From here on, ω is omitted unless in definitions, and θ in the steering vectors appears in subscripts where necessary. When no angle is shown, d = d ( ω , θ d ) is assumed to be the desired signal steering vector for conciseness.

Beamformer Metrics

Beampattern B , as a function of beamformer h and direction θ (through steering vector d θ ) is given by
B ( h , d θ ) = h H d θ .
Given the desired signal steering vector d , the white noise gain (WNG), desired signal distortion index (DSDI), and directivity factor (DF) are, respectively,
(6a) W ( h , d ) = | h H d | 2 h H h , (6b) υ d ( h , d ) = | h H d 1 | 2 , (6c) D ( h , d ) = | h H d | 2 h H Γ h ,
where Γ ( ω ) is the spherical isotropic noise field coherence matrix [27]. Using the DSDI, the distortionless constraint can also be written as υ d ( h , d ) = 0 .

3. Conventional Beamformers

This section briefly overviews the beamforming techniques used to construct the proposed beamformer, as well as the different methods for beamformer synthesis.

3.1. LCMV Beamformer

Assuming the presence of undesired sources in known directions, a linearly constrained minimum variance (LCMV) beamformer [17,18,19] is useful to position the nulls of the beampattern in those undesired directions. To this end, we assume the existence of N (with N < M ) uncorrelated interfering sources in the far field, each coming from a (different) direction θ i ( i { 1 ,   ,   N } ) that we want to cancel. We write v as
v = n = 1 N d θ n ν n + u ,
where ν n is the noise signal for the n-th undesired direction, and u is the portion of the noise signal not coming from the N undesired directions, also accounting for acoustic uncorrelated noise. We assume that E [ u H ( d θ n ν n ) ] = 0 . We then use N + 1 linear constraints, representing the distortionless constraint for the desired signal plus the canceling of the N undesired directions.
The LCMV constraint is written in matrix form as
(8a) C = [ d ,   d θ 1 ,   ,   d θ N ] , (8b) q = [ 1 ,   0 ,   ,   0 ] T , (8c) C H h = q ,
where C ( ω ) is an M × ( N + 1 ) matrix, q is an ( N + 1 ) × 1 vector, and the superscript T denotes the transpose operator. The LCMV beamformer is obtained by minimizing the variance of residual noise v rn in the beamformer output (from Equation (4)), given the constraints in (8), which translates to a minimization problem given as
h LCMV = argmin h   h H Φ v h   s . t .   C H h = q ,
where Φ v ( ω ) = E [ v ( ω ) v H ( ω ) ] is the correlation matrix of v , assumed to be a full-rank invertible matrix. As the undesired directions will be canceled (assuming an anechoic environment [28]), this minimization process minimizes u , the noise portion that is uncorrelated with the N undesired directions. The solution to this minimization problem is
h LCMV = Φ v 1 C [ C H Φ v 1 C ] 1 q .
To ensure the existence of a solution, the number of sensors should be larger than or equal to the number of constraints, i.e., M N + 1 . For N = 0 , only the distortionless constraint remains, and the LCMV beamformer reduces to a minimum variance distortionless response (MVDR) beamformer [29], which is given by
h MVDR = Φ v 1 d d H Φ v 1 d .
It is possible to show that the LCMV and MVDR beamformers are also defined in terms of the observed signal correlation matrix Φ y [24], in which case
(12a) h LCMV = Φ y 1 C [ C H Φ y 1 C ] 1 q (12b) h MVDR = Φ y 1 d d H Φ y 1 d .
This formulation depends only on the statistics of the observed signal, which are easier to compute than those of the noise signal.

3.2. CB Beamformer

Constant-beamwidth (CB) beamformers guarantee a certain beamwidth around the desired direction that is constant over frequency. This is important to ensure the correct receiving of the desired signal, even if θ d is not precisely calibrated.
We define θ B as the first-null beamwidth (FNBW), such that | B ( h , d ) | > 0 if | θ θ d | < θ B / 2 . That is, θ B / 2 is the first angle in which a null of the beampattern occurs. A constant-beamwidth beamformer can be achieved using a window-based design technique [22]. Here, the Kaiser window is used [30], which can be written as
[ w ] m = J 0 ( β 1 [ 2 m M 1 1 ] 2 ) J 0 ( β ) ,
where J 0 ( · ) is the zero-order modified Bessel function of the first kind, and β ( ω ) ([22], Section 3.3.2) is frequency-dependent to maintain the θ B constant. This technique requires that the desired source signal impinges on the array from the broadside direction [22]. To satisfy the distortionless constraint, we normalize w , obtaining
h CB = w m = 0 M 1 [ w ] m .

3.3. SD and DS Beamformers

Superdirective (SD) and delay-and-sum (DS) beamformers are obtained by maximizing the DF and the WNG, respectively [29,31], both subject to the distortionless constraint. The solutions to these minimization problems are respectively given by
(15a) h SD = G 1 d d H G 1 d , (15b) h DS = d M .

3.4. Kronecker Product Beamforming

Designing a single beamformer with different features is highly important, allowing different effects to be employed over a single array of sensors. Two methods to accomplish such a task are the Kronecker product (KP) [25,26] and linear convolutional Kronecker product (LCKP) methods [23]. In these processes, the sensor array is split into subarrays for which we design separate beamformers and use the chosen technique to synthesize the whole sensor array’s (or full array’s) beamformer.
The KP beamforming process is as follows: given steering vector d θ , we decompose it into two parts (namely, d 1 ; θ and d 2 ; θ ) satisfying the relation d θ = d 1 ; θ d 2 ; θ , where ⊗ represents the Kronecker product. By designing beamformers h 1 and h 2 for d 1 and d 2 , respectively, we obtain the beamformer for the full array as h = h 1 h 2 [32].
LCKP beamforming is achieved similarly: given a uniform linear array’s (ULA) steering vector d θ of length M, we define d 1 ; θ as the M 1 -th first elements of d θ and, similarly, d 2 ; θ with M 2 elements, respecting M 1 + M 2 1 = M . By designing beamformers h 1 and h 2 for each subarray, the full array’s beamformer is h = h 1 h 2 [23], where ∗ denotes the linear convolution.
For both the KP and LCKP methods, we have the following properties:
(16a) B ( h , d θ ) = B ( h 1 , d 1 ; θ ) B ( h 2 , d 2 ; θ ) , (16b) υ d ( h , d ) [ 1 + υ d ( h 1 , d 1 ) ] [ 1 + υ d ( h 2 , d 2 ) ] 1 .
The first one shows that the beampattern of h is the combination of the beampatterns of the subarrays. Through the second one, we can see that, if h 1 and h 2 are distortionless beamformers, h also will be. From these properties, we can see that the beampattern properties (including the distortionless feature) from h 1 and h 2 are maintained for h .

4. Constant-Beamwidth LCMV Beamformer with Rectangular Arrays

Both the KP and LCKP beamforming methods have advantages and disadvantages. While the LCKP is only usable over linear arrays, beamformers achieved through it have (virtually) more sensors than there are available in the physical array, generally leading to better performance when combining different techniques. Meanwhile, the KP method can be applied to rectangular arrays, which on its own is beneficial, but it does not have the virtual utilization of more sensors. To benefit from the lack of symmetry of the rectangular arrays and to also be able to have more sensors available for each subarray’s beamformer, we propose a generalization of the LCKP to URAs to take advantage of both methods, enabling virtual sensor augmentation while exploiting the rectangular array’s symmetry.

4.1. Rectangular Convolutional Kronecker Product Beamforming

Let S 1 be a subarray of S, including the M 1 , x × M 1 , y first sensors of S, with steering vector d 1 , and similarly for S 2 . These arrays sizes are such that M 1 , x + M 2 , x 1 = M x , and M 1 , y + M 2 , y 1 = M y . By designing beamformers h 1 and h 2 for S 1 and S 2 , respectively, we show in Appendix A that it is possible to synthesize a beamformer for full array S through
h   =   V ( V 1 ( h 1 ;   M 1 , y )     V 1 ( h 2 ;   M 2 , y ) ) ,
with ⊛ representing the 2-D convolution operation. We call this the rectangular convolutional Kronecker product (RCKP) method. A simple implementation of the proposed RCKP method is presented in Algorithm A1 in Appendix B, which is written in a Python-like pseudolanguage.
In the same way as with the KP and LCKP methods, with the RCKP, one can design beamformers for subarrays S 1 and S 2 , and one can synthesize the URA’s beamformer h through the 2-D convolution (and vectorization processes). Accordingly, h is an M × 1 vector; however, in this scenario we have a total of M = M 1 , x M 1 , y M 2 , x M 2 , y virtual sensors. Assuming that all Ms are ≥1, it is trivial to see that M M (equality happens if S 1 and S 2 are perpendicular ULAs, or if one of them has only one sensor). Therefore, we have more sensors virtually than if S was split into two VAs through the KP method.
It is easy to verify that the properties in (16) are also valid for the proposed RCKP. With this, like for the KP and LCKP, the full-array beamformer obtained through the RCKP inherits the beampattern and distortionless features from the subarray beamformers. Also, since convolution is commutative and associative, one can split the S array into more than two virtual arrays, design a beamformer for each subarray, and combine all the beamformers through the RCKP without the loss of generalization. In this case, assuming that we are using K beamformers, each with size M k , x × M k , y , then their dimensions must be such that
(18a) k = 1 K M k , x ( K 1 ) = M x , (18b) k = 1 K M k , y ( K 1 ) = M y .
This is easily verifiable by repeating the synthesis operation in Equation (17) K 1 times.

4.2. CB-LCMV Beamformer with RCKP

Here, we propose the use of the RCKP method as derived previously to construct a CB-LCMV beamformer with an increase in white noise gain and directivity measures. For this, full array S is separated into four subarrays: S 1 , S 2 ; x , S 2 ; x , and S 2 ; y . Each subarray is used to design one of the desired beamformers presented in Section 3.1, Section 3.2 and Section 3.3.
S 1 is used to design the LCMV beamformer, following the steps detailed in Section 3.1. Since the LCMV beamformer does not require the array to be linear, we use S 1 as a rectangular array of size M 1 , x × M 1 , y . As explained previously, this choice has the advantage of the rectangular array’s lesser symmetry compared to a linear array. We have M 1 = M 1 , x M 1 , y virtual sensors and, at most, M 1 1 nulls to be placed. If N = 0 , this same array is used for the MVDR beamformer instead.
The subarray S 2 ; y is used to construct the CB beamformer. Given that the CB is achieved through the window technique (as per Section 3.2), this subarray must be a ULA, and the desired source direction should be on its broadside direction. Here, this implies that θ d = 0 o .
The SD and DS beamformers are built from the S 2 ; x and S 2 ; x subarrays, respectively, based on Section 3.3. We assume that they are constructed from linear arrays, but this condition is flexible and can be changed in other implementations.
Once all the subarray beamformers are designed, we use the RCKP method to combine these beamformers into a full-array beamformer. Thus, we construct a beamformer with many desired features (null placement + constant beamwidth + white noise gain + directivity factor gain) that exploit the symmetry of the rectangular array, which implies more spatial information and more performance. Algebraically, the full-array beamformer is given by
h = V ( V 1 ( h LCMV ; M 1 , y ) V 1 ( h SD ; 1 ) V 1 ( h DS ; 1 ) V 1 ( h CB ; M 2 , y ) ) .
An implementation of the proposed CB-LCMV beamformer is shown in Algorithm A2 in Appendix B. The procedure calculates beamformer h for a single frequency.
For the CB beamformer to be effective, its condition must be valid for all beamformers. That is, | B ( h , d ) | > 0 if | θ θ d | < θ B / 2 , where h [ h LCMV , h CB , h SD , h DS ] . This is true by definition for the CB beamformer, and by setting N = M 1 1 for the LCMV, no nulls are free to end up inside the main beam. For the SD and DS beamformers, we assume that their FNBWs are sufficient to satisfy the condition.

5. Experimental Results

In this section, we present the simulations performed to verify the proposed method’s implementation and to compare its performance to that of existing methods. The code used for these simulations is available at https://github.com/VCurtarelli/py-cb-lcmv-rect (accessed on 27 July 2023). We test different combinations of the number of sensors for each subarray, as established in Table 1. The table presents the dimensions for each subarray, for the full array of sensors, and the virtual number of sensors being employed. [A, B] are obtained through the RCKP, [C, D] through the LCKP, and [E, F] through the KP + LCKP, with the KP and LCKP methods being as defined in Section 3.4. [A, C, E] utilize the SD and DS beamformers, while [B, D, F] “disable” them in favor of the CB beamformer. For [E, F], the LCMV beamformer is more spaced out.
In all situations, S has a total of M = 36 sensors. The sensors are assumed to be ideal omnidirectional sensors with a plain frequency response over the spectrum. The intersensor distances are δ x = 0.5 cm and δ y = 3.0 cm . We assume that θ d = 0 o and θ B = 40 o . Since the LCMV has four sensors, we use three interfering sources (for all situations), with directions θ i [ 90 o , 60 o , 130 o ] , each with ϕ v i = 1 , and we also assume v to be uncorrelated Gaussian white noise with unit variance (that is, ϕ v is the identity matrix). The variance of the desired signal is ϕ X = E [ | X | 2 ] = 5 . The simulations are performed for the range f [ 4 , 8 ] kHz . This range is chosen to satisfy the conditions from ([22], Equation (29)) for condition [A].
Figure 1 and Figure 2a–c show the simulation results for B , measured θ B , W , and D . θ B was measured as the first angle in which B ( h , d θ B / 2 ) 0.05 . All methods manage to achieve the nulls in the desired directions from the LCMV. Observing Figure 2a, all but [E] are capable of maintaining a reasonable beamwidth across all frequencies, with [E] not being able to maintain the beamwidth because of too few sensors in the CB beamformer. The “naïve” combination of KP and LCKP [E, F] leads to worse performance for all metrics; thus, their results are not further compared.
Both methods (the RCKP and LCKP) lead to very akin results for θ B , even though the LCKP simulations [C, D] have many more sensors in the CB beamformer than the RCKP ones [A, B]. This indicates that increasing the CB array size does not result in a better performance in either FNBW or directivity. However, in Figure 1, it can be seen that it generates a more focused beam. Both the RCKP and LCKP have similar WNG results of around 10 dB. The proposed method (especially [A]) leads to a better performance in terms of DF for all frequencies compared to the LCKP. The LCKP is marginally better for WNG in lower frequencies and does not suffer from performance loss at approximately 6.4   kHz , caused by using rectangular arrays.
Comparing [A] and [B], the former has better results for both WNG and DF, while both have a similar result in maintaining a constant θ B . This is a direct result of using the SD and DS beamformers in [A], and although [A] has less than half the number of sensors in the CB beamformer than [B], the FNBWs are similar. The beampattern of [B] is more focused than that of [A], but this does not cause a better performance in either beamwidth or directivity.
Comparing [C] and [D], both lead to almost identical results for all metrics. This can also be seen by comparing Figure 1c,d, where their beampatterns are almost indistinguishable. This indicates that the SD and DS beamformers are obfuscated by the CB in [C, D], caused by the CB beamformer having too many sensors compared to the other beamformers.

6. Conclusions

We have introduced a novel approach for designing a constant-beamwidth beamformer with null-direction constraints, utilizing uniform rectangular arrays and a generalized convolutional Kronecker product beamforming technique. By synthesizing beamformers for virtual arrays and applying our proposed technique to a full array, we successfully achieved the desired features in terms of beamwidth and null placement. Moreover, by using virtual arrays and beamformers, we demonstrated the ability to enhance the signal quality in terms of white noise gain and the directivity factor without compromising the beamwidth. The experimental results using simulated sensor arrays demonstrate that the performance of our method surpassed the performance obtained using only known techniques based on the Kronecker product for beamforming synthesis when assessing beamwidth, white noise gain, and directivity.

Author Contributions

Conceptualization, I.C. and V.P.C.; Methodology, V.P.C.; Software, V.P.C.; Writing—original draft: V.P.C.; Writing—review and editing, I.C. and V.P.C.; Supervision, V.P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Pazy Research Foundation and the Israel Science Foundation (grant no. 1449/23).

Data Availability Statement

The source code for the simulations developed herein is available at https://github.com/VCurtarelli/py-cb-lcmv-rect (accessed on 27 July 2023).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CBConstant beamwidth
CKPConvolutional KP
DFDirectivity factor
DSDelay and sum
DSDIDesired signal distortion index
FNBWFirst-null beamwidth
KPKronecker product
LCKPLinear CKP
LCMVLinearly constrained minimum variance
MVDRMinimum variance distortionless response
RCKPRectangular CKP
SDSuperdirective
ULAUniform linear array
URAUniform rectangular array
WNGWhite noise gain

Appendix A. Proof of CKP Beamforming for URAs

We assume a rectangular uniform sensor array of size M x × M y , with steering vector d , and two subarrays of sizes M 1 , x × M 1 , y and M 2 , x × M 2 , y , with steering vectors d 1 and d 2 , such that M 1 , x + M 2 , x 1 = M x and M 1 , y + M 2 , y 1 = M y . We define
(A1a) h ~   =   h 1     h 2 , (A1b) d ~   =   d 1     d 2 .
From this, it easily follows that
h ~ H d ~   =   ( h 1 H d 1 ) ( h 2 H d 2 ) .
Expanding both,
h H d = ( m 1 = 0 M 1 1 h 1 * [ m 1 ] d [ m 1 ] ) ( m 2 = 0 M 2 1 h 2 * [ m 2 ] d [ m 2 ] ) .
We define H 1 = V 1 ( h 1 ; M 1 , y ) as the inverse vectorization of h 1 ; similarly, D = V 1 ( d ; M y ) is the inverse vectorization of d , and H 2 = V 1 ( h 2 ; M 2 , y ) is the inverse vectorization of h 2 . Using the inverse vectorization on the sum over m 1 , we have
m 1 = 0 M 1 1 h 1 * [ m 1 ] d [ m 1 ] = m 1 , x = 0 M 1 , x 1 m 1 , y = 0 M 1 , y 1 H 1 * [ m 1 , x , m 1 , y ] D [ m 1 , x , m 1 , y ] .
We also assume that D [ m 1 , x , m 1 , y ] = d x [ m 1 , x ] d y [ m 1 , y ] , by the definition of the steering vector. Applying a similar process to the sum over m 2 from (A3),
h ~ H d ~ = ( m 1 , y = 0 m 1 , x = 0 M 1 , x 1 M 1 , y 1 H 1 * [ m 1 , x , m 1 , y ] d x [ m 1 , x ] d y [ m 1 , y ] ) × ( m 2 , y = 0 m 2 , x = 0 M 2 , x 1 M 2 , y 1 H 1 * [ m 2 , x , m 2 , y ] d x [ m 2 , x ] d y [ m 2 , y ] ) ,
By applying the Cauchy product [33] to the second sum twice,
h ~ H d ~ = m y = 0 m x = 0 M 1 , x + M 2 , x 2 M 1 , y + M 2 , y 2 n y = k 1 , y n x = k 1 , x k 2 , x k 2 , y ( H 1 * [ n x , n y ] d x [ n x ] d y [ n y ] ) × ( H 2 * [ m x n x , m y n y ] d x [ m x n x ] d y [ m y n y ] )       = m y = 0 m x = 0 M x 1 M y 1 n y = k 1 , y n x = k 1 , x k 2 , x k 2 , y ( H 1 * [ n x , n y ] H 2 * [ m x n x , m y n y ] ) × ( d x [ n x ] d y [ n y ] d x [ m x n x ] d y [ m y n y ] ) .
By the definition of the steering vectors, d x [ a ] d x [ b ] = d x [ a + b ] (also for d y ), and, therefore, the d term becomes d x [ m x ] d y [ m y ] = D [ m x , m y ] . By noting that the sum over H 1 and H 2 is the 2-D convolution between them (at the element [ m x , m y ] ), and with this defining H = H 1 H 2 (where ⊛ denotes the 2-D convolution), then
h ~ H d ~ = m y = 0 m x = 0 M x 1 M y 1 H * [ m x , m y ] D [ m x , m y ] .
Finally, by vectorizing H and D into h and d ,
h ~ H d ~ = m y = 0 m x = 0 M x 1 M y 1 h * [ M y m x + m y ] d [ M y m x + m y ] = m = 0 M 1 h * [ m ] d [ m ] = h H d ,
which concludes the proof.

Appendix B. Pseudocode Algorithms

These algorithms are written in a Python-like pseudolanguage. Brackets can denote both vector definition and vector slicing/indexing. For example, a = [ 1 , 2 , 3 , 4 , 5 ] denotes a 1 × 5 vector, while b = a [ 0 : 3 ] denotes a 1 × 3 vector, such that b = [ 1 , 2 , 3 ] (last exclusive slicing, zero indexing). Comments are added where necessary to clarify the steps.
Algorithm A1 RCKP beamforming algorithm
Algorithms 16 00385 i001
Algorithm A2 CB-LCMV beamformer algorithm
Algorithms 16 00385 i002

References

  1. Viswanath, P.; Tse, D.; Laroia, R. Opportunistic beamforming using dumb antennas. IEEE Trans. Inf. Theory 2002, 48, 1277–1294. [Google Scholar] [CrossRef] [Green Version]
  2. Herbordt, W.; Nakamura, S.; Kellermann, W. Joint Optimization of LCMV Beamforming and Acoustic Echo Cancellation for Automatic Speech Recognition. In Proceedings of the (ICASSP ’05), IEEE International Conference on Acoustics, Speech, and Signal Processing, Philadelphia, PA, USA, 23–23 March 2005; Volume 3, pp. 77–80. [Google Scholar] [CrossRef]
  3. Chiariotti, P.; Martarelli, M.; Castellini, P. Acoustic beamforming for noise source localization—Reviews, methodology and applications. Mech. Syst. Signal Process. 2019, 120, 422–448. [Google Scholar] [CrossRef]
  4. Haykin, S.; Liu, K.J.R. Handbook on Array Processing and Sensor Networks; Wiley-IEEE Press: New York, NY, USA, 2009. [Google Scholar]
  5. Van Veen, B.; Buckley, K. Beamforming: A versatile approach to spatial filtering. IEEE ASSP Mag. 1988, 5, 4–24. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, W.; Weiss, S. Wideband Beamforming Concepts and Techniques; Wiley: New York, NY, USA, 2010. [Google Scholar]
  7. Huang, Q.; Lin, M.; Wang, J.B.; Tsiftsis, T.A.; Wang, J. Energy Efficient Beamforming Schemes for Satellite-Aerial-Terrestrial Networks. IEEE Trans. Commun. 2020, 68, 3863–3875. [Google Scholar] [CrossRef]
  8. Elbir, A.M.; Mishra, K.V.; Vorobyov, S.A.; Heath, R.W. Twenty-Five Years of Advances in Beamforming: From convex and nonconvex optimization to learning techniques. IEEE Signal Process. Mag. 2023, 40, 118–131. [Google Scholar] [CrossRef]
  9. Gu, P.; Wang, G.; Fan, Z.; Chen, R. An Efficient Approach for the Synthesis of Large Sparse Planar Array. IEEE Trans. Antennas Propag. 2019, 67, 7320–7330. [Google Scholar] [CrossRef]
  10. Zhang, X.; Zheng, W.; Chen, W.; Shi, Z. Two-dimensional DOA estimation for generalized coprime planar arrays: A fast-convergence trilinear decomposition approach. Multidimens. Syst. Signal Process. 2019, 30, 239–256. [Google Scholar] [CrossRef]
  11. Lin, Z.; Lin, M.; Champagne, B.; Zhu, W.P.; Al-Dhahir, N. Secrecy-Energy Efficient Hybrid Beamforming for Satellite-Terrestrial Integrated Networks. IEEE Trans. Commun. 2021, 69, 6345–6360. [Google Scholar] [CrossRef]
  12. Heidenreich, P.; Zoubir, A.M.; Rubsamen, M. Joint 2-D DOA Estimation and Phase Calibration for Uniform Rectangular Arrays. IEEE Trans. Signal Process. 2012, 60, 4683–4693. [Google Scholar] [CrossRef]
  13. Ioannides, P.; Balanis, C. Uniform circular and rectangular arrays for adaptive beamforming applications. IEEE Antennas Wirel. Propag. Lett. 2005, 4, 351–354. [Google Scholar] [CrossRef]
  14. Singh, S. Minimal Redundancy Linear Array and Uniform Linear Arrays Beamforming Applications in 5G Smart Devices. Emerg. Sci. J. 2021, 4, 70–84. [Google Scholar] [CrossRef]
  15. Goodwin, M.; Elko, G. Constant beamwidth beamforming. In Proceedings of the IEEE International Conference on Acoustics Speech and Signal Processing, Minneapolis, MN, USA, 27–30 April 1993; Volume 1, pp. 169–172. [Google Scholar] [CrossRef]
  16. Zarifi, K.; Affes, S.; Ghrayeb, A. Collaborative Null-Steering Beamforming for Uniformly Distributed Wireless Sensor Networks. IEEE Trans. Signal Process. 2010, 58, 1889–1903. [Google Scholar] [CrossRef]
  17. Frost, O. An algorithm for linearly constrained adaptive array processing. Proc. IEEE 1972, 60, 926–935. [Google Scholar] [CrossRef]
  18. Buckley, K. Spatial/Spectral filtering with linearly constrained minimum variance beamformers. IEEE Trans. Acoust. Speech Signal Process. 1987, 35, 249–266. [Google Scholar] [CrossRef]
  19. Souden, M.; Benesty, J.; Affes, S. A Study of the LCMV and MVDR Noise Reduction Filters. IEEE Trans. Signal Process. 2010, 58, 4925–4935. [Google Scholar] [CrossRef]
  20. Hixson, E.L.; Au, K.T. Wide-Bandwidth Constant-Beamwidth Acoustic Array. J. Acoust. Soc. Am. 1970, 48, 117. [Google Scholar] [CrossRef]
  21. Wang, Z.; Li, J.; Stoica, P.; Nishida, T.; Sheplak, M. Constant-beamwidth and constant-powerwidth wideband robust Capon beamformers for acoustic imaging. J. Acoust. Soc. Am. 2004, 116, 1621–1631. [Google Scholar] [CrossRef]
  22. Long, T.; Cohen, I.; Berdugo, B.; Yang, Y.; Chen, J. Window-Based Constant Beamwidth Beamformer. Sensors 2019, 19, 2091. [Google Scholar] [CrossRef] [Green Version]
  23. Frank, A.; Ben-Kish, A.; Cohen, I. Constant-Beamwidth Linearly Constrained Minimum Variance Beamformer. In Proceedings of the 2022 30th European Signal Processing Conference (EUSIPCO), Belgrade, Serbia, 29 August–2 September 2022; pp. 50–54. [Google Scholar] [CrossRef]
  24. Benesty, J.; Chen, J.; Huang, Y. Microphone Array Signal Processing; Springer Topics in Signal Processing; Springer: Berlin/Heidelberg, Germany, 2008; Volume 1. [Google Scholar] [CrossRef]
  25. Abramovich, Y.I.; Frazer, G.J.; Johnson, B.A. Iterative Adaptive Kronecker MIMO Radar Beamformer: Description and Convergence Analysis. IEEE Trans. Signal Process. 2010, 58, 3681–3691. [Google Scholar] [CrossRef]
  26. Werner, K.; Jansson, M.; Stoica, P. On Estimation of Covariance Matrices With Kronecker Product Structure. IEEE Trans. Signal Process. 2008, 56, 478–491. [Google Scholar] [CrossRef]
  27. Habets, E.A.P.; Gannot, S. Generating sensor signals in isotropic noise fields. J. Acoust. Soc. Am. 2007, 122, 3464–3470. [Google Scholar] [CrossRef] [PubMed]
  28. Markovich-Golan, S.; Gannot, S.; Kellermann, W. Combined LCMV-TRINICON Beamforming for Separating Multiple Speech Sources in Noisy and Reverberant Environments. IEEE/ACM Trans. Audio Speech Lang. Process. 2017, 25, 320–332. [Google Scholar] [CrossRef]
  29. Erdogan, H.; Hershey, J.R.; Watanabe, S.; Mandel, M.I.; Roux, J.L. Improved MVDR Beamforming Using Single-Channel Mask Prediction Networks. In Proceedings of the Interspeech 2016, ISCA—17th Annual Conference of the International Speech Communication Association, San Francisco, CA, USA, 8–12 September 2016; pp. 1981–1985. [Google Scholar] [CrossRef] [Green Version]
  30. Kaiser, J.; Schafer, R. On the use of the I0-sinh window for spectrum analysis. IEEE Trans. Acoust. Speech Signal Process. 1980, 28, 105–107. [Google Scholar] [CrossRef]
  31. Brandstein, M.; Ward, D.; Lacroix, A.; Venetsanopoulos, A. (Eds.) Microphone Arrays: Signal Processing Techniques and Applications; Digital Signal Processing; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar] [CrossRef]
  32. Huang, G.; Benesty, J.; Chen, J.; Cohen, I. Robust and steerable kronecker product differential beamforming With rectangular microphone arrays. In Proceedings of the ICASSP 2020—2020 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Barcelona, Spain, 4–8 May 2020; pp. 211–215. [Google Scholar] [CrossRef]
  33. Hamahata, Y. Arithmetic functions and the Cauchy product. Arch. Math. 2020, 114, 41–50. [Google Scholar] [CrossRef]
Figure 1. Beampattern heatmap for all situations in Table 1. x -axis is direction of source (in o ), and y -axis is the frequency (in kHz).
Figure 1. Beampattern heatmap for all situations in Table 1. x -axis is direction of source (in o ), and y -axis is the frequency (in kHz).
Algorithms 16 00385 g001
Figure 2. (a) FNBW, (b) WNG, and (c) DF for the beamformers designed with the parameters in Table 1.
Figure 2. (a) FNBW, (b) WNG, and (c) DF for the beamformers designed with the parameters in Table 1.
Algorithms 16 00385 g002
Table 1. Number of sensors for each subarray, dimension of the full array (FA), and size of the virtual array in each simulation.
Table 1. Number of sensors for each subarray, dimension of the full array (FA), and size of the virtual array in each simulation.
Cond.LCMVSDDSCBFA M
A 2 × 2 2 × 1 2 × 1 1 × 8 4 × 9 128
B 2 × 2 1 × 1 1 × 1 1 × 17 2 × 18 68
C 1 × 4 1 × 2 1 × 2 1 × 31 1 × 36 496
D 1 × 4 1 × 1 1 × 1 1 × 33 1 × 36 132
E 2 × 2 2 × 1 2 × 1 1 × 3 6 × 6 48
F 2 × 2 1 × 1 1 × 1 1 × 9 1 × 18 36
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Curtarelli, V.P.; Cohen, I. Constant-Beamwidth LCMV Beamformer with Rectangular Arrays. Algorithms 2023, 16, 385. https://doi.org/10.3390/a16080385

AMA Style

Curtarelli VP, Cohen I. Constant-Beamwidth LCMV Beamformer with Rectangular Arrays. Algorithms. 2023; 16(8):385. https://doi.org/10.3390/a16080385

Chicago/Turabian Style

Curtarelli, Vitor Probst, and Israel Cohen. 2023. "Constant-Beamwidth LCMV Beamformer with Rectangular Arrays" Algorithms 16, no. 8: 385. https://doi.org/10.3390/a16080385

APA Style

Curtarelli, V. P., & Cohen, I. (2023). Constant-Beamwidth LCMV Beamformer with Rectangular Arrays. Algorithms, 16(8), 385. https://doi.org/10.3390/a16080385

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