#### Abstract

IEEE 802.16j is an emerging wireless broadband networking standard that integrates infrastructure base stations with multihop relay technology. Based on the idle mode operation in IEEE 802.16j, we propose a novel location management and paging scheme. It integrates the paging area-based and the timer-based location update mechanism. In paging area-based scheme, an idle mode mobile station updates when it moves to a new paging area. In timer-based scheme, an idle mode MS updates when the location update timer expires. In this work, we formulate the mathematical model to evaluate the performance of the proposed paging scheme. A new random walk mobility model that is suitable for modeling in multihop relay network is created. Optimization of location update timer is also investigated.

#### 1. Introduction

IEEE 802.16 standard [1] (or WIMAX) is an emerging broadband wireless access system to provide users with high-speed multimedia services. The IEEE 802.16e standard provides mobility support for WiMAX system. Mobile Stations (MSs) are usually powered by battery. Paging mechanism and MS idle mode operation are defined to save power in mobile IEEE 802.16e system. Recently, the IEEE 802.16j Multihop Relay (MR) standard is proposed to support for multihop relay communications with Relay station (RS) [2–4]. IEEE 802.16j standard provides better network coverage and enhance system throughput performance. In 802.16j network, the base station is called Multihop Relay BS (MR-BS). Relay Stations (RSs) relay signaling and data messages between the MR-BS and the MS.

In WiMAX system, MS enters idle mode to save power when there is no data to transmit or to receive. Whenever an incoming data message arrives, the network applies paging mechanism to wake up the dormant MS. During idle mode operation, MS still needs to update its location occasionally so that network only needs to perform broadcast paging in selected cells when a data message arrives. Tradeoff between signaling cost and location precision of idle mode MS is the main design issue in paging and location update protocol design.

Conventional cellular network paging and location management design could be categorized as follows Location-based paging area schemes [5]: users update when they move across the border between different paging areas. Paging area might be overlapping or nonoverlapping. Time-based schemes [6]: users update periodically when the update timer expires. Distance-based schemes [7–9]: users update when moving a fixed distance away from the last updating location. Movement-based schemes: users update based on the number of passing stations. Velocity-based schemes: users update based on the velocity. Profile-based schemes [10]: users update according to their behaviors. Some schemes apply an integrated approach to reduce the signaling cost [11]. Paging for microcell/macrocell overlay is also studied [12]. Pipeline paging technique could be applied to reduce the paging delay [13].

In this paper, we propose a novel paging and location update algorithm that integrates timer-based scheme and location-based paging area scheme for IEEE 802.16j system. For performance evaluation, we investigate a random walk mobility model that is suitable to evaluate the mobility issue in multihop relay cellular network like 802.16j, as base stations and relay stations operate differently but coexist in this type of network. The mobility model is described and validated in Section 2. The paging scheme design is presented in Section 3. In Section 4, we evaluate the system performance analytically. The optimization of location update period is presented in Section 5. Performance results are presented in Section 6. Finally, we conclude the paper in Section 7.

#### 2. Mobility Model

Random walk model is widely used for modelling mobility in cellular networks [6, 11, 14, 15]. Markov chain formulation is used to compute the probability that MS movement. Labelling and grouping cells based on geometric symmetry reduces the complexity of the model. Akyildiz et al. proposed a random walk model for MS mobility in cellular networks [14]. In this model, MSs move in the hexagonal cell. The probability that MS moves to an adjacent hexagonal cell is a system parameter. When the MS moves to an adjacent cell, it has the uniform probability to move to one of the 6 adjacent hexagonal cells. The cellular random walk model is no longer applicable in multihop relay network as some cells are base stations and some are relay stations.

In the proposed model, the probability of MS movement from arbitrary cell to arbitrary cell could be computed while computational complexity is limited. The goals of the proposed random walk mobility model for multihop relay networks are to uniquely identify the relay station cells and simplify the mathematical model based on the symmetric property.

An MR-BS (multihop relay base station) or an RS (relay station) is located in the center of a hexagonal cell. Random walk mobility model is applied to characterize the movement of mobile stations (MSs). The Absolute Geographical Location is applied to uniquely identify the hexagonal cells. The Relative Moving Distance is applied to reduce the complexity of the random walk mathematical model. Rules of mapping between Absolute Geographical Location and Relative Moving Distance will also be described in this section.

