Abstract

Low cost and high efficiency, defined as energy efficiency (EE) and spectral efficiency (SE), have raised more and more attention in the fifth generation (5G) communication systems due to steadily rising hardware cost, energy consumption, and mobile traffic. This paper studies the hybrid architecture of multiuser massive MIMO systems, where the digital domain utilizes the zero-forcing (ZF) precoding scheme and the analog domain uses discrete Fourier transform (DFT) processing that significantly reduces hardware cost and energy consumption. We derive analytical expressions on the total achievable SE and EE, as well as offering insight into some engineering parameters in the system performance. Our aim is to achieve low cost and high efficiency massive MIMO system, with constraints on the overall transmit power, the number of users, and the number of radio frequency (RF) chains. Results exhibit that the total achievable SE of the hybrid architectures with DFT precessing is inferior to the full digital architectures and hybrid architectures with the ideal phase shifters, but the performance attenuation can be compensated by providing the more input SNR and higher number of RF chains. Moreover, we find that the total achievable EE of hybrid architectures with DFT precessing outperforms other massive MIMO architectures that include a full digital implementation, ideal phase shifters, and a switched network.

1. Introduction

As environmental issue has attracted more and more concern from public, green radio communication has emerged as the most important concept in the fifth generation (5G) communication systems [14]. Along with improved traffic rate, transmission reliability, and quality of service, the upcoming 5G communication systems aim at improving the energy efficiency without compromising on the user experience [58]. Massive multiple-input multiple-output (MIMO), viewed as one of the important techniques for 5G wireless communication, has attracted a great deal of research interest in current years. For a conventional massive multiple-input multiple-output (MIMO) system with fully digital architecture, every antenna element requires to deploy a radio frequency (RF) chain that is equipped with a high-resolution analog-to-digital converter (ADC). The energy consumption of the ADCs grows exponentially with the number of ADC quantization bits and also increases with carrier frequency bandwidth and signal sampling rate. However, massive MIMO systems employing large-scale number antennas and high-resolution ADCs lead to unaffordable energy consumption and hardware cost [912].

Hybrid architectures massive MIMO systems, which consist of analog processing and digital precoding, are one such alternative with promising techniques in 5G communication systems [1315]. Such architectures not only achieve high total achievable SE, but also keep low cost and energy consumption. There has been a lot of work devoted to developing analog beamforming design. For example, the authors in [16, 17] proposed a novel algorithm for point-to-point MIMO systems, which aimed to reduce the computational complexity of system and maximize the total achievable spectral efficiency (SE) on different subcarriers. To handle the severe energy consumption and hardware cost, the authors in [18] proposed user scheduling algorithms and resource allocation schemes, which aimed to maximize the total achievable energy efficiency (EE) on hybrid massive system. However, all the aforementioned studies are based on the ideal phase shifters at the analog domain, which has a potential to enhance the complexity of the system, especially for a huge quantity of data. In order to simplify the complexity of system, the work of [19] proposed hybrid architecture with a switches network, where the ideal phase shifters are replaced by a switches network. Results reported that the switches network is a simpler solution but suffers from the significantly performance loss of the system.

Fortunately, the work of [20] found that using discrete Fourier transform (DFT) processing at analog domain has similar nature with electromagnetic lens and outperforms the fully connected phase shifters, which has low power consumption and is easy to implement. In light of these benefits, the DFT processing at the analog domain has gained great attention [21]. Results in [22] demonstrated that the DFT processing is carried out by the Butler circuit, which has low cost and high energy efficiency. Most importantly, the authors in [23] reported that the DFT processing is treated as the most promising technique used for signal transfer, which has been recommended by various industrial ports. Meanwhile, the authors in [24] investigated hybrid massive MIMO system based on DFT processing and further obtained the bounds on the achievable SE of such system, where the zero-forcing (ZF) precoding is adopted by the baseband processing. These studies confirm that the hybrid massive MIMO system based on DFT processing provides advantages for massive MIMO systems. However, the aforementioned work of [23, 24] did not derive an exact expression on the achievable SE and did not include analytical results on the achievable EE for hybrid massive MIMO systems.

