Multibeam Wideband Transmit Beamforming Using 2D Sparse FIR Trapezoidal Filters

: A low-complexity multibeam wideband transmit beamformer using a 2D sparse FIR filter design capable of multiple beams is proposed as a digital building block for fully digital beamformers. The 2D sparse FIR filter has multiple trapezoid-shaped passbands pertaining to wideband beams aimed at particular directions. The proposed multibeam digital beamformer drives a uniform linear array of wideband antenna elements to achieve the wideband multibeam transmit-mode signals desired by the communication system. The 2D sparse FIR filter is designed to be optimal in the minimax sense using convex optimization techniques. Full-wave electromagnetic simulations using real antenna models confirm that the proposed wideband transmit beamformer can achieve multibeam transmission in the 1.3–2.8 GHz frequency range, with more than 70% fractional bandwidth. Furthermore, the adoption of the wideband transmit multibeam beamformer leads to a significant reduction in digital arithmetic (computational) complexity compared with previously reported wideband transmit beamformers of similar transfer function type, without deteriorating beam directionality and causing increases in the side-lobe level. The proposed sparse 2D FIR multibeam beamformer architecture is well-suited for both sub-6 GHz (legacy) band transmit beamforming, frequency range three (FR3) beamforming up to 28 GHz, and mmWave operation for emerging 5G/6G applications.


Multibeam Multichannel Access Points
In massive-MIMO arrays, the number of elements N can be orders of magnitude greater than the number of mutually independent wireless connections served by an access point.Further, access points having large apertures must connect to mobile units having both line of sight (LOS) and non-LOS RF propagation channels, in the presence of both short-and long-term fading.Non-LOS channels arise naturally due to the presence of large reflectors in the environment.However, non-LOS channels are at times intentionally created by 6G systems using technologies such as large intelligent surfaces (LISs) and reconfigurable intelligent surfaces (RISs) [20][21][22][23].Non-LOS channels may consist of multiple ray-like paths that bridge an access point to mobile units.Typically, the number of such dominant non-LOS paths takes small integer values (usually fewer than five) in number for many practical situations.Therefore, non-LOS channels require energy from the transmit power amplifiers to be spatially directed to multiple directions using M beams pointed at said directions, such that the best received signal-to-noise ratio (SNR) can be realized while imposing the lowest possible harmful interference to unintended receivers.A steerable multibeam multiuser transmit beamformer with capabilities for both channel equalization, adaptive interference control via spatial nulling, and rapid beam reconfiguration at low complexity is therein a rather crucial part of 5G/6G array processing systems.

Wideband vs. Narrowband Systems in Modern Applications
A beamformer may feed an antenna array with delays in a manner such that the transmitted waves interfere constructively in the far-field region along one or more desired directions.Destructive interference occurs along other directions, as determined by the zeros of the array factor polynomial, in order to achieve directional energy transmission and minimum harmful interference system wide.An aperture consisting of a number of uniformly spaced antennas driven by RF power amplifiers (PAs) enables a number of independent plane waves to be launched from a single array.
Assuming all of the transmit antennas are identical to each other, the far-field beam pattern is the product of the array factor and the element pattern of the antennas.Wideband beamforming refers to the situation when the bandwidth of interest is at least 10% of the center frequency of the operational band.The number of elements used in a transmit array determines the sharpness and angular accuracy of the beam.The sharpest possible beams are achieved when all of the beamforming weights are unity within the frequency band, leading to a rectangular windowing operation with known side-lobe levels of about −13.6 dB for a large number of elements.Such beams are achieved using true-time delay beamformers.However, use cases in massive-MIMO for 5G/6G systems exist where multiple independent beams are required for illumination of different angular ranges by independent information-carrying signals.For serving non-LOS receivers, a single information-carrying signal needs to be directed to multiple directions with nulls placed at predetermined directions.

Beam Sharpness vs. Degrees of Freedom
The sharpness of a given beam can be increased only by increasing the size of the physical aperture; the larger the number of antennas used, the larger the physical aperture, and thus the sharper the beams.Let the number of independent plane-waves that must be launched from an aperture be K.Then, K ≤ N, where N is the number of elements used.When K = N, one obtains the maximum degrees of freedom, i.e., support for the maximum number of independent plane-waves from the array.