##### 2.1. Absolute Geographical Location Model

The Absolute Geographical Location is used to uniquely identify the geographical location of each hexagonal cell. Unlike the random mobility model described in [14], hexagonal cells have to be uniquely labelled to distinguish MR-BS and RS. As shown in Figure 1, we apply oblique coordinates with axis and to label the hexagonal cells. Each cell is uniquely identified as . The origin of the oblique coordinate is where MR-BS is usually located.

##### 2.2. Relative Moving Distance Model

As described previously, the Absolute Geographical Location indicates the geographical cell location. Due to the symmetric property of random walk mobility model, the probability of an MS moving from cell to new cell is the same as moving from cell to . Thus, in terms of moving probability between cells, we can model that the moving probability by considering the probability of an MS moves from the origin to in the Relative Moving Distance model. The MS moving probability in the Relative Moving Distance model is the same as and in the Absolute Geographical Location model.

The Relative Moving Distance model is consisted of tiers of hexagonal cells. A 5-tier Relative Moving Distance model is shown in Figure 2. In the boundary of the wireless network, an MS may enter an outer cell and does not come back to the network. In the Markov Chain models, those outer cells will be modeled as absorbing states. The outer cells are the fifth tier of the network, which is denoted as *out*, as shown in Figure 2.

There are three axes (, , and ) across the origin , and the network is divided into six regions. A hexagonal cell is labelled as . The tuple is labelled based on the oblique coordinate system with axes and . Note that, for the cells in the same tier, the sum of and is the same and is equal to the tier number . Cells in Relative Moving Distance model are symmetric. This model provides mobility information for Absolute Geographical Location. An MS movement from to in the Absolute Geographical Location model will be transformed to an MS movement from to in the Relative Moving Distance model. Notice that the in Relative Moving Distance model are all nonnegative integer; hence, absolute value operation is taken during the transformation.

##### 2.3. Simplified Moving Distance Model

Since the six regions in the Relative Moving Distance model shown in Figure 2 are symmetric in terms of MS moving probability, we could further simplify the moving distance model. Figure 3 illustrates the Simplified Moving Distance model, which is actually the Region I of the original Relative Moving Distance model. The cell in the Simplified Moving Distance model is denoted as , where are non-negative integers and

##### 2.4. Rules of Mapping

We will describe a set of mapping rules that transforms the relative moving distance to the absolute geographical location. Because of the Markov property, the future MS movement depends only on the current location state. In the Relative Moving Distance model, a mobile station always starts from as we proposed this *relative* mobility model for movement from the current location of the MS. The coordinate space is considered to be shifted so that the origin of the coordinate space is centered at the current MS location.

We observe the geometric property of the hexagonal topology to create 3 mapping rules to simplify the model. We classify the 6 regions in Figure 2 based on the geometric properties. Region I and IV will apply Mapping Rule I. In Regions I and IV, we find that is always true. Regions II and V will apply Mapping Rule II. In Regions II and V, we find that and is always true. Regions II and VI will apply Mapping Rule III. In Regions II and VI, we find that and is always true. Based on the geometric property, these 3 classifications of mapping rules will be discussed in Theorems 1, 2, and 3, respectively.

Moving from to in a given time interval is transformed to moving between and in the same time interval. If a user starts at and locates in after unit time, the probability is equal to that of moving from to after unit time. We define as the probability that an MS moves from to after unit time:

In the Relative Moving Distance model, three axes divide the network into six regions. As the Relative Moving Distance model applies an MS-centric view that considers relative movement from the starting location, the MS movement is always starting from . The MS movement in the original Absolute Geographical Location from to is equivalent to the transformed MS movement from to . The movement to could be classified based on values of and . The classification of the mapping rules also corresponds to the mobile movement in the six regions shown in Figure 2.

All cells in Regions I and IV have the property The relative movement vector can be denoted as a linear composition of two axes and with integer coefficients and

In the Simplified Moving Distance model , and are non-negative integers. We solve the above equation and derive the non-negative solution by taking absolute values , . Since in the Simplified Moving Distance model, as shown in Figure 4, is the larger one among and while the smaller one is .