In order to compensate for this gap in the study of achievable SE and EE, this paper presents a novel hybrid analog/digital architecture for massive MIMO systems, where analog processing is implemented using DFT processing. Moreover, we further investigate the EE and SE of the hybrid architecture massive MIMO system, which aims to effectively reduce the power consumption and hardware cost without compromising on the achievable SE. Our main contributions are as follows:(i)We present a novel hybrid architecture massive system based on the DFT processing and derive exact and approximated expressions on the total achievable SE of system. According to the derived analytical result, we further offer insight into the impact of several engineering parameters in the system performance, which constraints on the overall transmit power, the number of users, and the number of RF chains.(ii)By establishing the realistic power consumption model and employing the derived approximate expression, we evaluate the total EE for hybrid architectures with DFT processing. Closed-form expressions for the optimal the number of RF chains and the SNR by maximizing the achievable EE are derived.(iii)Results exhibit that the total achievable SE of the hybrid architectures with DFT precessing is inferior to the full digital architectures and hybrid architectures with the ideal phase shifters, but the performance loss can be compensated by offering more RF chains and more SNR. We also find that hybrid architectures with DFT precessing exhibit much better total achievable EE than the other hybrid architectures massive MIMO systems.

The reminder of the paper is listed as follows. Section 2 presents the signal model and hybrid architectures massive MIMO system with DFT precessing. In Section 3, we analyze the system’s performance that includes the achievable SE and EE. Based on the derived theoretical results, we attain the several engineering insights by considering the special cases. Section 4 provides the numerical results to confirm the obtained analytical results. Section 5 outlines the whole paper.

Notations. Matrix and vector are used by upper and lower case boldface, , , and stand for a matrix determinant, diagonal matrix, and factorial operation, respectively, denotes the -th diagonal elements of , accounts for rounding to the nearest integer, and and account for an identity matrix of dimension and an complex Gaussian matrix of dimension, respectively.

2. System Model

We firstly present the signal model associated with a hybrid analog/digital architecture in massive MIMO system, where the analog domain utilizes the low cost DFT processing, and then we describe channel model that consists of the small and large fading.

2.1. Signal Model

As depicted in Figure 1, we present a hybrid architecture of a multiuser massive MIMO system, where the BS comprises transmit antennas simultaneously serving users. Each user has a single-antenna element and the transmit antennas are connected to RF chains via a DFT processing with parameters restricted as and . Assume that denotes a transmitted signal vector of dimension for all users with . When the input signal is the Gaussian symbol [2527], the received signal after analog DFT processing technique and digital linear precoding scheme can be given byin which stands for the input SNR of the system, represents the digital baseband precoder matrix, denotes the analog processing matrix, denotes the additive white Gaussian noise (AWGN), and accounts for the channel matrix that includes the large-fading coefficient and fast-fading coefficient, which can be given bywhere accounts for the small-scale fading matrix, whose entries follow a complex Gaussian distribution, and denotes large-scale fading matrix that diagonal entries are given by and nondiagonal entries equal zero. In practical, large-scale fading coefficient contains the shadowing fading and path loss of the user , which is modeled as , where and account for a log-normal random variable of the user and the decay exponent of path loss, respectively, and accounts for the distance between the BS and the user .

2.2. DFT Processing Implemented

In the previous works about hybrid massive MIMO systems [13, 28], phase shifters or a network of switches requires the ideal CSI to be available at the BS. To overcome this drawback brought by phase shifters, this paper investigates a hybrid massive MIMO system as showed in Figure 1, in which the analog domain is implemented by using DFT processing. For hybrid architecture with DFT processing, all transmitted signals via a partial DFT matrix preprocessing before digital baseband precoding, in which the partial DFT matrix is derived by uniformly randomly selecting the rows of a DFT matrix. Therefore, the partial matrix is constructed bywhere .

Remark 1. From energy efficiency and economic perspective, the DFT processing is appealing for analog domain in massive MIMO systems since it has low energy consumption, less attenuation, and high performance, while the DFT processing is implemented by adopting a field programmable gate array (FPGA) [29] or an integrated Butler circuit [30], which is easy to implement and efficiently reduces hardware complexity. In addition, from the actual practical implementation, the DFT processing able to embrace all propagation signal from users and the DFT processing is implemented offline, which makes sure all RF chains always be activated.

3. Performance Analysis of System

In this section, our main task aims to derive the exact expression on the achievable SE. By using the obtained theoretical result, we offer the insight of engineering parameters on the achievable SE. And then the total achievable EE is investigated by considering the energy consumption of such a system. Finally a brief discussion is provided on how to achieve the EE maximization.

3.1. Hybrid Massive MIMO System

