Abstract

Massive multiple-input multiple-output (MIMO) systems based on millimeter-wave (mmWave) technology commonly utilize a combination of hybrid analog and digital signal processing to reduce the number of radiofrequency (RF) chains, cutting down the hardware costs and power consumption. The subconnected structure is one of the architectures that realize hybrid analog and digital processing, in which each RF chain is connected to a subset of antennas through phase shifters so as to reduce the implementation complexity. In this paper, we investigate the best subconnected structure design for mmWave massive MIMO systems. First, we formulate the optimization problem based on a simplified connection model. Next, we derive a suboptimal solution in closed form by leveraging the mmWave channel characteristics and discuss the impact of different connection parameters on the spectral efficiency. The spectral efficiency is closely related to the structure of the analog processing and the interantenna spacing. When the interantenna spacing is an even a multiple of half wavelength, the localized structure, i.e., achieves the maximum spectral efficiency. When the interantenna spacing is an odd multiple of half wavelength and the analog processing is accomplished by interleaved structure, the best performance is obtained. Numerical results demonstrate the correctness of the theoretical analysis.

1. Introduction

In the future generation of mobile communication systems, millimeter-wave (mmWave) combined with massive multiple-input multiple-output (MIMO) has been intensively envisaged as a promising technology [15], because of sufficient available spectrum resources which help to obtain higher data rates. Meanwhile, the high beamforming gains attributed to massive MIMO could overcome the extreme path loss of mmWave communications [6, 7]. However, the traditional full-digital structure that connects each antenna with a radio-frequency (RF) chain is not desirable for mmWave massive MIMO systems because of the high hardware cost and power consumption [8, 9].

By linking the antennas to far fewer RF chains via a phase shifter network, hybrid analog and digital signal processing provides a viable approach for mmWave massive MIMO communication, which is accomplished in the analog domain followed by baseband processing in the digital domain [10]. The system performance of the hybrid structure can be close to that of the full-digital structure [1113].

According to the connection topology between antennas and RF chains, the hybrid structure is divided into two types in general. For the fully connected structure, each RF chain is linked to all antennas. In contrast, for subconnected structures, an RF chain is only connected to a subarray of the total available antennas [14]. The performance of the subconnected structure is slightly lower than that of the fully connected structure, but it can be more energy-efficient by reducing the number of phase shifters [14, 15]. The subconnected structure achieves a trade-off between the system performance and hardware complexity. Based on the configuration of the antenna subarray, the subconnected structure is classified into two categories, localized structure and interleaved structure. For the localized structures, the elements in an antenna subarray are adjacent and consecutive [1417], while they are separated in the interleaved structure.

The interleaved structure is one of the most popular structures [8, 1828]. It is proved that the angle of arrival (AoA) estimation of interleaved structure could converge faster with lower complexity [2931]. The interleaved structure can eliminate the grating lobes as well as achieve multiple beams with high gain by adopting interleaved coding schemes [1820], [24], hence to reduce interference and achieve high signal to interference ratio (SIR). Moreover, the interleaved structure could also achieve high-speed communication even in densely populated areas [20].

The interleaved structure is first proposed in [29]. In [18, 19, 2125] and [2932], the authors studied the fixed structures interleaved with a constant number of elements. In [20, 26], arbitrary interleaved structures are proposed. In [27, 28], the antenna subarrays of the interleaved structures have the same size and are not allowed to overlap, which is different from the overlapped structure. In [18, 25], some prototype of interleaved structure in mmWave frequencies is realized through a phased array chip. In [31, 33], the spectral efficiency of the localized structure and the interleaved structure are compared. In [33], the performance of the interleaved structure with constant configuration of subarray is developed. The authors of [31] consider the influence of subarray configuration on the spectral efficiency. It shows that different configurations significantly affect the overall performance. However, existing works do not provide results on the optimal interleaved structure and the impact of connection parameters on the system performance.

In this paper, we consider the uplink of the mmWave massive MIMO system as shown in Figure 1, where the hybrid analog and digital signal processing is adopted. We compare the spectral efficiency of diverse subconnected structures, study the main parameters that affect the system performance, and derive the optimal structure in closed form. The main contributions of this paper are as follows: (i)We derive the suboptimal interleaved structure in closed form for channels with a single propagation path, and then the results are extended to the case with multiple propagation paths based on the properties of mmWave channel(ii)We propose an efficient algorithm to find the optimal interleaved structure that maximizes the spectral efficiency(iii)We study the impact of interleave factor () and interantenna spacing () on the spectral efficiency. The interleaved structure with achieves the best performance when the interantenna spacing is an odd multiple of half wavelength. The localized structure achieves the best performance when the interantenna spacing is an even a multiple of half wavelength