RF-Analog Phasing vs. Fully Digital Schemes
There are several techniques for achieving multiple steerable and possibly independent RF transmit beams from an aperture containing N elements.Analog-RF, digital, and RF-analog-digital hybrid approaches [24,25] can be utilized for achieving wideband multibeam directed energy beamforming.In an analog-RF approach, as shown in Figure 1, signals corresponding to each independent plane-wave are applied to a digital-to-analog converter (DAC) and transmit chain, followed by a phase shifting and feed network before being applied to a PA that drives each antenna [6,8].Note that each DAC output is bandpass filtered (BPF) and optionally upconverted to the carrier frequency, after which each channel is amplified and split into K outputs via an active or passive microwave splitter.The outputs of the power splitter are individually routed to N analog RF dot-product units, each having K phase-shifters and a K:1 passive microwave combiner whose output is applied to the PA of each of the antennas in the array.Furthermore, the N number of K-point dot-products is equivalent to a single RF-analog matrix-vector multiplication, where the input is a Kelement vector, and the matrix is K × N. Here, beamforming is achieved at RF using tunable delay lines.Wideband delays, however, are difficult to realize on integrated circuits making low size and weight integrated systems difficult to realize.Nevertheless, for narrowband systems, it is possible to realize tunable phase offsets for each PA input using programmable phasing circuits on chip.But the issue is that we need a dedicated set of phase-shifters for each independent beam, together with low-distortion microwave combiners and other associated electronics for signal combination prior to driving PA inputs.

Advantages of an All-Digital Approach
We can solve many such problems by adopting a fully digital transmit beamformer.In fully digital beamformers [26][27][28], beam personalities including main beams and null synthesis for each independent beam are achieved using discrete domain linear and timeinvariant filtering methods implemented in real time on a fast digital signal processor (DSP).Multibeam combination is also achieved in the digital domain [29,30].A dedicated DAC is used for each transmit chain driving the PAs so that maximum flexibility and reconfigurability are obtained.

Circuit Complexity vs. Flexibility
The overall circuit complexity and cost of such a solution is not scalable for even modest values of K. Fully digital approaches allow maximum flexibility on the number of independent beams, beam angles, beam sizes, and location of nulls.The PAs are directly driven by dedicated DACs at each antenna in the aperture, implying there is no need for RF-analog phase-shifters or true time-delay units to be used.Nevertheless, we note that fully digital approaches require up-link chains with the additional dynamic range due to the fact that K independent signals are multiplexed into the same N channels.Further, DAC resolution needs to increase by log 2 K bits compared with the reference case where phasing and combining occur at the RF stage.

Fully Digital Wideband FIR Beamformers
Finite impulse response (FIR) filters can be used for transmit beamforming in a fully digital array.Wideband digital FIR realizations can be achieved with a set of narrowband FIR realizations, called sub-band beamformers, where a wideband signal is decomposed into a set of narrowband channels, and each channel is processed with complex-valued weight-and-sum feeding networks [6,31].
In [32], a low-complexity wideband transmit beamformer has been proposed using two-dimensional (2-D) sparse finite-extent impulse response (FIR) trapezoidal filters, designed using the windowing method and hard thresholding.This beamformer can achieve wideband operation with 33% fractional bandwidth and provides a path for lower computational complexity compared with comparable discrete Fourier transform (DFT)-based techniques.However, the FIR beamformer in [32] is limited to a single beam (i.e., K = 1), and the 2D sparse FIR trapezoidal filter is suboptimal because the windowing technique has been employed to design the filter.In the context of multibeam beamforming using FIR filters, one key limitation is computational complexity, especially for wideband systems with variously shaped passbands and narrow transition bands.The accuracy of the designed filters for each beam direction in a multibeam system is required to be high, leading to higher-order FIR filters with increased computational complexity.Achieving a wide passband for wideband operation can be challenging while maintaining a sharp transition band for side-lobe suppression.Our work addresses these limitations by designing a 2D sparse FIR trapezoidal filter to be optimal in the minimax sense for multibeam beamforming.