For example, as shown in Figure 4, can be decomposed as the linear combination of and . Notice that the moving probability from to is the same as the moving probability from to

From observation, the Mapping Rule I maps the absolute geographical location to the relative moving distance model in Regions I and IV, as shown in Figure 2. Notice that the relative moving values and are both positive values (in Region I) or both negative values (in Region IV).

Theorem 1 (Mapping Rule I). *While *

If an MS moves to cells in Region II or V in Figure 2, the following two properties hold: and . The relative movement vector can be denoted as a linear combination of and . We can get , by solving the equation:

Theorem 2 (Mapping Rule II). * While and *

If an MS moves to Region III or VI, the following two properties hold: and . The relative movement vector can be denoted as a linear combination of and We can get , and by solving the equation:

Theorem 3 (Mapping Rule III). * While and *

An example of mapping movement to Region III is shown in Figure 5. The left part of the figure is the Absolute Geographical Location. An MS moves from to . The right part of figure is the equivalent Relative Moving Distance. Considering the starting point as the center of the map, the destination is in Region III. Applying Theorem 3 and setting and , we can obtain .

##### 2.5. Calculation of User Movement Probability

The user movement is modelled by the random walk mobility model. As described previously, the computation of MS movement could be simplified by the transformation and mapping to the Simplified Moving Distance Model . The mobile network model has tiers of cells. The value of must be large enough so that the probability of users moving outside is small. Depending on the requirements of modelling various mobility protocols, the value of should be selected accordingly.

To further simplify the notation, we map each cell in Simplified Moving Distance Model to a new state as shown in Figure 6. The states are relabelled from inner cells toward outer cells. For example, the origin is denoted as Likewise, the is denoted as and so forth. As we observed, the relabeling based on geometrical symmetry could be used to simplify the following mobility model formulation. A discrete-time Markov Chain model, as shown in Figure 7, is created to compute the MS movement probability. We denote the probability that an MS stays in the same cell in the next time slot as . The probability that an MS moves to a neighboring cell in the next time slot is thus which is denoted as In the random walk model, the MS has probability to stay in the same cell and to move to another adjacent cell (notice that there are 6 neighboring cells). By observation of the geometric properties of the hexagonal topology, the random walk mobility could be formulated as the Markov Chain shown in Figure 6.

We define the matrix to represent the probability that an MS is in state after unit time slots. The size of an -tier network is denoted as Hence, the size of is 1 by

In the relative moving model, the initial location of an MS is at the origin at time 0. The initial state is described as follows:

The probabilistic transition matrix in the Markov Chain model is denoted as It is an by matrix. As shown in Figure 2, the number of tiers in the hexagonal topology is symmetric. We derive the value of based on observing the geometric property of the hexagonal network topology:

By observing the mobility symmetry in Figure 6, the Markov Chain state transition diagram is drawn in Figure 7. Now, we will write down the state transition probability of the Markov Chain model of Figure 7 in matrix form. An element in is the probability of moving from one state to another state during one unit time in the Markov Chain model:

Based on the state diagram shown in Figure 7, the elements of can be obtained. From the definition, the Markovian state probability in time slot could be computed by iteratively multiply the current state probability with transition matrix. We can then calculate with iteratively:

The movement probability could be computed with (13). It multiplies by times. To reduce the computational complicity, we can diagonalize the matrix and derive matrix and is the diagonal matrix of eigenvalues. consists of the eigenvectors of . We can obtain the state probability quicker by applying (18):

##### 2.6. Validation of the Mobility Model

Similar to the previous work [14], we validate the mathematical model by simulation. The network tier is 3, and two mobility scenarios or are simulated. The movement probability values after time slots are shown in Table 1. Math 1 method is the result of computation based on (13). Math 2 method is the results of diagonalized computation based on(18)

The initial state probability matrix is 1 for the center cell and 0 for other cells:

The transition matrix is

The diagonal matrix is

The transpose of matrix is shown as the following:

