Channel assignment for cellular mobile networks with ... - IEEE Xplore

1 downloads 0 Views 274KB Size Report
algorithm. R. Ch!avez-Santiago, E. Gigi and V. Lyandres. Abstract: A heuristic sequential algorithm for fixed channel assignment in cellular mobile networks.
Channel assignment for cellular mobile networks with nonuniform cells – an improved heuristic algorithm ! R. Chavez-Santiago, E. Gigi and V. Lyandres Abstract: A heuristic sequential algorithm for fixed channel assignment in cellular mobile networks with nonuniform cells, i.e., cells of different sizes and shapes, is presented. The channel assignment is treated as a nonbinary constraints satisfaction problem. The algorithm uses both the co-channel and adjacent channel interference levels and the number of channel requirements per cell to define the corresponding sequence of assignment. The effects of intercell overlapping and complex propagation mechanisms on the channel assignment results are investigated too.

1

Introduction

The electromagnetic spectrum is one of the most important natural resources and must be administrated as efficiently as possible. Frequency reuse within a wireless communications network can offer considerable reduction of the required spectrum; however, it may also lead to a loss of quality of communication links owing to interference. The solution for such a problem balances the economies of frequency reuse and the loss of quality in the network. The channel assignment problem (CAP) has commonly been represented as a constraint satisfaction problem (CSP) [1], assuming a regular hexagonal network layout and binary constraints between pairs of transmitters (TXs). These constraints have the following form   fi  fj 4g 8g  0 ð1Þ where fi and fj are the frequencies assigned to transmitters i and j, respectively, and g is a positive integer. Based on this model, the co-channel, adjacent channel, and co-site constraints are defined [2, 3]. However, it has been demonstrated that modelling the CAP with only binary constraints might lose important elements of the original problem [1, 4], which produces solutions that are not necessarily the best ones. The use of nonbinary constraints is an alternative representation of the CAP that consists of finding an assignment of frequencies for all TXs in a mobile network in such a way that an acceptable carrier-tointerference ratio (CIR) over all points of the coverage area is obtained. This approach was used to solve the CAP with the help of simulated annealing (SA) in [3, 5]. Although SA provides very good solutions, it requires very long processing time to reach optimality, which makes SA unpractical for quasidynamic channel assignment schemes, in which the assignment is modified according to the traffic pattern at different hours of the day. r IEE, 2006 IEE Proceedings online no. 20045331 doi:10.1049/ip-com:20045331 Paper first received 12th December 2004 and in final revised form 28th June 2005 The authors are with Department of Electrical and Computer Engineering, Ben-Gurion University of the Negev, P.O.B. 653, Beer-Sheva 84105, Israel E-mail: [email protected] IEE Proc.-Commun., Vol. 153, No. 1, February 2006

A heuristic sequential algorithm with nonbinary constraints (CIR constraints) that takes into account the co-channel interference level to compute a difficulty degree (DD) that defines the sequence of assignment was presented in [6]; this CIR constraint algorithm is based on the representation of the CAP as a network colouring problem [7], and it solves the CAP fast enough for quasydynamic channel assignment applications. The algorithm solves the CAP in cellular mobile networks in which the coverage area is divided into a number of circular cells with different sizes. Conventional channel assignment (CA) techniques [2, 3, 8] cannot be directly applied to this kind of networks that do not have a regular frequency reuse pattern. Unfortunately, the algorithm considered in [6] is based on several assumptions that simplify the problem but, in most of the cases, do not capture properly the elements involved with real networks, reducing significantly its practical application. In this paper we present an improved algorithm that takes into account every base station’s (BS) transmit power and the adjacent channel interference to compute the CIR constraint. We also incorporate a more accurate urban propagation model. 2

CIR constraint algorithm

Let Iij be the co-channel interference from cell j to cell i in a cellular mobile network, and bi the set of all the cells (excluding cell i) that use the same channel as cell i. The CIR at the input of a mobile station (MS) located at the contour (assuming that the cells have circular shape) of cell i can be expressed in its more general form as Ci CIRi ¼ P Iij

ð2Þ

j2bi

where Ci is the signal power from the BS in cell i as received at the contour of the cell. Using the simplest model that expresses the propagation loss L proportionally to the path length d as L  d 4 in an urban environment, and assuming that all the BS antennas in the network radiate the same power, we can 61

3

rewrite (2) as follows CIRi ¼