Proposed 2-D FIR Transmit Beamformer
We propose a fully digital transmit beamforming method based on 2-D sparse FIR trapezoidal filters and a uniform linear array (ULA) of wideband Vivaldi antennas [33], as shown in Figure 2. The proposed transmit beamformer is different from and improves on the standard FIR filter-and-sum transmit beamformers [34][35][36], where separate onedimensional (1-D) FIR filters are employed for each antenna to implement the fractional time delays to achieve wideband operation.In particular, the proposed 2-D sparse FIR is optimal in the minimax sense and is designed using convex optimization techniques.Fullwave electromagnetic simulations using real antenna models confirm that the proposed wideband transmit beamformer can achieve multibeam transmission in the frequency range of 1.3-2.8GHz, with more than 70% fractional bandwidth.Importantly, the proposed transmit beamformer provides considerable reduction in the computational complexity compared with nonsparse filter-and-sum transmit beamformers [34][35][36], 2-D nonsparse FIR filters [37][38][39][40][41][42], and 2-D sparse FIR filter [32], without deteriorating beam directionality and causing increases in the side-lobe level.Furthermore, the computational efficiency of our sparse FIR filter design makes them suitable for real-time beamforming applications in wideband systems, while the linear phase response within the passband of the filter contributes to maintaining signal fidelity during beam steering.
The proposed sparse 2-D FIR multibeam beamformer architecture is well-suited for both sub-6 GHz (legacy) band transmit beamforming, frequency range three (FR3) beamforming up to 28 GHz, as well as mmWave operation for emerging 5G/6G applications.Multibeam beamforming systems with sparse trapezoidal FIR filters offer several advantages, including flexibility in beam control, improved signal quality due to sharp transition bands, and the potential for real-time implementation with the sparse filter design.These advantages make them suitable for various applications.In the context of wireless communications, our work on multibeam beamforming could be beneficial for cellular networks to improve network capacity and user experience, especially in dense urban environments [43,44].Additionally, for mmWave communication systems, the ability to form narrow beams with low side-lobes using 2D trapezoidal FIR filters is crucial for establishing high-data-rate links [29,30].Multibeam beamforming allows radars to track multiple targets simultaneously or scan a wider area electronically [45][46][47].Two-dimensional trapezoidal FIR filters can help shape the radar beam for specific applications like target tracking or clutter suppression.We note that the underlying concepts are equally applicable for receive mode as well, even though we present only the transmit mode in this paper.Therefore, in remote sensing, multibeam systems can be used to acquire data from different observation angles [48,49].Furthermore, 2D trapezoidal FIR filters can help achieve the desired beam patterns for specific features of interest in the scene.

Review of Spectra of 2D Plane Waves Received by Uniform Linear Arrays
We briefly review the spatiotemporal modeling of 2D plane waves received by a ULA and their spectra in this section.We consider the receive mode for the review because most of the related works on wideband beamforming using 2-D and three-dimensional filters consider receive-mode beamforming.A plane wave received by a ULA having N + 1 wideband antennas placed along the x axis in the three-dimensional space (x, y, z) ∈ R 3 is shown in Figure 3a.Such a plane wave can be modeled as a 2-D discrete-domain spatiotemporal signal, after synchronous analog to digital conversion, which can be expressed as [50,51] where (n x , n ct ) ∈ Z 2 , θ (the zenith angle) denotes the direction of arrival (DOA) of the 2-D plane wave, ∆ x and ∆ ct are the sampling intervals in the x and ct domains, respectively, c is the constant speed of propagation of an electromagnetic wave, and w(n ct ) specifies the 1-D temporal signal carried by the 2-D plane wave.
The ideal discrete-domain spectrum PW(ω x , ω ct ) of pw(n x , n ct ) can be obtained as [50,51] where which lies on a straight line going through the origin of the 2-D discrete frequency domain, as shown in Figure 3b.The angle α between the ROS of the plane wave and the ω ct axis is given by [39,50] In the case of K (∈ Z + ) plane waves having different DOAs θ k and w k (ct), k = 1, 2, . . ., K, received by a ULA, the ideal discrete-domian spectrum PW K (ω x , ω ct ) is given by In this case, the spectral ROS R K lies on K straight lines going through the origin inside the principal Nyquist square of the 2-D discrete frequency domain with α k = tan −1 (sin(θ k )).
In practice, the true ROS of the generated discrete-domain 2-D plane wave is a narrow parallelogram when the number of antennas is finite.In the case of a 2-D broadband bandpass plane wave having a temporal bandwidth B ct with center frequency Ω C ct , as shown in the Figure 4, the desired 2-D plane wave (PW) (signal of interest) can be extracted by using a 2D filter having narrow-angle double-trapezoidal-shaped passbands [39,40].This will attenuate noise and radio frequency interference (RFI) having DOAs that are different from the DOA of the signal of interest.The same analysis can be applied to the transmit beamforming, where 2-D wideband plane waves are emitted by a ULA in K different directions, i.e, with direction of departures (DODs) θ k ; k = 1, 2, . . ., K. In this case, the spectral ROS of the required 2-D digital filter should consist of K narrow-angle double-trapezoidal-shaped passbands, as shown in Figure 4.In our proposed system, we use 2-D sparse FIR filters with trapezoidal passbands, as shown in Figure 4.A single trapezoidal passband is specified by the temporal bandwidth of the broadband bandpass plane wave, as well as the angle α k ; k = 1, 2, . . ., K which depends on the DOD and the half-angle δ k .The operation of the multibeam transmit beamformer is explained in detail in the next section.

