Abstract

Following the recent advances in the Internet of Things (IoT), it is drawing lots of attention to design distributed algorithms for various network optimization problems under the SINR (Signal-to-Interference-and-Noise-Ratio) interference model, such as spanner construction. Since a spanner can maintain a linear number of links while still preserving efficient routes for any pair of nodes in wireless networks, it is important to design distributed algorithms for spanners. Given a constant as the required stretch factor, the problem of our concern is to design an efficient distributed algorithm to construct a -spanner of the communication graph under SINR such that the delay for the task completion is minimized, where the delay is the time interval between the time slot that the first node commences its operation to the time slot that all the nodes finish their task of constructing the -spanner. Our main contributions include four aspects. First, we propose a proximity range and proximity independent set (PISet) to increase the number of nodes transmitting successfully at the same time in order to reduce the delay. Second, we develop a distributed randomized algorithm SINR-Spanner to construct a required -spanner with high probability. Third, the approximation ratio of SINR-Spanner is proven to be a constant. Finally, extensive simulations are carried out to verify the effectiveness and efficiency of our proposed algorithm.

1. Introduction

The Internet of Things (IoT) has attracted great attention in recent years, owing to its potential military and civilian applications [1, 2]. Such a network generally consists of a large number of autonomous network nodes, in which algorithms are usually distributed since these algorithms have to work without global information and coordinated central control. Hence, there is an imperative need to design efficient distributed algorithms for various network optimization problems in the IoT.

Constructing a -spanner with a minimum number of edges is one of the fundamental network optimization problems since the spanner property is a critical requirement of topology control in the IoT [3]. The IoT is commonly modeled as a graph , in which is the set of wireless nodes and represents the set of communication links (edges) connecting the nodes in . A spanning subgraph of is called a -spanner, for if for all pairs of nodes , the length of the shortest path from to in is at most times of that in . Here, is called the stretch factor. A spanner can not only decrease the number of links and maintain connectivity but also ensure that the length of a path between any pair of communication nodes is within some constant factor from the shortest possible one. Therefore, constructing a spanner of the communication graph is enormously helpful for topology control, geographic routing, and compact routing in the IoT [3].

Spanners have been extensively studied in computational geometry [4], in which is generally the complete Euclidean graph of the node set in the Euclidean plane. However, even in the field of computational geometry, computing a minimum stretch factor spanner using not more than a given number of edges is NP-hard [5]. In the IoT, current spanner construction algorithms either do not consider interference, such as [6], or only handle it under the protocol interference model [7].

However, when multiple nodes send messages at the same time, a node may be unable to receive the message from its given sender owing to the interference caused by simultaneous transmissions. The protocol interference model and the SINR (Signal-to-Interference-and-Noise-Ratio) interference model are the commonly used interference models. SINR can take into account cumulative interference of wireless communications and is more realistic [8], thus has been widely adopted now. However, designing and analyzing algorithms are challenging under the SINR model since each given receiver should compute accumulated interference generated by all other senders at the same time.

1.1. Outline of the Problem

In this paper, we consider a general case of constructing a -spanner under SINR, namely -spanner-SINR. That is, given a constant as the required stretch factor, our objective is to design an efficient distributed algorithm for constructing a -spanner of the communication graph to minimize the delay of constructing the -spanner. The delay of constructing a -spanner is defined as the time interval from the start time-slot that nodes start to work to the last time-slot that all the nodes finish their task of constructing the -spanner.

A large scale IoT discussed in this work consists of sensor nodes, with uniform transmission powers, deployed randomly and uniformly in the two-dimensional Euclidean space. Nodes act in synchronous rounds; in every communication round, a node can transmit a message and attempt to receive a message. Each node initially knows only its own unique ID and its own coordinates. Since we adopt the SINR interference model when nodes communicate, the predesigned receivers can successfully decode the messages if and only if SINR constraints are satisfied. Owing to the accumulation and uncertainty of the SINR model, it is challenging to design a -spanner distributed algorithm based on SINR in wireless networks. How to make more nodes transmit simultaneously and meanwhile make their given neighbor nodes successfully decode the messages is crucial to the performance of the algorithm under SINR.

The authors in [9] proposed a distributed algorithm SINR-Undirected-YG under SINR and proved that the resultant graph is a -spanner. Zhang et al. [9] claim to be able to construct a spanner at time-slots, but the running time of its algorithm is very large during the simulation. In this paper, we study a general case of constructing a -spaner under SINR and try to reduce the delay.