R4 P i 4 dij j2bi

a

ð3Þ

where Ri is the radius of cell i, a is a specified CIR threshold, and dij is the worst case distance between the interfering cell j and cell i, given by qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi  2  2 xi  xj þ yi  yj  Ri ð4Þ dij ¼   where ðxi ; yi Þ and xj ; yj are the Cartesian coordinates of the BSs of cells i and j, respectively. Defining g ¼ 1=a, (3) can be written as follows X dij 4 g ð5Þ Ri j2b i

The last expression defines the co-channel interference constraint (hereinafter referred to as CIR constraint) that is used in the assignment algorithm as it will be explained in the next Section.

2.1

Algorithm description

The heuristic algorithm presented in [6] establishes a sequence of assignment in descending order according to the DD of every BS. Ten different DDs were defined in terms of interference from neighbouring cells and the number of channel requirements per cell, denoted by mi. The first and simplest DD is given by  4 N X dij DDi ð1Þ ¼ mi ð6Þ Ri j¼1 j6¼i

whereas the last and most complex one is given by "   4 # N X dij 4 dji DDi ð10Þ ¼ xij mi þ Ri Rj j¼1

ð7Þ

j6¼i

where N is the total number of cells in the network and xij is an influence factor whose value is 1 if the cells i and j use the same channel simultaneously, and 0 otherwise. The assignment process basically consists of a requirement–exhaustive strategy [6, 8, 9]. First, the algorithm arranges the cells in descending order in terms of their DD. Then, the first channel is taken and assigned to the first cell in the sequence, i.e., the cell that has the largest DD; after this, all the rest cells in the sequence are inspected in descending order according to their corresponding DD to determine whether the channel can be reused or not. Only if the CIR constraint given by (5) is satisfied, the channel reuse is allowed. After all the cells have been inspected, the sequence of assignment is reordered and then the second channel is considered. This process is repeated iteratively until all channel requirements have been satisfied. This algorithm belongs to the so-called exhaustive search techniques, which can prove optimality and are quicker than other heuristics and metaheuristics (e.g., genetic algorithms, simulated annealing, tabu search, etc.) in providing a solution. These facts were corroborated by a channel assignment example in a network with uniform cells [6]; the algorithm reached the theoretical lower bound of 275 channels predicted by [10]. Detailed analyses of the computational complexity of the algorithm can be found in [9, 11], which show that the algorithm is faster than any metaheuristic for problems of the same size. These characteristics justify further development and improvement of this channel assignment technique. 62

Improved CIR constraint algorithm

Although the algorithm described in the previous Section has provided satisfactory results in some practical examples [6], there are a number of inconveniences that limit its application to realistic scenarios:  The algorithm was designed for channel assignment in networks where cells have different sizes; therefore, it cannot be assumed that all the BS antennas radiate the same power because the size of a cell is determined, among other things, by its antenna’s transmit power.  The algorithm considers only co-channel interference. Although neighbouring cells using the same channel are the main source of interference, adjacent channels also contribute to the deterioration of the CIR and affect the assignment results. Considering only co-channel interference provides assignment plans that economise spectrum, but that in most of the cases maintain a considerable level of harmful interference.  The effects of complex radio propagation mechanisms also affect the final assignment results. Thus, a more sophisticated propagation model than L  d 4 must be considered in some special cases.  It is unrealistic to assume that all the cells have a regular shape, e.g., perfect hexagons or circles. An algorithm for more realistic scenarios must take into account the fact that most cells might be irregular-polygon-shaped.

3.1

Extended CIR constraint

The points listed above motivated the modification of the algorithm considered in [6] in order to make it practical for assignments in more realistic scenarios. The first improvement we introduced to the algorithm is the extension of the CIR constraint to include every BS’s transmit power and the adjacent channel interference. Adjacent channel interference is reduced by filter selectivity and its harmful effects decrease as the spectrum separation between channels increases. A simple model to compute the filter attenuation on an adjacent channel is given by [1] adj factorg ¼ ð1 þ log2 gÞ 8g40

ð8Þ

where g is the spectrum separation (in number of channels) between the adjacent channel frequency and the filter’s central frequency and L is an attenuation constant. A typical value for L is 18 dB; therefore, for g ¼ 1, an adjacent channel is attenuated by a factor equal to 0.015. Let bi;g define the set of all the cells (excluding cell i) that use an adjacent channel with a spectrum separation g from a channel used in cell i; thus (3) can be extended, including transmit power and adjacent channel interference, to the following form CIRi ¼

