Next Article in Journal
IoT–Blockchain: Harnessing the Power of Internet of Thing and Blockchain for Smart Supply Chain
Next Article in Special Issue
Dynamic Divide Grouping Non-Orthogonal Multiple Access in Terrestrial-Satellite Integrated Network
Previous Article in Journal
Icing Condition Predictions Using FBGS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Low Complexity Hybrid Precoding Designs for Multiuser mmWave/THz Ultra Massive MIMO Systems

1
Department of Information Science and Technology, ISCTE-Instituto Universitário de Lisboa, 1649-026 Lisboa, Portugal
2
Instituto de Telecomunicações, 1049-001 Lisboa, Portugal
3
Department of Electrical and Computer Engineering, Faculty of Science and Technology, Universidade Nova de Lisboa, 1099-085 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(18), 6054; https://doi.org/10.3390/s21186054
Submission received: 14 July 2021 / Revised: 25 August 2021 / Accepted: 8 September 2021 / Published: 9 September 2021
(This article belongs to the Special Issue Transmission Techniques for Future Global Mobile Systems)

Abstract

:
Millimeter-wave and terahertz technologies have been attracting attention from the wireless research community since they can offer large underutilized bandwidths which can enable the support of ultra-high-speed connections in future wireless communication systems. While the high signal attenuation occurring at these frequencies requires the adoption of very large (or the so-called ultra-massive) antenna arrays, in order to accomplish low complexity and low power consumption, hybrid analog/digital designs must be adopted. In this paper we present a hybrid design algorithm suitable for both mmWave and THz multiuser multiple-input multiple-output (MIMO) systems, which comprises separate computation steps for the digital precoder, analog precoder and multiuser interference mitigation. The design can also incorporate different analog architectures such as phase shifters, switches and inverters, antenna selection and so on. Furthermore, it is also applicable for different structures, namely fully-connected structures, arrays of subarrays (AoSA) and dynamic arrays of subarrays (DAoSA), making it suitable for the support of ultra-massive MIMO (UM-MIMO) in severely hardware constrained THz systems. We will show that, by using the proposed approach, it is possible to achieve good trade-offs between spectral efficiency and simplified implementation, even as the number of users and data streams increases.

1. Introduction

Over the last few years, significant advances have been made to provide higher-speed connections to users in wireless networks, with several novel technologies being proposed to achieve this objective. However, future generations of communication systems will have to fulfil more demanding requirements that cannot be met by the methods adopted in today’s communications systems. This motivates the exploration of other candidate technologies, like the millimeter wave (mmWave) and Terahertz (THz) bands, where many applications that require ultra-high data rates can be designed. These bands offer great underutilized bandwidths and also allow for a simplified implementation of large antenna arrays, which are crucial to combat the severe signal attenuation and path losses that occurs at these frequencies [1,2,3,4]. The first applications of the THz band were limited to imaging and sensing due to the unavailability of efficient devices that can work on these frequencies. However, recent advances in the field of THz devices give us strong indications that THz communications will be feasible in the near future [3]. These technologies (THz systems in particular), are expected to ease the spectrum limitations of today’s systems. They face several issues, such as the reflection and scattering losses through the transmission path, the high dependency between distance and frequency of channels at the THz band and the need for controllable time-delay phase shifters, since the phase shift will vary with frequencies based on the signal traveling time, which will also affect the system performance. These limitations require not only the proper system design, but also the definition of a set of strategies to enable communications [5,6].
The exploration of the potentialities of millimeter and sub-millimeter wavelengths is closely related to the paradigm of using very large arrays of antennas in beamforming architectures. This gives rise to so-called ultra-massive multiple-input multiple-output (UM-MIMO) systems. Still, to achieve the maximum potential of these systems, it is necessary to consider the requirements and the challenges related not only to the channel characteristics but also to the hardware component, especially regarding THz circuits [5,7,8]. Considering that high complexity and power usage are pointed out as major constraints of large-antenna systems, it is unfeasible to implement UM-MIMO schemes with a dedicated RF chain per antenna element in the mmWave and THz bands. Therefore, instead of fully-digital precoders and combiners, it becomes crucial to adopt hybrid digital/analog architectures as these require a reduced number of dedicated RF chains. By adopting this type of design, the signal processing is split into two separate parts: a low-dimensional digital part and a complementary analog part. This approach can enable a substantial reduction in the overall circuit complexity and power consumption [9]. By adopting a proper problem formulation, the analog design part can then be reduced to a simple projection operation in a flexible precoding or combining algorithm that can cope with different architectures, as we proposed in [10,11]. Despite the ultra-wide bandwidths available in mmWave and THz bands, and besides considering the problem of distance limitation, MIMO systems should take into account the operation in frequency selective channels [12]. To make the development of hybrid schemes for these systems a reality, it is necessary to handle the fading caused by multiple propagation paths typical in these types of channels [13]. Therefore, solutions inspired on multi-carrier schemes, such as orthogonal frequency division multiplexing (OFDM), are often adopted to address such problems [14].
Spectral Efficiency (SE) of point-to-point transmissions is a major concern in SingleUser (SU) and MultiUser (MU) systems. To achieve good performances, it is necessary to develop algorithms that are specially tailored to the architecture of these systems. Several hybrid precoding schemes have been proposed in the literature [15,16,17,18,19]. The authors of [15] proposed two algorithms for low-complexity hybrid precoding and beamforming for MU mmWave systems. Even though they assume only one stream per user, i.e., the number of data streams (Ns) is equal to the number of users (Nu), it is shown that the algorithms achieve interesting results when compared to the fully-digital solution. The concept of precoding based on adaptive RF-chain-to-antenna was only introduced in [16] for SU scenarios, but showed promising results. In [17], a nonlinear hybrid transceiver design relying on Tomlinson–Harashima precoding was proposed. Their approach only considers Fully-Connected (FC) architectures but can achieve a performance close to the fully-digital transceiver. A Kalman-based Hybrid Precoding method was proposed for MU scenarios in [18]. While designed for systems with only one stream per user and based on fully-connected structures, the performance of the algorithm is competitive with other existing solutions. A hybrid MMSE-based precoder and combiner design with low complexity was proposed in [19]. The algorithm is designed for MU-MIMO systems in narrowband channels, and it presents lower complexity and better results when compared to Kalman’s precoding. Most of the hybrid solutions for mmWave systems aim to achieve near-optimal performance using FC structures, resorting to phase shifters or switches. However, the difficulty of handling the hardware constraint imposed by the analog phase shifters or by switches in the THz band is an issue that limits the expected performance in terms of SE.
Array-of-Subarrays (AoSAs) structures have gained particular attention over the last few years as a more practical alternative to FC structures, especially for the THz band. In contrast to FC structures, in which every RF chain is connected to all antennas via an individual group of phase shifters (which is prohibitive for higher frequencies), the AoSA approach allows us to have each RF chain connected to only a reduced subset of antennas. The adoption of a disjoint structure with fewer phase shifters reduces the system complexity, the power consumption and the signal power loss. Moreover, all of the signal processing can be easily carried out at the subarray level by using an adequate number of antennas [6].
Following the AoSA approach, it was shown in [20] that, to balance SE and power consumption in THz communications, adaption and dynamic control capabilities should be included in the hybrid precoding design. Therefore, Dynamic Arrays-of-Subarrays (DAoSAs) architectures could be adopted. The same authors proposed a DAoSA hybrid precoding architecture which can intelligently adjust the connections between RF chains and subarrays through a network of switches. Their results showed that it is possible to achieve a good trade-off for the balancing between the SE and power consumption.
Within the context of multiuser downlink scenarios, the authors of [21] studied some precoding schemes considering THz massive MIMO systems for Beyond 5th Generation (B5G) networks. Besides showing the impact on EE and SE performance, carrier frequency, bandwidth and antenna gains, three different precoding schemes were evaluated and compared. It was observed that the hybrid precoding approach with baseband Zero Forcing for multiuser interference mitigation (HYB-ZF) achieved much better results than an Analog-only Beamsteering (AN-BST) scheme with no baseband precoder. In fact, this approach was capable of better approaching the upper bound defined by the singular value decomposition precoder (SVD-UB). The other relevant conclusion is that the design of precoding algorithms should be adapted to the communication schemes. While considering all the specific constraints may allow the maximization of the system performance of the system, formulating and solving the corresponding optimization problem may not be so simple. Motivated by the work above, in this paper we developed an algorithm for hybrid precoding design which can accommodate different low-complexity architectures suitable for both mmWave and THz MU-MIMO systems. It is based on the idea of accomplishing a near-optimal approximation of the fully-digital precoder for any configuration of antennas, RF chains and data streams through the application of the alternating direction method of multipliers (ADMM) [22]. ADMM is a well-known and effective method for solving convex optimization problems but can also be a powerful heuristic for several non-convex problems [22,23]. To use it effectively within the context of MU-MIMO, THE proper formulation of the hybrid design problem as a multiple constrained matrix factorization problem is first presented. Using the proposed formulation, an iterative algorithm comprising several reduced complexity steps is obtained.
The main contributions of this paper can be summarized as follows:
  • We propose a hybrid design algorithm with near fully-digital performance, where the digital precoder, analog precoder and multiuser interference mitigation are computed separately through simple closed-form solutions. Even though the hybrid design algorithm is developed independently of a specific channel or antenna configuration, it is particularly suitable for mmWave and THz systems where, on the one hand, very large antenna arrays are required to overcome distance limitations but, on the other hand, current hardware constraints in terms of cost and power consumption make the adoption fully-digital precoders/combiners with one dedicated RF chain per antenna element unviable. Whereas our previous work [10] also proposed a hybrid design algorithm for mmWave, it did not address multiuser systems, and in particular the MIMO broadcast channel. Therefore, it does not include any step for inter-user interference mitigation within its design. As we show here, for this multiuser channel, the hybrid design method must also deal with the residual inter-user interference as it can degrade system performance, particularly at high Signal Noise Ratios (SNRs);
  • Due to the separability of the different steps (analog precoder, digital precoder and interference suppression), the proposed algorithm can incorporate different architectures, making it suitable for supporting UM-MIMO in severely hardware-constrained systems typical in the THz band. Unlike [10], where we only considered the adoption of phase shifters, in this paper we present explicit solutions for some of the most common architectures, namely FC, AoSA and DAoSA structures based in either Unquantized Phase Shifters (UPS), Quantized Phase Shifters (QPS), Switches (Swi), Switches and Inverters (SI), Antenna Selection (AS) or Double Phase Shifters (DPS);
  • To cope with the large bandwidths available in mmWave/THz bands, where practical MIMO systems likely have to operate in frequency selective channels, the proposed hybrid design considers the application in a multicarrier context, where the same analog precoder is applied at different frequencies;
  • We explicitly show how the proposed design can be applied to a DAoSAs approach where a reduced number of switches are inserted at each AoSA panel, which allows the connections to the RF chains to be dynamically adjusted. Through extensive simulations, it is shown that our proposed solution is capable of achieving good trade-offs between spectral efficiency, hardware complexity and power consumption, proving to be a suitable solution for the deployment of UM-MIMO, especially in hardware-constrained THz systems.