We implement the Monte Carlo simulation in C++ to model the random walk mobility model in the hexagonal topology. Each MS has probability to stay in the same cell and probability to move to any adjacent hexagonal cell. Totally 1000000 simulation runs are conducted. The uniformly random walk mobility simulation results are compared with the Markov Chain analysis results. As shown in Table 1, the differences between the mathematical models and simulation results are always less than 1%. In addition, we observe that the diagonalized method effectively reduces the computation time.

#### 3. IEEE 802.16j Multihop Paging

##### 3.1. IEEE 802.16j Idle Mode

Idle mode operation reduces control signaling cost and MS energy consumption. An MS in idle mode periodically listens to the downlink broadcasting paging messages without registering to a specific BS. RSs relay all paging messages between MS and MR-BS. In this paper, we consider nontransparent mode operation in 802.16j system. Idle mode and paging operations are illustrated in Figure 8.

###### 3.1.1. Entering Idle Mode

Before entering idle mode, an MS sends Deregistration message (DREG-REQ) to the MR-BS. Then the MR-BS replies De/Reregister Command message (DREG-CMD) to MS. These two signaling messages are used to synchronize the paging listening time. For an MS serving by the relay stations, the access RS will relay all deregistration messages and paging messages between the MR-BS and the MS. Notice that the control signaling cost is multiplied by the number of relay hops in this scenario.

###### 3.1.2. Idle Mode Operation

As shown in Figure 9, there are two types of time intervals in idle mode operation: MS Paging Unavailable Interval and MS Paging Listening Interval. During MS Paging Unavailable Interval, an MS turns down radio interface to save power. In MS Paging Listening Interval, an MS listens to the downlink broadcast of paging advertisement messages (MOB_PAG-ADV). The listening interval has a period of *PAGING_CYCLE.* The *PAGING_OFFSET* parameter is used to separate MSs in different paging groups. An MS is synchronized to the periodic listening intervals based on the *PAGING_CYCLE* and *PAGING_OFFSET* given in a message.

###### 3.1.3. Termination of Idle Mode

At the end of MS listening interval, an MS must decide whether to leave idle mode or not. If an MS would like to transmit data, it must leave idle mode and enter active mode for normal operation. When an MS decides to terminate the idle mode, it will start the network reentry process by first sending Ranging Request (RNG-REQ) message to MR-BS. Then MR-BS will reply with Ranging Response (RNG-RSP) message to the MS. Then the MS can send the location update message and start the normal active mode operation. Relay stations will forward signaling messages, such as RNG-REQ and RNG-RSP, between MS and MR-BS when needed.

The paging operation is initiated when the system wants to find an MS. For example, a new data packet is arrived and is to be delivered to the MS. The network will check the paging information database that records the associated paging group of the to-be-paged MS. All MR-BS and access relay stations in the paging group will send broadcast paging message with the MS's MAC address. Once the MS receives the broadcast paging message, it will terminate the idle mode and go back to normal mode. The broadcasting is initiated from the MR-BS and is forwarded through relay stations.

##### 3.2. Paging Methods

In the network topology, MR-BS and RS are assumed to be located at the center of hexagonal cells. A cell is consisted of 1 MR-BS and 6 RSs as shown in Figure 10. Packets are either directly transmitted from MR-BS to MS, if an MS is located in the central cell, or forwarded through two-hop-relay transmission. When the network is going to page an MS, the paging message is forward from the MR-BS to the six RSs. Then the MR-BS and the 6 RSs will broadcast paging messages to MS (i.e., 7 transmissions are needed). Thus, the total signaling cost in one paging event is

Our paging scheme includes both paging area-based update mechanism and timer-based update mechanism. Several cells are grouped into one paging area. An MS roams between different paging areas and sends an update when it moves across the border. If a message arrives, the network only broadcasts the messages in one paging area to find the user. For example, the paging areas can be allocated as shown in Figure 11. There are totally 14 paging areas shown in this figure. In this example, one cell includes one base station and six relay stations, as shown in Figure 10. Notice that a hexagon that labelled with number has a base station, and other neighboring hexagons without number have relay stations, as shown in Figure 11.