Pi R4 i P j2bi

Pj dij4

þ

n P g¼1

adj factorg

P k2bi;g

!a Pk dik4 ð9Þ

where Pi is the transmit power of the BS antenna in cell i and n is the maximum number of spectrum separation channels to be considered for the computation of adjacent channel interference. In a straightforward way, (5) can be IEE Proc.-Commun., Vol. 153, No. 1, February 2006

rewritten to define an extended CIR constraint as 0 1 n X Pk d 4 X Pj dij4 X ik Ag þ adj factorg @ P R4 g ¼ 1 Pi R4 i j2b i i k2b i

3.2

ð10Þ

i;g

Propagation model

Propagation loss in urban environments is usually predicted by empirical and semideterministic models. The inclusion of one of these models in (9) and (10) is relatively easy; the propagation loss L computed with an adequate model should replace R4 and d 4 in all the cases. Most propagation models consist of a formula that defines L as a function of the propagation path d and the operation frequency f, i.e., Lðd; f Þ. Thus, the CIR constraint can be written in a general form as follows 0 1   n X X Pk Lðdik ; f Þ X Pj L dij ; f Ag þ adj factorg @ P i Lð R i ; f Þ g ¼ 1 Pi LðRi ; f Þ j2b k2b i

i;g

ð11Þ The Walfisch-Ikegami model (WIM) [12] is a semideterministic propagation model for a network with medium/large size cells in built-up areas, adopted within the framework of the Cooperation in the Field of Scientific and Technical Research (COST 231). It demonstrates a good fit to measured propagation data for frequencies in the range of 800–2000 MHz and propagation paths in the range of 0.02–5 km. The WIM provides formulas for line-of-sight (LOS) and non-line-of-sight (NLOS) propagation situations. In a LOS situation the WIM gives the propagation loss in dB by ð12Þ LLOS ¼ 42:64 þ 26 log10 d½km þ 20 log10 f½MHz For NLOS path situations, the WIM uses the parameters illustrated in Fig. 1. Such parameters are: hb BS antenna height over street level in metres (4–50 m) hm MS antenna height in metres (1–3 m) hB nominal height of building roofs in metres Dhb ¼ hb  hB height of BS antenna above rooftops in metres Dhm ¼ hB  hm height of MS antenna below rooftops in metres b building separation in metres (use 20–50 m if no data provided)

w width of street (use b/2 if no data provided) f angle of incident wave with respect to the street (use 901 if no data provided). The propagation loss in dB for NLOS conditions is given by  Lfs þ Lrts þ Lmds ; Lrts þ Lmds  0 ð13Þ LNLOS ¼ Lfs ; Lrts þ Lmds o0 where Lrts is the roof-to-street diffraction and scatter loss, Lmsd is the multi-screen diffraction loss, and Lfs is the free space loss given by Lfs ¼ 32:45 þ 20 log10 d½km þ 20 log10 f½MHz

Lrts and Lmsd are functions of the NLOS parameters listed above; a comprehensive description of their computation can be found in [12].

3.3

Irregular-polygon-shaped cells

For the case of networks with irregular-polygon-shaped cells, we define a matrix of ‘radiuses’ Ri and a matrix of distances dij for the worst-case interference point between each and every pair of cells. The matrices of ‘radiuses’ and distances have the following form, respectively 3 2 2 3 

6 R21 R ¼ 6 4 .. . Ri1

R12  .. . R12

... ... .. . ...

Rij R2j 7 .. 7 5; . 

 6 d21 6 D ¼ 6 . 4 .. di1

. . . dij . . . d2j 7 7 . 7 .. . .. 5 ... 

These matrices are used for the CIR computation between pairs of TXs. Let us consider an example to illustrate this point. Figure 2 depicts a network with three irregularpolygon-shaped cells, where the worst-case interference points between every pair of cell are indicated by black dots and the ‘radius’ and distance for every case are labelled. Thus, the co-channel CIRin cell  1 causedby 4 4 ¼ P R P2 d12 þ cell 2 and 3 is given by CIR 1 1 12   4  4 P3 d13 , using only the co-channel interference P1 R13 part of (10); in a straightforward way, the co-channel 3 caused by cell 2 is equal to CIR3 ¼ CIR 4incell  4 P2 d32 . P3 R32 cell 1 ( P1)