1.2. Summary of Contributions

The summary of contributions of this paper is as follows: (1)We identify the general case of constructing a -spanner under SINR (-spanner-SINR problem), and we design an efficient distributed algorithm SINR-Spanner to construct a required -spanner under SINR with high probability, i.e., with a probability of at least , where is the total number of nodes in the network. Moreover, the resultant -spanner has edges(2)Our distributed algorithm SINR-Spanner is also a local algorithm, in which the topology can be locally and self-adaptively maintained based on the information from the neighbor nodes without affecting the whole network. We reasonably utilize one kind of proximity graph—Yao graph (YG)—to construct spanner under SINR. YGs divide the surrounding area of each node into sectors of equal angles and add edges only to the nearest neighbor within each sector [10]. If there are two or more nearest neighbors in a sector, one can choose the first neighbor receiving the message. In our design, each node is capable of independently performing successful local broadcasts to collect its neighborhood information within a certain range, such that it can get the nearest neighbor in each sector and the resultant -spanner is a special YG(3)We introduce the definition of proximity range and proximity-independent set (PISet) to increase the number of nodes transmitting successfully at the same time and to reduce the delay in the SINR-Spanner algorithm. The approximation ratio of SINR-Spanner is proven to be a constant(4)Extensive simulations are carried out to verify the effectiveness and efficiency of our proposed distributed and randomized algorithm

The rest of this paper is organized as follows. Section 2 reports the most related work. Section 3 precisely defines the formulation of the problem and introduces relevant models and notations. The spanner construction algorithm SINR-Spanner is presented in Section 4. Section 5 gives a theoretical analysis of the algorithm. In Section 6, we evaluate the performance of the algorithm via simulation. Section 7 concludes the paper with suggestions for future work.

In this section, we first investigate the spanner algorithms in computational geometry and in the IoT, then discuss the method of applying randomized and distributed solutions under SINR.

2.1. -Spanner

The book [4] by Narasimhan and Smid is a comprehensive overview of geometric spanners. For geometric spanners, several structures and methods have been proposed, such as the Greedy method [11], Well-Separated Pair Decomposition method, Delaunay triangulation, -graphs, and YGs. Constructing YGs is one of the simplest ways of constructing -spanners. Yao [10] used YGs to simplify the computation of the Euclidean minimum spanning tree. Althöfer et al. [11] firstly proved that YGs are the -spanners for the corresponding complete graph. For the corresponding complete graph, YGs are -spanners with [12].

In wireless networks, the spanner property was first discussed by Li et al. in [6]. They modeled the network as a unit disk graph (UDG) and analyzed the energy stretch factor of several common subgraphs of a UDG: for the relative neighborhood graph (RNG), for the Gabriel graph (GG), and for YG. And these proximity graphs have been widely used in spanner construction as subgraphs of a UDG. There also exist spanner construction mechanisms for quasi-unit disk graphs, disk graphs, and unit ball graphs [13]. In [14], Kothapalli et al. proposed a local-control protocol for establishing a constant density spanner among a set of mobile stations. The LISE (low interference spanner establisher) algorithm was presented to establish a spanner in [7], where the interference definition is based on how many nodes are affected by the communication over a certain link. However, since the above spanner algorithms for wireless networks are all studied without considering interference or handling interference under the protocol interference model, they cannot deal with interference effectively under the SINR interference model. Zhang et al. [9] first consider spanner construction under SINR, and this work improves it.

Constructing spanners under a computational geometry field greatly promotes the study under the wireless network setting. Meanwhile, wireless network requirements, which generally need to efficiently satisfy various topological characteristics, encourage the development of geometric spanner construction. Recent results on sparse geometric spanners focused on satisfying one or multiple topological characteristics such as lightness, small degree [15], fault tolerance, no central agent [16], and multiple characteristics [17].

2.2. SINR Model

In wireless networks, the SINR model received increasing attention [8]. Despite the vast amount of researches in the design and analysis of centralized algorithms under SINR [18], few results are known about distributed solutions in this model, especially for global communication tasks, owing to the accumulation and uncertainty of interference.

