Next Article in Journal
PC5-Based Cellular-V2X Evolution and Deployment
Next Article in Special Issue
Physical-Layer Security Improvement with Reconfigurable Intelligent Surfaces for 6G Wireless Communication Systems
Previous Article in Journal
Recognition of Abnormal Chest Compression Depth Using One-Dimensional Convolutional Neural Networks
Previous Article in Special Issue
Design of mmWave Directional Antenna for Enhanced 5G Broadcasting Coverage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low-Cost Unattended Design of Miniaturized 4 × 4 Butler Matrices with Nonstandard Phase Differences

by
Adrian Bekasiewicz
1,* and
Slawomir Koziel
1,2
1
Faculty of Electronics, Telecommunications and Informatics, Gdansk University of Technology, 80-233 Gdansk, Poland
2
Engineering Optimization & Modeling Center, Reykjavik University, 101 Reykjavik, Iceland
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(3), 851; https://doi.org/10.3390/s21030851
Submission received: 28 December 2020 / Revised: 22 January 2021 / Accepted: 25 January 2021 / Published: 27 January 2021
(This article belongs to the Special Issue Advanced Antenna Techniques for IoT and 5G Applications)

Abstract

:
Design of Butler matrices dedicated to Internet of Things and 5th generation (5G) mobile systems—where small size and high performance are of primary concern—is a challenging task that often exceeds capabilities of conventional techniques. Lack of appropriate, unified design approaches is a serious bottleneck for the development of Butler structures for contemporary applications. In this work, a low-cost bottom-up procedure for rigorous and unattended design of miniaturized 4 × 4 Butler matrices is proposed. The presented approach exploits numerical algorithms (governed by a set of suitable objective functions) to control synthesis, implementation, optimization, and fine-tuning of the structure and its individual building blocks. The framework is demonstrated using two miniaturized matrices with nonstandard output-port phase differences. Numerical results indicate that the computational cost of the design process using the presented framework is over 80% lower compared to the conventional approach. The footprints of optimized matrices are only 696 and 767 mm2, respectively. Small size and operation frequency of around 2.6 GHz make the circuits of potential use for mobile devices dedicated to work within a sub-6 GHz 5G spectrum. Both structures have been benchmarked against the state-of-the-art designs from the literature in terms of performance and size. Measurements of the fabricated Butler matrix prototype are also provided.

1. Introduction

Antenna arrays are the key components of modern communication devices. Their popular applications include long-term evolution and 5th generation (5G) cellular technology, where good performance is crucial for sustaining high data-transfer rates. Apart from the radiators, feeding network is an integral component of the antenna array. Its role is to provide appropriate excitation of the radiators (both magnitude- and phase-wise) so as to ensure the desired beamforming capability of the array [1,2,3,4,5,6].
Common array feeding realizations include variants of corporate [7,8,9] and series networks [10,11,12]. More complex structures are based on a combination of series–parallel feeds [13], networks with tunable power dividers [14], as well as structures dedicated to operating in mm-wave spectrum [15]. Butler matrices (BMs) and their derivatives belong to another class of feeding networks [16,17,18,19,20]. Their characteristic feature is the beamforming capability resulting from the ability to provide different phase shifts between the output ports depending on the selected input [17,20]. The conventional 4 × 4 (i.e., with four input and four output ports) BM—which offers output-port phase shifts of ±45° and ±135°, respectively—is a composite structure comprising four 90° hybrid couplers, with two crossovers and two 45° phase shifters (PSs) [21,22,23,24]. The circuit is characterized by large dimensions and lack of flexibility in terms of attainable phase differences. Consequently, the usefulness of standard BMs for modern mobile systems, including sensor networks and/or Internet of Things (IoT) devices interconnected using the 5G backbone, is limited. The following paragraphs contain more details concerning design of BM structures for contemporary applications.
In recent years, the design of Butler matrices for space-limited applications has gained significant attention of the research community [16,25,26,27,28,29,30,31,32,33]. Popular miniaturization techniques are oriented toward replacing BM building blocks with their miniaturized counterparts [29,30]. The latter are often implemented in the form of composite cells, representing appropriate combination of the high- and low-impedance transmission lines (TLs) [29,34]. However, the application of left/right-handed TLs for reduction of BM footprint has also been reported [28]. Miniaturization of conventional crossovers—realized as a composition of two branch line couplers (BLCs) interconnected through a 90° TL—is one of the leading approaches to design of small BMs [28,29,30]. Notwithstanding, substitution of conventional BM building blocks with their miniaturized counterparts results in modest miniaturization rates [28,30]. From this perspective, replacement of standard crossovers with alternative topologies featuring either small footprints or improved phase-related functionality seems to be an interesting and preferable alternative to compact BM design [19,35]. Other miniaturization methods depart from standard BM topologies. Instead, they eliminate certain matrix subcomponents or exploit multilayer substrates for structure implementation [26,27,32]. On the other hand, miniaturization resulting from application of the mentioned techniques is often insufficient to increase usefulness of BMs for space-limited systems [28,30]. Furthermore, the discussed compact topologies do not address the problem of improving BM performance.
The design of Butler matrices with increased flexibility in terms of achievable output-port phase differences has also been investigated in the literature [19,23,25,36]. The considered methods involve redesign of conventional structures, either through modification of their building blocks (i.e., couplers or delay lines) [19,23] or introduction of additional sub-circuits affecting BM performance [25,36]. In [23], enhanced control over beamforming capabilities has been achieved by replacing conventional 90° hybrid couplers by structures with adjustable phase shifts. In [19], similar effect has been achieved using PSs with unequal electrical lengths. The reasoning behind both methods is that improved control of phase can be achieved by increasing the number of relevant design parameters. An alternative technique, proposed in [36], boils down to enhancement of the conventional BM through connection of phase-reconfigurable TLs to its output ports. The mentioned methods proved to be useful for performance enhancement. However, they neglect the need of circuit miniaturization, which is important for application of matrices in contemporary systems.
The design solutions discussed above indicate that the difficulties related to improving BMs’ functionality and their miniaturization are perceived as separate problems. In other words, size-reduction strategies considered in the literature are dedicated to conventional topologies [28,29], whereas BM designs featuring improved performance suffer from large dimensions [23,36]. The development of BMs characterized by both small size and improved functionality is important for contemporary, space-limited applications such as IoT and sensor networks. Another problem is that the discussion of methods used to obtain specific solutions of miniaturized and/or performance-enhanced Butler matrices is often neglected in the literature. Instead, the emphasis is put on structure synthesis or analysis of the specific case study. In practice, however, synthesis (or generation of BM components) is just the beginning of the design procedure, which has to be followed by circuit assembly and careful adjustment of parameters oriented toward determination of the desired performance.
Due to high complexity, Butler matrix design is realized as a multistage process. The structure is first divided into individual building blocks, which are subject to a set of simplified design tasks oriented toward obtaining the desired performance [29]. Once the design parameters of all subcomponents are found, the matrix is assembled. The main benefit of such a strategy is that each step gradually approximates a satisfactory design solution. On the other hand, conventional design approaches involve time-consuming and repetitive analysis of the structure (and/or its building blocks). These often involve modifications of geometry parameters followed by visual inspection of structure responses [34,37]. Although the concept of semi-manual design proved to be fairly successful for circuits characterized by a relatively low number of variables and simple responses [37,38,39], its applicability to more complex problems, such as tuning of corporate feeds or BM topologies, is limited. One of the reasons is that the discussed structures are characterized by complex electrical responses that cannot be reliably tracked using manual techniques. From this perspective, Butler matrix design stages should include synthesis, development of BM components, and fine-tuning of the assembled circuit oriented toward maximization of performance. For the sake of reliability, each step should be controlled by appropriate algorithm and rigorously defined objective function. Yet another problem—rarely considered in the literature in terms of BM design [29]—is a large design cost associated with BM optimization. The latter stems from the necessity of using numerically expensive electromagnetic (EM) simulation models in order to obtain accurate responses of the matrix, and a large number of evaluations required by the optimization algorithm to converge.
Analysis of the available literature indicates that the design of BM structures for IoT and 5G applications is pertinent to several challenges. First of all, the development of portable electronics is associated with steadily increasing requirements concerning both miniaturization and performance. However, for BM structures, these problems are treated separately, which results in the development of either relatively small components with standard (or close-to standard) performance or bulky circuits with improved functionality [19,23,25,36]. Another problem involves lack of unified and well-established procedure—ranging from components synthesis to final tuning of the assembled circuit—dedicated to BM design. Moreover, conventional approaches to the development of BM structures are unreliable and inefficient due to a large number of figures representing matrix performance (all of which have to be accounted for at a time) [27,28,30]. Although numerical optimization seems to be appropriate for addressing the difficulties pertinent to parameter adjustments, it suffers from high cost related to a large number of CPU-heavy evaluations required by the algorithm to converge. Consequently, maintaining low computational budget—e.g., using surrogate methods [40,41,42]—is of high importance in increasing the usefulness of algorithm-based approaches to BM design. The motivation of this work is to address the discussed challenges pertinent to design of modern BM circuits.
In this work, a bottom-up procedure for low-cost unattended design of miniaturized 4 × 4 Butler matrices with nonstandard output-port phase differences has been presented. The proposed method consists of two main steps involving synthesis and optimization of BM building blocks, followed by two-stage EM-driven tuning of the assembled matrix. The main contributions of the work include: (i) development of an automated methodology for sequential design of BM matrix components; (ii) determination of rigorously defined objective functions that facilitate handling the components/matrix design using robust numerical optimization algorithms; (iii) integration of all BM design stages (i.e., components synthesis, implementation, integration, and EM-driven tuning of assembled matrix) into a design framework; and (iv) validation of the presented method using two compact BMs with nonstandard output-port phases. The considered case studies involve design of structures featuring phase shifts of {−30°, 150°, −120°, 60°} and {−20°, 160°, −110°, 70°}, respectively. The footprints of the optimized circuits are only 696 and 767 mm2. To the best of the authors’ knowledge, this is the first work that thoroughly discusses the problem concerning automated BM design oriented for both miniaturization and enhancement of electrical performance. According to benchmark results, the computational cost of circuit design using the presented framework is over 80% lower compared to more conventional strategy based on direct EM-driven optimization of the assembled matrix. The optimized structures have been compared against state-of-the-art circuits from the literature in terms of size and performance. A comparison of simulation and measurement results is also provided.