cell 2 ( P2)

d

d12

R13 cell 3 ( P3)

buildings

d13 P2

hb

d12  .. . d12

ð15Þ

R12

BS antenna

ð14Þ

d32

R32

P3

hB

w b

BS antenna

hm

worst-case interference point MS antenna

Fig. 2 ∆hb = hb − hB MS

incident wave

Fig. 1

Network with irregular-polygon-shaped cells

∆hm = hB − hm

NLOS parameters for the WIM

IEE Proc.-Commun., Vol. 153, No. 1, February 2006

4 direction of travel

Examples and results

We implemented the heuristic algorithm including all our proposed modifications, and investigated the impact that such modifications have on the channel assignment result in networks with both regular and irregular layouts. All the 63

4.1

Effects of adjacent channel interference

We investigated the effects of adjacent channel interference on the assignment results using the example presented in [6], which is based on the Philadelphia benchmark [8]. In this example, the cells are represented by circles of radius equal to 1 km (Fig. 3), contrarily to the original Philadelphia benchmark where the cells are hexagons. In Fig. 3, every cell is labelled with its corresponding number of channel requirements. The resulting assignment’s span and order (the difference between the minimum and the maximum used channel, and the total number of the used channels, respectively) are summarised in Table 1.

changing the BS locations, we varied the radiuses of the cells to form different network configurations ranging from nonoverlapping to highly overlapping, as depicted in Fig. 4. When the radiuses are equal to 1 km, we obtain a transition configuration where cells are adjacent to each other but without overlapping, as shown in Fig. 4b. The assignment results for every configuration are summarised in Table 2.

base station location and coverage area (radius 500 m) 10 9 8 7 Y, Km

results we present were obtained using the DD(10) given by (7), and a CIR threshold a equal to 18 dB, unless otherwise indicated.

6 5 4

1 km 8

25

8

8

3

8

2 15

18

52

77

28

13

15

1 2

31

15

36

57

28

4

6

8 X, Km

8

10

12

14

a 10

13

8

base station location coverage area (radius 1 Km) 10

Fig. 3

Network structure of the Philadelphia-like example

9 8 7

Cochannel interference

Cochannel and adjacent channel interference

n¼0

n¼1

Y, Km

Table 1: Philadelphia-like example results

6 5 4

n¼2

Span

Order

Span

Order

Span

Order

275

275

448

345

448

345

3 2 1 2

4.2

Effects of intercell overlapping

Although intercell overlapping might be desirable for handover processes, an excessive overlapping might increase interference and therefore the values of the assignment’s span and order. This problem does not exist in networks with regular hexagonal layouts, but unfortunately, such ideal layouts are inexistent in real systems. In order to evaluate the impact of intercell overlapping, we selected 21 BS location points to form a regular network structure similar to the Philadelphia example, where every location point is the centre of a circular cell. Without 64

6

8

10

12

14

X, Km b base station location and coverage area (radius 1.25 Km) 10 9 8 7 Y, Km

In the case of only co-channel interference, we obtained an assignment’s span and order equal to 275 channels, the same result reported in [6]. However, interesting results are obtained when we consider also adjacent channel interference. In this case, both span and order increase to 448 and 345 channels, respectively. The same values are obtained whether n ¼ 1 or n ¼ 2, which suggests that only the first adjacent channels’ effects are significant and therefore they must be taken into account, whereas the second (and further) adjacent channels’ can be neglected. The overoptimistic assignment obtained with n ¼ 0 does not guarantee an interference-free network environment.

4

6 5 4 3 2 1 2

4

6

8

10

12

14

16

X, Km c

Fig. 4

Various network configurations

a Nonoverlapping network configuration b Nonoverlapping configuration with adjacent cells; this is the transition configuration c Overlapping configuration IEE Proc.-Commun., Vol. 153, No. 1, February 2006

Table 2: Intercell overlapping results

Table 3: Urban propagation results with DD(1)

Radius, m

Model

Cochannel interference n¼0 Span

Cochannel and adjacent channel Interference n¼1

Order

Span

Radius, m

n¼2 Order

Span

Order

Cochannel interference

Cochannel and adjacent channel interference

n¼0

n¼1

Span LBd–4

Order

Span

n¼2 Order

Span

Order

200

229

77

229

77

229

77

100