Proposed Low-Complexity Multibeam Wideband Transmit Beamformer
In this section, we present the proposed multibeam wideband beamformer in detail.
To this end, we first present an overview of the architecture of the proposed multibeam wideband beamformer.

Proposed Sparse 2D FIR Multibeam Beamformer
The proposed transmit beamformer comprises a 2-D sparse FIR filter or a filter bank, DACs, RF front-end including power amplifiers (PAs), and a ULA of broadband antennas, as shown in Figure 2. The proposed transmit beamformer is fully digital and supports for K ≤ N independently steerable wideband beams in two modes: 1) one wideband signal is transmitted in K independent directions (denoted as single-signal multi-beam (SSMB)) and 2) K independent wideband signals are transmitted in K independent directions (denoted as multi-signal multi-beam (MSMB)).For the SSMB case, we employ a 2-D sparse FIR filter having K trapezoidal passbands, and for the MSMB case, we employ a filter bank consisting of 2-D sparse FIR filters having singe trapezoidal passband.Trapezoidal passbands are chosen because of their constant half-angle δ k throughout the passbands, as shown in Figure 4.This enables the multibeam beamformer to generate constant beamwidth in the main lobes throughout the wideband frequency range.We note that, in the MSMB case, a dedicated 2-D sparse FIR filter is used for each independent far-field plane-wave corresponding to an independent information-carrying signal with a digital combination of plane-wave components after the filter bank, allowing a single set of DACs and RF chain to support K independent beams using the ULA.
In order to explain the operation of the transmit beamformer, we consider the SSMB case, where the 1-D signal w(n ct ), n ct ∈ Z is required to transmit in the K directions θ 1 , θ 2 , . . . ,θ K .We call θ i (∈ [−90 • , 90 • ], i = 1, 2, . . ., K) as the DODs.Here, broadband antennas are fed with the 2-D mixed-domain signal s m (n x ∆ x , ct), where (n x , ct) ∈ Z × R, which consists of appropriately scaled and delayed continuous-time counterparts of w(n ct ) in order to achieve transmit beamforming.Note that the electromagnetic waves emanating from broadband antennas constructively interfere in the far field to form a beampattern with K beams along the directions of θ 1 , θ 2 , . . ., θ K .The 2-D sparse FIR filter H(z x , z ct ) with K trapezoidal passbands is employed to scale and delay w(n ct ) appropriately to realize K beams, where the ith (= 1, 2, . . ., K) trapezoidal passband provides the scaled and delayed w(n ct ) for the ith beam.This is analogous to the use of a 2-D filter with K trapezoidal passbands to realize K beams in the receive mode beamforming, as reviewed in Section 2. The DACs and the PAs in the RF chains are used for the discrete-time to continuous-time conversion and to realize the required power gain, respectively.Now, let us consider the processing of the 2-D sparse FIR filter.The 2-D discretedomain input signal s in (n x , n ct ), (n x , n ct ) ∈ Z 2 of H(z x , z ct ) is defined as H(z x , z ct ) should have trapezoidal-shaped passbands in the principal Nyquist square N of the 2-D discrete frequency domain (ω x , ω ct ) ∈ R 2 , as shown in Figure 4, to achieve the required scaled and delayed versions of the wideband input signal w(n ct ), in the 2-D discrete-domain output signal s out (n x , n ct ) of the digital filter.Here, N = {(ω x , ω ct ) ∈ R 2 , −π ≤ ω x , ω ct < π}.The output signal s out (n x , n ct ) is given by where I is the set containing the indices of the nonzero coefficients of the impulse response h s (i x , i ct ) of H(z x , z ct ).Here, s out (n x , n ct ) is computed for −(N x −1) /2 ≤ n x ≤ (N x −1) /2 and 0 ≤ n ct ≤ (N ct − 1), where N x and N ct are the number of broadband antennas in the ULA and the number of temporal samples in w(n ct ), respectively.The design of H(z x , z ct ) with sparse coefficients is presented next.