2. Design Problem and Models of Butler Matrix Components

The automated design process presented here is governed by numerical optimization algorithms. This section contains formulation of the design problem and definition of models used by the proposed framework. For the sake of consistency, a brief discussion of numerical algorithms used in the work is also included. The details concerning mechanisms embedded into the proposed design framework are provided in Section 3.

2.1. Design Problem

The design problem can be formulated as a nonlinear minimization task of the form [40]
x * = arg min x U ( R ( x ) )
where R(x) is the response of the structure (be it the BM or its sub-circuit) under design. The vectors x and x* represent adjustable parameters of the structure and the optimal design to be found, whereas U denotes an objective function. The goal of (1) is to find x* through minimization of U. The latter “translates” the structure response into a scalar value, which is used by numerical algorithm to govern the design process toward the optimal solution. In other words, the function U represents the structure performance—calculated based on responses of the circuit—at the given design.

2.2. Models of Butler Matrix Components

In the proposed bottom-up design framework, the determination of the final BM geometry is preceded by synthesis, implementation, and optimization of its individual components. The responses of BM sub-circuits can be obtained through evaluation of the ideal transmission-line RT(xe), equivalent-circuit RC(xg), or electromagnetic RE(xg) models, respectively. Here, xe denotes the electrical parameters (characteristic impedance and electrical length of the transmission line), whereas xg is the vector of geometry variables. Conceptual illustrations of discussed structure representations are shown in Figure 1. It should be emphasized that the type of structure model matches the complexity of design task and size of the search space in terms of computational cost and accuracy. In other words, the model fidelity increases along with narrowing the search space to the region of interest. Here, the synthesis and correction/refinement of BM components at the level of their electrical parameters is realized using RT model, whereas RC and RE are applied for miniaturization-oriented development of building-blocks topologies and to obtain accurate responses for optimization and fine-tuning of the components.

2.3. Butler Matrix Representations

The Butler matrix models used here are defined as follows. Let RBe(y) and RBc(z) denote the high-fidelity EM model of the assembled BM (i.e., the one where all BM subcomponents are implemented and physically interconnected in the form of a single full-wave EM design) and the composite representation of the structure, respectively. The variables y and z represent their design parameters. The composite model responses are obtained using transmission-line theory from the S-parameter characteristics of individual BM subcircuits [43,44]. A conceptual illustration of the RBe and RBc models is shown in Figure 2. It should be noted that evaluation of the RBc model is to be preceded by simulations of its subcomponents (cf. Figure 2b). Nonetheless, the composite model is characterized by much lower evaluation cost compared to RBe. The vector z can be set equal to y but can also represent a composition of electrical and/or geometry variables of the sub-circuits selected for optimization (cf. Figure 2b). Example responses obtained from evaluation of RBe and RBc models for the z = y are shown in Figure 3. Despite reduced accuracy w.r.t. RBe (especially phase-wise), resulting from neglecting the coupling between adjacent components and loss at their interconnection, the composite model is useful for narrowing the search space to the region of interest before the fine-tuning stage.

2.4. Optimization Algorithms

As already mentioned, the proposed design method is governed by numerical optimization algorithms. The framework presented here exploits two routines: a trust-region-based gradient method and an unconstrained variant of the bisection algorithm. To make the work self-contained, both routines are briefly discussed below.

2.4.1. Trust-Region-Based Optimization

The main optimization engine is a gradient algorithm embedded in a trust-region framework. The method generates a series of approximations (i = 1, 2, 3, …) to the final design by solving [45]
x ( i + 1 ) = arg min x : x x ( i ) r ( i ) U ( G ( i ) ( x ) )
where G(i) is the first-order Taylor surrogate constructed from the S-parameter responses of the structure at hand. The model is given as [45]
G ( i ) ( x ) = R ( x ( i ) ) + J ( x ( i ) ) ( x x ( i ) )
Here, R(x(i)) is the response of the structure at hand. The perturbations for generation of the Jacobian J are obtained using a large-step finite differentiation [45]. Note that the model G(i) may be constructed using a combination of the responses determined from evaluations of the low- and high-fidelity representations of the circuit under design [46]. The parameter r(i) represents the trust-region radius, which is iteratively updated based on the calculated gain ratio, i.e., the obtained versus predicted change of the objective function. The radius is updated using standard rules [45]. A more detailed discussion of the algorithm can be found in [45,46].

2.4.2. Bisection-Based Heuristic