229

86

229

86

229

86

500

229

145

229

146

260

189

200

229

188

229

208

229

208

1000

240

240

330

330

338

337

250

237

237

325

325

457

276

1250

280

280

449

344

449

345

300

279

279

457

329

457

329

1500

331

331

439

365

481

374

400

360

360

439

388

601

375

2000

409

397

825

438

865

449 100

229

91

229

91

229

91

200

229

189

229

211

229

214

250

241

241

342

342

474

318

300

281

281

457

329

457

329

400

360

360

435

397

599

394

WIM

The results show that the radiuses of the cells have no significant impact on the assignment’ span when there is no overlapping, and only the assignment’s order increases significantly as the radiuses of the cells increase. This can be seen in the configurations with radiuses of 200 m and 500 m for n ¼ 0, in which the same assignment’s span of 229 channels was obtained in both cases, contrasting with an assignment’s order of 77 and 145 channels, respectively. Moreover, adjacent channel interference can be neglected in nonoverlapping configurations. These facts change drastically in overlapping network configurations, where adjacent channel interference must be taken into account to guarantee an interference-free environment. In highly overlapping configurations (radiuses of 1.5 km and 2 km), even second adjacent channels (n ¼ 2) affect the assignment results. However, in moderate overlapping networks (radius of 1.25 km), the assignment results are practically the same if we take into account the first and second adjacent channels (n ¼ 2) or only the first ones (n ¼ 1). Since the vast majority of networks have a moderate intercell overlapping because nonoverlapping and highly overlapping configurations are practically inexistent in real systems, the results obtained reaffirm the conclusion of the previous Section: for practical purposes, co-channel and first adjacent channel interference should be taken into account for CIR computation, whereas further adjacent channels’ effects can be neglected.

4.3

Table 4: Urban propagation results with DD(10) Model

LBd–4

WIM

Radius, m

Cochannel interference

Cochannel and adjacent channel Interference

n¼0

n¼1

n¼2

Span

Order

Span

Order

Span

Order

100

230

105

230

105

230

105

200

229

188

229

209

229

209

250

238

238

328

328

459

278

300

278

278

457

329

457

329

400

360

360

439

388

601

375

100

231

113

231

113

231

113

200

229

189

229

216

229

216

250

242

242

341

341

467

327

300

281

281

457

329

457

329

400

360

360

435

397

599

394

Urban propagation

The propagation model L  d 4 is a fair approximation in urban macrocellular environments. Nevertheless, more precise estimation of interference might be necessary when the scarcity of available channels is critical and overestimation of interference leads to waste of spectrum. Alternatively, underestimation of interference might produce assignments in which harmful interference levels still remain; also in this case, more accurate estimation of interference might be necessary. We incorporated the WIM [12] to our heuristic algorithm and obtained the assignment plans for different network configurations using the structure of the previous Section, but reduced to a scale 1 : 4; e.g., the same configuration of Fig. 4b will be obtained here with a radius equal to 250 m. We varied the radiuses to obtain several configurations ranging from nonoverlapping to highly overlapping, and obtained the assignment plan for each case using both, the L  d 4 and the WIM models. The input data assumed for the WIM are: f ¼ 850 MHz, hb ¼ 40 m, hm ¼ 2 m, hB ¼ 20 m, b ¼ 30 m, w ¼ 15 m, and f ¼ 901. The results are summarised in Tables 3 and 4 for DD(1) and DD(10), respectively. IEE Proc.-Commun., Vol. 153, No. 1, February 2006

In this example, the L  d 4 model underestimated the interference with respect to the WIM’s predictions in most of the cases, which resulted in a slight difference of assignment’s span and a more significant difference of assignment’s order for the co-channel interference case in the overlapping configurations. This difference is more significant when adjacent channel interference is considered. For nonoverlapping configurations (radiuses of 100 m and 200 m), the same assignment results are obtained with both propagation models. Though the assignment’s span and order obtained with the L  d 4 model are smaller than those obtained with the WIM, it does not necessarily mean that they are better solutions. It is important to remember that the optimal solution of the CAP consists of the assignment that best balances the usage of spectrum and interference in the network to satisfy the traffic demand. Nevertheless, the use of a complex propagation model should be limited to the special cases mentioned above, as for practical cases the results obtained with the conventional power law model are very similar to those obtained with the WIM. 65