The paper is organized as follows: Section 2 presents the adopted system model. The adopted formulation of the hybrid design problem for the MU-MIMO scenario and the proposed algorithm are described in detail in Section 3, which includes the implementation of the algorithm for different analog architectures. Performance results are then presented in Section 4. Finally, the conclusions are outlined in Section 5.
Notation: Matrices and vectors are denoted by uppercase and lowercase boldface letters, respectively. The superscript ( . ) T and ( . ) H denote the transpose and conjugate transpose of a matrix/vector,     p is the p - norm of a vector,     0 is its cardinality (i.e., the number of non-zero elements in a vector which is sometimes referred to as the 0 - norm in the literature) and I n is the n × n identity matrix.

2. System Model

In this section, we present the system and channel models adopted for the design of the hybrid precoding algorithm. Let us consider the OFDM base system illustrated in Figure 1. In this case, we have a mmWave/THz hybrid multiuser MIMO system, where a base station (BS) is equipped with N t x antennas and transmits to N u users equipped with N r x antennas over F carriers, as can be seen in Figure 1. On each subcarrier, N s data streams are transmitted to each user which are represented as s k = [ s k , 1 T s k , N u T ] T , with s k , u N s × 1 . Instead of a fully-digital design which would require a dedicated RF chain per antenna element, both the precoder and combiner comprise separate digital and analog processing blocks. This approach allows for the use of reduced digital blocks with only a few RF chains, which are complemented by the analog blocks, that can be supported solely on networks of phase shifters and switches. Since the analog precoder (combiner) is located after (before) the IFFT (FFT) blocks, it is shared between the different subcarriers, as in [24,25]. Regarding the analog precoder and combiner, which are represented by matrices F R F N t x × N R F t x and W R F u N r x × N R F r x with u = 1 , , N u , it is assumed that N u N s N R F t x N t x and N s N R F r x N r x , where N R F t x and N R F r x are the number of RF chains at the BS and each user, respectively. The received signal model at subcarrier k after the combiner can be written as
Y k , u = ρ u W BB k , u H W RF u H H k , u F RF F BB k s k + W BB k , u H W RF H n k , u ,
where H k , u N r x × N t x is the frequency domain channel matrix (assumed to be perfectly known at the transmitter and receiver) between the base station and the uth receiver at subcarrier k. Vector n k , u N r x × 1  contains independent zero-mean circularly symmetric Gaussian noise samples with covariance σ n 2 I N r x and ρ u denotes the average received power. The digital baseband precoders and combiners are denoted by F BB k N R F t x × N u N s and W BB k , u N R F r x × N s , respectively. Regarding the channel model, it is important to note that even though the mmWave and THz bands share a few commonalities, the THz channel has several peculiarities that distinguish it from the mmWave channel. For example, the very high scattering and diffraction losses in the THz band will typically result in a much sparser channel in the angular domain with fewer multipath components (typically less than 10) [21]. Furthermore, the gap between the line of sight (LOS) and non-line of sight (NLOS) components tends to be very large, making it often LOS-dominant with NLOS-assisted [25]. An additional aspect relies on the much larger bandwidth of THz signals which can suffer performance degradation due to the so-called beam split effect, where the transmission paths squint into different spatial directions depending on the subcarrier frequency [20]. In light of this, in this paper we consider a clustered wideband geometric channel, which is commonly adopted both in the mmWave [15] and THz literature [20,25,26,27,28] However, it should be noted that the hybrid precoding/combining approach proposed in this paper is independent of a specific MIMO channel. In this case, the frequency domain channel matrices can be characterized as
H k , u = γ ( α u L O S a r ( ϕ u r , L O S , θ u r , L O S ) a t ( ϕ u t , L O S , θ u t , L O S ) H + i = 1 N c l l = 1 N r a y α i , l , u a r ( ϕ i , l , u r , θ i , l , u r ) a t ( ϕ i , l , u t , θ i , l , u t ) H ) e j 2 π τ i , u f k ,
where N c l denotes the scattering clusters with each cluster i having a time delay of τ i , u and N r a y propagations paths. α u L O S α i , l , u are the complex gains of the LOS component and of the lth ray from cluster i. Index u is the user ( u = 1 , , N u ), f k = f c + B F ( k 1 F 1 2 ) ( k = 1 , , F ) is the kth subcarrier frequency, B is the bandwidth, fc is the central frequency and γ is a normalizing factor such that Ε [ H k , u F 2 ] = N t x N r x . Vectors a t ( ϕ i , l , u t , θ i , l , u t ) and a r ( ϕ i , l , u r , θ i , l , u r ) represent the transmit and receive antenna array responses at the azimuth and elevation angles of ( ϕ i , l , u t , θ i , l , u t ) and ( ϕ i , l , u r , θ i , l , u r ) , respectively. Vectors a t ( ϕ u t , L O S , θ u t , L O S ) and a r ( ϕ u r , L O S , θ u r , L O S ) have similar meanings but refer to the LOS path angles ( ϕ u t , L O S , θ u t , L O S ) and ( ϕ u r , L O S , θ u r , L O S ) . By carefully selecting the parameters of the channel model we can make it depict a mmWave or a THz channel. Considering Gaussian signaling, the spectral efficiency achieved by the system for the transmission to MS-u in subcarrier k is [28]
R k , u = log 2 | I N R F r x + R u 1 W BB k , u H W RF u H H k , u F RF F BB k , u × F BB k , u H F RF H H k , u H W RF u W BB k , u | ,
where R k , u is the covariance matrix of the total inter-user interference plus noise at MS-u, which is characterized by
R k , u = W BB k , u H W RF u H ( H k , u j u N u F RF F BB k , j F BB k , j H F RF H H k , u H + σ 2 I N r x ) W RF u W BB k , u .

3. Proposed Hybrid Design Algorithm

In this section, we will introduce the algorithm for the hybrid precoding problem and show how it can be adapted to different architectures. Although we will focus on the precoder design, a similar approach can be adopted for the combiner. However, since our design assumes that inter-user interference suppression is applied at the transmitter, only single-user detection is required at the receiver and therefore the algorithm reduces to the one described in [10].

3.1. Main Algorithm