The remainder of the paper is organized as follows. In Section 2, we give the system model and formulate the optimization problem for finding the optimal subconnected structure. In Section 3, we obtain the optimal interleaved structure by leveraging the characteristics of the mmWave channel. In Section 4, numerical results are provided. Finally, Section 5 draws conclusions.

Notations: vectors and matrices are denoted by boldface lower and uppercase letters, respectively. represents a set. denotes the conjugate transpose, and denotes the inverse of a matrix. Expectation and determinant operators are denoted by and , respectively. The diagonal matrix is denoted by with diagonal entries represented by , . denotes an identity matrix. denotes a complex Gaussian random vector, where the mean is and the covariance matrix is .

2. System Model and Problem Formulation

In this section, we describe the signal model of the mmWave massive MIMO system and the connection between RF chains and antennas. Then, the optimization problem of the subconnected structure is formulated. Consider the uplink of a multiuser mmWave massive MIMO system which comprises of a base station (BS) with a uniform linear array (ULA) of antennas and single antennas utilizing resources at the same time and frequency. The considered hybrid processing subconnected structure at the BS is illustrated in Figure 1, which has two successive subprocedures: analog processing and baseband processing. The antennas and the RF chains are connected arbitrarily through a phase shifter network which enables analog processing. The output of the RF chains is sent to the baseband for digital processing. To ensure the multiple data stream communication, we assume that the number of RF chains is equal to users, i.e., . The received signal vector at the BS can be presented as follows:

where denotes the channel matrix from users to the BS; , with being the transmit power of the -th user, ; is the transmit signal vector, and is the additive white Gaussian noise (AWGN) vector, which follows the distribution of . Consider the large-scale fading, and the channel is modeled as

where , and is the large-scale fading coefficient of the -th user to the BS; is composed of the small-scale fading channel coefficients of all users, and it can be presented as

where is the channel vector of the -th user to the BS. Different from the low-frequency band, the number of scatters in the mmWave propagation environment is limited, and the correlation between different antennas is high. Therefore, the expanded Saleh-Valenzuela channel model in [9] for the mmWave propagation environment is adopted. Thus, can be modeled as

where is the number of valid transmission paths from the -th user to the BS; and are the complex gain and the AoA of the -th propagation path of the -th user, respectively. For the ULA, the steering vector is given by

where is the distance between adjacent antennas and , and is the wavelength. Because of the near-far effect, the power of different users at the BS is different. To obtain the consistent received signal power , a power control scheme is adopted in the uplink [34]. Hence, (1) can be rewritten as

To recover at the BS with low power consumption and hardware cost, we focus on the hybrid processing scheme which has gained great research interest [10, 1417]. As shown in Figure 1, the received signal vector after analog processing can be presented as

where is the analog processing matrix, which is realized by the phase shifter network and is the equivalent RF channel. Denote as the baseband processing matrix, and the received signal after baseband processing is given by

The optimal connection between the RF links and the antenna array that maximizes the sum rate of all users is very difficult to derive, due to its complicated expression. Besides, the connections only affect the analog processing. To simplify the problem at hand, we denote the mutual information of the link as . Based on information theory, we have

where is the channel capacity for channel matrix . is the upper bound of the sum rate . Therefore, the optimization problem for is altered to optimize , which is given by

where is the noise covariance matrix after analog processing.

From (10), it is apparent that different choices of greatly affect the system performance. Note that the elements of are related to the way how antennas are connected to RF chains by the phase shifter network. To find the optimal connection, we will focus on the analog processing matrix . For the sake of describing the case of connection between the antenna and the RF chains, inspired by Zhu et al. [35], the -th element of can be modeled as

where and . It is worth mentioning that in general, the elements of can be any complex value with unit norm (i.e., ). To concentrate on the connection structure, they are assumed to be binary here. To describe the arbitrary connection between the antenna array and the RF connections, we use the connection model (11). For example, when the whole elements in are equal to one, it is actually a fully connected structure. When is a block diagonal matrix, it represents the localized structure [1417]. In this work, we aim to optimize the interleaved structure which has gained great research interest and is shown to provide superior performance [20, 29, 31].

As shown in Figure 2, it is the interleaved structure. antennas are uniformly parted into antenna subarrays, so that antennas belong to one subarray and connect to the same RF chain. We denote as the interleaved factor. Limited to the structure, ought to be a divisor of . There are antennas between two consecutive antennas of the same subarray. Therefore, the distance of them is . To elaborate on the interleaved structure from another perspective, RF chains and the antennas connected to them form a group. Therefore, the antenna arrays are parted into groups. In Figure 2, is the index of the first RF chain of corresponding to the last group, and the value of can be represented as . From Figure 2, we know that each subarray’s elements are dispersed through the antenna array group. When , it means that two consecutive antennas in the same subarray are separated by antenna element, which is actually the localized structure [31].

