Abstract
Public Safety Network (PSN) is the network for critical communication when disaster occurs. As a key technology in 5G, CloudRadio Access Network (CRAN) can play an important role in PSN instead of LTEbased RAN. This paper firstly introduces CRAN based PSN architecture and models the OFDM resource allocation problem in CRAN based PSN as an integer quadratic programming, which allows the tradeoff between expected bitrates and allocating fairness of PSN Service User (PSU). However, CRAN based PSN needs to improve the efficiency of allocating algorithm because of a mass of PSURRH associations when disaster occurs. To deal with it, the resources allocating problem with integer variables is relaxed into one with continuous variables in the first step and an algorithm based on Generalized Bender’s Decomposition (GBD) is proposed to solve it. Then we use Feasible Pump (FP) method to get a feasible integer solution on the original OFDM resources allocation problem. The final experiments show the total throughput achieved by CRAN based PSN is at most higher by 19.17% than the LTEbased one. And the average computational time of the proposed GBD and FP algorithm is at most lower than Barrier by 51.5% and GBD with no relaxation by 30.1%, respectively.
1. Introduction
With the rapid upgrading of emergency application, Public Safety Network (PSN) is getting increasingly popular. It is being recognized that wirelessbased PSN is a key to a successful response to emergency situations. Article [1] provides a good introduction to the field by surveying public safety use cases, current status of the related standards activities within the 3rd Generation Partnership Project (3GPP), and challenges ahead. Several works [2–4] have pointed out a safety network is a basic guarantee in different industries. Article [5] analyzes the infrastructure of LTEbased PSN and examines the capacity of LTE to meet the requirements of critical communications. The researchers and standards organizations both firmly believe PSN can be further enhanced by 5G emerging technologies, such as D2D, massiveMIMO, and CRAN [1, 6, 7].
There are many PSN issues that needs to be settled by new B4G or 5G technologies. Article [8] uses D2D communication to extend the coverage for PSN communications. Although D2D improves the efficiency of wireless communication, the key or the basic guarantee in PSN is the utilization on the bandwidth spectrum. Article [9] proposes a spectrum allocation scheme to avoid the interference with commercial network. And the American FirstNet is assigned Band 14 on 700 MHz to run the PSN service. Actually, as the increasing demand for richcontentbased emergency applications, PSN needs not only a clean sharing spectrum, but also a flexible centralized bandwidth allocation scheme that can maximally satisfy the datarate required by all PSN users. CloudRadio Access Network (CRAN), as an important access architecture in 5G, provides a centralized manner of bandwidth management [10]. In CRAN structure, the spectrum resource of a user associated with a certain RRH is processed by a centralized resource allocator. This integrated resource allocation mechanism enables the users to utilize the service resources more efficiently and makes the system allocate the user’s demanded bandwidth flexibly [11]. In addition, RRHs are always placed high above the ground, which means they are more difficult to be destroyed in disaster than the conventional BBU placed on the ground. In a consequence of that CRAN lets BBU together to form a BBU pool in a position which is far away from RRHs; CRAN is safer than LTE basestation infrastructure. It can be seen that these features of CRAN are more suitable for the requirement of PSN service than LTE EUTRAN.
However, the PSN Service Users (PSUs) will produce a large amount of traffic in times of disaster. Hence, the ability of bandwidth allocation algorithm is a key to assure the critical communication. Article [12] compares the performance of different scheme disciplines for PSN service. But the work is only designed to LTEbased PSN which is in a distributed management scheme [13]. Therefore it is necessary to propose a more higheffective centralized bandwidth allocation algorithm if we introduce a CRAN based PSN. The previous works concerning the resource allocation in CRAN networkmainly involve the flexible RRH configuration [14, 15] and the application of the spatial multiplexing technologies in CRAN [16, 17]. Although spectrum allocation in CRAN is also discussed in some of them, it is merely treated as a simple constraint. In addition, compared with the normal CRAN, CRAN based PSN maybe serve more users when disaster occurs. That means more efficient algorithm to solve the resource allocation problem is required. On the basis of the aforesaid works, in particular for their demerits, this paper proposes a Generalized Benders Decomposition (GBD) based OFDM resource allocation algorithm in CRAN based PSN to find out the solution on the system performance. This algorithm perfects the optimization through primal relaxation and the GBD algorithm that sufficiently enhances the efficiency of resource allocation.
The contributions of this paper include two key aspects. The first is we introduce CRAN based PSN architecture and propose the OFDM resources allocation model in this network. To the best of our knowledge, CRAN based PSN has not been investigated in the previous works. We formulate OFDM resources allocation problem in CRAN based PSN as a quadratic programming (QP) with integer variables of OFDM resource block that aims to maximize the emergency system performance. In the formulated problem, the objective system performance is affected by the expected bitrates revenue and the allocation fairness among different PSUs. By analyzing the wireless channel model, we propose an objective function of the CRAN based PSN system performance that allows tradeoff between the expected bitrates revenue and the PSU’s fairness when allocating resources.
Second, we present an efficient solution for the problem facing the centralized resource allocation in CRAN based PSN. We firstly relax the integer variable in this QP problem into the continuous one. And a GBDbased algorithm is proposed to address this relaxed problem. The traditional Bender’s Decomposition (BD) algorithm is used to solve the MixedInteger Nonlinear Programming (MINP) problem by decomposing the original problem into a master and a subproblem [18–21]. GBD extends the BD approach to a more general class of problems by adopting nonlinear duality theory, and some nonlinear problems are thereby brought into range, including QP we want to solve in this paper. Because of its high computational efficiency for solving the relaxed QP, the proposed algorithm is suitable in CRAN based PSN which has many PSU associations. Proved by theory, the proposed algorithm is guaranteed to obtain the optimal solution on the relaxed problem. Then we use the Feasible Pump (FP) [22] method to guarantee we can find out a feasible integer solution on the original QP from the optimal solution on the relaxed problem with continuous variables. We also analyze the error tolerance as the algorithm terminated condition to allow the balance tradeoff between the nice convergence properties and the system performance.
The rest of the paper is organized as follows: In Section 2, we describe the CRAN based PSN architecture and formulate the OFDM resource allocation problem. In Section 3, we solve the formulated QP of integer variables by problem relaxation, Generalized Bender’s Decomposition, and Feasible Pump. The specific algorithm is also proposed and analyzed in this section. In Section 4, we firstly investigate the performance of CRAN based PSN. And then the convergence and computing time of the proposed algorithm are analyzed. Finally, we drew some concluding remarks.
2. System Model
2.1. CRAN Based PSN Architecture
Figure 1 shows a schematic diagram of our introduced CRAN based PSN architecture for public safety applications. Each PSU runs different PSN service covered by a RRH. CRAN based PSN adopts a centralized processing manner, namely, the BBU processing pool. The spectrum resource of a PSU associated with a certain RRH is processed by a centralized resource allocator in CRAN. This integrated resource allocation mechanism enables the PSU to utilize the dedicated PSN service bandwidth resources more effectively and ondemand.
Our CRAN based PSN architecture allocates PSN spectrum resources centrally to improve the whole network performance. But to meet the critical communication requirement, this raises one important issue that needs to be considered: the assignment algorithm efficiency, especially for that there are a lot of PSURRH associations when disaster occurs. This is the key problem we want to solve in this paper. Since the OFDM is the dominant resourcemultiplexing technology currently, in this paper, we focus on the OFDM resource allocation in CRAN based PSN with a view to maximizing the system performance by optimally scheduling the centralized OFDM bandwidth resources.
2.2. PSURRH Association Model
Firstly, it needs to establish an association model between PSU and RRH. Assume that in a CRAN based PSN system, expresses the PSUs set and the RRHs set. We consider the network is an OFDM system whose resource element is called Resource Block (RB). Let denote the bandwidth of each RB and the total bandwidth. Hence the total number of RBs is . The RB scheduling matrix is denoted by , where is the number of RBs allocated to PSU by RRH . Particularly, means PSU is not served by RRH , which also stands for no association between PSU and RRH .
To facilitate analysis, we leave the reusable resource case to our further study, and temporarily in this section, the resource cannot be reused in the same CRAN based PSN system. Then we get one of the constraints of that . To make the model more convenient to solve, the RB allocation matrix should be transformed into a vector, the definition of which is given below.
Definition 1 (RB allocation vector). Let , where , if Then is called RB allocation vector for RB allocation matrix , where stands for the amount of RBs allocated to PSU from RRH . Particularly, means PSU is not served by RRH . A simple illustration of the transformation from a RB scheduling matrix to a corresponding RB allocation vector is shown in Figure 2.
In Figure 2, 2 RRHs and 3 PSUs make up the network, and the total number of RBs amounting to 100 is allocated by a CRAN allocator uniformly. It can be seen from the RB allocation matrix in the above example that PSU 1 and PSU 3 are served by RRH 1 with 23 RBs and RRH 2 with 24 RBs bandwidth, respectively. As multiple RRHs can perform joint transmission such as CoMP [23], PSU 2 is served by RRH1 and RRH2 simultaneously, and the spectrum resources allocated to PSU 2 totals at 53 RBs. The matrix is transformed into a column vector with the same 6 elements, following the rules illustrated in (1) that the columns in are rearranged into one column. The column vector is the corresponding RB allocation vector.
2.3. System BitRate Revenue
Bitrate capacity is a key feature to the performance of future wideband PSN richcontent services. The typical aim of the resource allocation scheme is to maximize the total bitrate capacity in the traditional cellular. Hence, in this paper we use the achievable bitrate capacity of all PSUs as the CRAN based PSN system’s revenue. This section analyzes the characteristics of the wireless fading channel to derive the outage probability of the system, based on which the achievable bitrate capacity of a PSU in unit bandwidth can be resolved. Finally, the achievable PSU bitrate per unit bandwidth is defined to represent the system revenue.
Let denote the transmit power of RRH , the antenna gain, and the signal power received by from RRH . The fading channel is characterized by a matrix , which is composed of 3 parts: path loss , slow fading , and fast fading . can be expressed by where , is the physical distance between PSU and RRH , and is the path loss exponent. Slow fading , where ; that is, follows lognormal distribution, and the Probability Distribution Function (PDF) is given in
Assume that the PSU uses Maximal Ratio Combining (MRC) scheme to achieve full diversity gain, and the fast fading can be calculated by where is the fast fading diversity gain from the th antenna of PSU and can be modeled as a random variable following standard normal distribution, that is, . is the number of antennas in the PSU . If , , and the PDF is given by
then the received power of UE from RRH , is given in
Supposing the threshold of the received power guaranteeing the basic communication service is , the outage probability of PSU to access RRH at a distance can be written as
In this paper, we take the impacts of two different scales of fading into consideration, which results in an outage probability resolving problem in elaborately modeled fading channel. Let , because and are independently and identically distributed; then
The cumulative distribution function of a gamma distributed variable with scale parameter being 1 can be written as , so the integral term of the slow fading can be written as
Substituting (9) into (7), the outage probability of PSU to access RRH at a distance can be rewritten as
Then we can use the GaussHermite quadrature to solve the numerical result of (10). Next, we are to derive the achievable bandwidth efficiency of PSU served by RRH . Firstly, according to (10), the PDF of SNR is
As is used to define CDF and its infinity integral is convergent, according to the derivation rule of the uncertain limit integral, we have
The PDF of SNR in (12) can be rewritten as
Let random variable denote the bitrate capacity per unit bandwidth of PSU served by RRH , in which is the average power of Gauss white noise. The expectation of , can be written as
Then (14) can be rewritten as
Then we can adopt the GaussLaguerre quadrature to solve the numerical result of (15). Based on (15), we give the definition of the expected revenue vector in CRAN based PSN.
Definition 2 (expected revenue vector). Let . If (16) holdswhere , then is the average bitrate capacity of PSU in unit bandwidth served by RRH , and is called the expected revenue vector of the system.
In addition, the variance of denoted as is
The correlation coefficients of and , , is
As complete joint transmission technology is adopted in CRAN and a PSU can be served by multiple RRHs simultaneously, in the case of and , and in other cases, .
2.4. Resource Allocation Problem Formulation
This paper takes into account both expected revenue of the system and PSU’s fairness when allocating the RBs in CRAN based PSN. Hence, for the sake of maximizing the system performance, the following quadratic programming (QP) is derived:where is a matrix relative to the covariance of the system revenue. The element , where and . It is noticed that the optimization problem aims at maximizing the system performance, where stands for the expectation of the bitrate capacity, that is, the average revenue of the radio access system with a RB allocation scheme . Based on the definition of , the second term stands for the variance of the system performance. If the parameter , the model is inclined to allocate more resources to PSUs contributing more revenue to the system, which means poor fairness of PSUs. If increases, the model will care more about the allocation fairness of PSUs.
Because of the integer decision variables of PSU, when having a large amount of PSURRH associations when disaster occurs, the time complexity of solving this NP problem will increase distinctly. Hence, this paper proposes a corresponding Generalized Bender’s Decomposition method on the basis of integer variable continuation, which satisfies the demand on efficiency of algorithm in PSN service.
3. GBD Based Resource Allocation Algorithm
3.1. Relaxation and Generalized Bender’s Decomposition
Firstly, relax problem (19) into
It can be seen that (20) is a typical QP problem with continuous variables after relaxing. Based on the definition of , it is easy to know is positive definite matrix. Therefore, the relaxed problem is a convex quadratic programming problem. However, this kind of problem is also hard to solve directly where there are a large number of quadratic variables. So in this paper, an algorithm based on Generalized Bender’s Decomposition is adopted to resolve the complex quadratic variables, resulting in a much simplified twostage solving process.
Let in be defined as . Then the other can be taken as the main variable. As is taken as an auxiliary variable, it can be regarded as a constant parameter in resolving main variable . And (20) can be rewritten as
Problem in (21) is called the slaveproblem of GBD, and the dual of (21) can be written as
From (22), it can be seen that due to the limit of the total bandwidth resources in the CRAN based PSN, the main variable conforms an equality constraint. Therefore the dual problem turned out to be a simple variable , though (22) the optimal solution of the dual problem can be computed by simply finding the maximum element in vector . Because there are only a limited number of optimal solutions on the dual problem, according to the strong duality theorem [24], there are also a limited number of optimal solutions on (20) and we have
Let denote the current number of Bender’s Decomposition iteration steps. Equation (20) can be written as a twostage optimization problem, composed of a Benders Decomposition slaveproblem as (24) and a relaxed Benders Decomposition masterproblem as (25):where ( is the optimal solution in ()th iteration).
Lemma 3 will prove the convergence of the above algorithm by demonstrating that the gap between the upper bound and lower bound of the optimal objective function value in (20) can be tightened, where is the value of objective function in (25); is the maximum value of problem in (24) in all the iterations, that is, the maximum element in .
Lemma 3. Each iteration in Generalized Bender’s Decomposition will narrow the gap between the upper bound and the lower bound of the value of objective function , where at th iteration, is and is .
Proof. Firstly, we prove the upper bound of being . Based on the strong duality theorem, we know problem in (21) is equivalent to the problem in (26) according to (22):So with the optimal solution of problem in (20) being , the value of the objective function is . As (25) is the relaxation of (26), there should bewhere is the optimal solution to question in (25). Therefore, the upper bound of the optimal solution is .
Then, we prove the lower bound of is . Because of the asymmetry in constraints of two dual problems, can be bounded or unbounded. If , , . In this case, the proposition that the bound of is holds. If , let . Suppose , and according to the above assumptions and strong duality properties, there should be an that conforms the following constraint:Equation (28) indicates that makes the objective function get a larger value to the objective function, which conflicts with the precondition that is the optimal solution of (20). Therefore, is the lower bound of the value to the problem in (19).
According to Lemma 3, the system can find the optimal solution of problem (20), as . Hence, the relaxation problem can be solved by Algorithm 1, and the specific solution is as follows: firstly, in the first iteration, that is, , the total system performance is greater than 0. According to the duality property, let and , where indicates a great enough value. These preconditions will ensure (25) can be feasibly constrained and the solution on the masterproblem in (25) can be resolved. Then let and solve the slaveproblem in (24). As slaveproblem in (24) is a Linear Programming (LP), the result should be an extreme point in the feasible domain of the problem. Afterwards, let iterations plus 1, and go on with the iterations to obtain an , until , which means the gap between the upper bound and lower bound is narrowed to 0. The extreme point corresponds to the optimal solution to the problem. The Generalized Bender’s Decomposition algorithm needs to obtain the channel information from all PSUs to construct a system revenue matrix, so it is a centralized algorithm and it is suitable to perform the system like CRAN based PSN, which pools all resources in a computing center.