Another algorithm used in this work is a simple unconstrained variant of the bisection method [47]. Let x0.1 and x0.2 represent the starting points for the algorithm. Here, the vector x0.1 is obtained as a result of structure synthesis, whereas x0.2 represents a perturbation of all design parameters w.r.t. x0.1. The algorithm flow is as follows:
  • Set i = 0, x0.1, and x0.2; set Δx = |x0.1x0.2|; set x(0) = 0.5(x0.1 + x0.2).
  • Generate interval [x1(i), x2(i)] around the point x(i), where x1(i) = x(i) − αΔx and x2(i) = x(i) + αΔx.
  • If U(x1(i))U(x2(i)) ≥ 0, set x ( i + 1 ) = { x 1 ( i ) α Δ x ( i )     when   U ( x 1 ( i ) ) > 0 x 2 ( i ) + α Δ x ( i )     when   U ( x 1 ( i ) ) < 0 , i = i + 1 and go to step 2; otherwise go to step 4.
  • Set x(i+1) = 0.5(x1(i) + x2(i)), x1(i+1) = x1(i), x2(i+1) = x2(i), i = i + 1.
  • If U(x1(i))U(x(i)) < 0, set x2(i) = x(i) and go to step 7; otherwise go to step 6.
  • If U(x(i))U(x2(i)) < 0, set x1(i) = x(i) and go to step 7.
  • If |x(i)x(i−1)| ≤ ε, set x* = x(i) and END; otherwise go to step 4.
It should be noted that α and ε are user-defined parameters. Here, U(x) = U(R(x)) represents the objective function value calculated based on the model response at the design x. The discussed bisection-based algorithm is useful for approximating dimensions of the miniaturized BM components.

3. Methodology

The framework presented involves two main design stages: (i) synthesis and design of individual BM components and (ii) optimization of the composite BM model followed by fine-tuning of the assembled structure. Here, a detailed discussion of each design step and a summary of the presented methodology are provided. The numerical and experimental validation of the framework is considered in Section 4.3.

3.1. Synthesis of Butler Matrix and Its Components

A conceptual illustration of the considered 4 × 4 Butler matrix with nonstandard output-port phase differences is shown in Figure 4. The structure comprises two pairs of hybrid branch line couplers with adjustable phase, as well as two crossovers and four phase shifters. Excitation of the matrix through the selected input port Pj (j = 1, 2, 3, 4) allows for obtaining the phase differences Δθj between its output ports P5–8. The relation between output phase shifts and electrical lengths of the structure components is given as [23]
β 1 = 0.5 β 2 0.25 π
β 2 = 2 Δ θ 1
β 3 = 0.25 π
Here, β1 and β2 represent phase shifts introduced by the first and second pair of BLCs, whereas β3 is the electrical length of the phase shifter. As shown in Figure 4 and Figure 5a, the component βc represents the electrical length of crossovers. Based on (4)–(6), one can infer that phase differences at output ports of the BM are a function of β2. Consequently, Δθ1 = 0.5β2, Δθ2 = 0.5β2 + π, Δθ3 = 0.5β2 − 0.5π, and Δθ4 = 0.5β2 + 0.5π [23]. The feasible ranges of phase shifts β1 and β2—i.e., the ones for which realizable topologies of compact BLCs can be obtained—vary from −30° to −90°. Therefore, the attainable output-port phase differences are Δθ1 ∈ [−45°; −15°], Δθ2 ∈ [135°; 165°], Δθ3 ∈ [−135°; −105°], and Δθ4 ∈ [45°; 75°], respectively. The parameters β1–3 are used as the starting points for design of individual matrix components.
The illustration of the BLC structure capable of obtaining the desired phase and its comparison with conventional 90° hybrid is shown in Figure 5b,c. Given the phase difference βk (k = 1,2), kth coupler consists of a pair of 90° TL sections with characteristic impedance z1.k and equal electrical length of φ1.k, as well as a pair of TLs with impedance z2.k and electrical lengths φ2.k and φ3.k, respectively. Electrical parameters of the circuit can be obtained from the following equations [48]:
z 1 . k = z 2 . k 2 1 z 2 . k 2
z 2 . k = tan ( β k ) ( tan 2 ( 0.5 φ 2 . k ) 1 ) 2 tan ( 0.5 φ 2 . k )
φ 3 . k = 2 tan 1 ( 1 tan ( 0.5 φ 2 . k ) )
Here, (8) is solved for φ2.k with z2.k = 0.5 × 20.5 to provide closed-form estimation of BLC parameters:
φ 2 . k = π tan 1 ( 2 tan ( β k ) )
Although solving (7), (9), and (10) using the considered z2.k provides good estimation of the φ2.k and φ3.k w.r.t. the required phase shift, it does introduce power split error at the operating frequency. To address the problem, the ideal model of the coupler is optimized using algorithm of Section 2.4.1. The objective function is
U 1 = α 1 Δ C 2 + α 2 ( max ( M M 0 , 0 ) M 0 ) 2 + α 3 ( ( ϕ 1 ϕ 0.1 ) 2 + ( ϕ 2 ϕ 0.2 ) 2 )
Here, (11) is minimized based on the RT(xe) model responses, where xe = xe.k = [z1.k z2.k φ2.k φ3.k]T represents the vector of electrical parameters of kth BLC (note that φ1.k = 90°). The parameter ΔC = ||S31| − |S21|| denotes the power-split imbalance at the center frequency f0, whereas M = max(|S11|, |S41|) is an in-band performance of the BLC over the frequency range of interest defined around f0 [46]. The figures φ1 = ∠(S21/S31) and φ2 = ∠(S24/S34) are phase shifts at f0. The parameters M0 = −20 dB, φ0.1 = βk, and φ0.2 = βk − π represent target values. The scaling coefficients [α1 α2 α3] = [1000 500 1000] are determined so as to ensure balanced contribution of the design requirements to the aggregated objective function (11), [49,50]. In other words, they maintain similar relative importance of sub-elements in (11) during the optimization process. The selected values are appropriate for BM circuit components synthesized using (7)–(10).

3.2. Sequential Design of BM Components

The design of individual BM components is realized as a sequential process that involves determination of BLCs and crossovers geometries, as well as by tuning of the phase shifters. For the given center frequency f0, the crossover—see Figure 5a for illustration—can be considered as a “static” component of the matrix. In other words, its dimensions do not change when the BM is redesigned for another set of output phase differences. Consequently, the same design can be used to realize a range of Δθj. The crossover dimensions are adjusted in two steps. First, the selected geometry is optimized for minimization of U2.1 = max(|S11|, |S33|) at the center frequency f0, where |S11|, |S33| represent reflection of the crossed TLs. Next, the design is oriented toward ensuring that electrical lengths of the crossed transmission lines βc.1 and βc.2 are equal at f0. This is achieved by minimization of the objective function U2.2 = (βc.1βc.2)2. In each design step, the optimization is carried out using the algorithm of Section 2.4.1 and the RE(xg.c) model responses. Here, the vector xg.c represents geometry parameters of the crossover. The final design xg.c* structure provides equal length lines (βc = βc.1βc.2) with low reflection and high isolation levels (note that the term isolation refers to attenuation of the signal and is expressed as an absolute value of the transmission—in dB—between the selected pair of ports), all important for high BM performance.
The design stage involves development of compact BLCs. For each coupler, the initial design is synthesized as described in Section 3.1. The design of kth BLC (cf. Section 3.1) can be summarized as follows:
  • Decompose ideal BLC model to individual TLs.
  • Use electrical parameters of TLs as the reference for development of miniaturized BLC sections.
  • Optimize the BLC sections to match electrical parameters of the reference TLs.
  • Construct miniaturized BLC using the optimized cells and define the vector of its design parameters xg.k.
  • Optimize compact BLC using objective function (11) and algorithm of Section 2.4.1.
Note that, in step 3, each section of the miniaturized BLC must be optimized to ensure geometrical consistency of the structure, as well as to provide sufficient flexibility for the tuning of phase shifts. In steps 2, 4, and 5, the BLC optimization is carried out using the algorithm of Section 2.4.1, whereas in step 3 the routine of Section 2.4.2 is used. It should be emphasized that RE model evaluations are used only in the last stage of the BLC development, whereas models RT and RC are used in stage 1 and stages 2–4, respectively. Shifting the optimization burden to the simplified models is important for maintaining low cost of coupler development. For more detailed discussion on design of miniaturized BLCs, see [29,34,46,51]. The optimized high-fidelity BLC designs xg.1* and xg.2* are used as the starting point for BM tuning.
The final stage of the sequential design process involves development of phase shifters. Conventional BM structures comprise equal-length PSs. Here, however, unequal-length shifters are used to increase flexibility of the BM in terms of control over the output-port phase differences. The electrical parameters of PSs are optimized using the composite BM model. In this step, the model integrates EM responses of optimized BLCs and crossovers (which remain fixed in the optimization process). The phase shifters are implemented in the form of ideal TLs. The composite model is optimized to minimize the following function:
U 3 = α 1 ( max ( M B M B 0 , 0 ) M B 0 ) 2 + α 2 ( max ( Δ C B Δ C 0 , 0 ) Δ C 0 ) 2 +               + α 3 ( max ( P B P 0 , 0 ) P 0 ) 2 + α 4 ( max ( M B f 0 M f 0 , 0 ) M f 0 ) 2
where
Δ C B = 1 N 1 N | max ( Δ C E . j ) min ( Δ C E . j ) |
P B = 1 N 1 N | max ( P E . j ) min ( P E . j ) |
Here, MB and MBf0 represent the maximum value of the reflection and isolation responses between BM input ports (expressed in dB) within the frequency range of interest and at f0, respectively. The figure ΔCE.j denotes the magnitude (in dB) of transmission to the output ports at f0 when the structure is fed through the jth port (N = 4). Similarly, PE.j is a normalized phase difference at the output ports for excitation through jth port. The parameters MB0 = −15 dB, ΔC0 = 0.2 dB, P0 = 1.5°, and Mf0 = −30 dB represent the target values for matrix design. The weights [α1 α2 α3 α4] = [400 1 1 10] are determined based on numerical studies.
The optimized electrical lengths of PSs are used as the target for determination of each shifter’s physical dimensions. The phase shifters are implemented in the form of meandered TLs. The optimization of each PS is realized separately and is oriented toward minimization of (12). The initial dimensions of the meander lines are determined based on the transmission-line theory. The design objective for refinement of pth (p = 1, 2, 3, 4) phase shifter geometry is U4 = (e0.pep)2, where ep is the electrical length of the structure under design and e0.p represents the target length. Due to low evaluation cost and quick convergence of the algorithm (2), each meandered PS is implemented only in the form of high-fidelity model RE. The Butler matrix components determined in this stage of the design process are used for further optimization and tuning of the assembled structure.

3.3. Butler Matrix Optimization and Fine-Tuning

For the sake of low computational cost, the BM optimization is realized using the composite model RBc. Here, the parameters of all BLCs and PSs are enabled for adjustment. The optimization is oriented toward minimization of (12) using algorithm of Section 2.4.1. The main goal of this step is to further narrow down the search space to the region of interest so as to reduce the number of RBe evaluations required to find the final design. It should be reiterated that the evaluation cost of RBe is much higher (at least 3-fold) compared to the simulation cost of the composite structure.
Fine-tuning of the structure is again realized through minimization of (12). Here, the low cost of the process is maintained using a modified Taylor-expansion model (3), which exploits Jacobian J constructed based on simulations of the composite model, whereas the high-fidelity model simulations are performed only at the center design (i.e., R(x(i)) = RBe(y(i))) [46]. Consequently, each iteration of the tuning process requires only two EM simulations of each subcomponent (i.e., a total of four and eight simulations for couplers and phase shifters, respectively) and single evaluation of the assembled BM. The design y* obtained after the fine-tuning stage is the final solution of the presented design process.

3.4. Summary of the Design Framework

The proposed framework for design of miniaturized Butler matrices with unconventional phase differences can be summarized as follows:
  • Define the desired performance of the Butler matrix.
  • Perform synthesis of the BM structure using (4)–(6).
  • Perform fine-tuning of the crossover to match length of its TLs and extract βc (cf. Section 3.2).
  • Use (7), (9), (10) to synthesize the BLCs.
  • Adjust performance of synthesized BLCs through minimization of (11).
  • Set k = 1.
  • Perform topology development of kth miniaturized BLC and optimize its EM model.
  • If k = 2, go to step 9; otherwise set k = k + 1 and go to step 7.
  • Optimize ideal models of phase shifters through minimization of (12), set p = 1.
  • Generate initial dimensions of pth PS and optimize its EM model (cf. Section 3.2).
  • If p = 4, go to step 12; otherwise set p = p + 1 and go to step 10.
  • Optimize composite model of the BM.
  • Perform fine-tuning of the BM.
It should be noted that the design process described here is automated. Consequently, once the models of components and assembled BM are prepared, the role of the user is reduced only to definition of the performance requirements. The design bounds for constrained optimization stages (i.e., the ones governed by the TR algorithm) are defined as ±30% around the starting point for each step. The computational cost of the design process realized using the proposed method is comparable to around a dozen of RBe model simulations. Typically, each step that involves EM model evaluations requires no more than 10 iterations of the algorithm (2) to converge. The cost associated with synthesis of BM, BLCs, and PSs is negligible as it only requires evaluations of the transmission line or equivalent-circuit models. It is worth mentioning that conventional matrix with phase shifts of ±45° and ±135° is just a special case for the presented methodology. Consequently, the proposed framework is applicable to the design of eight port BMs, for which the requirements concern size reduction, performance enhancement, or combination of thereof. From this perspective, the methodology represents a generalized approach to design of 4 × 4 Butler matrices discussed in Section 3.

4. Numerical Results and Experiment

In this section, the proposed bottom-up design framework is demonstrated based on two examples of a compact 4 × 4 Butler matrices. Both structures are implemented on a dielectric substrate with εr = 3.48, h = 0.168 mm, and tan δ = 0.0037. The center frequency for the considered BMs is set to f0 = 2.6 GHz, whereas the desired operational bandwidth is from 2.5 to 2.7 GHz. The considered range covers sub-6 GHz bands used by the 5G technology. The presented framework has been benchmarked against conventional approach to BM design. Furthermore, the considered matrices have been compared against the state-of-the-art designs from the literature. The measurement results obtained for one of the matrices have also been included and discussed.
The presented framework is implemented in MATLAB. The latter controls synthesis of BM components, integration of circuits, optimization, as well as bidirectional communication with external simulation packages. Evaluations of the equivalent-circuit models are performed using Keysight ADS software, whereas simulations of the EM models are handled using CST Studio packages.

4.1. Butler Matrix 1

The first design example is a Butler matrix featuring an output-port phase shifts of Δθ1–4 = {−30°, 150°, −120°, 60°}. The electrical parameters of the structure components (cf. Section 3.1), i.e., β1 = −75°, β2 = −60°, and β3 = −45°, have been synthesized from (4)(6).
In order to maintain small BM dimensions, the crossover used here departs from the standard designs composed as a cascade connection of hybrid couplers [28,29,33]. Instead, it is based on a microstrip-to-coplanar-waveguide (CPW) transition. The design is characterized by small dimensions, as well as good isolation between crossed TLs [52]. The crossover topology is shown in Figure 6a. The structure introduces a meander to one of the transmission lines which can be adjusted to match electrical length of the crossed TLs. The design parameters of the component are xg.c = [l1 l2 w1 s1 s2 d1 d2 d3]T. The relative dimensions are w2 = 2d1 + r, w3 = 2d2 + r, l0.cross = 2l1 + w0 + 2(d1 + 2d2 + 3r + 2s2), l3 = w0 + 2(4r + 2s2 + 2d1 +2d2), l4 = 0.5w1 + s1 + 2r + 2d2 − 2s2, whereas r = 0.3. The parameter w0 = 0.35 to ensure 50 Ω input impedance. Note that the unit for all dimensions is “mm.” The initial design xg.c(0) = [1.5 0 0.3 0.2 0.2 0.1 0.15 0.6]T has been obtained in the course of structure development. The final design xg.c* = [1.5 0.2 0.49 0.2 0.2 0.1 0.2 0.59]T has been found through minimization of the objective function defined in Section 3.2. The optimized structure is characterized by a small footprint of only 6.05 mm × 6.05 mm = 36.6 mm2. Furthermore, it offers reflection below −30.5 dB, as well as transmission and isolation of around −0.03 and 39.5 dB (all at the center frequency), respectively. Owing to the use of meander lines, the difference of electrical length between crossed lines is kept below 0.02°, whereas the phase shift introduced by the structure amounts to βc = −61.2° (cf. Section 3.2). The magnitude and phase characteristics of the optimized component are shown in Figure 6b.
Next, the design of individual BLC structures has been performed. The initial parameters of the first BLC Z1.1 = 35.35 Ω, φ2.1 = 100.73°, φ3.1 = 79.27° and the second coupler Z1.2 = 35.35 Ω, φ2.2 = 112.21°, φ3.2 = 67.79° have been obtained from (7), (9), and (10). Note that Z1–2.k = Z0·z1–2.k (here, z2.k = 0.5·20.5 and φ1–2.k = 90°). The corrected parameters, obtained through optimization of ideal models (cf. Section 3.1), are Z1.1 = 44.65 Ω, Z2.1 = 32.16 Ω, φ2.1 = 100.56°, φ3.1 = 79.44° for BLC1 and Z1.2 = 42.59 Ω, Z2.2 = 32.20 Ω, φ2.2 = 110.71°, φ3.2 = 69.30° for BLC2. The obtained electrical parameters have been used as the initial designs for microstrip-line-based implementations of the miniaturized BLCs. The geometry of the compact BLC structure with highlight on its individual transmission line sections is shown in Figure 7. The vector of kth coupler design parameters is xg.k = [w1.k w2.k c1.k c2.k l1.k l2.k l3.k]T. The dimensions of each section have been recalculated from electrical parameters using the transmission line theory and adjusted using the bisection-based algorithm of Section 2.4.2 assuming c1.k = c2.k = 0.3 mm. Then, the obtained initial designs xg.1(0) = [0.39 0.58 0.3 0.3 2.44 2.29 1.55]T and xg.2(0) = [0.41 0.59 0.3 0.3 2.35 2.61 1.23]T of both BLCs have been optimized (at the EM model level) using the algorithm of Section 2.4.1. The final designs xg.1* = [0.35 0.53 0.42 0.4 2.4 2.25 1.52]T and xg.2* = [0.38 0.54 0.23 0.41 2.63 2.68 1.27]T have been obtained after seven and nine iterations of (2), respectively. Magnitude and phase responses of both BLCs at xg.k(0) and xg.k* are shown in Figure 8. The optimized designs are characterized by small footprints of 9 mm × 6.7 mm = 60.2 mm2 for BLC1 and 8.6 mm × 6.9 mm = 58.8 mm2 for BLC2, which corresponds to 79% miniaturization compared to couplers based on conventional TL sections [29].
The last step of sequential design involves optimization of phase shifters. The initial vector of adjustable electrical parameters is xe.p(0) = [β3 + βc β3 + βc βc βc]T = [106.2 106.2 61.2 61.2]T degrees (note that impedance of shifters is fixed to 50 Ω). The optimized values xe.p* = [105.9 106.7 60.9 61.3]T degrees are obtained through optimization of shifters embedded into the composite BM model as described in Section 3.2. The PS geometry is shown in Figure 9a. The adjustable parameter of pth structure (p = 1, 2, 3, 4) is xg.p = lp, whereas the relative dimensions that allow for maintaining geometrical consistency of the BM are lc = 0.2(l0.cross − 4w0), lp.2 = lp − Δl, and Δl = hBLC.2 − hBLC.1 (here, hBLC.k = 2w0 + l2.k + 2w2.k + c2.k + 4w1.k + 4c1.k + l3.k, k = 1, 2); the width of meander lines is fixed to w0. Note that for p = 3 and p = 4, Δl = 0 and thus lp.2 = lp. The optimized lengths xg.1–4 = l1–4 = {3.56, 3.6, 1.34, 1.36} mm of each meander section have been obtained sequentially as explained in Section 3.2.
The geometry of the assembled Butler matrix is shown in Figure 9b. The design parameters of the composite model RBc used for optimization are z = [xg.1 xg.2 l1 l2 l3 l4]T. The starting point is set to z(0) = [0.35 0.53 0.42 0.4 2.4 2.25 1.52 0.38 0.54 0.23 0.41 2.63 2.68 1.27 3.56 3.6 1.34 1.36]T. The optimal design z* = [0.33 0.53 0.43 0.42 2.36 2.21 1.51 0.4 0.58 0.22 0.4 2.62 2.66 1.23 3.7 3.69 1.36 1.36]T has been found after six iterations of (2). It should be emphasized that, owing to the use of composite model for optimization of z, each successful iteration of (2) required only two EM model evaluations of each BM component. The design z* has been used as the initial point for fine-tuning of the matrix. The final design resulting from the tuning process of the RBe model y* = [0.35 0.54 0.46 0.40 2.39 2.26 1.48 0.43 0.61 0.22 0.35 2.56 2.72 1.22 3.75 3.66 1.37 1.32]T has been obtained in eight iterations of (3). The dimensions of the optimized matrix are only 24 mm × 29 mm with overall footprint of 696 mm2. A comparison of the BM responses obtained through simulation of the RBe model at the designs z(0), z*, and y* is shown in Figure 10.
The simulation results indicate that optimization of the composite model is important for improving transmission responses of the matrix, whereas fine-tuning provides correction of the output-port phase differences. At the final design y*, the RBe model response features in-band matching and isolation above the level of 15.5 dB. Furthermore, at the center frequency it offers insertion loss imbalance below 0.5 dB and phase shift errors below 2.6°. It should be noted that although the optimized design slightly violates the target values of (12)—defined in Section 3.2—this is justified as the objective function comprises a composition of design requirements, balanced by the user-defined weighting factors. Table 1 provides more information on performance of the matrix in terms of reflection RBW, isolation IBW, transmission imbalance ΔMBW, and phase imbalance ΔPBW within 2.5 to 2.7 GHz band, as well as transmission ΔMf0 and phase ΔPf0 imbalances at the center frequency. The quantities RBW and IBW refer to the worst-case in-band performance (across all considered responses), whereas imbalance represents maximum difference between the considered groups of characteristics, either within the band of interest or at the center frequency.
The proposed design procedure has been benchmarked against the method where the fine-tuning step governed by algorithm (2) involves only evaluations of the RBe model [45]. For fair comparison, it is assumed that the design z(0)—obtained through individual optimization of BM components—is used as a starting point for adjustment of the structure topology. The assumption seems justified considering that it follows the industry-wide divide-and-conquer strategy to design of complex circuits [27,36,53,54]. The tests have been performed on an Intel Xeon machine with 32 GB RAM. For the given machine, the average cost of EM simulation of BLC, PS, and assembled BM amounts to 4.5 min, 36 s, and 41 min, respectively. The results shown in Table 2 indicate that the proposed procedure is capable of yielding a design with competitive performance compared to the method that does not blend RBe and RBc models, yet at a fraction of its computational cost (13.2 h for the presented approach vs. 91.6 h for conventional method). The number of model evaluations in Table 2 refers to all simulations of particular models required to find the final design.

4.2. Butler Matrix 2

The second design example is the Butler matrix structure featuring phase shifts Δθ1–4 = {−20°, 160°, −110°, 70°}. The electrical properties of the matrix components required to obtain desired Δθj values are β1 = −65°, β2 = −40°, and β3 = −45° (cf. Section 3.1). The BM topology is shown in Figure 9b. The crossover dimensions remain the same as for the one used in Section 4.1. The initial dimensions of both BLCs xg.1(0) = [0.41 0.59 0.3 0.3 2.31 2.48 1.39]T and xg.2(0) = [0.59 0.72 0.3 0.3 1.5 2.79 0.54]T have been determined as described in Section 3.2 and Section 4.1 The optimized dimensions xg.1* = [0.37 0.53 0.36 0.26 3.06 2.89 1.6]T and xg.2* = [0.57 0.67 0.33 0.34 1.48 2.95 0.58]T have been obtained after seven and nine iterations of (2), respectively.
The initial design for Butler matrix optimization (composite model) is z(0) = [0.37 0.53 0.36 0.26 3.06 2.89 1.6 0.57 0.67 0.33 0.34 1.48 2.95 0.58 3.56 3.6 1.34 1.36]T. The optimized design (used as a starting point for fine-tuning of the BM) z* = [0.36 0.52 0.38 0.26 3.09 2.89 1.59 0.58 0.69 0.33 0.34 1.49 2.95 0.59 3.57 3.6 1.35 1.37]T has been obtained after seven iterations of (2). The final design y* = [0.35 0.51 0.42 0.25 3.11 2.90 1.56 0.64 0.73 0.34 0.27 1.53 2.95 0.52 3.57 3.56 1.39 1.36]T has been found after nine iterations of (2). The optimized BM is characterized by the dimensions of 25.9 mm × 29.6 mm and overall footprint of only 766.6 mm2. The response characteristics of the matrix are shown in Figure 11. The BM features in-band reflection below −14 dB and isolation higher than 20 dB, respectively. Moreover, at f0 the structure offers insertion loss imbalance below 0.55 dB and phase shift imbalance below ±2.5° w.r.t. the target values. It should be noted that the slightly worsened in-band reflection compared to the structure of Section 4.1 results from the narrow bandwidth of the BLC required for the realization of the −45° phase shift. The performance characteristics of the optimized structure are gathered in Table 3 (for explanation of used performance metrics, see Section 4.2).

4.3. Numerical and Experimental Validation

Simulation-based validation of the optimization results has been performed through evaluation of the optimized designs using a solver based on the method of moments (Keysight ADS). A comparison of the performance characteristics obtained using CST Studio and Keysight ADS is shown in Figure 12 and Figure 13, respectively. The simulation results for BMs feature reflection below −12 dB and isolation above 13 dB within the bandwidth of interest. Regardless of the selected excitation port, the first design offers the insertion-loss imbalance below 0.9 dB and phase shift error below 8°. For the second matrix, the values are 1.5 dB and 9°, respectively. The obtained results indicate that the EM simulation results are valid.
For further verification of the results, the structure of Section 4.2 has been fabricated and measured. The photograph of a manufactured prototype is shown in Figure 14. To facilitate the measurement process, feed lines of the structure have been modified to make space for connectors. A comparison of reflection/isolation characteristics obtained from simulations and measurements is shown in Figure 15. The misalignment between the responses (most notably the frequency shift) is the consequence of the systematic error caused by differences between the relative permittivity and thickness of the substrate used for simulations and measurements. As can be seen from Figure 16, accounting for systematic errors significantly improves alignment between the responses. Other factors affecting the discrepancies include manufacturing tolerances, as well as assembly- (manual positioning/soldering of the connectors) and measurement-related errors (also resulting from tolerances of the matching loads) [55]. It should be emphasized that connectors alter electrical properties of the measured BM, which is mostly due to discontinuities at their interface with microstrip lines. The effect contributes to degradation of measured performance as compared to simulations. Nevertheless, the agreement between the simulations and measurements can be considered acceptable. A summary of the measured structure performance is provided in Table 4. In many cases, resemblance between EM simulations and measurements can be further increased by excluding the effects of the fixture (here, connectors along with introduced feeding lines) on the device under test [56]. Unfortunately, high sensitivity to precision of the assembly makes the method unsuitable for the discussed setup.

4.4. Comparison with Benchmark Structures

The optimized designs from Section 4.1 and Section 4.2 have been compared in terms of size and performance with other planar BMs from the literature [16,21,23,27,28,29,30]. Whenever possible (all benchmark designs except the one reported in [28]), the performance figures have been obtained based on the EM simulation results. The considered figures include magnitude imbalance ΔMf0 and phase shift ΔPf0 at the center frequency, as well as bandwidth BW. The latter is expressed in percent and calculated as a ratio of the difference between the upper and lower frequencies (i.e., the ones for which isolation is above 15 dB and reflection below −15 dB) to the specified f0. All of the considered figures represent the worst-case scenario for a series of analyses w.r.t. all input ports. For fair comparison of size, the dimensions of all matrices are expressed in terms of a guided wavelength λg (defined for the given center frequency and electrical parameters of the substrate used to implement the circuit). The results shown in Table 5 indicate that the designs obtained using the proposed procedure are characterized by competitive performance (with particular emphasis on unconventional output-port phase differences). Moreover, at roughly 10-fold smaller size w.r.t. conventional BM, the obtained designs outperform other structures in terms of miniaturization. It should be emphasized that compact dimensions, competitive magnitude/phase performance, and the center frequency of around 2.6 GHz make the optimized BMs of potential use for IoT devices interconnected through a 5G backbone [57].

5. Conclusions

In this work, a bottom-up framework for low-cost automated design of 4 × 4 Butler matrices with nonstandard output-port phase shifts has been presented. The technique involves sequential development of BM components followed by optimization of the composite structure representation and fine-tuning of the assembled matrix. Each design step is governed by optimization algorithm. The proposed approach has been demonstrated using two compact BM structures designed to operate at 2.6 GHz frequency. The first circuit, designed to provide phase shifts of {−30°, 150°, −120°, 60°}, offers simulation-based reflection below −15 dB within the 2.5–2.7 GHz band, along with low transmission and phase shift imbalance of 0.5 dB and ±2.5°. The second structure realizes phase shifts of {−20°, 160°, −110°, 70°} with in-band reflection below −14 dB, as well as phase and magnitude imbalance below 0.55 dB and ±2.5°, respectively.
The computational cost of BM design using the proposed strategy is over 80% lower compared to more conventional approach that does not exploit composite models. The structures have been compared against state-of-the-art BMs from the literature in terms of performance and size. With the footprints of only 696 and 767 mm2, the optimized circuits outclassed BMs from the literature in terms of size reduction while providing competitive performance. The unique feature of the presented designs—which makes them of potential use for IoT systems interconnected through the 5G communication network—is that they address requirements concerning both size reduction and enhanced performance, whereas the benchmark designs address only one of these figures at a time. The measurements of the fabricated BM prototype are also provided.
Future work will focus on increasing the flexibility of BMs in terms of attainable output-port phase shift as well as development of design methods that support algorithm-driven integration of feeding networks and radiators into a complete antenna array.

Author Contributions

A.B. realized the conceptual work, developed the design method, implemented the framework, and performed numerical tests. S.K. performed measurements and provided interpretation of the obtained results. A.B. and S.K. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Science Centre of Poland Grant 2017/27/B/ST7/00563 and by the National Centre for Research and Development Grant NOR/POLNOR/HAPADS/0049/2019-00.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Song, C.; Pan, L.; Jiao, Y.; Jia, J. A High-Performance Transmitarray Antenna with Thin Metasurface for 5G Communication Based on PSO (Particle Swarm Optimization). Sensors 2020, 20, 4460. [Google Scholar] [CrossRef] [PubMed]
  2. Naqvi, S.I.; Hussain, N.; Iqbal, A.; Rahman, M.; Forsat, M.; Mirjavadi, S.S.; Amin, Y. Integrated LTE and Millimeter-Wave 5G MIMO Antenna System for 4G/5G Wireless Terminals. Sensors 2020, 20, 3926. [Google Scholar] [CrossRef] [PubMed]
  3. Lian, R.; Wang, Z.; Yin, Y.; Wu, J.; Song, X. Design of a Low-Profile Dual-Polarized Stepped Slot Antenna Array for Base Station. IEEE Antennas Wirel. Propag. Lett. 2016, 15, 362–365. [Google Scholar] [CrossRef]
  4. Alhalabi, R.A.; Rebeiz, G.M. High-Efficiency Angled-Dipole Antennas for Millimeter-Wave Phased Array Applications. IEEE Trans. Antennas Propag. 2008, 56, 3136–3142. [Google Scholar] [CrossRef]
  5. Rafii, V.; Nourinia, J.; Ghobadi, C.; Pourahmadazar, J.; Virdee, B.S. Broadband Circularly Polarized Slot Antenna Array Using Sequentially Rotated Technique for CC-Band Applications. IEEE Antennas Wirel. Propag. Lett. 2013, 12, 128–131. [Google Scholar] [CrossRef]
  6. Ojaroudiparchin, N.; Shen, M.; Zhang, S.; Pedersen, G.F. A Switchable 3-D-Coverage-Phased Array Antenna Package for 5G Mobile Terminals. IEEE Antennas Wirel. Propag. Lett. 2016, 15, 1747–1750. [Google Scholar] [CrossRef] [Green Version]
  7. Hansen, R.C. Phased Array Antennas, 2nd ed.; Wiley: Hoboken, NJ, USA, 2009. [Google Scholar]
  8. Ogurtsov, S.; Koziel, S. Systematic Approach To Sidelobe Reduction In Linear Antenna Arrays Through Corporate-Feed-Controlled Excitation. IET Microw. Antennas Propag. 2017, 11, 779–786. [Google Scholar] [CrossRef]
  9. Ogurtsov, S.; Koziel, S. On Alternative Approaches to Design of Corporate Feeds for Low-Sidelobe Microstrip Linear Arrays. IEEE Trans. Antennas Propag. 2018, 66, 3781–3786. [Google Scholar] [CrossRef]
  10. Ogurtsov, S.; Koziel, S. A Conformal Circularly Polarized Series-Fed Microstrip Antenna Array Design. IEEE Trans. Antennas Propag. 2019, 68, 873–881. [Google Scholar] [CrossRef]
  11. Chu, H.; Li, P.; Guo, Y.-X. A Beam-Shaping Feeding Network in Series Configuration for Antenna Array With Cosecant-Square Pattern and Low Sidelobes. IEEE Antennas Wirel. Propag. Lett. 2019, 18, 742–746. [Google Scholar] [CrossRef]
  12. Slomian, I.; Piekarz, I.; Wincza, K.; Gruszczynski, S. Microstrip Antenna Array With Series Feeding Network Designed With the Use of Slot-Coupled Three-Way Power Divider. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 667–670. [Google Scholar] [CrossRef]
  13. Wincza, K.; Gruszczynski, S. Microstrip Antenna Arrays Fed by a Series-Parallel Slot-Coupled Feeding Network. IEEE Antennas Wirel. Propag. Lett. 2011, 10, 991–994. [Google Scholar] [CrossRef]
  14. Al Abbas, E.; Abbosh, A.M.; Bialkowski, K. Tunable In-Phase Power Divider for 5G Cellular Networks. IEEE Microw. Wirel. Compon. Lett. 2017, 27, 551–553. [Google Scholar] [CrossRef]
  15. Huang, Y.; Wu, K.-L.; Fang, D.-G.; Ehlert, M. An Integrated LTCC Millimeter-Wave Planar Array Antenna With Low-Loss Feeding Network. IEEE Trans. Antennas Propag. 2005, 53, 1232–1234. [Google Scholar] [CrossRef] [Green Version]
  16. Kim, S.; Yoon, S.; Lee, Y.; Shin, H. A Miniaturized Butler Matrix Based Switched Beamforming Antenna System in a Two-Layer Hybrid Stackup Substrate for 5G Applications. Electronics 2019, 8, 1232. [Google Scholar] [CrossRef] [Green Version]
  17. Park, Y.; Bang, J.; Choi, J. Dual-Circularly Polarized 60 GHz Beam-Steerable Antenna Array with 8 × 8 Butler Matrix. Appl. Sci. 2020, 10, 2413. [Google Scholar] [CrossRef] [Green Version]
  18. Kim, Y.-J.; Kim, Y.-B.; Dong, H.-J.; Cho, Y.S.; Lee, H.L. Compact Switched-Beam Array Antenna with a Butler Matrix and a Folded Ground Structure. Electronics 2019, 9, 2. [Google Scholar] [CrossRef] [Green Version]
  19. Tajik, A. Shafiei Alavijehm A.; Fakharzadeh, M. Asymmetrical 4 × 4 Butler Matrix and Its Application For Single Layer 8 x 8 Butler Matrix. IEEE Trans. Ant. Prop. 2019, 67, 5372–5379. [Google Scholar] [CrossRef]
  20. Liu, C.; Xiao, S.; Guo, Y.-X.; Tang, M.-C.; Bai, Y.-Y.; Wang, B.-Z. Circularly Polarized Beam-Steering Antenna Array With Butler Matrix Network. IEEE Antennas Wirel. Propag. Lett. 2011, 10, 1278–1281. [Google Scholar] [CrossRef]
  21. Balanis, C.A. Modern Antenna Handbook; John Wiley & Sons: New York, NY, USA, 2008. [Google Scholar]
  22. Tseng, C.-H.; Chen, C.-J.; Chu, T.-H. A Low-Cost 60-GHz Switched-Beam Patch Antenna Array with Butler Matrix Network. IEEE Antennas Wirel. Propag. Lett. 2008, 7, 432–435. [Google Scholar] [CrossRef]
  23. Ren, H.; Arigong, B.; Zhou, M.; Ding, J.; Zhang, H. A Novel Design of 4×4 Butler Matrix with Relatively Flexible Phase Differences. IEEE Antennas Wirel. Propag. Lett. 2015, 15, 1277–1280. [Google Scholar] [CrossRef]
  24. Moody, H. The Systematic Design of The Butler Matrix. IRE Trans. Antennas Propag. 2004, 12, 786–788. [Google Scholar] [CrossRef]
  25. Der, E.T.; Jones, T.R.; Daneshmand, M. Miniaturized 4 × 4 Butler Matrix and Tunable Phase Shifter Using Ridged Half-Mode Substrate Integrated Waveguide. IEEE Trans. Microw. Theory Tech. 2020, 68, 3379–3388. [Google Scholar] [CrossRef]
  26. Ali, A.A.M.; Fonseca, N.J.G.; Coccetti, F.; Aubert, H. Design and Implementation of Two-Layer Compact Wideband Butler Matrices in SIW Technology for Ku-Band Applications. IEEE Trans. Antennas Propag. 2010, 59, 503–512. [Google Scholar] [CrossRef]
  27. Tian, G.; Yang, J.-P.; Wu, W. A Novel Compact Butler Matrix Without Phase Shifter. IEEE Microw. Wirel. Compon. Lett. 2014, 24, 306–308. [Google Scholar] [CrossRef]
  28. Xu, H.-X.; Wang, G.-M.; Wang, X. Compact Butler Matrix Using Composite Right/Left Handed Transmission Line. Electron. Lett. 2011, 47, 1081. [Google Scholar] [CrossRef]
  29. Koziel, S.; Kurgan, P. Low-Cost Optimization Of Compact Branch-Line Couplers And Its Application To Miniaturized Butler Matrix Design. In Proceedings of the 44th European Microwave Conference, Rome, Italy, 6–9 October 2014; pp. 227–230. [Google Scholar] [CrossRef]
  30. Letavin, D.A.; Shabunin, S.N. Construction of a 4x4 Butler Microstrip Matrix with Reduced Dimensions. In Proceedings of the 2020 Ural Symposium on Biomedical Engineering, Radioelectronics and Information Technology (USBEREIT), Yekaterinburg, Russia, 14–15 May 2020; pp. 257–260. [Google Scholar] [CrossRef]
  31. Karamzadeh, S.; Rafii, V.; Kartal, M.; Virdee, B.S. Compact and Broadband 4 × 4 SIW Butler Matrix With Phase and Magnitude Error Reduction. IEEE Microw. Wirel. Compon. Lett. 2015, 25, 772–774. [Google Scholar] [CrossRef]
  32. Lian, J.-W.; Ban, Y.-L.; Xiao, C.; Yu, Z.-F. Compact Substrate-Integrated 4 × 8 Butler Matrix With Sidelobe Suppression for Millimeter-Wave Multibeam Application. IEEE Antennas Wirel. Propag. Lett. 2018, 17, 928–932. [Google Scholar] [CrossRef]
  33. Bhowmik, P.; Moyra, T. Modelling and Validation of a Compact Planar Butler Matrix by Removing Crossover. Wirel. Pers. Commun. 2017, 27, 788–5132. [Google Scholar] [CrossRef]
  34. Tseng, C.-H.; Chang, C.-L. A Rigorous Design Methodology for Compact Planar Branch-Line and Rat-Race Couplers With Asymmetrical T-Structures. IEEE Trans. Microw. Theory Tech. 2012, 60, 2085–2092. [Google Scholar] [CrossRef]
  35. Ren, H.; Zhou, M.; Zhang, H.; Arigong, B. A Novel Dual-Band Zero-Phase True Crossover With Arbitrary Port Impedance. IEEE Microw. Wirel. Components Lett. 2019, 29, 29–31. [Google Scholar] [CrossRef]
  36. Chu, H.N.; Ma, T.-G. An Extended 4×4 Butler Matrix With Enhanced Beam Controllability and Widened Spatial Coverage. IEEE Trans. Microw. Theory Tech. 2018, 66, 1301–1311. [Google Scholar] [CrossRef]
  37. Tseng, C.-H.; Chen, H.-J. Compact Rat-Race Coupler Using Shunt-Stub-Based Artificial Transmission Lines. IEEE Microw. Wirel. Compon. Lett. 2008, 18, 734–736. [Google Scholar] [CrossRef]
  38. Ahmed, O.; Sebak, A.; Denidni, T. Compact UWB Printed Monopole Loaded With Dielectric Resonator Antenna. Electron. Lett. 2011, 47, 7–8. [Google Scholar] [CrossRef]
  39. Lin, C.-C. Compact Bow-Tie Quasi-Self-Complementary Antenna for UWB Applications. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 987–989. [Google Scholar] [CrossRef]
  40. Koziel, S.; Ciaurri, D.E.; Leifsson, L. Surrogate-Based Methods. In Computational Optimization, Methods and Algorithms. Studies in Computational Intelligence; Koziel, S., Yang, X.S., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; Volume 356, pp. 33–59. [Google Scholar]
  41. Leal-Romo, F.D.J.; Rayas-Sanchez, J.E.; Chavez-Hurtado, J.L. Surrogate-Based Analysis and Design Optimization of Power Delivery Networks. IEEE Trans. Electromagn. Compat. 2020, 62, 2528–2537. [Google Scholar] [CrossRef]
  42. Feng, F.; Na, W.; Liu, W.; Yan, S.; Zhu, L.; Ma, J.; Zhang, Q.-J. Multifeature-Assisted Neuro-transfer Function Surrogate-Based EM Optimization Exploiting Trust-Region Algorithms for Microwave Filter Design. IEEE Trans. Microw. Theory Tech. 2020, 68, 531–542. [Google Scholar] [CrossRef]
  43. Pozar, D.M. Microwave Engineering, 4th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  44. Hong, J.; Lancaster, M.J. Microstrip Filters for RF/Microwave Applications; Wiley: New York, NY, USA, 2001. [Google Scholar]
  45. Conn, A.R.; Gould, N.I.M.; Toint, P.L. Trust Region Methods; Society for Industrial & Applied Mathematics (SIAM): Philadelphia, PA, USA, 2000. [Google Scholar]
  46. Bekasiewicz, A. Low-Cost Automated Design of Compact Branch-Line Couplers. Sensors 2020, 20, 3562. [Google Scholar] [CrossRef]
  47. Arfken, G. Mathematical Methods for Physicists, 3rd ed.; Academic Press: Orlando, FL, USA, 1985. [Google Scholar]
  48. Wong, Y.S.; Zheng, S.Y.; Chan, W.S. Quasi-arbitrary phase-difference hybrid coupler. IEEE Trans. Microw. Theory Tech. 2012, 60, 1530–1539. [Google Scholar] [CrossRef]
  49. Deb, K. Multi-Objective Optimization Using Evolutionary Algorithms; John Wiley Sons, Ltd.: Hoboken, NJ, USA, 2001; p. 497. [Google Scholar]
  50. Coello, C.A.C.; Van Veldhuizen, D.A.; Lamont, G.B. Evolutionary Algorithms for Solving Multi-Objective Problems; Springer Science and Business Media LLC: New York, NY, USA, 2002. [Google Scholar]
  51. Kurgan, P.; Filipcewicz, J.; Kitlinski, M. Development of a Compact Microstrip Resonant Cell Aimed at Efficient Microwave Component Size Reduction. IET Microw. Antennas Propag. 2012, 6, 1291. [Google Scholar] [CrossRef]
  52. Tajik, A.; Fakharzadeh, M.; Mehrany, K. DC to 40-GHz Compact Single-Layer Crossover. IEEE Microw. Wirel. Compon. Lett. 2018, 28, 642–644. [Google Scholar] [CrossRef]
  53. Ding, K.; Kishk, A.A. 2-D Butler Matrix and Phase-Shifter Group. IEEE Trans. Microw. Theory Tech. 2018, 66, 5554–5562. [Google Scholar] [CrossRef]
  54. Hassan, A.T.; Kishk, A.A. Efficient Procedure to Design Large Finite Array and Its Feeding Network With Examples of ME-Dipole Array and Microstrip Ridge Gap Waveguide Feed. IEEE Trans. Antennas Propag. 2020, 68, 4560–4570. [Google Scholar] [CrossRef]
  55. Mroczka, J. The Cognitive Process in Metrology. Measurement 2013, 46, 2896–2907. [Google Scholar] [CrossRef]
  56. Chen, B.; He, J.; Guo, Y.; Pan, S.; Ye, X.; Fan, J. Multi-Ports ([2n) 2×-Thru De-Embedding: Theory, Validation, and Mode Conversion Characterization. IEEE Trans. Electromagn. Compat. 2019, 61, 1261–1270. [Google Scholar] [CrossRef]
  57. Ji, B.; Chen, Z.; Chen, S.; Zhou, B.; Li, C.; Wen, H. Joint Optimization for Ambient Backscatter Communication System With Energy Harvesting for IoT. Mech. Syst. Signal Process. 2020, 135, 106412. [Google Scholar] [CrossRef]
Figure 1. Illustration of structure representations (here, a phase shifter) supported by the design method: (a) ideal transmission line (TL) model RT(xe) with xe = [Z θ]T, (b) equivalent-circuit model RC(xg), and (c) high-fidelity electromagnetic (EM) model RE(xg). The vector of geometry parameters xg = [w l1 l2]T is the same for RC and RE models.
Figure 1. Illustration of structure representations (here, a phase shifter) supported by the design method: (a) ideal transmission line (TL) model RT(xe) with xe = [Z θ]T, (b) equivalent-circuit model RC(xg), and (c) high-fidelity electromagnetic (EM) model RE(xg). The vector of geometry parameters xg = [w l1 l2]T is the same for RC and RE models.
Sensors 21 00851 g001
Figure 2. A conceptual illustration of the Butler matrix (BM) models: (a) high-fidelity EM model of the assembled matrix where (i)–(iii) denote branch line coupler (BLC), crossover, and phase-shifter, respectively, and (b) a composite representation of the structure. Frequency characteristics of the composite BM are obtained from the responses of individual components. In this illustration, RBc is calculated using EM model responses of crossover (design xg.c) and BLCs (two designs: xg.1 and xg.2), as well as ideal TL model of phase shifters (two designs: xe.1 and xe.2).
Figure 2. A conceptual illustration of the Butler matrix (BM) models: (a) high-fidelity EM model of the assembled matrix where (i)–(iii) denote branch line coupler (BLC), crossover, and phase-shifter, respectively, and (b) a composite representation of the structure. Frequency characteristics of the composite BM are obtained from the responses of individual components. In this illustration, RBc is calculated using EM model responses of crossover (design xg.c) and BLCs (two designs: xg.1 and xg.2), as well as ideal TL model of phase shifters (two designs: xe.1 and xe.2).
Sensors 21 00851 g002aSensors 21 00851 g002b
Figure 3. Frequency characteristics of the composite (solid gray) and assembled (dashed black) BM models at some design: (a) reflection and isolation, (b) transmission, and (c) phase difference between output ports. The responses are obtained for excitation of both structures through the first input port.
Figure 3. Frequency characteristics of the composite (solid gray) and assembled (dashed black) BM models at some design: (a) reflection and isolation, (b) transmission, and (c) phase difference between output ports. The responses are obtained for excitation of both structures through the first input port.
Sensors 21 00851 g003
Figure 4. Butler matrix with nonstandard output-port phase differences: a conceptual illustration.
Figure 4. Butler matrix with nonstandard output-port phase differences: a conceptual illustration.
Sensors 21 00851 g004
Figure 5. Crossover and hybrid BLCs configurations: (a) crossover structure with highlight on phase shifts between crossed TLs (desired phase shift of the crossover is βc = βc.1 = βc.2), (b) conventional BLC with z1 = 0.707, z2 = 1, and φ = 90° (phase shift: β = 90°), and (c) modified BLC circuit with extra degrees of freedom φ2.k, φ3.k w.r.t. conventional design that enables adjustment of phase shift βk (the parameters z1 and z2 are normalized w.r.t. Z0 = 50 Ω).
Figure 5. Crossover and hybrid BLCs configurations: (a) crossover structure with highlight on phase shifts between crossed TLs (desired phase shift of the crossover is βc = βc.1 = βc.2), (b) conventional BLC with z1 = 0.707, z2 = 1, and φ = 90° (phase shift: β = 90°), and (c) modified BLC circuit with extra degrees of freedom φ2.k, φ3.k w.r.t. conventional design that enables adjustment of phase shift βk (the parameters z1 and z2 are normalized w.r.t. Z0 = 50 Ω).
Sensors 21 00851 g005
Figure 6. Compact crossover: (a) geometry of the structure, as well as (b) magnitude and phase characteristics of the circuit obtained based on EM simulations. Note that the electrical lengths of crossed lines (solid gray and dashed black) are virtually the same.
Figure 6. Compact crossover: (a) geometry of the structure, as well as (b) magnitude and phase characteristics of the circuit obtained based on EM simulations. Note that the electrical lengths of crossed lines (solid gray and dashed black) are virtually the same.
Sensors 21 00851 g006
Figure 7. Compact BLC with flexible phase difference: (a) individual sections of the structure and (b) geometry of the assembled circuit.
Figure 7. Compact BLC with flexible phase difference: (a) individual sections of the structure and (b) geometry of the assembled circuit.
Sensors 21 00851 g007
Figure 8. EM model responses of miniaturized BLC structures – magnitude (gray) and phase (black) at the initial (····) and optimized (––) designs: (a) BLC1 (realized phase difference β1 = −75°) and (b) BLC2 (realized phase difference β2 = −60°).
Figure 8. EM model responses of miniaturized BLC structures – magnitude (gray) and phase (black) at the initial (····) and optimized (––) designs: (a) BLC1 (realized phase difference β1 = −75°) and (b) BLC2 (realized phase difference β2 = −60°).
Sensors 21 00851 g008
Figure 9. Geometry of (a) phase shifter with highlight on design parameters and (b) assembled structure of a compact Butler matrix.
Figure 9. Geometry of (a) phase shifter with highlight on design parameters and (b) assembled structure of a compact Butler matrix.
Sensors 21 00851 g009
Figure 10. Butler matrix 1: structure characteristics at different stages of the design process. The responses are obtained from evaluations of the composite model RBe at (a) the initial design z(0), (b) the final design z* (used as a starting point for fine-tuning), and (c) the optimized design y* found based on simulations of the assembled BM model RBe. Each row shows responses of the structure “seen” from different input port (P1…4). Gray and black lines represent magnitude and phase responses, respectively.
Figure 10. Butler matrix 1: structure characteristics at different stages of the design process. The responses are obtained from evaluations of the composite model RBe at (a) the initial design z(0), (b) the final design z* (used as a starting point for fine-tuning), and (c) the optimized design y* found based on simulations of the assembled BM model RBe. Each row shows responses of the structure “seen” from different input port (P1…4). Gray and black lines represent magnitude and phase responses, respectively.
Sensors 21 00851 g010
Figure 11. Butler matrix 2: EM-model-based magnitude (gray) and phase (black) responses of the optimized structure.
Figure 11. Butler matrix 2: EM-model-based magnitude (gray) and phase (black) responses of the optimized structure.
Sensors 21 00851 g011
Figure 12. Butler matrix 1: comparison of the EM model responses obtained using CST Studio (gray) and Keysight ADS (black). Solid and dotted lines represent magnitude and phase responses, respectively.
Figure 12. Butler matrix 1: comparison of the EM model responses obtained using CST Studio (gray) and Keysight ADS (black). Solid and dotted lines represent magnitude and phase responses, respectively.
Sensors 21 00851 g012
Figure 13. Butler matrix 2: comparison of the EM model responses obtained using CST Studio (gray) and Keysight ADS (black). Solid and dotted lines represent magnitude and phase responses, respectively.
Figure 13. Butler matrix 2: comparison of the EM model responses obtained using CST Studio (gray) and Keysight ADS (black). Solid and dotted lines represent magnitude and phase responses, respectively.
Sensors 21 00851 g013
Figure 14. BM of Section 4.2: photograph of the fabricated prototype.
Figure 14. BM of Section 4.2: photograph of the fabricated prototype.
Sensors 21 00851 g014
Figure 15. Butler matrix 2: comparison of measurements (––) with EM simulations (···) in terms of reflection/isolation. Here, the EM simulation model does not account for systematic errors.
Figure 15. Butler matrix 2: comparison of measurements (––) with EM simulations (···) in terms of reflection/isolation. Here, the EM simulation model does not account for systematic errors.
Sensors 21 00851 g015
Figure 16. Measurements (––) vs. simulations (···) obtained using the EM model that account for the systematic errors: (a) reflection/isolation, (b) transmission, and (c) output-port phase differences.
Figure 16. Measurements (––) vs. simulations (···) obtained using the EM model that account for the systematic errors: (a) reflection/isolation, (b) transmission, and (c) output-port phase differences.
Sensors 21 00851 g016aSensors 21 00851 g016b
Table 1. Butler matrix 1: performance figures of the optimized structure (EM simulation).
Table 1. Butler matrix 1: performance figures of the optimized structure (EM simulation).
ExcitationPort P1Port P2Port P3Port P4
RBW (dB)−21.2−15.5−17.3−16.1
IBW (dB)18.417.617.618.1
ΔMBW (dB)0.790.750.740.45
ΔMf0 (dB)0.510.450.410.14
Target phase (°)−30150−12060
ΔPBW (°)±4.9±4.7±4.0±4.3
ΔPf0 (°)±1.6±1.3±2.5±2.6
Table 2. Butler matrix 1: numerical benchmark of the proposed design framework.
Table 2. Butler matrix 1: numerical benchmark of the proposed design framework.
MethodEM Model TypeNo. of Model EvaluationsCostTotal CostPerformance Figures *
(RBe)(h)(RBe)(h)RBW (dB)ΔMf0 (dB)ΔPf0 (°)
DirectRBe13413491.613491.6−15.20.59±3.2
This workBLC #606.64.519.413.2−15.50.51±2.6
PS #1201.81.2
RBe11117.5
* Worst-case performance obtained for excitation of the BM using ports P1–4; # Different subcomponents of the BM are represented using the same EM model.
Table 3. Butler matrix 2: performance figures of the optimized structure (EM simulation).
Table 3. Butler matrix 2: performance figures of the optimized structure (EM simulation).
ExcitationPort P1Port P2Port P3Port P4
RBW (dB)−14.9−17.1−18.3−14.5
IBW (dB)20.120.720.720.1
ΔMBW (dB)0.680.930.680.50
ΔMf0 (dB)0.380.520.280.10
Target phase (°)−30150−12060
ΔPBW (°)±3.6±4.3±3.8±3.1
ΔPf0 (°)±1.3±1.0±1.6±2.1
Table 4. Butler matrix 2: measured performance figures of the fabricated prototype.
Table 4. Butler matrix 2: measured performance figures of the fabricated prototype.
Excitation *Port P1Port P2Port P3Port P4
RBW (dB)−17.6−10.8−14.5−9.1
IBW (dB)15.814.510.310.9
ΔMBW (dB)2.012.753.991.05
ΔMf0 (dB)1.341.912.590.85
Target phase (°)−30150−12060
ΔPBW (°)±14.3±9.06±18.4±9.23
ΔPf0 (°)±7.48±7.83±10.9±8.79
* Due to frequency shift, the measured values are obtained for f0 = 2.45 GHz and BW from 2.35 to 2.55 GHz.
Table 5. Compact BM: comparison with state-of-the-art structures from the literature.
Table 5. Compact BM: comparison with state-of-the-art structures from the literature.
Designs Performance Figures *SizeRelative Size
(%) $
f0
(GHz)
BW &
(%)
ΔMf0
(dB)
Δθ1–4
(°)
ΔPf0
(°)
Dimensions
(mm × mm)
Dimensions
(λg × λg)
Conventional [21]1.07.00.60−45, 135, −135, 45±4.0233 × 1911.29 × 1.06100
[23]5.87.30.45−30, 150, −120, 60±3.871.4 × 1191.89 × 3.15437
[27]6.07.20.40−45, 135, −135, 45±0.958.9 × 57.91.96 × 1.93279
[16]283.24.70−45, 135, −135, 45±1616.8 × 14.9N/A #N/A #
[28] !1.85.52.40−45, 135, −135, 45±5.999.5 × 1270.82 × 1.0462.8
[30]1.05.92.30−45, 135, −135, 45±1.8164 × 1211.00 × 0.7454.1
[29]1.0N/A1.20−45, 135, −135, 45±1.087.8 × 82.40.49 × 0.4616.3
This work (BM1)2.68.40.51−30, 150, −120, 60±2.624.0 × 29.00.34 × 0.4210.5
This work (BM2)2.67.60.52−20, 160, −110, 70±2.125.9 × 29.60.37 × 0.4211.6
* Worst-case performance obtained for excitation of the BM using ports P1–4; # Multilayer structure implemented using substrates with different electrical parameters; $ With respect to conventional Butler matrix; ! Simulation results not available—data obtained from measurements; & Calculated for reflection below −15 dB and isolation above 15 dB.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bekasiewicz, A.; Koziel, S. Low-Cost Unattended Design of Miniaturized 4 × 4 Butler Matrices with Nonstandard Phase Differences. Sensors 2021, 21, 851. https://doi.org/10.3390/s21030851

AMA Style

Bekasiewicz A, Koziel S. Low-Cost Unattended Design of Miniaturized 4 × 4 Butler Matrices with Nonstandard Phase Differences. Sensors. 2021; 21(3):851. https://doi.org/10.3390/s21030851

Chicago/Turabian Style

Bekasiewicz, Adrian, and Slawomir Koziel. 2021. "Low-Cost Unattended Design of Miniaturized 4 × 4 Butler Matrices with Nonstandard Phase Differences" Sensors 21, no. 3: 851. https://doi.org/10.3390/s21030851

APA Style

Bekasiewicz, A., & Koziel, S. (2021). Low-Cost Unattended Design of Miniaturized 4 × 4 Butler Matrices with Nonstandard Phase Differences. Sensors, 21(3), 851. https://doi.org/10.3390/s21030851

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