Although there are several problem formulations for the hybrid design proposed in the literature, one of the most effective relies on the minimization of the Frobenius norm of the difference between the fully-digital precoder and the hybrid precoder [22,29,30,31]. In this paper we follow this matrix approximation-based approach. First, we compute fully-digital precoders, which we assume to be designed so as to enforce zero inter-user interference, using, for example, the block-diagonalization (BD) approach described in [32]. Using the BD procedure, we obtain one different digital precoder matrix for each subcarrier, F opt k (with k = 1 , , F , and F opt k = [ F opt k , 1 , , F opt k , N u ] ), which satisfy H k , u F k , u = 0 for all u u ( u , u = 1 , , N u ), thus guaranteeing no inter-user interference. Using these digital precoder matrices, we design the hybrid digital/analog precoder by solving a matrix approximation problem formulated as
min F RF , F BB k k = 1 F F o p t k F RF F BB k F 2
subject   to   F R F C N t x × N R F t x
F RF F BB k F 2 = N u N s
where Equation (7) enforces the transmitter’s total power constraint and C N t x × N R F is the set of feasible analog precoding matrices, which is defined according to the adopted RF architecture (it will be formally defined for several different architectures in the next subsection). Matrix F opt k denotes the fully-digital precoder. Even if F opt k is selected in order to cancel all interference between users, the hybrid design resulting as a solution of Equations (5)–(7) will correspond to an approximation and, as such, residual inter-user interference will remain. To avoid the performance degradation that will result from this, an additional constraint can be added to the problem formulation, namely
u = 1 u u N u H k , u F RF F BB k , u = 0 ,   k = 1 , , F ,   u = 1 , , N u
where F BB k , u = F BB k [ : , ( u 1 ) N s + 1 : u N s ] . This restriction is equivalent to enforcing F RF F BB k , u to lie in the null space of H ¯ k , u ( N u 1 ) N r x × N t x ( H ¯ k , u is a matrix corresponding to H k with the N r x lines of user u removed) which we denote as N ( H ¯ k , u ) . The overall optimization problem can be then expressed as
min F RF , F BB k k = 1 F F opt k F RF F BB k F 2
subject   to   F RF C N t x × N R F t x
F RF F BB k F 2 = N u N s
F RF F BB k , u N ( H ¯ k , u ) ,   k = 1 , , F ,   u = 1 , , N u .
To derive a hybrid precoder/design algorithm that can cope with the different RF architectures, we can integrate the RF constraint directly into the objective function of the optimization problem. This can be accomplished through the addition of an auxiliary variable, R, combined with the use of the indicator function. The indicator function for a generic set A is defined as I A ( x ) , returning 0 if x A and +∞ otherwise. A similar approach can be adopted for integrating the other constraints, Equations (11) and (12), also into the objective function. The optimization problem can then be rewritten as
min F RF , F BB k , R , B k , F aprox k , u k = 1 F F opt k F RF F BB k F 2 + I C ( R ) + k = 1 F I . F 2 = N u N s ( B k ) + k = 1 F u = 1 N u I N ( H ¯ k , u ) ( F aprox k , u )
subject   to   R = F RF
B k = F RF F BB k
F aprox k , u = F RF F BB k ,
where F aprox k = [ F aprox k , 1 , , F aprox k , N u ] . The augmented Lagrangian function (ALF) for Equations (13)–(16) can be written as
L ρ , η , μ ( F RF , F BB , R , B , F aprox , Λ , Ψ , Γ ) = k = 1 F F opt k F RF F BB k F 2 + I C N t x × N R F t x ( R ) + k = 1 K I . F 2 = N u N s ( B k ) + k = 1 F u = 1 N u I N ( H ¯ k , u ) ( F aprox k , u ) + 2 Re { tr ( Λ H ( F RF R ) + tr ( Ψ k H k = 1 F ( B k + F RF F BB k ) ) + tr ( Γ k H k = 1 F ( F aprox k + F RF F BB k ) ) } + ρ F RF R F 2 + η k = 1 F B k + F RF F BB k F 2 + μ k = 1 F F aprox k + F RF F BB k F 2 ,
where Λ N t x × N R F , Ψ N t x × N R F and Γ N t x × N R F are dual variables and ρ , η , μ are penalty parameters. After some straightforward algebraic manipulation and working with scaled dual variables U = Λ / ρ , W k = Ψ / η and Z k = Γ k / μ we can rewrite the ALF as
L ρ , η , μ ( F RF , F BB , R , B , F aprox , U , W , Z ) = k = 1 K F opt k F RF F BB k F 2 + I C N t x × N R F t x ( R ) + k = 1 F I . F 2 = N u N s ( B k ) + k = 1 F u = 1 N u I N ( H ¯ k , u ) ( F aprox k , u ) + ρ F RF R + U F 2 ρ U F 2 + η k = 1 F B k + F RF F BB k + W k F 2 η k = 1 F W k F 2 + μ k = 1 F F aprox k + F RF F BB k + Z k F 2 μ k = 1 F Z k F 2
In the following, we apply ADMM [22] as a heuristic for solving problem Formulation (13)–(16). To accomplish this, we can apply the gradient ascent to the dual problem involving the ALF, which allows us to obtain an iterative precoding algorithm comprising the following sequence of steps. We start with the minimization of the ALF over F RF for iteration t + 1 defined as
F RF ( t + 1 ) = min F RF L ρ , η , μ ( F RF , F BB ( t ) , R ( t ) , B ( t ) , F aprox ( t ) , U ( t ) , W ( t ) , Z ( t ) ) ,
which can be obtained from
F RF H L ρ , η , μ ( F RF , F BB ( t ) , R ( t ) , B ( t ) , F aprox ( t ) , U ( t ) , W ( t ) , Z ( t ) ) = 0
leading to the closed form expression
F RF ( t + 1 ) = [ k = 0 F 1 [ F opt k + η ( B k ( t ) W k ( t ) ) + μ ( F aprox k ( t ) Z k ( t ) ) ] . F BB ( t ) k H + ρ ( R ( t ) U ( t ) ) ] × [ ( 1 + η + μ ) k = 0 F 1 F BB ( t ) k F BB k ( t ) H + ρ I N R F t x ] 1 ,   k = 0 , , F 1 .
After obtaining the expression for F RF , F BB ( t + 1 ) can be found by following the same methodology. In this case the minimization is expressed as
F BB ( t + 1 ) = min F BB L ρ , η , μ ( F RF ( t ) , F BB , R ( t ) , B ( t ) , F aprox ( t ) , U ( t ) , W ( t ) , Z ( t ) ) ,
which can be obtained from
F BB H L ρ , η , μ ( F RF ( t ) , F BB , R ( t ) , B ( t ) , F aprox ( t ) , U ( t ) , W ( t ) , Z ( t ) ) = 0
and leads to the closed form expression
F B B k ( t + 1 ) = ( 1 + η + μ ) ( F RF ( t + 1 ) H F RF ( t + 1 ) ) 1 F RF ( t + 1 ) H . ( F opt k + η ( B k ( t ) W k ( t ) ) + μ ( F aprox k ( t ) Z k ( t ) ) ) , k = 1 , , F .
The next steps consist of the minimization over R and B k . The minimization of Equation (18) with respect to R and B k can be written as
R ( t + 1 ) = min R { I C N t x × N R F ( R ) + ρ F R F ( t + 1 ) R + U ( t ) F 2 } = C N t x × N R F ( F RF ( t + 1 ) + U ( t ) ) ,
and
B k ( t + 1 ) = min R { I . F 2 = N u N s ( B k ) + η F RF ( t + 1 ) F BB k ( t + 1 ) B k + W ( t ) F 2 } = . F 2 = N u N s ( F RF ( t + 1 ) F BB k ( t + 1 ) + W ( t ) ) , k = 1 , , F ,
where C a × b ( ) and . F 2 = N u N s ( ) denote the projection onto set C a × b and onto the set of matrices whose squared Frobenius norm is N u N s , respectively. While the former projection depends on the adopted analog architecture and will be explained in the next subsection, the second projection is simply computed as
B k ( t + 1 ) = ( F RF ( t + 1 ) F BB k ( t + 1 ) + W ( t ) ) N u N s F RF ( t + 1 ) F BB k ( t + 1 ) + W ( t ) F 2 .
The minimization of (18) with respect of F aprox k , u can be written as
F aprox k , u ( t + 1 ) = min F aprox k , u { I N ( H ¯ k , u ) ( F aprox k , u ) + μ F RF ( t + 1 ) F BB k , ( u 1 ) N s + 1 : u N s ( t + 1 ) F aprox k , u + Z k , ( u 1 ) N s + 1 : u N s ( t ) F 2 } = N ( H ¯ k , u ) ( ( F RF ( t + 1 ) F BB k ( t + 1 ) ) k , ( u 1 ) N s + 1 : u N s + Z k , ( u 1 ) N s + 1 : u N s ( t ) ) ,
which also involves a projection, N ( H ¯ k , u ) ( ) , but in this case onto the null-space of H ¯ k , u . Let us use A to denote A = F RF ( t + 1 ) F BB k , u ( t + 1 ) + Z k , ( u 1 ) N s + 1 : u N s ( t ) . The procedure to compute the projection of matrix A onto the null-space of H ¯ k , u can be formulated as another optimization problem, which can be expressed as
min i = 1 N s A : , i X : , i F 2
subject   to   H ¯ k , u X : , i = 0 .
The general solution for this problem is presented in [29] corresponding to
X : , i = ( I N t x H ¯ k , u H ( H ¯ k , u H ¯ k , u H ) 1 H ¯ k , u ) A : , i , i = 1 , , N s .
Reordering the column vectors in the original matrix form results in the final expression which can be rewritten as
X = ( I N t x H ¯ k , u H ( H ¯ k , u H ¯ k , u H ) 1 H ¯ k , u ) A = ( I N t x V ¯ k , u ( 1 ) ( V ¯ k , u ( 1 ) ) H ) A .
In this expression, V ¯ k , u ( 1 ) denotes the matrix containing the right singular vectors corresponding to the non-zero singular values associated with the singular value decomposition (SVD) given by H ¯ k , u = U ¯ k , u Λ ¯ k , u [ V ¯ k , u ( 1 )   V ¯ k , u ( 0 ) ] H . Therefore, to compute matrix X, one can perform a single value decomposition of H ¯ k , u and then use this to remove the projection of A onto the row space of H ¯ k , u . Finally, the expressions for the update of dual variables U, W and Z are given by
U ( t + 1 ) = U ( t ) + F RF ( t + 1 ) R ( t + 1 ) ,
W k ( t + 1 ) = W k ( t ) + F RF ( t + 1 ) F BB k ( t + 1 ) B k ( t + 1 ) ,
Z k ( t + 1 ) = Z k ( t ) + F RF ( t + 1 ) F BB k ( t + 1 ) F aprox k ( t + 1 ) .
Appropriate values for the penalty parameters can be obtained in a heuristic manner by performing numerical simulations. Regarding the initialization and termination of the algorithm, the same approach described in [10] can be adopted. The whole algorithm is summarized in Table 1. In this table, Q denotes the maximum number of iterations. The projection operation is the only step specific to the implemented architecture, as will be explained in the next subsection.

3.2. Analog RF Precoder/Combiner Structure

The projection required for obtaining matrix R in step 5 of the precoding algorithm has to be implemented according to the specific analog beamformer [6,20,33,34,35,36,37]. This makes the proposed scheme very generic, allowing it to be easily adapted to different RF architectures. In the following, we will consider a broad range of architectures that can be adopted for the RF precoder for achieving reduced complexity and power consumption implementations. We will consider FC, AoSA and DAoSA structures as illustrated in Figure 2, where we assume single phase shifters (SPS). Besides SPS, we will also consider other alternative implementations for these structures, as illustrated in Figure 3 for AoSA. The different solutions either rely on selectors, switches, inverters or phase shifters, or combinations of these. The overall analog structure is defined as a combination of one of the architectures in Figure 2 with either SPS or one of the alternatives illustrated in Figure 3.
(1)
Unquantized Phase Shifters (UPS)
In the first case, we consider the use of infinite resolution phase shifters which, while being ideal, are often used as a reference benchmark. For this architecture the RF constraint set is given by
C a × b = { X a × b : | X i , j | = 1 }
and the corresponding projection can be performed simply using
R ( t + 1 ) = ( F RF ( t + 1 ) + W ( t ) ) | F RF ( t + 1 ) + W ( t ) | ,
where denotes the Hadamard (i.e., element-wise) division.
(2)
Quantized Phase Shifters (QPS)
The second case considers a more realistic scenario, in which phase shifters can be digitally controlled with N b bits. These devices allow the selection of 2 N b different quantized phases and the RF constraint set becomes
C a × b = { X a × b : X i , j = e 2 π k i / 2 N b ,   k = 0 , , 2 N b 1 } .
The implementation of the projection in line 5 of Table 1 can be obtained as the following element-wise quantization
R i , j ( t + 1 ) = e min k = 0 , , 2 N b 1 { angle ( F RF i , j ( t + 1 ) + W i , j ( t ) ) 2 π k / 2 N b } , i = 1 , , N t x , j = 1 , , N R F t x .
Phase shifters are typically one of the best solutions for analog processing blocks but also have a higher implementation cost and power consumption, especially when they have high resolutions.
(3)
Double Phase Shifters (DPS)
Another appealing architecture relies on the use of double phase shifters (DPS) since these remove the constant modulus restriction on the elements of F RF , following the idea in [37]. The main difference between SPS and DPS structures relies on the number of phase shifters in use to compose each connection from an RF chain to a connected antenna element, which in this case is doubled. Even though it increases the implementation complexity and power consumption, this solution can increase the spectral efficiency and approach the performance of the fully-digital one [37].
In this case, the projection can be implemented element-wise simply as
R i , j ( t + 1 ) = ( F RF i , j ( t + 1 ) + W i , j ( t ) ) e i angle ( F RF i , j ( t + 1 ) + W i , j ( t ) ) × max ( 0 , | F RF i , j ( t + 1 ) + W i , j ( t ) | 2 ) .
Similarly to other architectures, DPS can be used not only in the fully-connected approach but also in the AoSA and DAoSA cases, replacing the constant modulus setting operation.
(4)
Switches (Swi)
In an architecture based on switches, each of the variable phase shifters can be replaced by a switch which typically consumes less power [34]. This simplification results in a network of switches connecting each RF chain to the antennas. The RF constraint set can be represented as
C a × b = { X a × b : X i , j = 0   or   X i , j = 1 }
(set of matrices having solely ‘0′s or ‘1′s as elements) and the projection can be implemented element-wise as
R i , j ( t + 1 ) = 1 / 2 + 1 / 2 sign ( 2 Re [ F RF i , j ( t + 1 ) + W i , j ( t ) ] 1 ) .
(5)
Switches and Inverters (SI)
Assuming that N b = 1 , then each variable phase shifter of the QPS architecture can be replaced by a pair of switched lines, including also an inverter. The corresponding constraint set can be reduced to
C a × b = { X a × b : X i , j = ± 1 }
and the implementation of the projection simplifies to
R i , j ( t + 1 ) = sign ( Re [ F RF i , j ( t + 1 ) + W i , j ( t ) ] ) .
(6)
Antenna Selection (AS)
The simplest scenario that we can consider corresponds to an architecture, where each RF chain can be only connected to a single antenna (and vice-versa). Antenna selection is a low-cost low-complexity alternative where only a specified subset of antennas are active at any given time [34]. The RF constraint set will comprise a matrix with only one non-zero element per column and per row, i.e.,
C a × b = { X a × b : X i , j = 0   or   X i , j = 1 , X i , : 0 = 1 , X : , j 0 = 1 } .
In the definition,     0 represents the cardinality of a vector (sometimes referred to as the 0 - norm in the literature). Defining X = F RF ( t + 1 ) + W ( t ) , the projection can be approximately implemented by setting all the elements in X as 0 except for X t j , j = 1 , where t j is the row position with the highest real component in column j:
t j = arg max i = 1 , , N t x { Real [ X i , j ] } .
The computation of t j is performed for all columns j = 1, …, N R F t x , sorted by descending order in terms of highest real components. It should be noted that during this operation, the same row cannot be repeated.
(7)
Array-of-Subarrays (AoSAs)
Within the context of UM-MIMO, one of the most appealing architectures for keeping the complexity acceptable relies on the use of AoSA, where each RF chain is only connected to one or more subsets of antennas (subarrays). Denoting the number of subarrays as n S A , which is typically set as n S A = N R F , and the size of each subarray as N t x S A , then we have N t x S A = N t x n S A = N t x N R F . To limit the complexity of the architecture, each RF chain can connect to a maximum of Lmax consecutive subarrays. In this case, the RF constraint set comprises matrices where each column has a maximum of Lmax blocks of N t x S A constant modulus elements, with all the remaining elements being zero. Defining X = F RF ( t + 1 ) + W ( t ) , the projection can be implemented by setting all the elements in X as 0 except for the subblocks in each column j which fulfill
X { ( [ ( j 1 ) N t x S A + ( i 1 ) N t x S A + 1 : ( j 1 ) N t x S A + i N t x S A ] 1 ) mod N t x } + 1 , j 1 > N t x S A 2
with i = 1 , , Lmax and j = 1 , , N R F . In this case, the corresponding elements of R are set as R i , j ( t + 1 ) = ( X i , j ) / | X i , j | , assuming UPS in these connections. Clearly, the phase shifters can be replaced by any of the other alternatives presented previously.
(8)
Dynamic Array-of-Subarrays (DAoSAs)
As a variation of the previous AoSA architecture, we also consider an implementation where each subarray can be connected to a maximum of Lmax RF chains (which can be non-adjacent). In this case, the constraint set comprises matrices where each N t x S A × N R F component submatrix contains a maximum of Lmax columns with constant modulus elements. The rest of the matrix contains only zeros. In this case, starting with X = 0, the projection can be obtained by selecting the Lmax columns of
X { [ ( j 1 ) N t x S A + ( i 1 ) N t x S A + 1 : ( j 1 ) N t x S A + i N t x S A ] 1 ) mod N t x } + 1 , j 1 > N t x S A 2
where j = 1 , , n S A with the largest 1-norm and setting the corresponding elements of R as R i , j ( t + 1 ) = ( X i , j ) / | X i , j | , assuming the use of UPS. Care must be taken to guarantee that at least one sub-block will be active in every column of R. Similarly to the AoSA, the phase shifters can be replaced by any of the other presented alternatives.