A hybrid architecture massive MIMO system with DFT processing is displayed in Figure 1, which is split into DFT processing and digital precoding processing successively. Before DFT processing at the analog domain and digital precoding at the baseband domain, an equivalent channel matrix is equivalent to a composite matrix composed of the channel matrix and partial DFT matrix, i.e., . Thus, the digital precoding matrices are based on . In addition, we make an assumption that the accurate and instantaneous CSI is available. In practical setups, the CSI is efficiently derived by channel estimation [31] and shared with the users using limited feedback techniques such as given in [32]. The ZF precoding scheme is able to completely eliminate interference among users, which achieves the performance of the nonlinear precoding schemes, for example, dirty paper coding (DPC) [33]. Thus, we would like to employ zero-forcing (ZF) precoding at the digital domain; can be written aswhere accounts for the normalization matrix that satisfies the constraint condition, i.e., , which makes sure the precoding power constraint is satisfied. When the input symbol is a Gaussian random variable, the achievable SE of the user with a ZF precoding scheme is expressed asConsidering all the users, the total achievable ergodic SE of system is identical to a sum of the achievable ergodic SE per user, which is calculated asAccording to (5), we see that the expectation is carried out via ergodic of the equivalent channel matrix. In order to derive the exact expression of the achievable SE, we have to calculate the probability density function (p.d.f) of the SINR term. In the following section, we shall present the detailed process about the achievable SE.

3.2. Spectral Efficiency Analysis

Theorem 2. Consider hybrid architecture massive MIMO systems with DFT processing; the exact expression on the achievable ergodic SE of the user with ZF precoding scheme is calculated aswhere denotes the exponential integral function of order . For real nonzero values of and the integer , the exponential integral is defined as

Proof. See Appendix A.

From Theorem 2, it is shown that the achievable SE of the user is related to the number of RF chains, the input SNR, and the number of users. However, it is hard to get an intuition about the impact of the system parameters on the achievable SE, since it invokes the exponential integral function. Therefore, we try to calculate an approximate expression on the achievable SE, which enables us to directly observe the general trends of the SE as a function of the most determinant parameters.

Corollary 3. Considering a hybrid architecture massive MIMO system based on DFT processing, an approximate expression on the achievable SE can be given by

Proof. To begin with, we invoke the results in [34] and directly obtain the following approximate expression:According to the above approximate equation, we have the following result:Substituting (11) into (8), we can derive the desired result.

With the derived tractable expression in (9) at hand, we obtain clear insights into the SE behavior with respect to different system parameters. From (9), we observe that is concerned with the large-fading coefficient, the input SNR, the number of RF chains, and the number of users. It can be seen that is a monotonically increasing function of the number of RF chains and the input SNR. Additionally, assuming a fixed number of BS antennas and the input SNR, is a monotonically decreasing function of the number of users. However, as the number of users grows, the total achievable SE increases since it is the multiplication of the achievable ergodic per rate by . Besides, the total achievable SE is a concave function with respect to the number of users , and there exists a unique globally optimal number of users maximizing the total achievable SE according to convex optimization theory. In the following corollary, we keep the other parameters being fixed and concentrate on finding the best number of users to achieve the total achievable SE maximization.

Corollary 4. For hybrid architecture massive MIMO systems based on DFT processing, there exists a unique globally optimal number of users that achieve the total achievable SE maximization, which iswhere and represents the Lambert function.

Proof. See Appendix B.

Corollary 4 reveals an important conclusion that the achievable SE will not linearly increase with the number of users, but it is a convex function. This implies that the optimal number of users in (12) achieves the total achievable SE maximization. Additionally, it worth noting that Corollary 4 typically shows us a noninteger value with respect to ; however, the quasi-concavity of the solution in (12) indicates that is obtained at one of the two closest integers.

3.3. Energy Efficiency Analysis

Green radio communication has emerged as a new prominent indicator in the design of 5G communication systems [35, 36]. Thus, the total achievable EE has become a primary concern, which can be viewed as the total achievable SE divided by the total energy consumption, which is established aswhere stands for the available bandwidth, whose unit is MHz, accounts for the overall energy consumption, and was defined as the total achievable SE of system, which is provided in (6).