Regarding the performance of the algorithm, we can notice that the DD(1) provides better results than the DD(10) in this particular example. This suggests that in some cases it is not necessary to use complex difficulty degrees, since simpler ones can produce similar or even better results.

10 Y, km

4.4

15

Irregular network layouts

One of the important advantages of modelling the CAP with nonbinary constraints is the possibility of finding assignments in networks in which the cells have different sizes and shapes, and where no regular frequency reuse pattern exists. Conventional CA techniques cannot be applied directly to this kind of irregular network layouts. The work in [13] suggests the possibility of finding an optimal regular reuse pattern in irregular network layouts, which might enable the application of conventional CA algorithms; however, to the best of our knowledge, no further results in this direction were reported. We applied our heuristic algorithm to a practical example based on the structure of a real network in the city of W.urzburg, Germany [14–16], depicted in Fig. 5. In this example, the spatial traffic distribution is represented by discrete points called demand nodes. Such points are the centre of an area containing a quantum of demand from a traffic point of view, accounted in a fixed number of call requests per time unit. Demand nodes are generated by a partition clustering algorithm [14] using available population data; thus, demand nodes are densely distributed in traffic hotspots and sparse in areas of low demand, as shown in Fig. 6. The so-called set cover positioning algorithm (SCBPA) [15] selects the optimal positions of BSs that maximise the number of covered demand nodes. The lines indicating the convex hull around the set of demand nodes covered by each BS delimit the size and shape of the cells (Fig. 5). The transmit power necessary for every BS’s antenna to cover its corresponding cell area is also indicated in Fig. 5. We assumed three cases of traffic load conditions where each demand node takes a value of 1, 3, and 5 calls/h [17], respectively, and computed the number of channel requirements per cell in all the cases (Table 5). The network operates with the global system for mobile communications at 1800 MHz (GSM1800) standard, which means that a maximum of 374 frequency channels (200 kHz wide-carrier per channel) are available [18]. We defined the matrices R

5

0

0

5

10

15

X, km

Fig. 6

Spatial distribution of traffic demand nodes

. Table 5: Channel requirements of the Wurzburg-like example Cell

Number of channel requirements 1 call/h

3 call/h

5 call/h

1

10

22

33

2

9

19

29 13

3

5

9

4

25

61

97

5

5

9

13

6

6

11

16

7

14

31

47

8

7

15

22

9

5

9

13

and D as given by (16), and obtained the assignment plans for each case using the DD(1) and DD(10). The results are summarised in Table 6. This practical example once again reaffirms the conclusions elaborated upon in the previous Section regarding the adjacent channel interference. Moreover, we notice a slight better performance of the DD(10) for the case of only cochannel interference, but it is outperformed by the DD(1) when adjacent channel interference is also taken into

15 cell 1 (54 dBm)

. Table 6: Wurzburg-like example results

cell 3 (55 dBm)

cell 2 (56 dBm)

Cochannel interference n¼0

10 cell 5 (47 dBm)

Span

Y, km

cell 4 (61 dBm)

n¼1 Order

Span

n¼2 Order

Span

Order

DD(1) 1 call/h

cell 7 (61 dBm)

cell 6 (60 dBm)

5

cell 8 (59 dBm)

74

64

101

69

101

72

3 calls/h 182

144

241

157

241

159

5 calls/h 290

222

380

237

380

245

cell 9 (51 dBm)

DD(10) 1 call/h

0 0

5

10 X, km

Fig. 5 66

Cochannel and adjacent channel Interference

Network structure of the W.urzburg-like example

15

73

64

113

72

113

72

3 calls/h 181

144

270

157

270

157

5 calls/h 289

222

424

242

424

242

IEE Proc.-Commun., Vol. 153, No. 1, February 2006

