Research Article | Open Access
Lei Li, Hao Jin, Zhipeng Yan, Changqing Yang, Yong Wu, "Energy-Efficient User Association Strategy for Hyperdense Heterogeneous Networking in the Fifth Generation Systems", International Journal of Antennas and Propagation, vol. 2015, Article ID 686783, 9 pages, 2015. https://doi.org/10.1155/2015/686783
Energy-Efficient User Association Strategy for Hyperdense Heterogeneous Networking in the Fifth Generation Systems
Redesigning user association strategies to improve energy efficiency (EE) has been viewed as one of the promising shifting paradigms for the fifth generation (5G) cellular networks. In this paper, we investigate how to optimize users’ association to enhance EE for hyper dense heterogeneous networking in the 5G cellular networks, where the low-power node (LPN) much outnumbers the high-power node (HPN). To characterize that densely deployed LPNs would undertake a majority of high-rate services, while HPNs mainly support coverage; the EE metric is defined as average weighted EE of access nodes with the unit of bit per joule. Then, the EE optimization objective function is formulated and proved to be nonconvex. Two mathematical transformation techniques are presented to solve the nonconvex problem. In the first case, the original problem is reformulated as an equivalent problem involving the maximization of a biconcave function. In the second case, it is equivalent to a concave minimization problem. We focus on the solution of the biconcave framework, and, by exploiting the biconcave structure, a novel iterative algorithm based on dual theory is proposed, where a partially optimal solution can be achieved. Simulation results have verified the effectiveness of the proposed algorithm.
To provide universal high-data coverage and a seamless user experience, it is anticipated that the fifth generation (5G) cellular networks have extreme base station (BS) density and heterogeneity. Besides, the core networks can reach unprecedented levels of flexibility and intelligence , which allows the cellular strategies designed towards high efficiency with moderate complexity. Meanwhile the issue on improving energy efficiency (EE) has gained big momentum due to the increasing awareness of environmental protection and cost-efficiency. Energy-aware design and planning are motivated by the fact that wireless networks are responsible for a fraction between 0.2 and 0.4 percent of total carbon dioxide emissions , and this value is expected to grow due to the ever-increasing number of subscribers.
To curtail expenditures and improve EE performance, cloud radio access networks (C-RANs) are by now recognized as a promising system structure evolution for the 5G cellular networks . The densely deployed remote radio heads (RRHs) operate as soft relay by compressing and forwarding the received signals from UEs to the centralized baseband unit (BBU) pool. Then, the centralized large-scale cooperative processing, such as the joint decompression and decoding schemes, can proceed in the BBUs. Since RRHs are mainly deployed to provide high capacity in special zones, to guarantee backward compatibility with the existing cellular systems and support seamless coverage, high-power nodes (HPNs) are still critical in C-RANs. With the help of HPNs, the multiple heterogeneous radio networks can be converged, and all system control signals are delivered wherein. Consequently, HPNs should be incorporated into C-RANs, and thus heterogeneous cloud radio access networks (H-CRANs) are proposed in [4, 5], to take full advantage of both HetNets and C-RANs. The abovementioned advantages and challenges for cellular networks request a redesign of user association strategies . Therefore, this paper focuses on how to design user association strategies to improve EE for the hyperdense heterogeneous networking with a mass of low-power nodes (LPNs) and some HPNs.
Traditionally, users are associated with the BS providing the maximum received signal reference power (RSRP), which imposes a heavy burden on the tower-mounted macro BSs. To make the most of the dense low-power infrastructure, mobile users are actively pushed onto small BSs by using biases in . Considering the fairness among users, load balancing between macrocells and small cells is investigated in [8, 9]. In , the metric of cell selection is changed from the signal strength to the average throughput. Based on this metric, handover will happen once it brings positive gain of network throughput. In , the received signal to interference and noise ratios (SINRs) at users from BSs are multiplied by the designed factors to make small cells more attractive than macrocells. To control the bias towards improving throughput and enhancing users’ fairness during association, a -fairness network utility function is formulated and optimized in . The user association problem in conjunction with the almost blank subframe (ABSF) technique is considered in heterogeneous networks (HetNets) in , and the optimal ABSF density is proved to be the proportion of vulnerable users in total users. All aforementioned works have a network utility maximization objective that was adopted with log-utility to obtain network-wide proportional fairness , but few of them focus on the EE performance of cellular networks.
Recently, attentions have been also paid to the energy-efficient design of user association strategies [13, 14]. A distributed association strategy is developed to minimize the total power cost of heterogeneous cloud cellular networks in . An energy-efficient user association problem is studied from a population game-theoretic perspective to minimize power consumption in . However, the researches on energy-efficient user associations are still much limited and insufficient. One of the open issues is how to define EE for different networks. Ordinarily, EE metrics are mainly designed either to minimize the network power consumption under quality of service (QoS) constraints or to maximize the ratio of the network throughput to the area power consumption. However, for hyperdense heterogeneous networking in the 5G cellular networks, with the loads continuing to be transferred from macrocells to small cells, the roles of macro BSs and small BSs are gradually distinguished . High-power macrocells are mainly responsible for coverage, while a large number of low-power small cells undertake a majority of high-rate services. It indicates that power is consumed for both providing coverage and enhancing network capacity. Such distinction should be reflected in the EE metric. Therefore, the EE metric may need to be redesigned for above characteristics.
Another important issue is how to solve nonconvex EE optimization problems. In , the weighted EE is proved to be quasiconcave, and a bisection based resource allocation (RA) strategy is proposed accordingly. In , to monotonically increase EE, a quasi-distributed iterative RA algorithm is proposed for heterogeneous orthogonal frequency division multiplexing (OFDM) systems. In , an equivalent model for the nonconvex EE optimization object function is researched and solved by an iterative algorithm. These aforementioned works mainly focus on the design of energy-efficient RA strategies and suggest that there is no common solution for nonconvex optimization problems. In this paper, two mathematical transformation techniques will be presented to solve the nonconvex energy-efficient association strategy.
Contributions of this paper are summarized as follows. First, the EE performance of HPNs and LPNs is distinguished, and a new EE metric, that is, average weighted EE of access nodes, with the unit of bit per joule is proposed. Second, an EE optimization objective function is formulated to optimize user association strategy under the constraints of backhaul capacity, users’ data rate, and the maximum transmit power of BSs, which is proved to be neither convex nor quasiconvex. Third, two mathematical transformation techniques are introduced. In the first case, the nonconvex original problem is reformulated as a biconcave maximization framework. In the second case, the original problem is equivalent to a concave minimization framework. Furthermore, the relationship between two reformulated problems is presented. Fourth, an iterative algorithm based on dual theory and properties of the biconcavity is proposed to solve the equivalent biconcave maximization problem. Simulation results have verified the effectiveness of the iterative algorithm and suggest that transferring loads from HPNs to LPNs can improve the network EE.
The remainder of the paper is organized as follows. Section 2 gives the system model and formulates the problem. In Section 3, two equivalent transformation techniques are presented, respectively; then the optimization algorithm for the first equivalent problem is introduced. In Section 4, simulation results are presented, followed by some conclusions drawn in Section 5.
2. System Model and Problem Formulation
A general K-tier downlink radio network is considered, which can apply to a H-CRAN network or a dense HetNet. All deployed access nodes are called a BS for simplicity. Then, denote by the set of all BSs in K-tiers and the set of all users. As shown in Figure 1, users are allowed to be associated with multiple BSs at the same time. The association indicator for the th user and the th BS is represented as , where , , , and . Denote by the average spectral efficiency (SE) on the radio link between the th user and the th BS, and assume that the system has the knowledge of during association (similar assumptions can be found in [7–11], in which the received average SINRs at users from BSs for different radio resources are fixed during association; in addition, can be easily obtained in H-CRANs, wherein HPNs can achieve all system control signals in a centralized way). The data rate requirement of the th user is represented as , and the th BS is responsible for . Denote by and the total frequency band and maximum transmit power of the th BS, respectively. Then, the th BS needs to assign the user with frequency band, which accounts for power consumption. Based on above assumptions, it is easy to verify that if the requirement of the th user is satisfied, the achieved data rate is . The total transmit power and throughput in the th BS are written as and , respectively.
2.1. Energy Efficiency Metric
Then, how to define EE for the dense K-tier radio network is our focus in this part. As mentioned in the Introduction, for hyperdense heterogeneous networking, a large number of LPNs like pico BSs undertake a majority of high-rate services, while HPNs like macro BSs are mainly responsible for coverage. Though HPNs could contribute a small proportion of throughput compared with LPNs, they account for a significant proportion of power consumption due to high transmit, circuit, and cooling power. The existing EE metric defined as the ratio of area network throughput to area power consumption can not reflect the above characteristics. To distinguish the power consumption for ensuring coverage and providing high-rate services, we propose a new EE metric:which also has the unit bit per joule, where denotes the cardinality of the BS set and and are the predefined weight factor and the EE value of the th BS, respectively. A tradeoff between ensuring coverage and enhancing data services can be achieved by adjusting the weight factors of HPNs and LPNs.
2.2. Problem Formulation
Denote by the association matrix. Mathematically, we formulate the EE optimization objective function for the association aswhere and are the amplifier coefficient and the fixed power consumption in the th BS, respectively. The fixed power includes circuit power and cooling power (only macro BSs). The expressions (3), (4), and (5) constrain the transmit power of BSs, the backhaul throughput, and users’ data rate, respectively. In (4), represents the capacity limitation for the backhaul links of the th BS.
Remark 1. Problem (PX) is presented as the maximization of the sum of linear fractional functions over the domain constrained by (3)~(5). Though the domain can be proved to be convex, it is shown in  that problem (PX) is neither quasiconcave nor quasiconvex, even when and is set equal to 1, that is, even when problem (PX) is simplified as the sum of one linear term and one fractional term.
3. Two Mathematical Transformation Techniques and Association Optimization Methods
In this section, two mathematical transformation techniques for the nonconvex problem (PX) are introduced individually. First, problem (PX) is reformulated as an equivalent biconcave maximization problem. Second, problem (PX) is equivalent to a concave minimization problem. The relationship between the two equivalent problems will be introduced. Then, we focus on how to solve the biconcave maximization framework by employing dual theory and properties of biconcave functions.
3.1. The First Mathematical Transformation Technique
Since the numerator in (6) is differentiable and concave on and the denominator in (6) is positive and linear on , problem () is pseudoconcave on according to . Therefore, for each , any local maximum of problem () is also a global maximum of problem (). Furthermore, problem () can be solved by using any appropriate convex programming methods. For example, a simple roll polling algorithm can be employed here; that is, the th BS first serves the user with the largest , and other users can be served if and only if their association can enhance the EE of the th BS. Note that the association optimization for a single user to a single BS is convex. The detailed procedures for the roll polling algorithm are omitted here to save space, and we directly denote the maximum value of problem () by . Then, a vector parameter can be defined in the following set: where denotes -dimensional real space. Since problem (PX) is nonconvex, we introduce the following problem:
The objective of problem (PX1) is to achieve the optimal combination of that maximizes the network-wide utility function . Both and can be proved to be convex sets. Problems (PX) and (PX1) are equivalent in the aspects shown in Theorem 2 according to .
Theorem 2. Problem (PX1) is biconcave on . If is a globally optimal combination for problem (PX1), then is a globally optimal solution for the EE optimization problem (PX). Besides, it holds the relationship
Proof. First, and are convex sets. For all , letwhere and . Since, for each and ,to show that (PX1) is biconcave, it suffices to show that, for all , is biconcave. Given a fixed , is a quadratic function with respect to . Since the coefficient is always positive, is concave to on . In addition, given a fixed , the power function can be easily proved to be concave on . Besides, is linear. Therefore, is concave on the vector , where . The notation represents the transpose operation. The biconcavity has been proved.
Suppose that is the globally optimal solution to problem (PX1). can be directly obtained based on the biconcavity of problem (PX1) by solvingIt is easy to verify thatTo prove the optimality of , we first deny that is the globally optimal solution for the EE optimization problem (PX). Then, for some ,Also, because is biconcave on , the globally optimal solution to the subproblemcan be calculated asSince must be smaller than according to subproblem (), the combination is a feasible point in . Substituting (16) into problem (PX1), the utility value of problem (PX1) can be obtained, which equals . Similarly, substituting (13) into problem (PX1), the result is . Since is the global optimal solution of problem (PX1) in , it requires that which contrasts (14). Therefore, must be a globally optimal association matrix for the EE optimization problem (PX).
According to Theorem 2, the nonconvex EE optimization problem in (2) can be reformulated as an equivalent optimization problem (PX1) involving the maximization of a biconcave function over , where is a hyperrectangle.
3.2. The Second Mathematical Transformation Technique
Based on above discussions, the nonconvex problem (PX) can be reformulated as an equivalent biconcave problem. In this part, we will present another transformation technique which reformulates problem (PX) as a concave minimization framework.
For each , let and be defined as before. For a given , a new vector is introduced, which is constrained as
Then, a utility function is formulated on as
Notice that is a convex set. Furthermore, for each , is the minimum value of a set of linear functions of evaluated at . Since linear functions are also concave functions, is a concave function on according to . Based on the definition of function, we formulate a utility minimization problem as follows:
Problem (PX2) involves the minimization of a concave function on the convex set , and it is an equivalent problem of the EE optimization problem (PX) in the aspects shown in Theorem 3.
Theorem 3. If is a globally optimal solution for problem , then is a globally optimal solution for problem , where is any association matrix that achieves the minimum in (18) with .
Proof. First, suppose that is a globally optimal solution for problem (PX2) and achieves the minimum in (18) with . Since in cellular networks the expressions and are always positive for each , to achieve the minimum of , the variable should be assigned with a minimum value in . Therefore, based on the definition of , it holds that , and this implies that globally minimizes . It is equivalent to saying that is the globally optimal solution for problem (PX1), which maximizes . Therefore, is the globally optimal solution for problem (PX) according to Theorem 2.
Based on above discussion, the nonconvex EE optimization problem (PX) can be reformulated as a concave minimization problem (PX2). This new property provides a new way to solve the original EE optimization problem. In addition, as proved in Theorem 3, problem (PX2) can be further transformed to problem (PX1) by setting ; hence our next focus is how to solve problem (PX1) by taking advantage of its biconcavity.
3.3. The Energy-Efficient Association Strategy
In this section, the biconcave problem (PX1) is solved by using the Alternate Convex Search (ACS) algorithm as in . First, given a feasible association matrix , the subproblemis concave on and has a globally optimal solution represented as , where
The solution is a feasible point in for problem (PX1). Similarly, given the , the subproblemis concave on . To obtain the globally optimal solution of problem (S2), the dual theory is applied here. The duality gap is zeros since problem (S2) is concave. The Lagrange dual function for problem (S2) is formulated aswhere is the Lagrangian function and , , and are the dual variables for constraints (3), (4), and (5), respectively. The notation is an element-wide operation. can be easily proved to be concave on , and the optimal solution to (23) is represented aswhere .
Considering that is concave on , the optimal association matrix for problem (S2) can be achieved by using the iterative coordinate search method in . The explicit procedures are omitted here to save space. For simplicity of presentation, the solution of problem (S2) with input is represented as . The dual problem is given as
Obviously, is linear to . Hence, is the maximum of linear function and the dual problem (D-S2) is always convex. Furthermore, the subgradient of is presented as
Then, problem (S2) can be solved by updating via the ellipsoid method , of which the details are omitted here for brevity. Finally, the ACS algorithm is shown in Algorithm 1 at the next page. For simplicity, denote by and the solutions to problem (S1) and problem (S2) at the th iteration, respectively.
The complexity of Algorithm 1 is related with the value of . In each iteration, problem (S1) can be directly and optimally solved, while problem (S2) needs an inner loop to attain the optimal for a fixed . The main operation expenditure of Algorithm 1 is equivalent to that of solving problem (S2), which is further equivalent to the complexity of solving () linear equations. Therefore, the complexity of Algorithm 1 is .
Theorem 4. The obtained sequence generated by Algorithm 1 monotonically converges. Further, it holds that . ( denotes the set of nonnegative integers, and the notation denotes the Euclidean norm of a matrix.)
Proof. Since the sequence generated by Algorithm 1 is monotonically increasing and the function is bounded from above, the sequence converges to a finite value. Furthermore, suppose a feasible point . Since both problem (S1) with fixed and problem (S2) with fixed are convex and have a globally optimal solution, holds according to Theorem 4.9 in .
Although it has been proved in Theorem 3 that both the sequence and the sequence converge in finite iterations, the obtained by Algorithm 1 is a partially optimal solution to problem (PX) instead of a globally optimal one.
In this section, simulations are carried out to verify the effectiveness of Algorithm 1. A macro-pico HetNet is considered, and the parameters of BSs are set as in , where, for pico BSs, W, W, and ; for macro BSs, W, W, and . , Hz. In addition, . The initialed value of the association matrix is set as , that is, an average connection. We focus on the EE and the utility performances of a macro BS and several pico BSs.
Figure 2 shows how the network EE and the utility value change with respect to the iteration number in Algorithm 1. First, for both scenarios with ideal and nonideal backhaul, the EE lines and the utility lines monotonically increase with respect to the iteration number and match with each other in the end. The blue line represents the maximum EE value for the ideal backhaul scenario, where the optimal association matrix is obtained by a brute-force search. The solid lines of EE and the utility can gradually approach the optimal value as the iteration number increases, which verifies the effectiveness of Algorithm 1. In addition, as expected, nonideal backhaul can reduce the performance of network EE compared with the scenario with ideal backhaul.
Figure 3 shows how the EE and the utility value of individual pico BSs and macro BSs change with respect to the iteration number in Algorithm 1. First, as with Figure 2, all the dashed utility lines and solid EE lines of pico BSs and macro BSs match well when the iteration ends. Second, pico BSs perform much better than macro BSs in terms of EE, and their EE values are about 10 times larger than that of the macro BS. Third, it can be seen that, during the association, loads are transferred from macro BSs to pico BSs and among pico BSs. From the iteration to the iteration, the EE values of the above three pico BSs’ lines almost do not change, while the EE of the macro BS decreases, and the EE of the below pico BS’s line increases. It indicates that the loads of the macro BS are transferred to the pico BS. From the iteration to the iteration, the EE of the macro BS remains the same, while the blue pico BS’s line continues to increase; the other pico BSs’ lines monotonically decrease. It indicates that loads are being transferred among pico BSs. Since the network EE can monotonically increase, the above results suggest that one simple energy-efficient association strategy is to transfer the loads of macro BSs to pico BSs as much as possible.
Figure 4 presents the impact of the weight factors of macro BSs on the performance of the network EE. For both scenarios with ideal and nonideal backhaul, the network EE monotonically decreases as the weight factor of the macro BS increases. It further suggests that transferring loads from macro BSs to pico BSs is a promising approach to improve EE. For the blue line, the EE value at the point of is about 7 times larger than the EE value at the point of . Therefore, a proper design of the weight factors can achieve a large gain of network EE.
In this paper, we investigate the design of energy-efficient association strategy for hyperdense heterogeneous networking in the 5G cellular networks. The network EE is defined as average weighted EE of access nodes to characterize the feature that high-power macro BSs are mainly responsible for coverage while low-power small BSs undertake a majority of high-data rate services. The formulated EE optimization objective function is proved to be nonconvex. Two mathematical transformation techniques are presented to equivalently reformulate the original nonconvex problem as a biconcave maximization and a concave minimization problems, respectively. Furthermore, an iterative algorithm based on the properties of biconcave functions and dual theory is proposed to solve the biconcave problem. Simulation results indicate that transferring loads from macro BSs to small BSs can improve the network EE.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported in part by the National Natural Science Foundation of China (no. 61471062), the National Basic Research Program of China (973 Program) (Grant no. 2013CB336600), the Beijing Natural Science Foundation (Grant no. 4131003), the Key Program of Science and Technology Development Project of Beijing Municipal Education Commission of China (KZ201511232036), and the State Major Science and Technology Special Projects (Grant no. 2012ZX03001037-004).
- J. G. Andrews, S. Buzzi, W. Choi et al., “What will 5G be?” IEEE Journal on Selected Areas in Communications, vol. 32, no. 6, pp. 1065–1082, 2014.
- D. Feng, C. Jiang, G. Lim, L. J. Cimini Jr., G. Feng, and G. Y. Li, “A survey of energy-efficient wireless communications,” IEEE Communications Surveys & Tutorials, vol. 15, no. 1, pp. 167–178, 2013.
- M. Peng, C. Wang, V. Lau, and H. V. Poor, “Fronthaul-constrained cloud radio access networks: insights and challenges,” IEEE Wireless Communications, vol. 22, no. 2, pp. 152–160, 2015.
- M. Peng, Y. Li, Z. Zhao, and C. Wang, “System architecture and key technologies for 5G heterogeneous cloud radio access networks,” IEEE Network, vol. 29, no. 2, pp. 6–14, 2015.
- M. Peng, Y. Li, J. Jiang, and J. Li, “Heterogeneous cloud radio access networks: a new perspective for enhancing spectral and energy efficiencies,” IEEE Wireless Communications, vol. 21, no. 6, pp. 126–135, 2014.
- Z. Hasan, H. Boostanimehr, and V. K. Bhargava, “Green cellular networks: a survey, some research issues and challenges,” IEEE Communications Surveys & Tutorials, vol. 13, no. 4, pp. 524–540, 2011.
- S. Singh, H. S. Dhillon, and J. G. Andrews, “Offloading in heterogeneous networks: modeling, analysis, and design insights,” IEEE Transactions on Wireless Communications, vol. 12, no. 5, pp. 2484–2497, 2013.
- K. Son, S. Chong, and G. De Veciana, “Dynamic association for load balancing and interference avoidance in multi-cell networks,” IEEE Transactions on Wireless Communications, vol. 8, no. 7, pp. 3566–3576, 2009.
- Q. Ye, B. Rong, Y. Chen, M. Al-Shalash, C. Caramanis, and J. G. Andrews, “User association for load balancing in heterogeneous cellular networks,” IEEE Transactions on Wireless Communications, vol. 12, no. 6, pp. 2706–2716, 2013.
- C. Guo, M. Sheng, X. Wang, and Y. Zhang, “Joint scheduling and association for α-fairness network utility maximization in cellular networks,” in Proceedings of the 24th Annual International Symposium on Personal, Indoor, and Mobile Radio Communications (PIMRC '13), pp. 1769–1773, IEEE, London, UK, September 2013.
- Y. Jin and L. Qiu, “Joint user association and interference coordination in heterogeneous cellular networks,” IEEE Communications Letters, vol. 17, no. 12, pp. 2296–2299, 2013.
- T. Bu, L. Li, and R. Ramjee, “Generalized proportional fair scheduling in third generation wireless data networks,” in Proceedings of the 25th IEEE International Conference on Computer Communications (INFOCOM '06), pp. 1–12, April 2006.
- S. Moon, Y. Yi, and H. Kim, “Energy-efficient user association in cellular networks: a population game approach,” in Proceedings of the 11th International Symposium and Workshops on Modeling and Optimization in Mobile, Ad Hoc and Wireless Networks (WiOpt' 13), pp. 388–395, May 2013.
- H. Zhu, S. Wang, and D. Chen, “Energy-efficient user association for heterogenous cloud cellular networks,” in Proceedings of the IEEE Globecom Workshops (GC Wkshps '12), pp. 273–278, Anaheim, Calif, USA, December 2012.
- M. Peng, C. Wang, J. Li et al., “Recent advances in underlay heterogeneous networks: interference control, resource allocation, and self-organization,” IEEE Communications Surveys & Tutorials. In press.
- C. Xiong, G. Y. Li, S. Zhang, Y. Chen, and S. Xu, “Energy-efficient resource allocation in OFDMA networks,” IEEE Transactions on Communications, vol. 60, no. 12, pp. 3767–3778, 2012.
- K. Zhang, J. Jiang, M. Peng, and L. Li, “Energy-efficient power allocation algorithm for heterogeneous OFDM downlink systems,” in Proceedings of the International Conference on Wireless Communications and Signal Processing (WCSP '13), pp. 1–5, Hangzhou, China, October 2013.
- M. Peng, K. Zhang, J. Jiang, J. Wang, and W. Wang, “Energy-efficient resource assignment and power allocation in heterogeneous cloud radio access networks,” IEEE Transactions on Vehicular Technology, vol. PP, no. 99, p. 1, 2014.
- H. P. Benson, “On the global optimization of sums of linear fractional functions over a convex set,” Journal of Optimization Theory and Applications, vol. 121, no. 1, pp. 19–39, 2004.
- J. Gorski, F. Pfeuffer, and K. Klamroth, “Biconvex sets and optimization with biconvex functions—a survey and extensions,” Computers & Mathematics with Applications, vol. 15, no. 1, pp. 167–178, 2012.
- S. Boyd and L. Vandenberghe, Convex optimization, Cambridge University Press, Cambridge, UK, 2004.
Copyright © 2015 Lei Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.