We now present a practical energy consumption study of the system. According to the prior works of [3739], the main energy consumption comes in several varieties: circuit energy , signal processing energy , and signal transmission energy . Therefore, the total energy consumption in our system can then be calculated asNext, we shall introduce a detailed explanation of each energy consumption part in the following. (1)Circuit energy: the circuit loss of system mainly comes from three parts: which includes the LNA module, DFT processor, and one RF chain. Thus, the circuit energy is calculated aswhere accounts for the energy consumption from a , denotes the maximum energy consumption from DFT processor at the analog domain, and is the energy consumption for a single RF chain.(2)Signal processing energy: the total energy consumption mainly comes from the signal processing, which is given by [39]where , and stand for the energy consumption caused by the coding and decoding symbols, respectively. The second term can be understood as the energy computation of precoding due to the matrix inversion of the ZF scheme. The values of , , and stand for the coherence time, the computational efficiency, and the number of blocks per sec, respectively. The third term can be understood as the energy computation of information symbols during signal transmission.(3)Signal transmission energy: the energy consumption mainly comes from the loss of signal transmission in the circuits and transmission, which is function of the input SNR and is equivalent to , where accounts for the effective coefficient of the input SNR.

Under the above consideration about the energy consumption, we shall derive a closed-form expression of total achievable EE. By plugging (14) into (13), the total achievable EE of system is expressed aswhere . Based on the above optimization problem, it is shown that the total achievable EE in (17) is related to other fixed energy consumption, the number of BS antennas, the number of RF chains, the input SNR, and the number of users. Maintaining the number of BS antennas, the number of users, the input SNR, and the number of RF chains constant, we will focus on finding the best input SNR that achieves the total achievable EE maximization.

Corollary 5. Keeping the number of BS antennas, the number of RF chains, the number of users, and the other parameters constant, the optimal SNR is calculated as (18), which can be found below.

Proof. To complete the proof, we consider the optimization function as follows:with constants , and . With the aid of the conclusion provided in [40], we infer that the objective function is strictly quasi-concave and there exists an optimal that maximizes ; the optimal can be given bywhere is the natural number. By capitalizing on the above result and keep the system parameters constant except the SNR, we arrive at the desired result.

Corollary 5 provides an explicit guideline on how much SNR should be offered at the BS to maximize the EE. But, it is worth noting that the optimal SNR involves the Lambert function , which is a strictly increasing formula in accordance with the random variance [41, 42]. In the following, we shall present the numerical result to verify the derived analytical result.

Corollary 6. Keeping the number of BS antennas, the SNR, the number of users, and the other parameters as constant, the number of RF chains is calculated as where .

Proof. Using the similar method, the desired result is directly obtained by directly utilizing the results provided in (20) and setting the other parameters being constant except the number of RF chains.

4. Numerical Results

In this section, we provide the Monte-Carlo simulations to confirm the theoretical results on the achievable SE, which have been derived in Section 3, and then we compare the achievable EE for full digital architectures and hybrid architectures by considering realistic power consumption. In all simulations, the cell radius is 1000 meters, the guard zone ring region is set to 100 meters, the number of users is set to eight, the bandwidth of carrier frequency is = 10 MHz, the computational efficiency is = 12.8, the coherence time is set to ms, the decay exponent of large-scale fading is given by , and the standard deviation of large-scale fading is . The large-fading coefficients () are randomly created as follows: . Unless specified otherwise, the above simulation parameters are used. All the numerical results are averaged over 10,000 independent channel realizations.

Figure 2 depicts the total achievable SE versus the input SNR, where all curves include the exact expression of Theorem 2, the approximate expression of Corollary 3, and the simulated results. We can see that the curves for the simulated results, the approximate expression, and exact expression coincide perfectly in the entire SNR regime. This observation confirms the derived analytical expressions in Theorem 2 and Corollary 3. Meanwhile, Figure 2 also provides the total achievable SE for the full digital architectures [43], as well as the hybrid architectures with the ideal phase shifters and switch network [19]. It is shown that the total achievable SE of full digital architectures always outperforms the one of hybrid architectures, regardless of ideal phase shifters, switch network, and DFT precessing. However, the total achievable SE of the hybrid architectures and DFT precessing is superior to that of the hybrid architectures with the ideal phase shifters, but it is inferior to that of the hybrid architectures with the ideal phase shifters. As expected, there is a significant performance gap that is caused by inaccurate CSI and randomly DFT selection. According to (9), we know that the performance gap can be compensated for offering the more input SNR and the more number of RF chains.

In Figure 3, we depict the total achievable SE versus the number of RF chains. As expected, the total achievable SE increases with the number of RF chains, which keeps pace with theoretical results in Corollary 3. For comparison, the total achievable SE for the full digital architectures is also provided using [43, Table I], respectively. We can see that, with fixed number of users, number of BS antennas, and SNR (15 dB and 5 dB), the total achievable SE for full digital architectures maintains a constant value and is 41.54 and 20.07 bit/s/Hz, respectively, which keeps constant. This is because each antenna is connected to the RF chain. In addition, only when the number of available RF chains equals the number of BS antennas ( = = 100), the difference between hybrid architectures and full digital architectures tends to zero. Moreover, we observe that the total achievable SE with larger SNR regime always outperforms that one of the smaller SNR regime for the reason that the larger transmitted power boosts the total achievable SE.