account. However, the practical importance of the assignment plan of Table 6 for the heavy traffic-loaded case, i.e., assuming 5 calls/h, is very poor as the assignment span exceeds the maximum number of available channels. This situation generates a problem referred to as the fixed spectrum problem [3], which consists in finding a feasible assignment subjected to a limited number of available channels. In the case that an assignment without constraints violations is unfeasible, the solution of the fixed spectrum problem must admit some violations in order to satisfy the channel requirements of the network. For our particular way of modelling the CAP, the violation of the CIR constraint in order to obtain a lower assignment span consists in tolerating higher levels of interference. This is done by decreasing the CIR threshold a until a valid assignment is found. Thus, we reduced the CIR threshold a in steps of 1 dB and computed the assignment plan in each case until the span was equal or lower than 374 channels. The results using the DD(1) are summarised in Table 7. Notice that for a ¼ 17 dB the assignment plan is exactly the same as in Table 6, and only with a ¼ 16 dB the span in all the cases is below the number of available channels. A further reduction to a ¼ 15 dB leaves the previous results practically unchanged. This assignment meets perfectly the requirements of the GSM1800 standard, which recommends a CIR threshold of 15 dB for planning purposes [19]. This example showed the effectiveness of our algorithm at generating appropriate assignments with realistic constraints regarding number of available channels, cells size and shape, and transmit power. In this case, the transmit power was not the same for all BSs and this fact was taken into account by our algorithm. The impact of modelling different transmit powers into the channel assignment process can be shown by obtaining the channel assignment assuming the same transmit power for all BSs and comparing it to the results in Table 7. We obtained such an assignment for a ¼ 15 dB using the DD(1) (Table 8). We observe that assuming the same transmit power for all BSs produces significant lower assignment spans than those shown in Table 7 when both co-channel and adjacent channel interference are taken into account. For the case of only co-channel interference the results are practically the same. This means that the assignments in Table 8 might not

Table 7: Results for different CIR thresholds Cochannel interference

Cochannel and adjacent channel Interference

n¼0

n¼1

n¼2

Span

Order Span

Order Span

Order

a ¼ 17 dB 1 call/h

74

64

101

69

101

72

3 calls/h

182

144

241

157

241

159

5 calls/h

290

222

380

237

380

245

a ¼ 16 dB 1 call/h

74

64

83

65

83

69

3 calls/h

182

144

200

143

200

152

5 calls/h

290

221

317

221

317

233

a ¼ 15 dB 1 call/h

74

56

83

62

83

65

3 calls/h

181

124

200

138

200

145

5 calls/h

289

189

317

213

317

223

IEE Proc.-Commun., Vol. 153, No. 1, February 2006

Table 8: Results for equal transmit power in all BSs Cochannel interference

Cochannel and adjacent channel Interference

n¼0 Span

n¼1 Order

Span

n¼2 Order

Span

Order

1 call/h

75

63

83

64

83

64

3 calls/h

182

142

194

144

194

144

5 calls/h

290

219

307

222

307

222

guarantee an interference-free network environment. The difference between both results is due to a less inaccurate estimation of CIR values during the assignment process when the assumption of equal transmit power is made. 5

Conclusions

We have presented an improved heuristic sequential algorithm for fixed channel assignment that uses nonbinary constraints. The algorithm provides an acceptable carrierto-interference value at the worst case interference point of each and every cell using the same channel. We demonstrated that it is necessary to take into account adjacent channel interference for the computation of carrier-tointerference ratio, but for practical purposes we recommend to consider only the first adjacent channels since that is enough to obtain interference-free assignments. This algorithm solves the minimum span problem for a fixed carrier-to-interference ratio, i.e., the problem of finding an assignment to meet a given carrier-to-interference ratio threshold criteria. In the case that the algorithm produces an assignment with more channels than the limit assigned to a network operator, the given threshold should be reduced to tolerate more interference. We recommend reducing such a threshold in steps of 2 dB; lower values might not produce significant changes. The channel assignment must be computed using every new threshold value until a valid assignment is found. If the achieved carrier-to-interference ratio is below the standard value for a specific system, the network might not work properly and other techniques such as cell splitting should be applied. Further research in sequential algorithms for fixed channel assignment should not be devoted to the development of more elaborated difficulty degrees or different sequence arrange criteria, since we have demonstrated that simple difficulty degrees provide similar and even better results than complex ones in practical examples. With the development of metaheuristic algorithms that can solve assignment problems with binary constraints even to optimality in some cases, it seems that the only way to improve the usage of spectrum is by modelling the channel assignment problem using nonbinary constraints. However, the key to obtain maximum profit of the nonbinary constraints approach resides in the proper estimation of interference. Thus, research efforts must be concentrated in this direction. 6

Acknowledgments

The authors are grateful to the referees who reviewed the manuscript. Their comments were extremely helpful and improved the paper significantly. R. Cha! vez-Santiago thanks the National Council of Science and Technology (Conacyt–Mexico) for the financial support given under Graduate Scholarship number 130737. 67