There is a growing interest in developing randomized distributed solutions to local broadcast [19], which is defined as successfully transmitting a message to all neighbors in the corresponding reachable proximity of a node. Randomized distributed solutions to local broadcast are often used as a building block for global communication tasks, such as multiple-message broadcast [20], synchronization [21], and multiple channels broadcast [22]. In [20], selected leader nodes adopt local broadcast to collect the messages that arrive at their dominated nonleader nodes and then disseminate the received messages to the whole network. The algorithm in [21] starts from very low probabilities and increases them gradually until nodes can hear a reasonable number of messages and implement all nodes’ clock synchronization. In [22], the selected leader in different channel collects locally the messages of its dominated nonleader nodes in the same channel to speedup multiple channel broadcast. Note that most of the existing distributed and randomized works under SINR are still in the theoretical stage except [19] and Fuchs’s coloring study, such as [23]. Using local broadcasting as a basic unit, we propose an algorithm of spanner construction in this paper and perform simulations to verify the performance of the algorithm.

3. Model and Problem Formulations

Assume that a set of wireless network nodes, modeled as a graph , is deployed randomly and uniformly in a 2-dimensional geographic plane. Nodes act in synchronous rounds. Each node is conscious of its ID and coordinates and has the same transmission power (). Let the Euclidean distance for the two endpoints and , denoted by , be the length of an edge in . And let be the length of the shortest path between and in which is defined as the sum of the lengths of its edges.

A commonly assumed model for the propagation effect of wireless nodes is deterministic path loss, i.e., , where transmits a message to , is the received power at a receiver , and is the path loss exponent (typically, ). A deterministic path loss model is applied to the following interference model.

3.1. SINR (Signal-to-Interference-and-Noise-Ratio) Interference Model

In the SINR model, a transmission from node to node is successful iff the SINR condition holds: where is the set of transmitting nodes, is the path loss exponent depending on the network environment, is a hardware-defined threshold, and is the environmental noise.

The transmission range of a node is the maximum distance at which a node can receive a clear transmission from while no other node is transmitting at the same time, i.e., . The SINR condition (1) tells us that for the given power level . Note that is for only one node transmitting in the whole network at the time slot.

3.2. Local Broadcasting Range ()

We set the local broadcasting range where is a fixed model parameter.

We say a node transmits successfully in a time slot if it transmits a message, and this message is received by all its neighbors in a distance smaller or equal to in the time slot. In Section 4, we define the proximity range such that the nodes which are in a distance greater or equal to can transmit successfully at the same time slot.

We denote the region within of node as a local broadcasting region and the number of nodes in it as . Furthermore, let .

3.3. Communication Graph

The communication graph of a given network consists of all network nodes and edges such that . Since is for only one transmission in the whole network at the same time, we adopt a slightly smaller range as [22] which suffice for practical communication. In the communication graph , which is simply denoted by , a node is a neighbor of node if .

3.4. Yao Graph (YG)

The directed Yao graph with a fixed integer parameter is defined as follows. Any equally separated rays starting at the origin node define sectors. The orientation of the cut is identical for all nodes. Translate the sectors to each node . In each sector with a node , pick the shortest directed edge , if there is one, to . Ties are broken arbitrarily.

This implies that preserves the shortest outcoming edge in each sector. Accordingly, preserves the shortest incoming edge in each sector.

An undirected Yao graph, in which the edge directions are ignored, is denoted by . is the Yao graph in which only the edges whose lengths are no more than are preserved from . In other words, is the spanning subgraph of a complete Euclidean graph on node set and is the spanning subgraph of .

3.5. -Spanner

Let be a real number. A spanning subgraph of is said to be a -spanner of , if for any two nodes and in , the shortest path between and in , whose length is at most times that of the shortest path in , i.e.,

The constant is called the stretch factor of (w.r.t. ). Note that can be or a communication graph that is a spanning subgraph of .

3.6. Delay of Constructing a -Spanner

The delay of constructing a -spanner is defined as the time interval from the start time slot that the first node starts to work to the last time slot that all the nodes finish their task of constructing the -spanner.

Next, we present the definition of the problem, the -spanner under the SINR problem, which is our focus in this paper.

3.7. -Spanner under the SINR Problem (-Spanner-SINR)

Assume that a set of wireless network nodes deployed randomly and uniformly in a 2-dimensional geographic plane, in which each node is aware of its ID and coordinates, has the same transmission power () and a local broadcasting range ; given a constant , the goal is to design a distributed algorithm to find a -spanner of the corresponding communication graph under SINR, such that the delay of constructing a -spanner is minimized.