Figure 4 depicts the total achievable SE varying with the number of users, in which the curves only showed the hybrid architecture with DFT processing. To guarantee fairness among users, we assume that the large-scale fading of each user is set to be identical, namely, . Moreover, given the number of BS antennas and the input SNR ( and SNR = 10 dB), we observe that the total achievable SE increases first and then decreases as the number of users increases from 2 to 40. Meanwhile, we depict the total achievable SE for different numbers of RF chains ( = 40, 50); no matter what the number of RF chains is utilized, a globally optimal number of users always exists that achieves the total achievable SE optimum, which keeps pace with theoretical result in Corollary 4. In addition, keeping the parameters , , and constant and plugging these parameters into (B.1) and (12), the best number of users and corresponding total achievable SE of system is given by (22, 5.521) and (26, 8.370), respectively, which further support our theoretical result in Corollary 4.

In Figures 5 and 6, we compare the total achievable EE for a hybrid massive MIMO system based on DFT processing, ideal phase shifters, and switch network, as well as massive MIMO systems with full digital architecture. Before depicting the total achievable EE, we shall present the energy consumption for full digital architecture, hybrid architectures with ideal phase shifters, and switch network hybrid architectures because the energy consumption of hybrid architectures with DFT processing has been provided in Section 3. With the help of the energy consumption models that have been provided in [19, 37, 44], we know that, for full digital architecture, hybrid architectures with ideal phase shifters, and switch network hybrid architectures, the energy consumption of the signal processing part and the signal transmission part is the same as the hybrid architectures with DFT processing, and only the energy consumption of the circuit part has changed drastically. For full digital architecture, the energy consumption of the circuit part is provided by [44]where = 1 Watt, = 2 Watt, and = 0.3 Watt. For a hybrid architecture massive MIMO system based on ideal phase shifters and switch network, the energy consumption of circuit part is provided by [19]where the subscript represents two different cases of the ideal phase shifter and switch network. The detailed energy consumption is provided by Watt, Watt, Watt, Watt, Watt, and Watt, and assume that the efficiency of transmitter power is = 0.5 at the BS. These above parameters are obtained from the results of [19, 37] for hybrid architectures.

Figure 5 displays the total achievable EE as a function of the input SNR. These curves are shown for hybrid architectures based on DFT processing, ideal phase shifters, and switch network. It is shown that the total achievable EE of hybrid architectures based on DFT processing is always superior to full digital architectures, as well as hybrid architectures ideal phase shifters and switch network. This is because the DFT processing has low cost and uses a few number of RF chains. However, the total achievable EE of the architectures with ideal phase shifters is better than that of the architectures with switch network for the reason that the degradation of the achievable SE with switch network as compared with the ideal phase shifters case is significant, especially for switch network case. Furthermore, the total achievable EE of hybrid architectures with DFT processing invariably outperforms the other in the input SNR regime, which indicates that the DFT processing is the best one, because the energy consumption of DFT processing is very small. In addition, it is shown that the total achievable EE first improves the input SNR and then reduces while the input SNR increases. Indeed, there exists the best input SNR that makes the total achievable EE of system maximization, which is consistent with analytical results in Corollary 5.

Figure 6 depicts the total achievable EE versus the number of BS antennas and includes full digital architecture and hybrid architectures with DFT processing, ideal phase shifters, and switch network. It is shown that the total achievable EE of hybrid architectures with DFT processing decreases by increasing the number of BS antennas. This is because keeping the number of RF chains and the SNR fixed, the total achievable SE becomes a constant, but the total energy consumption increases by increasing the number of BS antennas. Finally, this leads to the total achievable EE of hybrid architectures with DFT processing decreases. Obviously, the total achievable EE with DFT processing outperforms that of other hybrid architectures. For the hybrid architectures with ideal phase shifters and switch network, it is shown that the best number of BS antennas exists that makes the total achievable EE optimum, which is consistent with the conclusion described in [45]. For full digital architectures, it is shown that as the number of BS antennas increases, the total achievable EE decreases. This is because each antenna needs to deploy a unique RF chain, which is augmented by increasing the number of BS antennas.