Two-Dimensional Sparse FIR Trapezoidal Filter
The 2-D sparse FIR trapezoidal filter H(z x , z ct ) is designed by using the FIR filter design method presented in [52].To this end, the required passband R PB of the 2D filter H(z x , z ct ) having K double trapezoidal passbands inside the principal Nyquist square N can be derived as where R k is the trapezoidal passband corresponding to the kth bandpass 2-D PW, which is given by where B k ct is the temporal bandwidth of the k th PW, Ω C k ct is the center frequency of the k th PW, and δ k is the half-angle of the double trapezoidal-shaped k th passband.The value of δ k affects the beamwidth of the multibeams produced by the transmit beamformer.
The ideal frequency response H d (e jω x , e jω ct ) of H(z x , z ct ) inside N is specified as Note that H d (e jω x , e jω ct ) has a zero group delay with respect to both ω x and ω ct .

Minimax Design of the 2D Sparse FIR Filter
The design of H(z x , z ct ) is presented in this subsection, adapting the 2-D sparse FIR filter design method proposed in [52].We choose the filter design method in the minimax sense to reduce the side-lobe levels.To this end, we express the frequency response of H(z x , z ct ) as where h(n x , n c ) is the impulse response of size The order of the filter depends on the number of broadband antennas in the ULA and the expected accuracy of the beams produced by the transmit beamformer.H(z x , z ct ) is designed as a zero-phase filter ([53], Chapter 3).The impulse response of the filter is centrosymmetric in this case, i.e., h(n x , n c ) = h(−n x , −n c ).Therefore, considering the centrosymmetric property, H(e jω x , e jω ct ) can be simplified as which is represented in vector form by where, and In [52], a two-phase design method for 2D FIR digital filters in weighted minimax design applicable for filters having quadrantally symmetric impulse responses was presented.With (13), this approach can be used to design H(z x , z ct ) despite its centro-symmetric impulse response.

Step 1
We obtain an intermediate sparse impulse response in the weighted minimax sense.The objective function to be minimized is expressed as where H(e jω x , e jω ct ) is given in ( 13), H d (e jω x , e jω ct ) is given in (10), and W(ω x , ω ct ) is a weighting function used to control the passband ripple and stopband ripple, which affects the side-lobes of the transmit beamformer across the bandwidth.Parameter µ is a small positive number (typically between 0.001 and 1 [54]), and F is the region corresponding to the passband and stopband, i.e., without the transition band [52].Due to the centrosymmetric impulse response of H(e jω x , e jω ct ), only the region [−π, π] × [0, π] in the 2-D frequency domain (ω x , ω ct ) is considered to define F .By introducing an upper bound β for W(ω x , ω ct )|H(e jω x , e jω ct ) − H d (e jω x , e jω ct )| over a finite M set of frequency grids 16) is reduced to a tractable constrained problem as To deal with the nondifferentiability of ∥h∥ 1 , taking an upper bound for each entry of h, i.e., with n = 1 + N x +N x N c +N c 2 , leads (17) to Treating the bounds β and d i (1 ≤ i ≤ n) as auxiliary design variables in (19), the optimization problem becomes a convex optimization problem, which can be converted as a linear programming (LP) problem [52] expressed as where with 1 , e jω (1) 2 ) . . .

−W(ω
The unique and globally optimal solution h of the linear programming problem is an approximately sparse [52].To obtain a sparse impulse response h s , hard thresholding is employed, i.e., where ϵ th (∈ [10 −4 , 10 −2 ] typically [54]) is the threshold value.

Step 2
We again optimize h s in the weighted minimax sense in order to obtain the optimal solution for the 2-D sparse FIR filter.This optimization problem is expressed as minimize where I ∞ is the set containing indices i for which h s i = 0. Again, by introducing an upper bound β for W(ω x , ω ct )|c(ω x , ω ct ) T h s − H d (e jω x , e jω ct )| over a finite M set of frequency grids 24) is reduced to a tractable constrained problem as minimize β By defining vectors ĥs and ĉ(ω) by deleting those entries of h s and c(ω) whose indices belong to set I ∞ , we have c(ω) T h s = ĉ(ω) T ĥs .Hence, (25) can be written as minimize β With β as an auxiliary variable, problem (26) becomes an LP problem of the form where 1 , e jω (1) 2 ) . . .
By solving (27), the impulse response of the required 2-D FIR filter h s can be obtained.