7

References

1 Dunkin, N.W., and Allen, S.M.: ‘Frequency assignment problems: representations and solutions’, Tech. Rep. CSD–TR–97–14, University of London, UK, June 1997 2 Aardal, K.I., Van Hoesel, S.P.M., Koster, A., Mannino, C., and Sassano, A.: ‘Models and solution techniques for frequency assignment problems’. Tech. Rep. 01–40, ZIB, Berlin, Germany, December 2001 3 Smith, D.H., Allen, S.M., Hurley, S., and Watkins, W.J.: ‘Frequency assignment: methods and algorithms’. Proc. NATO RTA SET/ISET Symp. Frequency Assignment, Sharing and Conservation in Systems (Aerospace), Aalborg, Denmark, October 1998 4 Dunkin, N.W., Bater, J.E., Jeavons, P.G., and Cohen, D.A.: ‘Towards high order constraint representations for the frequency assignment problem’, Tech. Rep. CSD-TR-98-05, University of London, UK, June 1998 5 Hurley, S., Whitaker, R.M., and Smith, D.H.: ‘Channel assignment in cellular networks without channel separation constraints’. Proc. IEEE Vehicular Technology Conf. Fall 2000, Boston, MA, 24–28 September, 2000 6 Wu, J.-L.C., and Wey, L.-Y.: ‘Channel assignment for cellular mobile networks with nonuniform cells’, IEE Proc. Commun., 1998, 145, pp. 451–456 7 Tamura, H., Sengoku, M., Shinoda, S., and Abe, T.: ‘Channel assignment problem in a cellular mobile system and a new coloring problem of networks’, IEICE Trans. Commun., 1991, E 74, (10), pp. 2983–2989 8 Hurley, S., Smith, D.H., and Thiel, S.U.: ‘FASoft: a system for discrete channel frequency assignment’, Radio Sci., 1997, 32, (5), pp. 1921–1939 9 Yeung, K.L., and Yum, T.-S.P.: ‘Fixed channel assignment optimization for cellular mobile networks’, IEICE Trans. Commun., 2000, E83-B, (8), pp. 1783–1791

68

10 Gamst, A.: ‘Some lower bounds for a class of frequency assignment problems’, IEEE Trans. Veh. Technol., 1986, VT-35, (1), pp. 8–14 11 Ch!avez-Santiago, R., Gigi, E., and Lyandres, V.: ‘The time complexity cost of adjacent channel interference in a heuristic for frequency assignment’. Proc. IEEE Int. Symp. Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications, Beijing, China, 8–12 August, 2005 12 Lee, J.S., and Miller, L.E.: ‘Walfisch-Ikegami formula’ in ‘CDMA systems engineering handbook’ (Artech House, Boston, 1998), pp. 190–199 13 Bater, J., Jeavons, P., and Cohen, D.: ‘Are there optimal reuse distance constraints for FAPs with random Tx placements?’. Tech. Rep. CSD–TR–98–01, University of London, UK, March 1998 14 Tutschku, K., Leskien, T., and Tran-Gia, P.: ‘Traffic estimation and characterization for the design of mobile communication networks’. Tech. Rep. no. 171, Institute of Computer Science, University of W.urzburg, Germany, April 1997 15 Tutschku, K.: ‘Demand-based radio network planning of cellular mobile communication systems’. Tech. Rep. no. 177, Institute of Computer Science, University of W.urzburg, Germany, July 1997 16 Tutschku, K., and Tran–Gia, P.: ‘Spatial traffic estimation and characterization for mobile communication network design’, IEEE J. Sel. Areas Commun., 1998, 16, pp. 804–811 17 Ch!avez-Santiago, R., Raymond, A., Lyandres, V., and Kontorovitch, V.: ‘Effective base stations location and frequency assignment in mobile radio networks’. Proc. IEEE Int. Symp. Electromagnetic Compatibility, Istanbul, Turkey, 11–16 May, 2003 18 Noll, A.M.: ‘Global system for mobile (GSM) communication’ in ‘Introduction to telephones and telephone systems’ (Artech House, Boston, 1998), pp. 233–237 19 Walke, B., Seidenberg, P., and Althoff, M.P.: ‘Interference and cochannel interference’ in ‘UMTS: the fundamentals’ (Wiley, 2003), pp. 16–17

IEE Proc.-Commun., Vol. 153, No. 1, February 2006