Before an MS enters idle mode, the serving base station exchanges DREG-REQ and DREG-CMD messages with it. The last serving cell will be denoted as paging areas 1 as shown in Figure 11. In idle mode, an MS still needs to listen to paging-related information periodically. During every MS Listening interval, the MS listens to broadcast messages, which contains paging-related information. From this information, if the MS detects that it moves to a different paging area, it must notify the network about the paging area change. We call this update the Paging Area Notification (PA Notification). Hence, when a data message arrives, the network knows the right paging area to find the idle mode MS. When an MS moves to a new paging area, the MS will always first send update to the RS and then forwarded the signaling message to BS. In PA Notification, there are totally signaling cost, which is defined by the number of signaling message transmitted weighed by the number of hops to be forwarded. In the 2-hop multihop cellular structure, as shown in Figure 10, the PA Notification signaling cost is

In the proposed paging scheme, the paging area topology is MS-centric. When an MS updates the exact cell location to the network, the system recomputes the paging area, and the current cell becomes the centralized cell in the paging area, which is labelled with as shown in Figure 11. Paging areas will only be reset in two circumstances: data message arrival (and system create paging message to locate the MS) or timer-based update (timer expires after ).

The first case occurs after the data messages arrive, and the network starts the broadcast paging procedure. All cells in the paging area, where the MS located, will send broadcast paging messages. The second case is timer-based location update. If no message arrives after time slots, the MS must update its location to avoid losing track of its location. After timer expires, the MS goes into active mode, updates its location, and resets the paging area (set the current cell as paging area 1) before it enters idle mode again.

#### 4. Paging Performance Analysis

Signaling cost in wireless network paging design is critical. In this section, we will investigate the signaling cost in the proposed paging scheme. The 802.16j paging cycle strucutre is shown in Figure 12. MS Paging unavailable interval and MS Paging Listening Interval appear alternatively. In MS Paging unavailable interval, the MS enters idle mode and does not receive packets from the network. In MS Paging Listening Interval, the MS listens to the paging channel to find whether paging messages are sent. The process could be modeled as discrete events including MS movement and paging arrival occur at the MS Paging Listening Interval. For performance evaluation, we compute the probabilities of the MS movement events and paging arrival events accordingly.

##### 4.1. Interrupted Versus Uninterrupted Idle Periods

We denote the overall time duration as . During this time, we could further categorize the time period into two types: interrupted idle period and uninterrupted idle period.

###### 4.1.1. Interrupted Idle Period

A paging message arrives and terminates an interrupted idle period. We calculate the the number of interrupted Idle periods during the total duration and the number of uninterrupted idle periods during The paging message arrival follows Poisson random process with rate Hence, the expected number of paging message arrival during time is The number of interrupted idle periods is

###### 4.1.2. Uninterrupted Idle Period

No paging messages arrive during an uninterrupted idle period. An uninterrupted idle period is terminated due to the timer-based forcing update. The mobile-centric location area is reset after timer-based paging area update period One additional cycle for active mode operation for the location area reset is needed. Thus, the length of an uninterrupted idle period is cycles. We denote the time duration from entering idle mode to the paging arrival time as The expected value of is denoted as

During total duration the expected interrupted time periods is cycles. So the number of uninterrupted idle periods is the remaining uninterrupted time during divided by the duration of an uninterrupted idle period. The expected number of uninterrupted idle periods is

In an interrupted idle period, the signaling messages include paging and location updating. In an uninterrupted idle period, the signaling messages only include location updating at the end of the period.

##### 4.2. Broadcast Paging

Broadcast paging event only occurs during an interrupted period. If the call arrives between and cycle, the system broadcasts a paging message to the paging area where the MS locates. We can derive the probability of the MS in a paging area from the probability computation in Section 2.

The total paging signaling cost of one MS at cycle is “the probability of the MS in paging area ” multiplied by “the signaling cost in paging area ” We have calculated the paging signaling cost in one multihop cell, in (21). Thus, the total paging signaling cost is (the probability of the MS in paging area ) (the number of multihop cells in paging area ).

Based on the mobility model described in Section 2, we can readily compute the probability of an MS in an paging area after time For example, the paging area shown in the center of Figure 11 has 1 multihop relay cell, which includes 1 BS hexagonal cell marked with 1 and 6 RS hexagonal cells surrounding the BS. The probability of an MS is located with the paging area after time is . The mobility matrix that corresponds to paging area is denoted as . Similarly, is the matrix corresponding to paging area . Notice that only depends on the paging area topology and is independent of . Considering the whole wireless networks, we have