3.3. Complexity

In the proposed algorithm, the F RF ( t + 1 ) and F BB ( t + 1 ) updates (steps 3 and 4 in Table 1) are defined using closed-form expressions that encompass several matrix multiplications, sums and an N R F × N R F matrix inverse (with an assumed complexity order of O ( N R F 3 ) ). These steps require a complexity order of O ( Q N u N s N R F N t x + F 1 Q N R F 2 N t x ) and O ( Q N u N s N R F N t x + Q N R F 2 N t x ) , respectively.
The R ( t + 1 ) update (step 5) involves simple element-wise division (assuming UPS) with O ( Q N R F N t x ) while variable B k ( t + 1 ) (step 6) comprises a Frobenius norm computation with O ( Q N u N s N R F N t x ) . Step 7, the F aprox k , u ( t + 1 ) update, has a complexity order of O ( Q N t x 2 N u N s + N u 3 N t x N r x 2 + N u 4 N r x 3 ) , whereas the dual variables updates (steps 8–10) have a complexity of O ( Q N u N s N R F N t x + Q N R F 2 N t x ) . Therefore, keeping only the dominant terms, the overall complexity order for the proposed algorithm is O ( Q ( N t x 2 N u N s + N R F 2 N t x ) + N u 3 N t x N r x 2 + N u 4 N r x 3 ) . Table 2 presents the total complexity order of the proposed method and compares it against other existing low-complexity alternatives, namely AM–Based [15], LASSO–Based Alt-Min (SPS and DPS) [14] and element-by-element (EBE) [20] algorithms. Taking into account that in UM-MIMO, Ntx will tend to be very large, it means the algorithms with higher complexity will typically be EBE and the one proposed in this paper due to the terms O ( Q N t x 2 ) and O ( Q N t x 2 N u N s ) . It is important to note, however, that while the computational complexity of these two design methods may be higher, both algorithms can be applied to simple AoSA/DAoSA architectures. In particular, the proposed approach directly supports structures with lower practical implementation complexity (and are more energy-efficient) such as those based on switches. Furthermore, in a single-user scenario, the interference cancellation step of the proposed algorithm is unnecessary, and the complexity reduces to O ( Q ( N u N s N R F N t x + N R F 2 N t x ) ) . Regarding the other algorithms, they have similar complexities. However, the AM-based algorithm is designed for single stream scenarios whereas the others consider multiuser multi-stream scenarios.