5. Conclusion

In this paper, low cost and high efficiency hybrid massive MIMO system has been studied, in which the analog domain uses DFT processing and the digital domain utilizes ZF precoding. We derived the exact expressions on the achievable SE. Results showed that compared with full digital architecture and hybrid architecture with ideal phase shifters, there is some degradation of achievable SE for hybrid architecture with DFT processing. Fortunately, the degraded achievable SE can be compensated by providing the more number of RF chains and the input SNR. Moreover, we also studied the total achievable EE and derived exact expression of the achievable EE. Compared to conventional massive MIMO system, regardless of full digital and hybrid architecture, results showcased that the total achievable EE of hybrid architectures with DFT processing invariably outperforms the full digital architectures and hybrid architectures with switch network and ideal phase shifters. We also found out that the achievable EE decreases as the number of antennas chains increases, while the best SNR exists that realizes the total achievable EE optimum. We believe that the hybrid massive MIMO systems with DFT processing have low cost and high efficiency for 5G communication system. In the future, we will incorporate some other wireless techniques such as [4649] to further enhance the performance of the massive MIMO system.

Appendix

A. Proof of Theorem 2

To calculate the exact expression on the achievable SE in (5), we shall recall that the definition of with . By plugging these results into (8), the achievable ergodic SE in (8) is expressed asFor convenience, we firstly define an auxiliary matrix , which is identical to . According to the properties of the DFT matrix, we easily know that the column vectors of the DFT matrix are mutually orthogonal, and then the auxiliary matrix is transformed by the eigen-decomposition, which can be written aswhere is the diagonal matrix, whose entries are one or zero. Since the partial matrix is a dimension of DFT, by applying the eigen-decomposition is a diagonal matrix having unit singular values, whose diagonal entries can be expressed as . By applying the unitary transformation, the matrix has the same distribution form as , where reduces to a complex Gaussian matrix of dimension and is equivalent to columns of . Thus, (A.1) can be written asWe defineThus, (A.3) can be rewritten asWe study the evaluation of according to the following expression:According to random matrix theory, when the entries of small-scale fading are i.i.d. Rayleigh random variances [5052], the p.d.f. of is given bySubstituting (A.7) into (A.6) and employing the following integration identity:By doing more basic mathematical manipulations, we finish the detailed proof.

B. Proof of Corollary 4

For fairness among users, we assume that the large-scale fading of each user is identical, namely, . According to the achievable SE of per user (9), the total achievable SE of system can be expressed asWe differentiate the total achievable SE in (B.1) in terms of random variable , the first-order partial derivative of is calculated asSince the second term is negative, it is hard to judge directly if the value of (B.2) is negative or positive, and we need to further solve the second order partial derivative of . The detailed computation and associated analysis of second-order partial can be found in the works [5355]. The second-order partial derivative of is calculated asWe found that due to the value of numerator is a negative; this shows that the function is concave with respect to . According to convex optimization theory, it is easy to know that a unique globally optimal number of users always exist, which enables us to achieve the total achievable SE maximization. To check the optimal value, let the first-order partial derivative be equal to zero that holdsPlugging (B.2) into (B.4), we attainIn order to solve the above equation in terms of random variable , we observe that the equation is so sophisticated because it involves the logarithmic function. We start by letting ; (B.5) can be simplified asBy applying the properties of the Lambert function [41] that holds , we can attain the following equation:Substituting into (B.7) and playing some basic mathematical transformations, we finish the proof.

Data Availability

The authors derived the writing material from different journals as provided in the references. A MATLAB tool has been utilized to simulate the concept.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported in part by the Natural Science Foundation of Guangdong Province of China under Grant 2018A030310338, the National Natural Science Foundation of China under Grants 61671144, 61801132, and 61871139, the Project of Educational Commission of Guangdong Province of China under Grant 2017KQNCX155, the Guangdong Natural Science Funds for Distinguished Young Scholar under Grant 2014A030306027, the Innovation Team Project of Guangdong Province University under Grant 2016KCXTD017, the Science and Technology Program of Guangzhou under Grant 201807010103, the Project of Fundamental Science and Frontier Technology Research Project of Chongqing under Grant cstc2017jcyjAX0380, the Guangdong Science and Technology Project under Grant 2017A050506058, the Major Research Project of Guangdong Provincial Department of Education under Grant 2016KZDXM035, the Innovation Group Project of Guangzhou under Grant 1201610010, and the Yangcheng Scholars Project of Guangzhou under Grant 1201561560.