The proposed GBD algorithm takes the difference between and as the criteria to judge the termination of the iteration. When the difference between and becomes 0, the whole algorithm terminates and the optimal solution to the problem in (19) is retrieved. The convergence of the algorithm will be proved by Theorem 4.
Theorem 4. The Generalized Bender’s Decompositionbased Algorithm 1 converges to an optimal solution after finite iterations.
Proof. As an optimization problem with continuous variables and linear constraints [24], there exists a polyhedral convex cone to make problems in (24) and (25) equivalent towhere the convex cone is composed of finite radials along with extreme directions. In each iteration, the set of in (25) will be tighten eventually, because a new radial will be added as a new extreme directions to refine the solution region; namely, . In addition, it is considered that is a convex hull composed of finite radials in extreme directions, and the newly introduced radial is different from the existing ones, so the , which is determined by a complete set of constraints, can be obtained after finite iterations. This suggests that the optimal resource allocation vector and the value of the objective function in (24), that is, the upper bound , can be obtained after finite iterations. So the optimal objective function value of relaxed problem (20) can also be obtained after finite iterations. Because the sequence is not increasing, can be obtained after is completely determined by the constraints , . Additionally, due to , the sequence is nondecreasing. After retrieving for the problem in (25), we have . According to the properties of the two sequences, terminated condition can be satisfied and the global optimal solution to the problem (20) can be obtained after finite iterations.
3.2. Feasible Pump (FP)
GBD algorithm can solve the solution of relaxation problem (20). Nevertheless, the vector solved by each iteration is not the integer vector, and it is not a feasible solution for problem (19). Therefore, it needs to find the actual feasible integer vector solution. There is a kind of simple method to find the integer vector solution; that is, the general vector solution is taken as the starting point to look for its nearest integer vector solution in the solution space. Article [22] proposes a kind of FP way which was taken as the standard to look for the integer vector solution. Specifically, the calculation is carried on each element in the optimal solution of th iteration in the relaxation problem (20), where signifies that the nearest integer is taken for the scalar. So the new integer vector is obtained. The function is defined as the sum of norm distance of various elements between a general vector and an integer vector; then
The principle of FP adopts a heuristic way to find the feasible integer solution . The determination standard of FP is that there exists a certain point in the general feasible region after relaxation and this point meets . In this case, is the integer vector and meets . Aiming at the relaxed problem (20), the FP standard is equivalent to the following minimum optimization problem:
It is clearly observed that if exists in the feasible region, then . This is equivalent to a fact that is a feasible solution. The entire execution flow of FP is as follows. Initially, optimal solution , which was obtained by the relaxed problem (20), carried on operator to obtain the initial integer solutions , and may not be the feasible integer vector solution. If is not the feasible solution, is fixed, and then the general vector solution closest to the norm distance of is looked for in the feasible region of relaxation. This process is called as a Pumping Cycling. Afterwards, FP standard (i.e., whether the norm distance is 0) was used to verify this solution. If the norm distance is equal to 0, the solution is . Otherwise, operator is carried on the newly obtained general vector solution to update and the Pumping Cycling iteration is continually carried on until is finally found.
3.3. OFDM RBs Allocation Algorithm Based on GBD and FP
Thus, this paper proposes the allocation algorithm on the basis of taking the OFDM RB as basic allocated unit, which is shown in Algorithm 2. The characteristic of this algorithm is that, firstly, the original allocation problem with integer RB is simplified by relaxing into one with continuous variable. And GBD method is used to solve this relaxed problem in each iteration. Then, the FPbased distance criterion is used to find the feasible integer solution. Afterwards, terminated condition is used to judge whether this solution is the acceptable approximate optimal solution.