4. Algorithm

4.1. Algorithm Outline

The main idea of our algorithm is to construct for the -spanner-SINR problem. The reason is that the graph is a -spanner [9]. Consequently, given a constant , we compute a and construct under SINR in our distributed algorithm to get the required -spanner.

To construct , each node needs the node ’s information which is the closest to in the sector belongs. Accordingly, each node should locally broadcast its ID and coordinates to its neighbors. However, if all the nodes transmit together under SINR, no node will receive any message owing to the interference. To avoid collision, each node could locally broadcast one by one. However, in a distributed algorithm, there is not a centralized coordination for arranging nodes to transmit in sequence. Accordingly, each node can only transmit with the probability . Furthermore, in order to reduce the delay, it has an obligation to have as many nodes as possible to transmit simultaneously. Therefore, each node will make a range as the radius of its neighborhood circle region and next take the inverse of the number of nodes in its neighborhood circle region as the sending probability.

First, we will compute the range and give the node set in which all nodes can transmit simultaneously.

4.2. Proximity Range and Proximity Independent Set (PISet)

How to make more nodes transmit simultaneously and meanwhile make their neighbor nodes successfully decode the messages is crucial for the performance of the algorithms under SINR. Thus, we define the proximity range such that any nodes and can transmit simultaneously if . Intuitively, a tiny implies a high degree of channel utilization. We examine how to set a proper to guarantee SINR threshold and meanwhile the highest channel utilization degree. For clarity, we define the proximity range and proximity independent set as follows.

4.2.1. Proximity Range and the Corresponding Proximity Independent Set

Proximity range is a length, and a is a subset of that satisfies for . A is maximal if and only if for and . How to design ? The basic idea underlying the design is to ensure the nodes in the same transmit simultaneously and all their neighbors receive the message successfully.

We refer to the region within of node as proximity region and the number of nodes in it as . Besides, let .

Before designing , we first give an example for and . and are illustrated in Figure 1 where 19 nodes represented by solid circles are in the same . Note that this is maximal. The distance between the given receiver with the sender and the nearest of other senders, which is in Figure 1, is at least . In other words, in Figure 1.

By observing the above example, we give the specific relationship between and in the following theorem.

Theorem 1. Suppose that where and , the nodes in the same can transmit successfully at the same time slot.

Proof. Let . We begin by estimating the smallest value of when transmits to successfully and transmits simultaneously with other nodes in the same . Further, we prove that the nodes in the same can transmit successfully at the same time slot with the above . We say a node transmits successfully in a time slot if it transmits a message and this message is received by all its neighbors in a distance smaller or equal to in the time slot.
In order for to be able to receive the message from , we require .
Thus, Since the equations (2) and (3), where is a fixed parameter.
Now We derive the lower bound of the above formula. First, since is the maximum local broadcasting range of a node. Furthermore, if a node transmit together with as shown in Figure 1, produces the largest interference when and the given receiver have the closest distance . If we represent a link as a node as shown in Figure 2(a), for the nodes in the , the densest packing of interfering links is the hexagon packing [24] with edge length as shown in Figure 2(b). There are at most six nodes in the first layer, and the distance is with respect to the abstracting node . Furthermore, the distance between and any node in the th layer is no less than with the th layer having at most nodes. Since , where is the Riemann zeta function, considering that , then . Thus, we have where .
Therefore, to make valid, it is sufficient to have Therefore,
Hence, , where and

Conversely, if have the value as shown in the above, the nodes in the same , which is maximal as shown in Figure 1 or not maximal, can transmit successfully at the same time slot since the receiver power is no less than in (9) and the cumulative interference is less than .

Figure 3 depicts visually the relation between proximity range and local broadcasting range with different , , and , which is helpful to pick these values during subsequent algorithm simulation. These parameters are reasonably set as follows: the path loss exponent , the threshold , and . From Figure 3, one can see that increases when decreases, increases, and decreases. So we can get that is as about 2.57 times at a minimum as when , , and , while is as about 6.34 times at a maximum as when , , and .

4.3. Algorithm

From the above Theorem 1, the nodes in the same can transmit successfully at the same time slot while each node takes the value of as Theorem 1, so the sending probability for each node is set to the inverse of the number of nodes in its proximity region . Thus, we guess all its neighbors can receive the message successfully with high probability, and we will give the proof in the next section. The pseudo-code for node is given in Algorithm 1, which implies that each node runs it independently.

