1. Introduction
Multi-hop relay systems have been studied as a key technology to increase the transmission rate and reliability in cellular communication systems, with a low infrastructure cost [
1,
2,
3]. The multi-hop relay system improves the channel quality between a base station (BS) and mobile stations (MSs) by deploying a relay station (RS) between the existing BS and MSs. RS may operate in either transparent or non-transparent mode [
4]. When an RS operates in transparent mode, the donor BS and the RS transmit the same signal to MSs simultaneously to increase the channel capacity, and the MSs do not recognize the existence of the RS. Meanwhile, when an RS operates in the non-transparent mode, it decodes the signal received from a donor BS and transmits the re-encoded signal to MSs, which eliminates shadow areas caused by signal attenuation. Moreover, since RSs can be installed at a lower infrastructure cost than BSs, the cell coverage of a conventional single-hop system can be extended more economically and efficiently.
In a multi-hop relay system, additional wireless resource should be allocated to the link between a donor BS and RSs. In setting aside the required resource for the communication between a donor BS and RSs, a time-division duplex (TDD) or a frequency-division duplex (FDD) scheme is applied to minimize the interference between the BS and RSs [
5]. Because whether an MS is registered to a BS or an RS determines wireless resource to be used for the MS, it is important for the MS to select an appropriate service node, i.e., BS or RS, to which it registers. If a certain BS or an RS is overloaded with many MSs, a resource shortage is inevitable, and the radio resource of other BSs or RSs would be under-used resulting in degradation of the system performance.
To prevent such performance degradation by overloading, methods of switching the service node of MSs were proposed, considering the traffic load condition [
6,
7]. However, these schemes do not guarantee the best channel quality for the MSs after switching the service nodes. Moreover, the scheme in [
6] cannot completely resolve the performance degradation because the service node of MSs changes after the overload has already occurred. Moreover, once installed, it is rather costly to move RSs. In [
8,
9], the authors proposed a method of increasing the amount of resource to be allocated to BS or RSs by adopting a frequency reuse scheme when resource shortage happens due to overloading. However, reusing the frequency resources may cause severe intracell interference. Because the above-mentioned schemes also commonly increase system complexity, it is important to minimize the possibility of overloading at the design level of a relay system.
The signal to interference plus noise ratio (SINR) distribution of MSs and the shadow areas in a cell are greatly affected by the RSs’ position. Accordingly, the installation location of each RS is a very important parameter determining the overall performance of a relay system. In [
10], a heuristic algorithm to determine the number and location of RS was studied by taking into account the probability distribution of MSs’ location; however, a detailed mechanism of obtaining achievable transmission rate is not specified. In [
11], an RS placement strategy with an advanced coding scheme was considered, and traffic demand was taken into account in determining the optimal RS location. In this scheme, a realistic user distribution is not considered. In [
12], the optimal BS and RS locations was selected from among a candidate set of sites, and an integer programming technique was adopted in solving this problem. In this work, the tradeoff between load-balancing and throughput enhancement was not considered. In [
13], the authors proposed a system capacity maximization scheme with multiple BSs and RSs for IEEE 802.16j networks, and experiments were conducted to confirm performance improvement through BS and RS placement in uniformly and non-uniformly distributed user scenarios. Several studies have assumed that the location of MSs in a cell follows a uniform distribution [
11,
14,
15,
16]; however, the distribution of MSs is greatly influenced by the existence of collective living infrastructure such as housing complexes, commercial districts, roads, large stadiums, etc. Therefore, the optimal RS location should be obtained considering unequal MS position distribution.
In this paper, we propose a method to find the sub-optimal location of RSs in a cell, based on the probability distribution of MSs’ location in a dual-hop relay system, where RSs operate in non-transparent mode. Specifically, the proposed scheme is constituted with the two stages of load-balancing and throughput enhancement. In the first stage, a rough position of RSs is determined by balancing overall traffic load in a cell site. In this stage, the computational complexity is very low, because the position of the RSs are adjusted by changing the distances and angles of RSs in a polar coordinate. In the second stage, the position of RSs are fined-tuned to enhance throughput, which results in the overall increase of spectral efficiency. This newly proposed two stage approach has an advantage in attaining the tradeoff between load-balancing and throughput enhancement. Moreover, in the load-balancing, unequal user distribution is considered for the purpose of reflecting a realistic user distribution scenario. The advantages of the proposed scheme are summarized as follows:
RS positioning with low computational complexity considering load-balancing over a cell site. Instead of intriguing global optimizations, which adjust all the combinations of RSs’ position in a cell site, the proposed scheme can reduce the complexity by adjusting the distance and angular position of RSs in sector-by-sector manner.
Reliable tradeoff between load-balancing and throughput enhancement through two stage approach. The fine-tuning of RSs’ position follows the first stage of load-balancing. Hence, the final result enhances the cell throughput with little sacrifice of load-balancing.
Algorithm framework considers unequal user distributions to reflect realist user distribution scenarios. Instead of simple equal user distribution, a user distribution with clusters is taken into account.
The simulation studies show that the proposed scheme has fairly good performance.
3. Proposed RS Positioning Scheme
In this section, we propose a method of finding the sub-optimal location for an RS. The proposed scheme achieves load-balancing to reduce the risk of overloading. In addition, a fine adjustment of the RS position is conducted to enhance throughput. This method prevents overloading and improves the overall throughput.
3.1. RS Positioning for Load Balancing
Please note that all the service nodes in a sector have the same amount of AZ resource, as shown in
Figure 2 and
Figure 3. If RSs are deployed in locations where service nodes can serve a similar number of MSs, the risk of overloading is reduced. Default RS locations are arranged considering the range of the sector. Because there are six RSs in a cell, the location coordinate of the
r-th RS in the
s-th sector can be denoted as
.
To adjust the location of each RS, the distance and the angle (degree) from the BS are required, and each position is expressed using the polar coordinate system as follows:
Here,
is a distance between the BS and each RS,
is an angle from the reference angle. The reference angle is
in Sector 1,
in Sector 2, and
in Sector 3.
Figure 4 shows the location of RSs using polar coordinates in Sector 1.
To find the RS installation location, we use the newly generated MS location MAP. From
of the generated MAP and the expected SINR at spot
n, we can estimate the number of MSs to be registered to each service node. According to the service node determined for the MSs at spot
n, the expected the number of MSs to be registered to each service node
,
, and
are determined. For instance, if the service node of MSs at spot
n is BS,
,
, and
are satisfied. Similarly, if the service node of MSs at spot
n is RS1 or RS2,
and
or
will be satisfied. Let
denote the new coordinate of an RS considering load-balancing among the service nodes. The process of obtaining
is conducted in a sector-by-sector manner by updating the default coordinate
L. First, the MS share for the BS (i.e., the ratio of the number of MSs registered to the BS to the total number of MSs in the sector) is adjusted, and this ratio can be expressed as follows:
In this step, and are adjusted so that should reside in the range .
If
,
and
are increased by
, the RSs are moved away from the BS until
is satisfied.
When the RSs move away from the BS, a large number of the spots between the BS and the RSs become areas served by the BS, resulting in the increment of . At this step, if the distance from the RSs to the BS becomes too long, interference from the RSs with the neighboring cells becomes too severe. Hence, the maximum RS distance is defined, and if either or is to exceed this upper bound, the current step is terminated.
In the case of
,
and
are decreased by
, moving the RSs toward the BS until
is satisfied.
As the RSs approach the BS, is decreased. Likewise, the minimum RS distance is set, and if either or is to exceed this lower bound, the current step is also terminated.
The next step is to adjust the MS shares of the RSs (i.e., the ratios of the number of the MSs registered to a specific RS, to the number of MSs registered to the RSs). The RS with the higher MS share is denoted as
, and the other is denoted as
. The MS shares for the two RSs are represented by
and
, respectively, and these can be calculated as follows:
In this step, , which is the angle of , is adjusted so that the lies in the range , where .
When
should be decreased by
to move
toward
until
is satisfied.
Then, the service nodes of many spots located between and are changed from to . At this time, if the two RSs are too close to each other, cell coverage cannot be efficiently extended. Therefore, the maximum angle needed to prevent this inefficient positioning is set. If is satisfied, the algorithm is terminated. However, if is not satisfied even after changing to the maximum angle, we change the angle of so that is satisfied.
To satisfy
(the angle of
), is increased by
, moving
toward the boundary of the sector.
The maximum RS angle is also defined in the same way to prevent severe intracell interference between different sectors. If is satisfied, the algorithm is terminated. The location of is adjusted before because moves toward the center of the sector, while moves toward the boundary of the sector. The latter has a higher risk of incurring severe intracell interference between two different sectors. If the is not satisfied even after moving and to the maximum change angle, the distance of is decreased by and the angles of the RSs are initialized. Then, this algorithm goes back to the first step. The whole process is detailed in Algorithm 1.
3.2. Fine Adjustment of RS Position for Throughput Enhancement
The RSs located at
can prevent overloading by adequately distributing the traffic load among the service nodes; however, additional fine adjustment of the RS locations is needed to increase the cell throughput. Cartesian coordinates are adopted in updating the RSs’ location to improve the cell throughput.
is the Cartesian coordinate representing
, obtained by Algorithm 1, and can be expressed as follows:
Algorithm 1 RS location for load balancing |
- 1:
Initialize: - 2:
- 3:
Input: - 4:
- 5:
step 1: adjust BS share by changing RSs’ distance - 6:
while ( and ) - 7:
- 8:
endwhile - 9:
while ( and ) - 10:
- 11:
endwhile - 12:
step 2: adjust RS share by changing RSs’ angle - 13:
while ( and ) - 14:
- 15:
endwhile - 16:
while ( and ) - 17:
- 18:
endwhile - 19:
if () - 20:
- 21:
- 22:
goto step 1: - 23:
endif - 24:
Output:
|
First, from the RSs’ position of Algorithm 1, the the total spectral efficiency of a sector is calculated. In this calculation, each sector should be divided by small spots, and the modulation and coding scheme (MCS)-levels and the average number of MSs for each spot should be taken into account. This is given by
where
n is a spot index,
is the selected service node at spot
n,
is the spectral efficiency at spot
n with the SINR
from service node
, where spectral efficiency the number of bits can be transmitted over 1 Hz [
24] and various spectral efficiency can be achieved by combining modulation schemes and channel coding rates. This can be obtained by referring to
Table 1.
From the current RS position, an RS could move in eight directions with angles equally spaced by
. The candidate RS positions (including the current position) are denoted as
, which can be expressed as:
where
is a unit distance of movement,
. For example,
is adhering to the current position,
is moving to upper right direction,
is moving to upper direction, etc. The total spectral efficiency with respect to each candidate position is denoted as
. If
is out of the bound of the sector,
is set to 0. Now, we select a direction
j that has the highest
as follows:
At this time, if
is not 0, the RS position is updated as
, and these processes are repeated until
becomes 0, because
means that the performance will be degraded when the RS is moved further. Therefore, when
, this algorithm terminates and
becomes the final RS installation location. The above procedure is detailed in Algorithm 2.
Algorithm 2 RS location for throughput enhancement |
- 1:
Initialize: - 2:
- 3:
Input: - 4:
- 5:
while () - 6:
- 7:
if () - 8:
- 9:
endif - 10:
endwhile - 11:
Output:
|
The proposed algorithm requires the MAP to use the information of the average number of users on each spot; however, it is not dependent on the MAP generation method. Accordingly, the proposed algorithm is separated from the MAP generation.
4. Performance Evaluation
To analyze the performance of the proposed schemes, we performed a simulation of a dual-hop relay system. There were six cells around a center cell, and the sub-optimal positions for RSs in the center cell were obtained considering the interference from neighboring cells. The radius
R of each cell was assumed to be 800 m. As shown in
Figure 2, the cell was divided into three sectors, and each sector had two RSs. The main parameters considered in the simulation are shown in
Table 2.
The proposed algorithms were run after determining appropriate input parameters. Preliminary tests were conducted to obtain these parameters, which are listed in
Table 3. In determining RS moving step size
, the expected total spectral efficiencies with different step size
is depicted in
Figure 5.
In this figure, the last point of each line indicates the number of iterations before the termination of Algorithm 2 and the achieved spectral efficiency. As shown in this figure, when the step size is small, it takes many iterations before Algorithm 2 terminates by finding the sub-optimal spectral efficiency. On the contrary, if the step size is big, the total spectral efficiency increases sharply at early iterations; however, it fails to achieve the highest sub-optimal spectral efficiency. One important thing should be noticed is that the smallest step size 10 m fails to achieve the highest sup-optimal spectral efficiency. Since we adopt the unequal user distribution and log-normal shadowing, there are many local maximums. If the step size is too small, Algorithm 2 may stop at one of the local maxima. Other test cases showed that the unit distance of movement and angle domain step size are rather insensitive to throughput performance, because Algorithm 1 focuses on load-balancing. Moreover, the ranges and also need to be carefully selected, because it affects the performance of Algorithm 1. In this paper, and were selected to achieve the lowest standard deviation of the number of MSs registered to the service nodes over the cell, considering the load-balancing capability of Algorithm 1.
Figure 6 shows
in a generated sample MAP. As mentioned above,
represents the expected number of MSs located at spot
n and
is equivalent to the expected number of MSs in the cell. The number of dense areas
C is 15, and the weight of dense area at the cell center is 50 and the remaining dense area is 20, so that an average of 330 MSs can be generated in the cell. The term
result is 296.71 on the generated MAP, which is smaller than 330 because the spots outside of the cell were excluded.
The generated MAP cannot be directly used in simulation, because it describes the average number of users on each spot in a real value. Hence, the integer number of users for spot
n over dense area
c is acquired through Bernoulli trials with probability
conducted
times.
Figure 7 shows a MS distribution sample, where user thinning is applied after being drawn from a generated MS location MAP. Before running the simulation, 10 different MS distribution samples were created with the sample MS location MAP. The default distance
for RSs was set to 1/2 of the cell radius, and the default RS location angle
was set to
. The simulation was performed for each distribution over 200 downlink frames.
In
Figure 8 each color point represents an MS and its SINR, and the black points represent RSs. As we can see in this figure, the SINRs of the MSs change as the locations of the RSs change.
Figure 8a is the SINR distribution with the default RSs locations,
Figure 8b is the SINR distribution with the RSs positions obtained by Algorithm 1. In
Figure 8b, compared to
Figure 8a, since the RSs are moved toward the edge of the cell, the SINRs of the MSs located at cell edges were greatly improved by Algorithm 1. Meanwhile, the SINRs of the MSs located around default RS locations were somewhat decreased. In
Figure 8c which represents the SINR distribution with the final RSs locations. As shown in
Figure 8c, the SINRs of MSs were further improved by moving the RSs toward to the centers of the clusters. Hence, the enhanced throughput is expected in
Figure 8c.
Figure 9 shows the cell throughputs for the ten samples of MS distribution and a mean throughput averaged over the ten samples. Please note that the purpose of Algorithm 1 is to balance the load among the service nodes and that the spectral efficiencies for the MSs are not considered in running Algorithm 1. Accordingly, the cell throughput may be reduced after Algorithm 1, because load-balancing may have some disadvantage in throughput maximization. After adjusting the RS positions using Algorithm 1, the average cell throughput decreased 8.1% compared to the default position. Using Algorithm 2, we found the sub-optimal RS location considering the expected number of MSs and their channel quality. By locating the RSs at the final positions obtained using Algorithm 2, the average cell throughput was increased by 16.4% compared to Algorithm 1, and increased by 14.9% compared to the default location. This confirms that Algorithm 2 greatly improves the average spectral efficiency, and that the purpose of Algorithm 2 was achieved. This improvement of spectral efficiency comes from the increment of SINR as clearly stated in Shannon capacity formula
[
26]. Accordingly, if the amount of bandwidth given in
B Hz, the achievable throughput will be
.
In addition, for the purpose of proving the optimality of the proposed algorithm, a semi-exhaustive search-based sub-optimal scheme was compared with the proposed algorithm. Finding the sub-optimal RSs position through an exhaustive full search is an extremely difficult task. For instance, if 800 m cell radius and 10 m positioning resolution are assumed, each sector can have approximately
different RS sites. If six RSs are in a cell, as in this paper,
test cases should be examined, and it is computationally intractable. In the semi-exhaustive search, 35 m RS positioning resolution was adopted, and the sub-optimal RS position was searched per sector. This semi-exhaustive search and the proposed scheme achieve the similar performance with each other as shown in the
Figure 9. However, note the computational complexity of semi-exhaustive search is still too high to be practically used.
Finally, the simulation result evaluating the load-balancing capability of the proposed scheme is presented. Because the cell was divided into three sectors, there were nine service nodes in the cell. We evaluated the state of load-balancing through the standard deviation of the number of MSs registered to each service node. If all the service nodes have the similar number of MSs registered to oneself, the standard deviation taken over these numbers of MSs will be small; however, if some service node have many MSs, while others have small number of MSs, the standard deviation will be relatively large. Accordingly, standard deviation can be adopted as a metric representing the level of load-balancing. On average, there were 296.71 MSs on the generated MAP, and the actual number of MSs was reduced to 267.3 after user thinning. Accordingly, the average number of MSs per service node was 29.7. The standard deviation of each sample was calculated over nine service nodes, and the results are depicted in
Figure 10.
As we can see in this figure, the overall standard deviation with the default RSs position is 13.7, which is a relatively high value. It means that some service nodes have too many MSs to serve while others have a small number of MSs to serve, which may incur a risk of performance degradation due to overload. However, Algorithm 1 decreased the standard deviation to 8.71, which is about 36.4% lower than for the default location, which makes the load more balanced. After running Algorithm 2, the standard deviation became 12.86, which is higher that Algorithm 1, but still lower than default position. Even though Algorithm 2 has increased the standard deviation by 4.15 compared to the result of Algorithm 1, it should be noted that the loss in the load-balancing is low while the throughput has increased greatly. As a result, we can improve both the load-balancing and the throughput by the proposed scheme, where an average of 6% standard deviation is decreased, and the throughput is increased by 14.9%.