Diverse interleaved structures could be described by the aforementioned model. Solving the following problem can yield the optimal structure, i.e., the optimal .

where , and (a) indicates that there are elements equal to one in each column of . In (b), represents the index of the -th -valued element in the -th column, and . Constraint b) means that two continuous -valued elements in a column are separated by elements.

The optimization objective involves expectation, (12) is discrete, resulting in a nonconvex optimization problem, and thus the globally optimal solution of (12) is very difficult.

3. Suboptimal Design of the Interleaved Structure

In this section, we derive a closed form suboptimal solution to (12) by leveraging the properties of the mmWave channels.

According to (11), we have , and thus the mutual information can be written as

Applying Jensen’s inequality results in

Problem (12) is changed to optimize the upper bound and we have.

by taking advantage of the characteristics of the mmWave channel, the -th element of in (14) is given by

where and can be any integer from to .Then, we write (16)

where denotes the largest integer not more than , and denotes to the remainder of . Since (16) is very complicated and contains multiple propagation paths, it is not easy to solve (15) directly. Therefore, we first consider , which means that there is only one effective propagation path available for each user, and then the results are extended to cases with .

3.1. Solution for

Assuming and from (16), we have

By some mathematical manipulations, can be rewritten as

where and the -th element is given by

In (19), is a diagonal matrix, and the -th diagonal element is provided by

Substituting (19) into (14) yields

where and follow . Note that generally. However, a characteristic of large random matrices states that the columns of become asymptotically orthogonal when is large, which results in . The average correlation coefficient of two different columns of decreases as the amount of increases, as in Figure 3. Hence, can be further simplified with (19) as

where , and is given by

where follows

and is obtained by using

Substituting (24) into (23), is given by

in which is simplified to be a function of and . For given and , (15) can be simplified into the following problem:

Then, (28) is equivalent to

Since contains a summation of weighted Bessel functions, the recursive formula of zero-order Bessel function is used to find the optimal value of . The derivative of is obtained as

Denote as the solution to , and we obtain the following theorem.

Theorem 1. When the system has a large number of users and each user has only one valid propagation path , the optimal interleaved structure is satisfied .

provides a sufficient condition for the optimal interleaved structure. However, the solution of is not easy to obtain in general. Therefore, to gain more insights, we analyze some special cases in the following corollaries.

Corollary 2. In a particular case, when each subarray contains only two antenna elements, i.e.,, supposing that , the interleaved structure that maximizes the spectral efficiency requires .

Substituting into (30) and setting it to zero, we have

Based on the asymptotic expansion of the first order Bessel function , the zero point is , where Since , the solution to (31) is approximated by

When is even, is the locally maximum point of the function . Because the first-order Bessel function varies and disappears, achieves the maximum when , i.e., . However, considering the actual situation, is the divisor of user number, the interleaved structure of can be realized, and the value of is the largest at this time.

Corollary 3. In more general cases, for that is an integer greater than , when , the interleaved structures that maximizes the spectral efficiency should satisfy .

Although the first-order Bessel function is not periodic, the zero points of the first-order Bessel function tend to be periodic, i.e., where is the nonnegative zero point of , When is small, the distance between two adjacent zero points is close to . Therefore, its period is infinitely close to . Assuming , the period of is approximately equal to . Note that is a summation of multiple weighted first-order Bessel functions. According to the addition principle of the periodic function, the period of the sum is equal to the least common multiple of the period of each term in the summation. Therefore, the approximate period of is . Besides, according to the expansion of the Bessel function, could make each term, i.e., in equals zero. In other words, is the solution of . The locally maximum point of is obtained. Limited to the interleaved structure, ought to be .

For the case of multiple effective propagation paths, more results are provided in the following subsection.

3.2. Solution for

Consider that each user’s channel contains multiple valid propagation paths, i.e., . According to (16), can be expressed as

where represents the component corresponding to the -th path in . Similar to the case when , is decomposed into

where and the -th element is defined as

Moreover, is a diagonal matrix, and the -th element of is

To simplify the derivation, in (14) can be further deduced to

According to (34) and (35), can be rewritten as

where

where the first part is the component of the same path in the equivalent channel , and the second part is the mutual interference between different paths. Since different paths are independent of each other, when the number of antennas is large, the mutual interference between different paths is approximately to be zero. Then, is approximated by