3.4. TradeOff between Performance and Convergence Time
It can be seen that is used as the terminated condition of Algorithm 2. Therefore the value of will determine the accuracy of approximatelyoptimal solution. And it is not surprising that when is small the accuracy of the solution is higher, however with a greater computing time. Hence we use Theorem 5 to illustrate the proper value of as a reference to the implemented system.
Theorem 5. As for problem (19) and Algorithm 2, the system performance is at most less than the maximum system performance when (the operator signifies that the nearest integer is taken downward.), and as the difference between the system performance of Algorithm 2 and the maximum one is where is a constant that is greater than and can be expressed aswhere , and , are all integers
Proof. For problem (19), there exist and ; hence, the difference between the maximum system performance and the actual system performance of Algorithm 2 isBecause is an integer that is greater than 0, thenSince , , and are all the integers, there existsBecause the operator signifies that the nearest integer is taken downward, when , namely, , then . In this case, the result obtained from Algorithm 2 is the original approximatelyoptimal solution on the integer programming (19).
4. Simulation Results
This paper evaluates the performance from both CRAN based PSN architecture and proposed GBD algorithm. The throughput and PSU’s allocation fairness in CRAN based PSN are investigated in Section 4.1. The convergence properties and computational time of the proposed GBD algorithm are investigated in Sections 4.2 and 4.3, respectively. In Section 4.4, the tradeoff between the performance and the terminated parameter is discussed. We select CPLEX as the experimental solver and write the problem model and GBD and FP algorithm script into the CPLEX solver. The computational time evaluation is based on a Lenovo desktop with 3.7 GHz dualcore processors and 4 Gb of RAM. All corresponding results are averaged over 500 different simulation runs. The basic channel and radio parameters are shown in Table 1.
4.1. Performance Evaluation on CRAN Based PSN
Our experiment assumes a CRAN based PSN architecture of 10 RRHs of 20 W transmitting power with a fixed position placed in a 2000 m by 2000 m area. Each RRH serves 20 PSUs. These PSUs are randomly placed in each run and satisfies that 50% of them have average SNRs lower than 3 dB. In the simulation a variable bandwidth value per RRH is assumed (5, 10, 15, 20) MHz in both downlink and uplink. These bandwidths per RRH are equivalent to 25, 50, 75, and 100 OFDM RBs, respectively. The fairness factor of Algorithm 2 and is selected. On the base of the previous analysis, when is small, the allocator pursues to a high network performance revenue rather than PSU’s fairness of allocation. It is opposite when is large. This paper assumes all the RBs (Num. of . of RBs per RRH) can be allocated centrally without interference in the CRAN based architecture because all the channel information can be collected by the allocator; the interference pair of RBs can be arranged to the users who are far away enough to avoid the signal interference. We choose LTEbased PSN architecture as comparison with a typical allocation algorithm Proportional Fair (PF) [25] and an improved one called PS which is dedicatedly proposed for the PSN service [12]. The interference between PSUs in LTEbased PSN can also be ignored by the interference notice signaling message that communicated on the LTE X2 interfaces. We compare CRAN based PSN and LTEbased PSN from two perspectives: system total throughput and Jain’s Fairness Index (<1) [26]. Jain’s Fairness Index is used to measure the fairness among the PSUs. The ideal case of fairness is achieved when the index is equal to 1 reflecting all users receive equal resources while decreasing fairness is reflected with decreasing value of the index.
From Figures 3 and 4, we can see in CRAN based PSN architecture the system total throughput with the case of is better than while Jain’s Fairness Index with the case of is higher, which is consistent with the aim of the proposed algorithm discussed in Section 2.3. In the LTEbased PSN, PF can obtain a better throughput performance than PS while PS reveals better PSU’s allocation fairness. The results in LTEbased PSN are similar with the work [12]. This is because PS, which prioritizes the users with bad channel condition, is dedicatedly designed for the PSN situation with a considerable ratio of bad SNR users. PS scheduler provides a compromise in allocating the RBs among PSUs with good and bad channel conditions, which results in that the system throughput achieved using PS scheduler is lower than PF. From these two figures, it can be seen that the performance of CRAN based PSN is commonly better than LTE based. Particularly, if considering the PSU’s allocating fairness less, the system total throughput in CRAN based PSN of GBD scheduler with is higher than in LTEbased PSN of PF scheduler by 19.17% with 20 MHz 10 RRHs. Meanwhile, average Jain’s Fairness Index in CRAN based PSN of GBD scheduler with is also higher than in LTEbased PSN of PF scheduler by 0.045. If considering the PSU’s allocating fairness more, the system total throughput in CRAN based PSN of GBD scheduler with is at most higher than in LTEbased PSN of PS scheduler by 12.79% and average Jain’s Fairness Index is also higher by 0.03. We find out that CRAN based PSN has a better throughput performance than LTEbased without a reduction of Jain’s Fairness Index if having the same consideration on the PSU’s allocating fairness. The reason is that CRAN based architecture owns the capability of scheduling available RBs globally, which brings out a more obvious advantage of userondemand allocating. Meanwhile, since introducing the fairness factor as an important part of objective function model, CRAN based PSN with GBD algorithm can also obtain a good PSU’s allocating fairness. However, it can also be seen that Jain’s Fairness Index in CRAN based PSN of GBD scheduler with is lower than LTEbased with PS although the throughput is much better. This implies a suitable value can allow a wellbalanced tradeoff between the maximal throughput and the allocation fairness among each PSU.
4.2. Convergence of GBD and FPBased RB Allocation Algorithm
In Section 4.1, the performance of CRAN based PSN with GBD scheduler has been investigated. Then we move on to the algorithm’s convergence and efficiency. Figure 5 shows the convergence of proposed Generalized Bender’s Decomposition and FPbased algorithm compared with one only on the basis of Bender’s Decomposition with no relaxation. We assume a CRAN based PSN consisting of 10 RRHs and 60 PSUs. Firstly, this figure implies that difference between the upper bound and lower bound is tightened by Bender’s Decomposition so that an approximatelyoptimal solution of problem (19) can be obtained after finite iterations. This change is consistent with Lemma 3 and Theorem 4 in Section 3. And secondly, it can also be seen that the proposed algorithm needs 14 iterations to obtain the approximately optimal solution of the system, which is less than the other one of only Bender’s Decomposition. The acceleration is brought by the relaxation to the original problem (19). We use continuation to deal with the integer variable and then find the feasible integer solution by FP whose converged speed is faster than the method that directly search for the integer solution in the regional feasible set constrained by in each iteration.
Figure 6 illustrates the convergence of proposed Generalized Bender’s Decomposition algorithm with differentsized network consisting of 10 RRHs × 50 PSUs, 10 RRHs × 60 PSUs, and 10 RRHs × 70 PSUs. It can be seen that the algorithms in any assumptions converge to the approximatelyoptimal solution quickly. The figure also implies that the iterations in largesized network are less than a smallsized one. The reason is that when the RRH’s position and simulation area are assumed and PSUs are in large quantity, it is more likely to place more center users who are nearby the RRHs. Hence, the system performance in largesized network is generally higher than it in smallsized one. Judging from the same upper bound the algorithm started from and analysis in Section 3, it can be concluded that the iterations to a great value of approximatelyoptimal solution in a largesized network are less than to a small one in a smallsized network. One can conclude that the iteration of the proposed GBD and FP basedalgorithm is more relevant to the achievable optimal solution rather than the size of the network. This result on convergence features the applicability of the proposed algorithm in the centralized scheduling system like CRAN based PSN in this paper.
4.3. Computational Time of GBD and FPBased RB Allocation Algorithm
The average computational time of proposed GBD algorithm is expressed against the differentsized network in Figure 7. We select Simplex, Barrier and Bender’s Decomposition with no relaxation for comparison. Simplex is a general method to search for the (approximately) optimal solution of programming problem. Barrier is a commonly used way to solve QP problem. It is not surprising that the computational time needed in proposed GBD and FPbased algorithm is much less than that in Simplex and Barrier. And the time of Barrier and Simplex algorithms tends to increase in a nearly linear and exponential fashion, respectively, along with the network size while that of GBD algorithm has a more flatted one. Since the analysis on Figure 3, iterations of the proposed algorithm are less than that of one based on Bender’s Decomposition with no relaxation for an assumed network size. So the average computational of the proposed algorithm is least. It is reasonable that the average computational time of Simplex tends to increase in a nearly exponential fashion along with the network size. However the computational time of the proposed GBD algorithm increases in a more smooth fashion. The reason is that the computational time of Bender’s Decomposition depends on two factors: total iterations and the computing time in each iteration by the experimental desktop. The iterations decrease as the size of network increases as shown in Figure 6. But the experimental observation shows that the computing time in each iteration increases almost in a linear way along with the increased network size. Associated with the iteration steps, the average computational time of the proposed GBD and FPbased algorithm has a smoothly increasing fashion along the network size.
The similar results are obtained from the perspective of the jitter of the computational time which is shown in Figure 8. GBD converts the relaxed master QP problem into a dual problem as analyzed in Section 3.1 so that the jitter mainly depends on the experimental desktop performance rather than the solutionsearching method in each iteration on which the Simplex and Barrier depends.
4.4. TradeOff between Objective Function Value and Computational Time
Algorithm 2 is also run with different terminated parameter τ. We assume a CRAN based PSN network consisting of 10 RRHs and 60 PSUs. In Figure 9, the objective function value of the proposed algorithm is plotted against τ. It is not surprising to see that the objective function value degrades as the terminated parameter decreases. Moreover, the objective function value obtained is not sensitive to the value of τ when it is small enough, namely, smaller than 100 in our instance, which is consistent with Theorem 5 in Section 3. And the degradations of the objective function value with each τ are also not greater than the value proposed by Theorem 5. It is easy to know the greater τ is, the less computation is timeconsuming. Therefore we can select ondemand value of error tolerance τ that considers the tradeoff between the computational time and the objective function value in the CRAN based PSN system.
5. Conclusion
This paper researches on the CRAN based PSN architecture and OFDM resources allocation problem in it. Firstly, we design the CRAN based PSN and then by formulating the OFDM RB resources allocation problem into a QP with integer variable, we work out a solution considering both the expected system revenue and the PSU’s allocation fairness. To solve this complicated QP with integer variable, the continuation of integer variable is used to relax the original problem into a general QP problem. The relaxed problem is optimally solved by a proposed Generalized Bender’s Decomposition algorithm after finite iterations. Then FP standard is adopted to obtain the actual integer solution on the original problem. The tradeoff between the system degradation and the error tolerance of the proposed algorithm is given theoretically. The numerical results finally show the CRAN based PSN obtains a good throughput performance without a cost of PSU’s allocation fairness and the proposed GBD and FP based algorithm has the fine convergence and low computational time.
One achievement of this paper is to propose a PSN architecture based on CRAN and model the resource allocation problem in this centralized allocation system. In addition, OFDM resource allocation problem in CRAN based PSN is complicated because of the integer variables and a large amount of PSURRH associations when disaster occurs. However, the allocating efficiency is more important in PSN than commercial network. Hence, another achievement is that we use a corresponding GBD and FP method to solve resources allocating problem efficiently in the proposed CRAN based PSN, which is referential to the future works on the similar problem.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This paper is supported by the 863 Program (2015AA01A705) and NSFC (61271187).