1: Initialize the stretch factor ;
2: Initialize ;
3: Compute the number of nodes in its proximity region, i.e., ;
4: Initialize ;
5: Initialize for ;
6: for to time-slots do
7: Send a message containing its ID and coordinates with probability , and remains listening with probability ;
8: while receiving a message from some node do
9:   the index of the sector to which belongs;
10:  if or the distance then
11:   
12   end if
13:    end while
14: end for
15: ;
16: for to do
17:  ifthen
18:  ;
19:  end if
20: end for
21: ;
22: for to time-slots do
23: Broadcast the incoming neighbor set with probability ;
24: while receiving a message from some node do
25:  if is the incoming neighbor of and then
26:   ;
27:  end if
28: end while
29: end for

Now, we describe the full operation of our distributed algorithm SINR-Spanner. Each node carries out the same operations and has its local memory. The algorithm consists of three parts. In part 1 (line 1-line 5), each node performs initialization work. Each node first computes the required number of sectors depending on the given stretch factor . Then, each node takes the inverse of as the sending probability. Thus how to obtain a priori is critical to our algorithm design. Some existing works, such as [19, 23], assume the availability of to facilitate the algorithm design and analysis. Such requirements are common under SINR in order to enable initial communication. Therefore, in our algorithm, we also assume that is available to simplify the presentation. In part 2 (line 6-line 20), each node obtains the incoming neighbor in each sector by receiving the neighbors’ messages, and thus, a directed Yao graph forms. Figure 4 shows an intermediate result after part 2 in which there are 6 nodes; the edge length is at most and . Here, , , , and so on. In part 3 (line 21-line 29), each node sends an acknowledgement message back to its incoming neighbor, and thus, an undirected Yao graph is constructed which is a -spanner we require. An example of is presented in Figure 5, in which edge directions are ignored from the graph shown in Figure 4. Here, , , , and so on. has local data , and has local data ; thus, the undirected edge is known by two endpoints.

Our distributed and randomized algorithm SINR-Spanner, which is different from the algorithms in [9], can solve the -spanner-SINR problem. In part 1, the algorithm first initializes the stretch factor to attain a -spanner. The sending probability for each node is set to the inverse of to reduce the delay; then, repeats randomized transmission for time slots, respectively, in part 2 and part 3.

5. Performance Analysis of Algorithm SINR-Spanner

5.1. The Delay of SINR-Spanner

In order to obtain the main result of this section in Theorem 8, we first prove that a graph is constructed with high probability after the algorithm SINR-Spanner terminates in time slots in Theorem 3; next, analyze why the resultant is the required -spanner of the communication graph in Theorem 3.

Now, we give a form of Chernoff bounds which can be found in [22] and some advanced textbooks, such as [25], for the proof of Theorem 3.

Lemma 2 [22, 25] (Chernoff bounds). Let and be independent Bernoulli random variables, and let and . Then, for any , it holds that And for ,

Theorem 3. A graph is constructed with high probability after the algorithm SINR-Spanner terminates in time-slots, where is the maximum of the node number in one proximity region and is a constant for the number of time slots for the initialization work.

Proof. Since is as about 2.57 times to 6.34 times as from Theorem 1 and the resultant graph is connected, for any node . The probability that transmits as the only transmitting node in its proximity region is The last inequality holds, since increases when increases, and it obtains the minimum value when is 2.
Since the nodes can transmit successfully if they belong to the same owing to Theorem 1, there are at least nodes transmitting together in one time slot. Consequently, there are at least nodes transmitting successfully together in one time slot.
Therefore, after time slots, there are at least nodes transmitting successfully in expectation. By Chernoff bound (Lemma 2), after time slots, the probability Therefore, the probability for all nodes transmitting successfully is Thus, after time slots in part 2 of the algorithm SINR-Spanner, all the node transmit its ID and coordinates successfully with high probability, i.e., each node obtains the “nearest” incoming neighbor in its proper sectors. Then, after a further time slots in part 3, all the node transmit acknowledgement messages successfully with high probability. As a result, a graph is constructed with high probability.
In addition, the number of time slots for the initialization work in part 1 is a constant, which is denoted by . So the algorithm SINR-Spanner terminates in time slots.
Obviously, the number of nodes in the proximity region is upper bound by . If , a graph is constructed with high probability in time slots.
Next, in order to prove that the resultant of the algorithm is a -spanner of the communication graph , we apply some conclusions about where is constructed from .