Since , it yields

It is worth pointing out that the approximation from (41) to (42) is valid when , goes to infinity, and it becomes less accurate for the limited value of and . However, this approximation simplifies the analysis, and the obtained optimal interleaved structure works well under moderate values of and . According to (37) and (42), for given , is rewritten as

where , and is defined as

in which is given by

It can be seen that is independent of , while it is related to . (44) is rewritten as

Note that in (46) is equal to in (24). With (43) and (46), the objective function is eventually simplified to

It is observed that is related to , and then the optimization problem (15) is simplified as

Because the optimization objective is a nonconvex function of , the optimization problem is complicated to solve. In this work, we propose to solve (48) by Algorithm 1. It calculates and compares the values of to obtain the optimal parameter that maximizes .

Input:, , , and ;
1: Calculate the number of antennas in one Subarray
;
2:Calculate the divisors of and denote as a set.
, where is the number of the
 divisors of ;
3: For do:
4:  Calculate ;
5:  ;
6: End;
7:Find ;
8: the index of ;
9: The optimal interleaved structure parameter is
Output:.

4. Numerical Results

In this section, the numerical results are presented to demonstrate the performance of various subconnected structures in terms of the spectral efficiency. In simulations, the number of BS antennas is , and the number of RF chains is or . Two cases of and are considered, respectively. The AoAs are assumed to be uniformly distributed in . is referred to as signal-to-noise ratio (SNR).

Under the condition that , when the RF chains is , the achievable spectral efficiency performance of diverse interleaved structures with variational SNR for and is illustrated in Figure 4(a) and 4(b), respectively. As shown in Figure 4(a), the interleaved structure obtains the maximum spectral efficiency when , which is followed by in descending order of spectral efficiency. This verified the analysis in Section 3.1. The interleaved structure with is actually the localized structure, which performs the worst. The same results can be found in Figure 5 where the number of RF chain in 8. Besides, the gap between interleaved and localized structures becomes larger as the SNR increases. In Figure 4(b), the spectral efficiency of various interleaved structures also increases as the SNR increases, and compared with Figure 4(a), the spectral efficiency is higher when due to the increased order of freedom. When the SNR is fixed, various interleaved structures are ranked as in descending order of spectral efficiency. The interleaved structure achieves the maximum spectral efficiency when , which is in accordance with the output of Algorithm 1 with a linear complexity of . Similarly as , the spectral efficiency is the smallest when , i.e., the localized structure. The same results can be found in Figure 5 where the number of RF chain in 8.

Figure 6 shows a comparison of the spectral efficiency for different subconnected structures in regard to the interantenna spacing in the case of and , respectively. It is shown that the spectral efficiency of the interleaved structure with fluctuates and attenuates as increases, and the performance of is in descending order. When is an odd multiple of the half wavelength, i.e., , , the maximum spectral efficiency of the interleaved structure is obtained when , and the minimum is got when . However, when the interantenna distance is an even multiple of the half wavelength, i.e., , we have in descending order of spectral efficiency. The same results as show in Figure 7 where the number of RF chain is 8. Figures 6(a) and 6(b) illustrate the connection between spectral efficiency and antenna array configuration, which provide help for the disposition of the mmWave massive MIMO antenna arrays.

Figure 8 shows the impact of the parameter of interleaved structure on the spectral efficiency for and , respectively. The interantenna spacing is . According to the above theoretical analysis, the optimal interleaved structure is related to the number of users. The parameter of the interleaved structure is a divisor of the number of users. For example, when the number of users is , there are five available interleaved structures, . It can be drawn from the figure that the structure with the largest spectral efficiency requires . Algorithm 1 gives the same results.

5. Conclusion

This paper has studied the connection between the antennas and the RF chains for the hybrid structure in mmWave massive MIMO systems. We investigated the best design of an interleaved subconnected hybrid structure in two situations using the features of the mmWave channel. One is the case with only one propagation path from each user to the BS, and the other is with multiple propagation paths. With the goal of maximizing the spectrum efficiency, the suboptimal interleaved structure has been derived, and numerical results were verified during the analysis. When the interantenna spacing equals even multiples of half wavelength, the best spectral efficiency is obtained by the localized structure. When the interantenna spacing equals an odd multiple of half wavelength, the interleaved structure with achieves the best performance.

Data Availability

The data used to support the funding of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Nature Science Foundation of China (Nos. 61701198 and 62171204), Nature Science Foundation of Jiangsu Province (No. BK20170557), and the Key Research & Development Plan of Jiangsu Province (No. BE2018108).