Analysis and Selection of Parameters
Finding an optimal balance between the performance and the computational complexity of the sparse FIR filter is a challenge.When solving the optimization problem in the design of the filter, the parameters µ and ϵ th were chosen such that the sparsity of the filter impulse response was significant without substantially degrading the filter performance.For a particular filter design, the value for the parameter µ is varied in the range [0.001, 1], fixing the number of nonzero values in the impulse response, and the filter performance is observed in terms of the optimal value of step 2 of the optimization problem, which is the weighted maximum error between the desired response and the actual response, E opt .The value that gave the best performance out of the tested values is chosen for the design of the particular filter.After choosing a suitable value for µ, to choose a suitable threshold ϵ th , the root-mean-square-error (RMSE) E, and E opt , are considered.E is given by When the number of nonzero coefficients is decreased by using a suitable ϵ th , the E opt increases rapidly after some point, as shown in Figure 5. Therefore, ϵ th is chosen such that it yields the required number of nonzero values without increasing the E opt substantially.As a rule of thumb, ϵ th is chosen such that E opt does not exceed more than 0.05 of its initial value and E < 10%.As an example, Figure 5 shows the variation in E opt vs. the nonzero values in the impulse response for the case of the single-beam filter design problem explained in Section 4.2.In that case, the chosen ϵ th is 0.007125, which corresponds to 143 nonzero values in the impulse response.The value E was 7.095%, and E opt was 0.1364, which had an initial value of 0.08795.To have the same amount of side-lobes as the passband ripple, the passband and the stopband weights are taken as 1.