According to the random walk mobility model, the MS location state probability is . For each paging event, the signaling cost is . The cost of paging signaling during total time duration is

Similar to [9], we will compute The Poisson arrival is

The number of arrived paging message is denoted as . If it implies that no message arrives. If it implies that at least one message arrives. The probability that paging message arrival time falls between and as shown in Figure 12, is

Then, when , we calculate , the expected value of a message arrival time that falls between and [16]:

The average value is

Thus, from (26), the overall paging signaling cost is

##### 4.3. Paging Area Notification (PA Notification)

If the MS moves across the border between two different paging areas, the MS must notify the network about the PA change. The MS update signaling cost of each PA Notification event is denoted as The corresponding PA notification probability between cycle and is the summation of the probability across the paging area border, according to the previously described random walk mobility model and the paging area topology. There are totally signaling message transmissions when an MS updates.

Similar to the formulation, the mobility matrix for PA notification event, in which an MS moves away from paging area is denoted as Similarly, when we consider the whole network, we have as follows:

###### 4.3.1. Uninterrupted Idle Period

The update signaling during time to is In an uninterrupted idle period, there are totally MS Paging listening intervals, since an uninterrupted idle period is terminated by the timer-based update after time . In each MS Paging listening interval, the MS checks if PA changes. The expected PA Notification signaling cost in one uninterrupted idle period is: During the total duration the number of uninterrupted idle periods is

The total update singling cost in all uninterrupted idle periods will be

###### 4.3.2. Interrupted Idle Period

In an interrupted idle period, there are totally cycles, since an interrupted idle period is terminated by message arrival at time . The expected PA Notification signaling cost in an interrupted idle period is During the total time duration the number of interrupted idle period is . The total update singling cost in all interrupted idle periods will be

##### 4.4. Timer-Based Paging Area Update

Timer-based paging area update (Timer-Based PA Update) occurs when the update timer expires. The system recomputes the MS-centric paging area, as shown in Figure 11. In addition, the same MS-centric paging area recomputation occurs when an MS goes into active mode, which happens after a data message arrives. During the expected data message arrival is . As the signaling message flow is the same in the timer-based PA update and the paging due to data arrival, we will lump together the signaling cost into one term in this subsection.

The number of total PA update, which includes both Timer-Based PA Update and PA update due to data arrival, is . For each PA update, the signaling cost is denoted as Notice that the timer-based PA update signaling messages are the same as the PA notification signaling messages since similar signaling message flow is applied.

So the total timer-based PA update signaling cost is

From (23), (31), (34), (35), (37), and (39) the total signaling cost is

#### 5. Optimized Timer-Based Location Update

In the previous section, we derive the signaling cost given parameters , , , and . In this section, we will optimize the timer-based update period to minimize the overall signaling cost. The total time duration which is just an observation time period, does not affect the optimization results. We will normalize the formulation by defining After normalization of (40), we have

The matrix is an by matrix. After diagonizing the matrix, the matrix is composed of eiganvalues to and some constant values. To simplify the notation, we define and as follows: Notice that the parameters and are constants, for all . Then, the normalized signaling cost is

After substituting (28) and (30) for and and some computation, we could obtain

To find the optimized , we take the first-order derivatives: where denotes

By solving we will get the optimal paging area update timer

#### 6. Performance Evaluation

The PA Notification signaling cost decreases as increases because, in our paging area topology, the size of paging area near the center is smaller than the size of paging area away from the center. As expected, the timer-based PA update signaling cost decreases as increases. As increases, the low PA update frequency reduces the signaling cost; however, the location tracking of MS becomes coarser. The broadcast paging signaling cost depends on the data message arrival rate In addition, if an MS goes to *outside* state of the paging area, mostly due to infrequent paging area update, the network needs to broadcast the whole network to locate the MS. Tradeoffs between frequency of paging area update and the broadcasting cost could be observed in the figures.

##### 6.1. Finding Optimized Location Update Timer