Lemma 4 ([11]). If is a set of points in the plane, and the integer , then the graph contains at most edges.

Lemma 5 ([12]). Let for the integer , is a -spanner of , where is the number of nodes in .
In the next, we show that is the required -spanner of the communication graph if is connected in the following lemma of our previous work [9]. Zhang et al. [9] first give the condition that is connected if and only if the longest edge in Euclidean minimum spanning tree of the node set is at most . Note that if the nearest neighbor of each node in every sector is within the local broadcast region, the graph is a and it is also a -spanner of .

Lemma 6 ([9]). Let . If is connected and the integer , is a -spanner of .
Now, we show that the resultant of the algorithm is a -spanner of the communication graph .

Theorem 7. Given a constant , setting , the resultant graph is the required -spanner with at most edges after the algorithm SINR-Spanner terminates.

Proof. As , the sector angle is not larger than . Hence, is a monotone decreasing function of , and its value is less than . As a result, increases as increases when and approaches as , i.e., is an injective and increasing function of . Hence, its inverse function is a decreasing function of . Therefore, given a constant , we will find a such that the Yao graph with the number of sectors of no less than is the required -spanner.
Since the number of sectors is an integer, and the number of edges becomes larger as the number of sectors increases from Lemma 4, is assigned to . Therefore, the resultant graph is the required -spanner. Furthermore, have edges by Lemma 4.

Based on Theorems 3 and 7, we can state the main conclusion of this section in Theorem 8.

Theorem 8. The distributed algorithm SINR-Spanner constructs the required -spanner with high probability, and the delay is time slots, where is the maximum of the node number in one proximity region and is a constant for the number of time slots for the initialization work.

Lastly, Table 1 illustrates the relation between the stretch factor and the number of sectors in theory. In SINR-Spanner, is the value in the second row given the corresponding in the first row. Note that is the upper theory bound of the stretch factor for with respect to so far, maybe a smaller is sufficient in practical.

5.2. The Approximation Ratio of SINR-Spanner Algorithm

The goal of -spanner-SINR problem is to find a suitable -spanner under SINR and to minimize the delay in the construction. Now, we try to give the approximation ratio of the SINR-Spanner algorithm. The basic idea of our algorithm is that each node should locally broadcast its ID and coordinates to its neighbors successfully under SINR. In order to reduce the delay, there should be as many nodes as possible to transmit simultaneously. As we see, in the SINR-Spanner algorithm, the nodes in different proximity region can transmit simultaneously. Now, we consider whether the nodes in a different local broadcasting region can transmit simultaneously, then discuss the approximation ratio of SINR-Spanner.

Theorem 9. The approximation ratio of the algorithm SINR-Spanner is bounded by .

Proof. The probability that any node transmits as the only transmitting node in its local broad region is Let be the total area of wireless nodes distribution. Assume that a node as the only transmitting node in its local broadcasting region can transmit successfully, the number of nodes transmitting successfully in one time slot is In fact, the assumption cannot be guaranteed by theory and the following simulation, i.e., even though a node is the only transmitting node in its local broadcasting region, the node may not transmit successfully owing to cumulative interference producing by other simultaneously transmitting nodes. So the solution adopted the above assumption is the low bound for -spanner-SINR problem.
From the analysis of SINR-Spanner algorithm, the delay of constructing a -spanner under SINR is mainly and inversely proportional to the number of nodes transmitting successfully in one time slot. Hence, the approximation ratio of the algorithm SINR-Spanner is The inequality holds owing to the following: for connectivity and decreases with increases when . When is large enough, replacing by does not cause much error and .
Since is at about 2.57 to 6.34 times as from Theorem 1 and Figure 3, the approximation ratio of SINR-Spanner is a constant.

6. Simulation

In the previous section, we theoretically prove that the algorithm SINR-Spanner performs well in the worst cases and the resultant graph is the required spanner. In this section, we conduct simulations to investigate the average performance of our algorithm.