Full-Wave Electromagnetic Simulation Results
Full-wave electromagnetic simulations, conducted in CST Microwave Studio © (CST-MWS) (Release Version 2020.01) for a ULA of 15 broadband antipodal Vivaldi antennas, are used to obtain the far-field transmit beam patterns produced by the proposed transmit beamformer.Here, we consider both multibeam and single-beam beamforming.The transmit signal frequency range of interest is GHz, which has a fractional bandwidth of more than B f = 70%.The ULA of Vivaldi antennas is shown in Figure 6a, and the S 11 curve of a single antenna element is shown in Figure 6b, verifying the wideband element response within the frequency range of interest.The interelement spacing is selected as λ/2 such that ∆x = ∆ct, where λ is the wavelength corresponding to the temporal sampling frequency F s = 6 GHz.The passband of the 2-D FIR trapezoidal filter is taken such that its temporal frequency is in this frequency range of interest, i.e., [0.43π, 0.93π] rad/sample.The order and the maximum transition band of the filter, and other parameters of the trapezoid which vary with different cases of beamforming are mentioned in the corresponding places in Sections 4.1 and 4.2.These include angles For the single-beam case, the results of the proposed method are compared with similar approaches adapted from [32,[39][40][41] such that the filters have the same specifications.In particular, we present the average of the DODs, their standard deviation (SD), average gain, and side-lobe levels in the passband frequency range 1.3-2.8GHz.Side-lobe levels are considered in the range θ ∈ [−90 • , 90 • ], and the peak and average side-lobe levels at a particular temporal frequency are computed with respect to the average gain of the main beams.Additionally, 2-D FIR filter performance in terms of the passband ripple and the stopband attenuation, and the DSP adder and multiplier complexity are presented.

Multibeam Beamforming
Here, beamforming is achieved in K ∈ Z + directions, where K ≥ 2. The DODs are  7e and Figure 7f, respectively.Furthermore, Table 1 provides a comparison of the wideband beam pattern and DSP complexity between the proposed 2D sparse and nonsparse FIR filter-based transmit beamformers.The effect of sparsification is minimal on the passband ripple and side-lobe levels, while saving more than 85% in DSP complexity.

Same Signal Transmission in Two DODs
As an example of two-direction beamforming of the same signal, consider the case where θ 1 = −26 • and θ 2 = 28 • , with corresponding α 1 = −23.67• and α 2 = 25.15 • , respectively.The magnitude frequency response |H(e jω x , e jω ct )| of the 2-D FIR filters of order 14 × 64, with a nonsparse impulse response and sparse impulse response (with only 161 nonzero elements), is shown in Figure 8a and Figure 8b, respectively.Values for the other parameters are µ = 0.1 and ϵ th = 0.00701.The wideband transmit beam patterns produced by the nonsparse and sparse filters are shown in Figure 8e and Figure 8f, respectively.Furthermore, Table 2 provides a comparison of the wideband beam pattern and DSP complexity between the proposed 2D sparse and a nonsparse FIR filter-based transmit beamformers.We note that the effect of sparsification is minimal, while saving more than 83% in DSP complexity.

Same Signal Transmission in Three DODs
As an example of three-direction beamforming, consider the case where θ 1 = −45 • , θ 2 = 0 • , and θ 3 = 45 • , with corresponding α 1 = −35.26• , α 2 = 0 • , and α 3 = 35.26• , respectively.The magnitude frequency response |H(e jω x , e jω ct )| of the 2-D FIR filters of order 14 × 64, with a nonsparse impulse response and a sparse impulse response (only with 161 nonzero elements), is shown in Figure 8c and Figure 8d, respectively.Values for the other parameters are µ = 0.02 and ϵ th = 0.010455.The wideband transmit beam patterns produced by the nonsparse and sparse filters are shown in Figure 8g and Figure 8h, respectively.Moreover, Table 2 provides a comparison of wideband beam pattern and DSP complexity between the proposed 2D sparse and nonsparse FIR filter-based transmit beamformers.The results clearly indicate the minimal effect of sparsification of the filter impulse response towards the performance of the filter with regard to beamforming, while saving more than 83% in DSP complexity.
As an example of single-direction beamforming, consider the case where the DOD is 30 • , where α = 26.57• .We design 2-D FIR trapezoidal filters using the windowing method and optimization techniques with nonsparse and sparse impulse responses.Here, the design of the filter with the windowing technique is adapted from [39,40] for the nonsparse impulse response whereas from [32] for the sparse impulse response.Furthermore, nonsparse filter designed using the optimization technique is adapted from [41] (without considering the mutual-coupling compensation for fair evaluation).For the windowing technique, the Dolph-Chebyshev window is used as the window function, because a Dolph-Chebyshev window provides an approximately equiripple frequency magnitude response ( [55], Chapter 9.4), and therefore an approximately minimax design.Furthermore, for the sparse impulse response [32], hard thresholding with a threshold of 2.544 × 10 −4 is employed.The order of all of the filters is selected as 14 × 64, and the half-angle of the trapezoidal filter δ is selected as 3 • .Furthermore, the maximum transition band is 0.1π rad/sample.The proposed filter is designed with the parameters µ = 0.05 and ϵ th = 0.007144.3 provides a comparison of the accuracy of the wideband beam patterns and DSP complexity for the proposed transmit beamformer and previously proposed filter designs.We observe that a significant reduction (more than 8 dB) in the passband ripple can be with the proposed transmit beamformer with similar stopband attenuation compared with the nonsparse [39,40] and sparse [32] filter designs based on the windowing technique.Furthermore, the average side-lobe level of the proposed beamformer is lower than 7 dB, and the peak side-lobe level is lower than 4 dB compared with these filters.Compared with the nonsparse optimization-based design [41], the proposed transmit beamformer provides a more than 85% reduction in DSP complexity without substantial degradation of the accuracy of the beam pattern.This confirms the superior performance of the proposed transmit beamformer compared with the transmit beamformers that employ 2-D FIR filters designed using previously proposed techniques [32,[39][40][41].[41] (c) windowing method [32] (sparse), and (d) proposed (sparse).Polar plots of the simulated beam patterns obtained using the filter made using (e) windowing method [39,40] (nonsparse), (f) optimization [41] (nonsparse), (g) windowing method [32] (sparse), and (h) proposed (sparse).
Table 3. Accuracy of the beam pattern and the DSP complexity of the proposed sparse transmit beamformer and the single-beam sparse transmit beamformer designed using the windowing method [32] for the DOD = 30

Conclusions and Future Work
We propose a digitally fed wideband transmit multibeam beamformer using lowcomplexity 2-D sparse FIR trapezoidal filters with multiple trapezoidal passbands.The 2-D sparse FIR trapezoidal filter transforms a given wideband temporal input signal into a 2-D spatiotemporal signal with appropriate delays so that directional transmission can be achieved with a ULA of wideband antennas.We adapt a two-step sparse filter design method to design the 2-D sparse FIR trapezoidal filter optimal in the minimax sense.Fullwave electromagnetic simulations verify accurate wideband transmit beam patterns with more than 70% fractional bandwidth.More importantly, our 2-D sparse FIR trapezoidal filter provides significantly lower computational complexity compared with nonsparse counterparts.Furthermore, the proposed beamformer provides significantly better beam patterns for single-beam beamforming compared with the previously proposed sparse beamformer designed using the windowing method.Future work may include an extension of the proposed transmit beamformer to planar and other array geometries, implementations using analog-domain multidimensional filters, and experimental validation.

KFigure 1 .
Figure 1.Conventional fully analog phasing matrices for achieving K beams: an N-element wideband antenna array is fed through N power amplifiers (PAs) using a digital intermediate frequency (IF) or baseband signal processor having N digital to analog converters (DACs).

Figure 2 .
Figure 2. Architecture of the proposed transmit beamformer for the MSMB case, consisting of a filter bank of 2D sparse FIR trapezoidal filters, DACs, and RF front end, as well as a uniform linear array of broadband antennas.
), and δ(•) is the 1-D continuous-domain impulse function.Note that ω x is the discrete-domain spatial frequency, and ω ct is the discrete-domain temporal frequency scaled by c.The region of support (ROS) R of PW(ω x , ω ct ) can be obtained as[50,51]

Figure 5 .
Figure 5. E opt variation with the number of nonzero values in the impulse response of the proposed 2-D FIR trapezoidal filter.
2, . . ., K and the half-angle δ k with respect to Figure 4. Matlab © (Version R2018a) software is used for generating the impulse response of the desired 2-D sparse FIR filter responsible for beamforming with the help of the CVX optimization toolbox and SeDuMi solver.It is then exported and fed into CST-MWS ULA model as port excitation signals for each antenna.To obtain the far-field transmit beam patterns in CST-MWS, farfield monitors are placed within the frequency range of interest.Transmit beam patterns produced by both the sparse and nonsparse 2-D FIR filters are then obtained and compared.
2, . . ., K. Results for different independent signal transmissions (i.e., MSMB mode) and the same signal transmissions in different desired directions (i.e., SSMB mode) are discussed below.In the following examples, all the trapezoidal filters' half-angles δ k are taken as 3 • and their maximum transition band is 0.1π rad/sample.The region [−π, π] × [0, π] in the 2-D frequency domain (ω x , ω ct ) is uniformly discretized into 201 × 101 points.The passband and stopband weights are taken as 1. 4.1.1.Two Independent Signals Transmission in Two DODs As an example of two-direction beamforming where two independent signals are to be transmitted in different directions, consider the case where θ 1 = −26 • and θ 2 = 28 • with corresponding α 1 = −23.67• and α 2 = 25.15 • , respectively.The magnitude frequency response |H(e jω x , e jω ct )| of the sum of 2-D FIR filters of order 14 × 64, with nonsparse impulse responses and with sparse impulse responses, is shown in Figure 7a and Figure 7b, respectively.Values for the other parameters are µ = 0.05 and ϵ th = 0.004997 for DOD = −26 • filter and µ = 0.05 and ϵ th = 0.007005 for DOD = 28 • filter.The wideband transmit beam patterns produced by the nonsparse filter and the proposed sparse filter outputs are shown in Figure

4. 1 . 2 .
Three Independent Signals Transmission in Three DODs As an example of three-direction beamforming with three independent signals, consider the case where θ 1 = −45 • , θ 2 = 0 • , and θ 3 = 45 • , with corresponding α 1 = −35.26• , α 2 = 0 • , and α 3 = 35.26• , respectively.The magnitude frequency response |H(e jω x , e jω ct )| of the sum of 2-D FIR filters of order 14 × 64, with nonsparse and with sparse impulse responses, is shown in Figure 8c and Figure 8d, respectively.Values for the other parameters are µ = 0.05 and ϵ th = 0.006998 for DOD = ±45 • filters and µ = 0.15 and ϵ th = 0.006094 for DOD = 0 • filter.The wideband beam patterns produced by the nonsparse and sparse filters are shown in Figure7gand Figure7h, respectively.Moreover, Table1provides a comparison of wideband beam pattern and DSP complexity between the proposed 2D sparse and nonsparse FIR filter-based transmit beamformers.The results clearly indicate the minimal effect of sparsification of the filter impulse response towards the performance of the filter with regard to beamforming, saving more than 86% in DSP complexity.
Furthermore, the passband and stopband weights are taken as 1, and the region [−π, π] × [0, π] in the 2-D frequency domain (ω x , ω ct ) is uniformly discretized into 201 × 101 points.The 2-D sparse FIR filters, i.e., those proposed and designed adapting the method in [41], have only 143 nonzero coefficients, whereas nonsparse filters have 975 coefficients.The magnitude frequency responses of the designed filters are shown in Figure 9a-d.Furthermore, Table

Table 1 .
Accuracy of the beam pattern and the DSP complexity of the proposed sparse transmit beamformer and a nonsparse counterpart for the MSMB mode.

Table 2 .
Accuracy of the beam pattern and the DSP complexity of the proposed sparse transmit beamformer and a nonsparse counterpart for the SSMB mode.