In Figures 13, 14, 15, and 16, we illustrate the signaling cost of the proposed IEEE 802.16j paging scheme in different mobility scenarios and paging arrival scenarios. In each figure, the three signaling cost components, PA Notification, Broadcast Paging, and Timer-Based PA Update, are shown, respectively. The optimal value of the total signaling cost is also labelled.

Figures 13 and 14 show the performance differences between a high message arrival rate () scenario and a low message arrival rate scenario. In the three signaling cost components, the broadcast paging cost changes the most. With small the signaling cost grows more steeply as increases. The reason is that the broadcast paging signaling cost becomes large when the MS location is updated infrequently. When an MS receives a message more frequently, it goes into active mode more frequently. When an MS goes into active mode and then reenters the idle mode, the paging area is updated. Consequently, the MS less likely goes to outside area.

Comparing Figures 13 and 15, the mobility parameter differs. Notice that a high indicates the low mobility scenario since defines the probability that an MS stays in the same cell during unit time. When MS mobility is high, the optimal is smaller to keep the needed precision of location tracking.

In Figure 16, as the MS mobility is low and the data message arrival rate is high, the probability that an MS stays in the central region of the paging area is high. The probability of an MS that moves out of the paging areas is low; hence, the cases of network-wide broadcasting to find MS rarely occur. The signaling cost of Broadcast Paging component is relative flat, compared with the other three figures. Hence, the overall signaling cost is related flat when is large.

##### 6.2. Comparing to Pure Timer-Based Scheme

In addition to update the paging area topology when an MS does not update its location for time , an MS notifies the network when an MS moves across the border of paging areas in the proposed paging scheme. On the contrary, a pure timer-based paging algorithm might only update the location of an MS only when the timer expires. In Figure 17, we compare the proposed scheme and the pure timer-based scheme. The proposed scheme has a lower signaling cost than the pure timer-based scheme as shown in the figure.

#### 7. Conclusion

In this work, we investigated the paging and location management scheme in the IEEE 802.16j multihop relay networks. The paging scheme is compatible with the idle mode operation in the IEEE 802.16j standard and integrates with the paging area design and timer-based location update mechanism scheme. We propose a generalized random walk mobility model that is suitable for investigating user mobility in multihop cellular relay system, for example, IEEE 802.16j. The analytical mobility model is shown to match the simulation results. We applied this random walk mobility model to analyze the proposed paging scheme. The proposed scheme performs well compared to naive timer-based scheme. In addition, the proposed paging area update optimization has been shown to minimize the signaling cost effectively. In the future, we plan to further invest age advanced paging and location update algorithms to further enhance the signaling cost and paging delay. Moreover, nonrandom-walk mobility model for IEEE 802.16j is an interesting future work item to study. Advanced paging and location update scheme over generalized user mobility model will play a critical role in optimization the IEEE 802.16j relay network.

#### Notations

Absolute Geographical Location | |

Transition matrix | |

Probability of MS stay in the same cell at next cycle | |

Total time of MS operation | |

Probability of MS stay in a different cell at next cycle | |

Relative Moving Distance | |

Index of a Relative Moving Distance cell | |

: | Number of tiers in Relative Moving Distance topology |

: | The number of states |

: | Probability at state |

Number of states at the th tier | |

: | Each state’s probability at time |

After time slots of the timer-based location update | |

The timer for timer-based location update | |

Diagonal matrix of | |

Eigenvector matrix of | |

Total paging signaling cost matrix | |

: | Paging signaling cost matrix in paging area |

Total update signaling cost matrix | |

: | Update signaling cost matrix in paging area |

Total signaling cost | |

: | Coefficients of eigenvalues in matrix |

: | Eigenvalue of matrix |

: | Coefficients of eigenvalues in matrix |

: | Number of interrupted idle period |

: | Number of uninterrupted idle period |

: | Signaling cost in 1 MS update operation |

: | Signaling cost in 1 broadcast paging operation |

: | Signaling cost in 1 paging area update operation |

: | Message arrival time |

: | Average message arrival time |

: | Message arrival number |

: | Total signaling cost during |

: | Normalized signaling cost in one time slot |

: | The optimal value of update timer . |

#### Acknowledgment

This work was partly supported by the Industrial Technology Research Institute (ITRI).