Our simulations are coded in the Sinalgo simulation framework [26], which is for testing and validating network algorithms and abstracts from the underlying layers. We consider a square area of 1000 by 1000 and deploy nodes within this network region randomly and uniformly, where . The local broadcasting range varies in . The ambient noise . Figure 3 has given the ranges of the path loss exponent , the threshold , and , and the relation between them and the proximity range . mainly affects the sending probability, which affects the delay of our algorithm. After testing various values, no matter what values , , and have in the ranges, the variety of the performance and the delay in our algorithm is similar. Therefore, we adopted , , and in our following reported result. Accordingly, the transmission power owing to the equation (2) and (3). The setting of simulation parameters refers to [9, 19]. With the proof of Theorem 8, we know that each node runs rounds and the algorithm could get the required spanner with high probability. However, the way we adopted was that the algorithm terminates when the resultant graph does not change in 50 continuous timeslots in the simulation, i.e., in 50 continuous time slots no node can find a nearer neighbor within each sector. Over 100 runs of the simulations have been made for each reported average result.

Now, we first explore the stretch factor of the resultant graph (real stretch factor) in Figure 6. In (a), we analyze the influence of the number of nodes. The local broadcasting range was set to , and the number of sectors was set to . The maximum and average real stretch factor slightly reduced with the number of nodes increasing. In (b), we investigate the impact of with the number of nodes and the number of sectors . The real stretch hardly changes with . In (c), we analyze the influence of the number of sectors with and . The maximum and average real stretch factor decrease with increases. When , the stretch factor approaches one. From all figures in Figure 6, it can be seen that both the maximum and average real stretch factor is much smaller than the theory value in Table 1 with the same . Hence, given the required , we can choose a smaller according to practical experience. In summary, the stretch performance of the algorithm SINR-Spanner is better than expected.

Due to randomization of SINR-Spanner, the constructed graphs may not be perfect ; for example, the connection link to the nearest neighbor in a sector might be missing. But the algorithm performance is guaranteed with high probability, i.e., the probability for all nodes transmitting successfully and all edges in being reserved is bounded by . To verify this, we evaluate the number of missing edges in the resultant graph from the corresponding perfect in Figure 7. With a different number of nodes, different values of , and different number of sectors, we compare the number of missing edges in three sending probability cases including , (SINR-Spanner), and (the algorithm SINR-Undirected-YG in previous work [9]). When the sending probability is , the number of missing edges is a little bit more than the other two cases; the reason is that the nodes transmitting at the same time are more and mutual interference is a little bigger. However, the number of missing edges in all three cases is no more than three in total 100 runs, though each run has no less than 3000 nodes and edges. Simulations indicate that the similarity between the resultant graph and the corresponding with no-missing edges is close to in all three cases.

We then consider the average delay needed by SINR-Spanner in Figure 8 from the influences of in (a), in (b), and in (c), respectively. We still compare the results with three sending probability cases including , (SINR-Spanner), and (the algorithm SINR-Undirected-YG in previous work [9]). Whether the sending probability is , , or in (a), the delay grows with the number of nodes increasing. In (b), the delay increases with increasing when the sending probability is and . However, when the sending probability is , the delay does not change in (b) since it is only related to the number of nodes. In (c), the delay does not vary with increasing. In summary, from Figure 8, the delay mainly changes with the change of the sending probability, while the sending probability of or mainly varies with and . Moreover, the average delay needed by the algorithm in [9] is close to the case where each node transmits with the probability , and the average delay needed by SINR-Spanner is much smaller than that of previous work [9]. Finally, in theory, is the delay upper bound of the algorithm SINR-Spanner from Theorem 8 and is the upper bound by , while in the simulation, the delay is much smaller than and the algorithm can achieve reliable performance when the algorithm terminates when the resultant graph does not change in 50 continuous time slots.

7. Summary

In this paper, we present a randomized and distributed algorithm SINR-Spanner to solve the -spanner-SINR problem using small delay in the IoT, which has the following characteristics: (1) being a distributed algorithm, (2) considering the SINR interference model, (3) applying the YG idea, and (4) theory and simulation guaranteed. In future research, the delay performance of the spanner construction algorithm under SINR may be able to be improved by adopting a smaller proximity region. Other methods for spanner construction except YG are also worthy of investigating.

Data Availability

Our simulations are coded in the Sinalgo simulation framework [26], which is for testing and validating network algorithms and abstracts from the underlying layers. We consider a square area of 1000 by 1000 and deploy nodes within this network region randomly and uniformly, where .

Conflicts of Interest

The author(s) declare(s) that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grants (12071478, 61972404) and the Natural Science Foundation of Shandong Province under Grants (ZR2019ZD10, F060505)