4. Numerical Results

In this section, the performance of the proposed algorithm will be evaluated and compared against other existing alternatives from the literature, considering multiuser MIMO systems. We consider that both the transmitter and receivers are equipped with uniform planar arrays (UPAs) with N t x × N t x antenna elements at the transmitter and N r x × N r x at the receiver. The respective array response vectors are given by
a t / r ( ϕ i , l , u t / r , θ i , l , u t / r ) = 1 N t x / r x × [ 1 , , e j 2 π λ d ( p sin ϕ i , l , u t / r sin θ i , l , u t / r + q cos θ i , l , u t / r ) , , e j 2 π λ d ( ( N t x / r x 1 ) sin ϕ i , l , u t / r sin θ i , l , u t / r + ( N t x / r x 1 ) cos θ i , l , u t / r ) ] T ,
where p , q = 0 , , N t x / r x 1 are the antenna indices, λ is the signal wavelength and d is the inter-element spacing, which we assume to be d = λ / 2 . We consider a sparse channel with limited scattering where N r a y = 4 and N c l = 6 . The angles of departure and arrival were selected according to a Gaussian distribution whose means are uniformly distributed in [ 0 , 2 π ] and whose angular spreads are 10 degrees. While we include a few results for a NLOS channel, which is often considered a possible scenario in mmWave communications [14,15,25,31], we also present results for a channel with a LOS component which is more realistic, especially in the THz band. In the scenarios with a LOS component, a ratio of E [ | α u L O S | 2 ] / i = 1 N c l l = 1 N r a y E [ | α i , l , u | 2 ] = 10 is assumed (in this case we are admitting very weak NLOS paths compared to LOS which is typical in the THz band [27]). A fully-digital combiner was considered at each receiver and all simulation results were computed with 5000 independent Monte Carlo runs.

4.1. Fully-Connected Structures

First, we evaluate the performance assuming a fully-connected structure. Simulation results for a scenario where a base station with N t x = 100 antennas transmits a single data stream ( N s = 1 ) to N u = 4 users with N r x = 4 antennas are shown in Figure 4 for F = 1 and Figure 5 for F = 64 . The number of RF chains in the transmitter ( N R F t x ) is equal to N u N s . Besides our proposed precoder, several alternative precoding schemes are compared against the fully-digital solution, namely the LASSO-Based Alt-Min, the AM-Based and ADMM-Based precoding [10,14,15].
It can be observed that when F = 1 , only the LASSO-Based Alt-Min with SPS and the ADMM-Based precoder from [10] (which does not remove the inter-user interference) lie far from fully-digital precoder. All the others achieve near optimum results and, in fact, can even match them when adopting DPS (proposed approach and LASSO-based Alt-Min). As explained in Section 2, whereas for F = 1 we have F BB and F RF designed for that specific carrier, when F = 64 , F RF has to be common to all subcarriers. While this reduces the implementation complexity, it also results in a more demanding restriction that makes the approximation of F opt k (problems (5)–(7)) to become worse. Additionally, when this approximation worsens, there can also be increased interference between users. Therefore, it can be observed in the results of Figure 5, that the gap between the fully-digital precoder and all the different hybrid algorithms is substantially wider. Still, the proposed precoder manages to achieve the best results.
Given the performances of the different approaches, it is important to keep in mind that the AM-based precoding algorithm has the lowest performance in wideband but also one of the lowest computational complexities (see Table 2 of Section 3.3). In general, the proposed precoding algorithm is the one that can achieve better results at the cost of some additional computational complexity. Later on, we will address strategies based on lower complexity architectures that will allow for reducing the power consumption associated with its complexity.
In Figure 6, we consider a scenario where the BS employs a larger array with N t x = 256 antennas to transmit N s = 2 simultaneous streams to each user, where N u = 2 . To better fit this scenario to a typical communication in the THz band, we consider the existence of a LOS component, a center frequency of f c = 300   GHz and a bandwidth of B = 15   GHz (it is important to note that the beam split effect is also considered in the channel model). The AM precoder from [15] requires a single stream per user and thus was not included in the figure. It is important to note that in this scenario the use of F = 64 with only 4 RF chains results in a more demanding restriction that makes the approximation to F opt k (Equation (5)) more difficult, thus widening the gap between all the schemes and the fully-digital curve. Still, the LASSO-based Alt-Min precoding schemes present a performance substantially lower when compared to the proposed approaches. Furthermore, the best performance is achieved with the use of double phase shifters, as expected. Once again, in comparing the curves of the proposed precoder against the ADMM-based precoder from [10], the advantage of adopting an interference-cancellation-based design over a simple matrix approximation one is clear.

4.2. Reduced Complexity Architectures

Next, we will focus on the adoption of different reduced complexity architectures according to the typologies presented in Section 3.2. The objective is to evaluate the performance degradation when simpler architectures are adopted.
Figure 7 considers a scenario in which we have more than one data stream ( N s = 2 ) being sent from the BS to each user ( N u = 4 ) in a system with N R F t x = N u N s , F = 1 , N t x = 256 and N r x = 4 . We considered the same penalty parameters configuration: ρ = 0.05 , μ = 1 and η = ρ . This figure is placed in a perspective of simplifying the implementation of the analog precoder but keeping a fully-connected structure. We can observe that the versions based on DPS and single UPS achieve the best results, as expected. Considering the more realistic QPS versions, the results can worsen but it is visible that it is not necessary to use high resolution phase shifters since, with only 3 bits resolution, the results are already very close to the UPS curve. It can also be observed that the simplest of the architectures, AS, results in the worst performance but the SE improves when the antenna selectors are replaced by a network of switches, or, even better, if branches with inverters are also included.
In Figure 8, we intend to simplify the implementation even further with the adoption of AoSAs. In this case we considered that the maximum number of subarrays that can be connected to an RF chain (Lmax) is only one. This imposes a very demanding restriction on matrix F R F since most of it will be filled with zeros, thus substantially deteriorating the approximation to F opt k (Equation (5)). The scenario is the same in Figure 7 but considers the existence of a LOS component with a few weak NLOS paths.
In fact, hereafter, the existence of an LOS component is assumed for the remaining figures of the paper in order to fit the AoSA/DAoSA results to a more typical scenario in the THz band. We can observe that for AoSA structures, the degradation of the SE is notorious, since all candidate versions present worse results when compared to the corresponding fully-connected design and are all far from the fully-digital solution. To reduce the large performance loss due to the adoption of a simple AoSA architecture, we can allow the dynamic connection of more subarrays to each RF chain by adopting a DAoSA structure, as introduced in Section 3.2.
In Figure 9, we study the effect of increasing the maximum number of subarrays that can be connected to an RF chain (Lmax) in the performance of these schemes. Each subarray has a size of 32 antennas (nt). Curves assuming the use of SPS as well as of DPS are included. It can be observed that the increase in the number of connections to subarrays, Lmax, has a dramatic effect on the performance, resulting in a huge improvement by simply going from Lmax = 1 to Lmax = 2. Increasing further to Lmax = 4, the results become close to the fully-connected case, showing that the DAoSA can be a very appealing approach for balancing spectral efficiency with hardware complexity and power consumption. Combining the increase in Lmax with the adoption of DPS can also improve the results but the gains become less pronounced for Lmax > 1. It is important to note that the penalty parameters can be fine-tuned for different system configurations.
One of the objectives of adopting these low-complexity solutions is to reduce overall power consumption. Based on [20], we can calculate the total power consumption of each precoding scheme using
P C = P B B N B B + ( P D A C + P O S + P M ) N R F t x + P P A N t x + P P C N t x + P P S N P S + P S W I N S W I + P t x ,
where PBB is the power of the baseband block (with NBB = 1), PDAC is the power of a DAC, POS is the power of an oscillator, PM is the power of a mixer, PPA is the power of a power amplifier, PPC is the power of a power combiner, PPS is the power of a phase shifter, PSWI is the power of a switch and Ptx denotes the transmit power. The Nx variable represents the number of elements of each device used in the precoder configuration.
Based on the values provided in [20,38] for the power consumption of individual devices in the 300 GHz band, we adopt the following values: PBB = 200 mW, PDAC = 110 mW, POS = 4 mW, PM = 22 mW, PPA = 60 mW, PPC = 6.6 mW, PSWI = 24 mW and PT = 100 mW.
Regarding the phase shifters, we assume values of PPS = 10, 20, 40 and 100 mW for 1, 2, 3 and 4 quantization bits. Considering the same configuration scenario as Figure 7, Figure 8 and Figure 9 with N u = 4 , N s = 2 , N R F t x = N u N s , F = 1 and N t x = 256 , we provide the values of power consumption for different precoder configurations in Table 3.
For the fully-connected structure with UPS, we assumed that PPS = 100 mW, which corresponds to quantized phase shifters with Nb = 4 bits [38]. For the remaining phase-shifter-based precoder structures, we assumed that PPS = 40 mW, which corresponds to quantized phase shifters with Nb = 3 bits, since with only 3 bits resolution the results are already very close to the UPS curve (see Figure 7). As can be seen from this table, the use of architectures based on DAoSAs allows us to reduce considerably the amount of power that is consumed by the precoder. In fact, we can reduce the amount of consumed power up to 55% if we consider a precoder scheme based on DAoSA with DPS and Lmax = 4 versus an FC structure precoder based on UPS, with only a small performance penalty (Figure 9). This reduction increases to 73% if the DPS structure is replaced by an SPS one.
In the particular case of architectures based on quantized phase shifters, we observed that by decreasing the number of quantization bits, it is possible to substantially reduce the power consumption without excessively compromising the complexity (as seen in Figure 7). This conclusion is corroborated by [20,38], since the architectures based on low resolution QPS, AoSAs and DAoSAs present superior energy efficiency when compared to the fully-connected structure with UPS.
In Figure 10 and Figure 11, we provide a comparison between our proposed precoder and the EBE precoder from [20], considering an architecture based on DAoSAs (with SPS) and a scenario configuration similar to Figure 9, i.e., with N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 . These figures present various curves where the maximum number of subarrays that can be connected to an RF chain, Lmax, is changed. Figure 10 refers to an SU scenario ( N u = 1 ) whereas Figure 11 corresponds to an MU scenario with N u = 4 . In the SU case, the proposed precoder achieves results very close to the fully-digital precoder, even with only Lmax = 2. Compared to the proposed algorithm, EBE shows a wider gap even though it has a smaller complexity (as presented in Table 2 of Section 3.3).
When we increase the number of users from N u = 1 to N u = 4 , we can clearly observe that the EBE algorithm suffers a substantial degradation compared to the proposed solution which can be explained by the lack of inter-user interference cancellation (it was not specifically designed for MU scenarios).
Even though a sub-6 GHz system often adopts fully-digital processing [39], where each antenna element has a dedicated RF chain, it is possible to apply the proposed hybrid design algorithm to a sub-6 GHz channel since it is independent of a specific MIMO channel (as are the other alternative algorithms that we used as benchmarks and which are targeted at solving the matrix approximation problem). To exemplify, Figure 12 presents the simulated results obtained for the same scenario of Figure 4 but considering an ideal uncorrelated channel which approximates a rich scattering environment that is typical in sub-6 GHz bands. It can be observed that the proposed approach displays similar behavior to the ones in the upper-bands channel, showing that it can also be used for this particular type of channel (even though it may require a higher number of RF chains to achieve a good approximation of the fully-digital solution in some scenarios, due to the channel not being sparse, as noted in [40]). It is important to highlight that even though the proposed approach can be applied to other channels, the algorithm was designed with the aim of dealing with architectures with a very large number of antennas and with large hardware constraints, making it especially interesting for mmWave and THz.
While we have shown how the proposed approach can deal with several relevant types of analog precoders/combiners, it is important to note that are other alternative structures that have been recently proposed in the literature. For example, some authors have considered precoding paradigms based on time-delayer structures for THz systems [27,41]. One of the most notorious is the Delay Phase Precoding (DPD), which consists in the use of a Time-Delay (TD) network between the RF chains and the traditional phase shifters network in order to convert phase-controlled analog precoding into delay-phase-controlled analog precoding. The main advantage related to this type of precoding is that the time delays in the TD network are carefully designed to generate frequency-dependent beams which are aligned with the spatial directions over the whole bandwidth [41]. While we do not address the adoption of time-delay structures in this paper, it should be possible to derive a projection algorithm that simultaneously takes into account the constraints imposed in both analog-precoding steps: time-delay networks and frequency-independent phase shifters.

5. Conclusions

In this paper, we proposed an iterative algorithm for hybrid precoding design which is suitable for multiuser UM-MIMO systems operating in mmWave and THz bands. The adopted approach replaces unfeasible fully-digital precoders/combiners relying on a dedicated RF chain per antenna with a hybrid architecture comprising low dimension digital blocks with only a few RF chains which are complemented by analog blocks supported solely on networks of phase shifters and switches. The proposed hybrid design algorithm is based on the approximation of the fully-digital approach whose problem formulation is split into a sequence of smaller subproblems with closed-form solutions and can work with a broad range of configuration of antennas, RF chains and data streams. The separability of the design process allows for the adaptability of the algorithm to different architectures, making it suitable for implementation with low-complexity AoSA and DAoSA structures, which are particularly relevant for the deployment of UM-MIMO in hardware-constrained THz systems. It was shown that good trade-offs between SE and hardware implementation complexity can in fact be achieved by the proposed algorithm for several different architectures. Numerical results showed that the use of architectures based on DAoSAs allows us to reduce considerably the amount of power that is consumed at the precoder. In fact, in a reference scenario, we showed that it was possible to reduce the amount of consumed power up to 55% if we consider a precoder scheme based on DAoSAs with DPS and Lmax = 4 versus an FC structure based on UPS, with only a small performance penalty. This reduction increases to 73% if the DPS structure is replaced by an SPS one.

Author Contributions

Conceptualization, J.P.P. and N.S.; methodology, J.P.P., N.S. and R.D.; software, J.P.P. and N.S.; validation, J.P.P., N.S., M.R. and R.D.; formal analysis, J.P.P., V.V., R.F. and N.S.; investigation, J.P.P., V.V., R.F. and N.S.; resources, N.S. and M.R.; data curation, J.P.P.; writing—original draft preparation, J.P.P. and N.S.; writing—review and editing, J.P.P., V.V., R.F., N.S., M.R., J.S. and R.D; visualization, J.P.P. and N.S.; supervision, N.S. and M.R.; project administration, N.S., M.R., J.S. and R.D.; funding acquisition, N.S. and R.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the FCT—Fundação para a Ciência e Tecnologia under the grant 2020.05621.BD. The authors also acknowledge the funding provided by FCT/MCTES through national funds and when applicable co-funded EU funds under the project UIDB/50008/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rappaport, T.S.; Sun, S.; Mayzus, R.; Zhao, H.; Azar, Y.; Wang, K.; Wong, G.; Schulz, J.K.; Samimi, M.; Gutierrez, F. Millimeter Wave Mobile Communications for 5G Cellular: It Will Work! IEEE Access 2013, 1, 335–349. [Google Scholar] [CrossRef]
  2. Uwaechia, A.N.; Mahyuddin, N.M. A Comprehensive Survey on Millimeter Wave Communications for Fifth-Generation Wireless Networks: Feasibility and Challenges. IEEE Access 2020, 8, 62367–62414. [Google Scholar] [CrossRef]
  3. Akyildiz, I.F.; Kak, A.; Nie, S. 6G and Beyond: The Future of Wireless Communications Systems. IEEE Access 2020, 8, 133995–134030. [Google Scholar] [CrossRef]
  4. Tan, J.; Dai, T. THz Precoding for 6G: Applications, Challenges, Solutions, and Opportunities. arXiv 2020, arXiv:2005.10752. Available online: https://arxiv.org/abs/2005.10752 (accessed on 2 July 2021).
  5. Sarieddeen, H.; Alouini, M.-S.; Al-Naffouri, T.Y. An Overview of Signal Processing Techniques for Terahertz Communications. Proc. IEEE 2021, 1–38. [Google Scholar] [CrossRef]
  6. Lin, C.; Li, G.Y.L. Terahertz Communications: An Array-of-Subarrays Solution. IEEE Commun. Mag. 2016, 54, 124–131. [Google Scholar] [CrossRef]
  7. Ahmed, I.; Khammari, H.; Shahid, A.; Musa, A.; Kim, K.S.; De Poorter, E.; Moerman, I. A Survey on Hybrid Beamforming Techniques in 5G: Architecture and System Model Perspectives. IEEE Commun. Surv. Tutor. 2018, 20, 3060–3097. [Google Scholar] [CrossRef] [Green Version]
  8. Chataut, R.; Akl, R. Massive MIMO Systems for 5G and beyond Networks—Overview, Recent Trends, Challenges, and Future Research Direction. Sensors 2020, 20, 2753. [Google Scholar] [CrossRef]
  9. Sohrabi, F.; Yu, W. Hybrid Digital and Analog Beamforming Design for Large-Scale Antenna Arrays. IEEE J. Sel. Top. Signal Process. 2016, 10, 501–513. [Google Scholar] [CrossRef] [Green Version]
  10. Souto, N.N.; Silva, J.; Pavia, J.; Ribeiro, M. An alternating direction algorithm for hybrid precoding and combining in millimeter wave MIMO systems. Phys. Commun. 2019, 34, 165–173. [Google Scholar] [CrossRef] [Green Version]
  11. Pavia, J.P.; Souto, N.; Ribeiro, M.; Silva, J.; Dinis, R. Hybrid Precoding and Combining Algorithm for Reduced Complexity and Power Consumption Architectures in mmWave Communications. In Proceedings of the 2020 IEEE 91st Vehicular Technology Conference: VTC2020-Spring, Antwerp, Belgium, 25–28 May 2020; pp. 1–5. [Google Scholar] [CrossRef]
  12. Guan, K.; Li, G.; Kuerner, T.; Molisch, A.F.; Peng, B.; He, R.; Hui, B.; Kim, J.; Zhong, Z. On Millimeter Wave and THz Mobile Radio Channel for Smart Rail Mobility. IEEE Trans. Veh. Technol. 2016, 66, 5658–5674. [Google Scholar] [CrossRef]
  13. Alkhateeb, A.; Mo, J.; Gonzalez-Prelcic, N.; Heath, R. MIMO Precoding and Combining Solutions for Millimeter-Wave Systems. IEEE Commun. Mag. 2014, 52, 122–131. [Google Scholar] [CrossRef]
  14. Yu, X.; Zhang, J.; Letaief, K.B. Alternating minimization for hybrid precoding in multiuser OFDM mmWave systems. In Proceedings of the 2016 50th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 6–9 November 2016; pp. 281–285. [Google Scholar] [CrossRef] [Green Version]
  15. Yuan, H.; An, J.; Yang, N.; Yang, K.; Duong, T.Q. Low Complexity Hybrid Precoding for Multiuser Millimeter Wave Systems Over Frequency Selective Channels. IEEE Trans. Veh. Technol. 2018, 68, 983–987. [Google Scholar] [CrossRef]
  16. Liu, F.; Kan, X.; Bai, X.; Du, R.; Liu, H.; Zhang, Y. Hybrid precoding based on adaptive RF-chain-to-antenna connection for millimeter wave MIMO systems. Phys. Commun. 2020, 39, 100997. [Google Scholar] [CrossRef]
  17. Xu, K.; Cai, Y.; Zhao, M.; Niu, Y.; Hanzo, L. MIMO-Aided Nonlinear Hybrid Transceiver Design for Multiuser Mmwave Systems Relying on Tomlinson-Harashima Precoding. IEEE Trans. Veh. Technol. 2021, 70, 6943–6957. [Google Scholar] [CrossRef]
  18. Vizziello, A.; Savazzi, P.; Chowdhury, K.R. A Kalman Based Hybrid Precoding for Multi-User Millimeter Wave MIMO Systems. IEEE Access 2018, 6, 55712–55722. [Google Scholar] [CrossRef]
  19. Elmagzoub, H.M. On the MMSE-based multiuser millimeter wave MIMO hybrid precoding design. Int. J. Commun. Syst. 2020, 33, e4409. [Google Scholar] [CrossRef]
  20. Yan, L.; Han, C.; Yuan, J. A Dynamic Array-of-Subarrays Architecture and Hybrid Precoding Algorithms for Terahertz Wireless Communications. IEEE J. Sel. Areas Commun. 2020, 38, 2041–2056. [Google Scholar] [CrossRef]
  21. Busari, S.; Huq, K.; Mumtaz, S.; Rodriguez, J. Terahertz Massive MIMO for Beyond-5G Wireless Communication. In Proceedings of the ICC 2019-2019 IEEE International Conference on Communications (ICC), Shanghai, China, 11–22 March 2019. [Google Scholar]
  22. Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers. Found. Trends Mach. Learn. 2010, 3, 1–122. [Google Scholar] [CrossRef]
  23. Wang, Y.; Yin, W.; Zeng, J. Global Convergence of ADMM in Nonconvex Nonsmooth Optimization. J. Sci. Comput. 2018, 78, 29–63. [Google Scholar] [CrossRef] [Green Version]
  24. Alkhateeb, A.; Heath, R.W. Frequency Selective Hybrid Precoding for Limited Feedback Millimeter Wave Systems. IEEE Trans. Commun. 2016, 64, 1801–1818. [Google Scholar] [CrossRef] [Green Version]
  25. Yuan, H.; Yang, N.; Yang, K.; Han, C.; An, J. Hybrid Beamforming for Terahertz Multi-Carrier Systems Over Frequency Selective Fading. IEEE Trans. Commun. 2020, 68, 6186–6199. [Google Scholar] [CrossRef]
  26. Tarboush, S.; Sarieddeen, H.; Chen, H.; Loukil, M.H.; Jemaa, H.; Alouini, M.S.; Al-Naffouri, T.Y. TeraMIMO: A Channel Simulator for Wideband Ultra-Massive MIMO Terahertz Communications. arXiv 2021, arXiv:2104.11054. Available online: https://arxiv.org/abs/2104.11054 (accessed on 2 July 2021).
  27. Lin, C.; Li, G.Y.; Wang, L. Subarray-Based Coordinated Beamforming Training for mmWave and Sub-THz Communications. IEEE J. Sel. Areas Commun. 2017, 35, 2115–2126. [Google Scholar] [CrossRef]
  28. Nguyen, D.D.; Le, L.; Le-Ngoc, T.; Heath, R. Hybrid MMSE Precoding and Combining Designs for mmWave Multiuser Systems. IEEE Access 2017, 5, 19167–19181. [Google Scholar] [CrossRef]
  29. Bertsekas, D. Nonlinear Programming; Athena Scientific: Belmont, MA, USA, 2016. [Google Scholar]
  30. El Ayach, O.; Rajagopal, S.; Abu-Surra, S.; Pi, Z.; Jr, R.W.H. Spatially Sparse Precoding in Millimeter Wave MIMO Systems. IEEE Trans. Wirel. Commun. 2014, 13, 1499–1513. [Google Scholar] [CrossRef] [Green Version]
  31. Yu, X.; Shen, J.-C.; Zhang, J.; Letaief, K.B. Alternating Minimization Algorithms for Hybrid Precoding in Millimeter Wave MIMO Systems. IEEE J. Sel. Top. Signal Process. 2016, 10, 485–500. [Google Scholar] [CrossRef] [Green Version]
  32. Spencer, Q.; Swindlehurst, A.L.; Haardt, M. Zero-Forcing Methods for Downlink Spatial Multiplexing in Multiuser MIMO Channels. IEEE Trans. Signal Process. 2004, 52, 461–471. [Google Scholar] [CrossRef]
  33. Mendez-Rial, R.; Rusu, C.; Gonzalez-Prelcic, N.; Alkhateeb, A.; Heath, R. Hybrid MIMO Architectures for Millimeter WaveCommunications: Phase Shifters or Switches? IEEE Access 2016, 4, 247–267. [Google Scholar] [CrossRef]
  34. Lee, J.; Lee, Y.H. AF relaying for millimeter wave communication systems with hybrid RF/baseband MIMO processing. In Proceedings of the 2014 IEEE International Conference on Communications (ICC), Sydney, Australia, 10–14 June 2014; pp. 5838–5842. [Google Scholar] [CrossRef]
  35. Payami, S.; Ghoraishi, M.; Dianati, M.; Sellathurai, M. Hybrid Beamforming with a Reduced Number of Phase Shifters for Massive MIMO Systems. IEEE Trans. Veh. Technol. 2018, 67, 4843–4851. [Google Scholar] [CrossRef] [Green Version]
  36. Tian, M.; Zhang, J.; Zhao, Y.; Yuan, L.; Yang, J.; Gui, G. Switch and Inverter Based Hybrid Precoding Algorithm for mmWave Massive MIMO System: Analysis on Sum-Rate and Energy-Efficiency. IEEE Access 2019, 7, 49448–49455. [Google Scholar] [CrossRef]
  37. Yu, X.; Zhang, J.; Letaief, K.B. Doubling Phase Shifters for Efficient Hybrid Precoder Design in Millimeter-Wave Communication Systems. J. Commun. Inform. Netw. 2019, 4, 51–67. [Google Scholar]
  38. Yan, L.; Han, C.; Yang, N.; Yuan, J. Dynamic-subarray with Quantized- and Fixed-phase Shifters for Terahertz Hybrid Beamforming. In Proceedings of the GLOBECOM 2020-2020 IEEE Global Communications Conference, Taipei, Taiwan, 7–11 December 2020; pp. 1–6. [Google Scholar] [CrossRef]
  39. Li, Z.; Zhang, C.; Lu, I.-T.; Jia, X. Hybrid Precoding Using Out-of-Band Spatial Information for Multi-User Multi-RF-Chain Millimeter Wave Systems. IEEE Access 2020, 8, 50872–50883. [Google Scholar] [CrossRef]
  40. Park, S.; Alkhateeb, A.; Heath, R.W., Jr. Dynamic Subarrays for Hybrid Precoding in Wideband mmWave MIMO Systems. IEEE Trans. Wirel. Commun. 2017, 16, 2907–2920. [Google Scholar] [CrossRef]
  41. Tan, J.; Dai, L. Delay-Phase Precoding for THz Massive MIMO with Beam Split. In Proceedings of the 2019 IEEE Global Communications Conference (GLOBECOM), Big Island, HI, USA, 9–13 December 2019; pp. 1–6. [Google Scholar] [CrossRef]
Figure 1. A multiuser OFDM mmWave/THz MIMO system with hybrid precoding.
Figure 1. A multiuser OFDM mmWave/THz MIMO system with hybrid precoding.
Sensors 21 06054 g001
Figure 2. Different precoder architectures for a mmWave/THz MIMO system based on phase shifters.
Figure 2. Different precoder architectures for a mmWave/THz MIMO system based on phase shifters.
Sensors 21 06054 g002
Figure 3. Alternative implementations to single phase shifters based on array-of-subarrays for a mmWave/THz MIMO system.
Figure 3. Alternative implementations to single phase shifters based on array-of-subarrays for a mmWave/THz MIMO system.
Sensors 21 06054 g003
Figure 4. Spectral efficiency versus SNR achieved by different methods with N u = 4 , N s = 1 , N R F t x = 4 , F = 1 , N t x = 100 and N r x = 4 (only NLOS).
Figure 4. Spectral efficiency versus SNR achieved by different methods with N u = 4 , N s = 1 , N R F t x = 4 , F = 1 , N t x = 100 and N r x = 4 (only NLOS).
Sensors 21 06054 g004
Figure 5. Spectral efficiency versus SNR achieved by different methods with N u = 4 , N s = 1 , N R F t x = 8 , F = 64 , N t x = 100 and N r x = 4 (only NLOS).
Figure 5. Spectral efficiency versus SNR achieved by different methods with N u = 4 , N s = 1 , N R F t x = 8 , F = 64 , N t x = 100 and N r x = 4 (only NLOS).
Sensors 21 06054 g005
Figure 6. Spectral efficiency versus SNR achieved by different methods with N u = 2 , N s = 2 , N R F t x = 4 , F = 64 , N t x = 256 and N r x = 4 (with LOS component).
Figure 6. Spectral efficiency versus SNR achieved by different methods with N u = 2 , N s = 2 , N R F t x = 4 , F = 64 , N t x = 256 and N r x = 4 (with LOS component).
Sensors 21 06054 g006
Figure 7. Spectral efficiency versus SNR achieved by the proposed precoder using different fully-connected architectures for N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (only NLOS).
Figure 7. Spectral efficiency versus SNR achieved by the proposed precoder using different fully-connected architectures for N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (only NLOS).
Sensors 21 06054 g007
Figure 8. Spectral efficiency versus SNR achieved by the proposed precoder using different AoSA architectures with Lmax = 1, N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Figure 8. Spectral efficiency versus SNR achieved by the proposed precoder using different AoSA architectures with Lmax = 1, N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Sensors 21 06054 g008
Figure 9. Spectral efficiency versus SNR achieved by the proposed precoder considering an architecture based on DAoSAs and the variation of the maximum number of subarrays that can be connected to an RF chain ( L m a x ) for N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Figure 9. Spectral efficiency versus SNR achieved by the proposed precoder considering an architecture based on DAoSAs and the variation of the maximum number of subarrays that can be connected to an RF chain ( L m a x ) for N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Sensors 21 06054 g009
Figure 10. Spectral efficiency versus SNR achieved by the proposed precoder and by the EBE algorithm considering an architecture based on DAoSAs and the variation of the maximum number of subarrays that can be connected to an RF chain (Lmax) for N u = 1 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Figure 10. Spectral efficiency versus SNR achieved by the proposed precoder and by the EBE algorithm considering an architecture based on DAoSAs and the variation of the maximum number of subarrays that can be connected to an RF chain (Lmax) for N u = 1 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Sensors 21 06054 g010
Figure 11. Spectral efficiency versus SNR achieved by the proposed precoder and by the EBE algorithm considering an architecture based on DAoSAs and the variation of the maximum number of subarrays that can be connected to an RF chain (Lmax) for a mmWave/THz system with N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Figure 11. Spectral efficiency versus SNR achieved by the proposed precoder and by the EBE algorithm considering an architecture based on DAoSAs and the variation of the maximum number of subarrays that can be connected to an RF chain (Lmax) for a mmWave/THz system with N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 and N r x = 4 (with LOS component).
Sensors 21 06054 g011
Figure 12. Spectral efficiency versus SNR achieved by different methods for a mmWave/THz MIMO-OFDM system with N u = 4 , N s = 1 , N R F t x = 4 , F = 1 , N t x = 100 and N r x = 4 considering an uncorrelated channel.
Figure 12. Spectral efficiency versus SNR achieved by different methods for a mmWave/THz MIMO-OFDM system with N u = 4 , N s = 1 , N R F t x = 4 , F = 1 , N t x = 100 and N r x = 4 considering an uncorrelated channel.
Sensors 21 06054 g012
Table 1. General Iterative Hybrid Design Algorithm.
Table 1. General Iterative Hybrid Design Algorithm.
1:  Input:  F opt k , F RF ( 0 ) , F BB k ( 0 ) , R ( 0 ) , B k ( 0 ) , F aprox k , u ( 0 ) , ρ, Q
2:  for t = 0, 1, …, Q − 1 do
3:  Compute F RF ( t + 1 ) using (21).
4:  Compute F BB k ( t + 1 ) using (24), for all k = 1, …, F.
5:  Compute R ( t + 1 ) using (25).
6:  Compute B k ( t + 1 ) using (26), for all k = 1, …, F.
7:  Compute F aprox k , u ( t + 1 ) using (28), for all k = 1, …, F and u = 1, …, Nu.
8:  Update U ( t + 1 ) using (33).
9:  Update W k ( t + 1 ) using (34), for all k = 1, …, F.
10:  Update Z k ( t + 1 ) using (35), for all k = 1, …, F.
11:  end for.
12:   F ^ RF R ( Q ) .
13:   F ^ BB k ( F ^ RF H F ^ RF ) 1 F ^ RF H F aprox k ( Q ) , for all k = 1, …, F.
14:   F ^ BB k N u N s F ^ BB k H F ^ BB k F 1 F ^ BB k .
15:  Output:  F ^ RF , F ^ BB .
Table 2. Overall Complexity of Different Hybrid Precoding Algorithms (per subcarrier).
Table 2. Overall Complexity of Different Hybrid Precoding Algorithms (per subcarrier).
AM—Based
OperationComplexity Order
Overall [15] O ( Q ( N u N s N R F N t x + N R F 2 N u N s + F 1 N R F 3 ) F 1 N R F 2 N t x + N u 3 N s 3 )
LASSO—Based Alt-Min (SPS)
OperationComplexity Order
Overall [14] O ( Q ( N u N s N R F N t x + N R F 2 N u N s + F 1 N R F 3 ) + N u 2 N s N R F N t x + N u 4 N s 3 )
ADMM
OperationComplexity Order
Overall [10] O ( Q ( N s N R F N t x + N R F 2 N t x ) )
EBE
OperationComplexity Order
Overall [20] O ( Q N t x 2 )
Proposed
OperationComplexity Order
F RF O ( Q N u N s N R F N t x + F 1 Q N R F 2 N t x )
F BB O ( Q N u N s N R F N t x + Q N R F 2 N t x )
R O ( Q N R F N t x )
B O ( Q N u N s N R F N t x )
F aprox O ( Q N t x 2 N u N s + N u 3 N t x N r x 2 + N u 4 N r x 3 )
U, W, Z O ( Q N u N s N R F N t x + Q N R F 2 N t x )
Overall O ( Q ( N t x 2 N u N s + N R F 2 N t x ) + N u 3 N t x N r x 2 + N u 4 N r x 3 )
Table 3. Power Consumption for Different Implementations of the Proposed Precoder for N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 .
Table 3. Power Consumption for Different Implementations of the Proposed Precoder for N u = 4 , N s = 2 , N R F t x = 8 , F = 1 , N t x = 256 .
PrecoderEstimated Power Consumption [W]
Fully-ConnectedDPS428.04
UPS223.24
QPS (Nb = 2)59.4
QPS (Nb = 3)100.36
SWI67.59
SI38.92
DAoSA SPSLmax = 128.87
Lmax = 239.30
Lmax = 349.73
Lmax = 460.17
DAoSA DPSLmax = 139.11
Lmax = 259.78
Lmax = 380.45
Lmax = 4101.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pavia, J.P.; Velez, V.; Ferreira, R.; Souto, N.; Ribeiro, M.; Silva, J.; Dinis, R. Low Complexity Hybrid Precoding Designs for Multiuser mmWave/THz Ultra Massive MIMO Systems. Sensors 2021, 21, 6054. https://doi.org/10.3390/s21186054

AMA Style

Pavia JP, Velez V, Ferreira R, Souto N, Ribeiro M, Silva J, Dinis R. Low Complexity Hybrid Precoding Designs for Multiuser mmWave/THz Ultra Massive MIMO Systems. Sensors. 2021; 21(18):6054. https://doi.org/10.3390/s21186054

Chicago/Turabian Style

Pavia, João Pedro, Vasco Velez, Renato Ferreira, Nuno Souto, Marco Ribeiro, João Silva, and Rui Dinis. 2021. "Low Complexity Hybrid Precoding Designs for Multiuser mmWave/THz Ultra Massive MIMO Systems" Sensors 21, no. 18: 6054. https://doi.org/10.3390/s21186054

APA Style

Pavia, J. P., Velez, V., Ferreira, R., Souto, N., Ribeiro, M., Silva, J., & Dinis, R. (2021). Low Complexity Hybrid Precoding Designs for Multiuser mmWave/THz Ultra Massive MIMO Systems. Sensors, 21(18), 6054. https://doi.org/10.3390/s21186054

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