Abstract
For realtime multicast services over wireless multihop networks, to minimize the energy of transmissions with satisfying the requirements of a fixed data rate and high reliabilities, we construct a conflict graph based framework by joint optimizing network coding and scheduling. Then, we propose a primaldual subgradient optimization algorithm by random sampling K maximal stable sets in a given conflict graph. This method transforms the NPhard scheduling subproblem into a normal linear programming problem to obtain an approximate solution. The proposed algorithm only needs to adopt centralized technique for solving the linear programming problem while all of the other computations can be distributed. The simulation results show that, comparing with the existing algorithm, this algorithm can not only achieve about 20% performance gain, but also have better performance in terms of convergence and robustness.
1. Introduction
Since Ahlswede et al. proposed network coding (NC) [1], many studies have shown that NC can not only increase the throughput significantly but also achieve better robustness. Actually, it can achieve network multicast capacity by using random network coding (RNC) over a multihop wireless network [2]. This has promoted many studies on RNC based distributed optimization algorithms. For example, Lun et al. proposed to decompose the multicast optimization problem into two subproblems [3]: one is to search for a NC subgraph with minimum cost by modelling as a linear or convex program problem, which can be solved by a distributed primaldual subgradient optimization algorithm; the other is to design a network coding scheme for the optimal subgraph obtained from the first subproblem; then a simple RNC solution can be employed for this subproblem. Similarly, Wu et al. proposed a distributed NC optimization algorithm over mobile ad hoc networks to minimize the energy for multicast services [4]. Lee and Vishwanath recently proposed a distributed algorithm for rate allocations to achieve the network capacity with minimizing the operation cost [5].
Another research direction is to introduce scheduling techniques in medium access control (MAC) layer into an optimization framework. The research in [3, 6] has shown that network performance can be improved significantly by optimizing a scheduling technique in MAC layer. Recently, using interference graph model, Jaramillo et al. have studied an optimization problem on resource scheduling when realtime and nonrealtime services coexist in wireless multihop networks [7]. In addition, they have also studied the optimal rate allocation problem under heterogeneous delay constraints [8]. Although these studies cover the resource allocation and scheduling optimization problem for realtime services, they do not introduce NC into the optimization framework.
Apparently, we must integrate the two problems mentioned above together for achieving the best overall performance. By taking both NC and scheduling into account jointly, Rajawat and Giannakis proposed a joint optimization technique to improve the wireless multicast throughput performance under strict delay constraints [9]. Using hyperarcs to model the natural properties of wireless multicasts, Traskov et al. proposed a conflict graph model to identify effective network settings for studying joint NC and scheduling optimization algorithms [10]. In fact, this conflict graph based framework can be used to build an interference model for those active nodes in wireless networks. In order to avoid interference, we can select an efficient scheduling policy by sampling stable sets in a conflict graph. The studies have shown that, comparing with the scheduling technique with simple orthogonal models, the joint scheduling and NC graph optimization algorithm can improve the multicast throughput significantly [10].
Unfortunately, the joint optimization framework proposed in [10] cannot be used for wireless realtime multicast services without changes, because realtime multicast services usually need networks to support a predefined fixed and qualified data rate. By addressing this issue, Lun et al. [3] and Wu et al. [4] proposed a NC subgraph based optimization framework with minimum cost for supporting fixed multicast rates. However, they do not integrate scheduling techniques into their optimization frameworks. Moreover, since the conflict graph based scheduling optimization problem built in [10] is a NPhard problem, they proposed a greedy algorithm with sampling the maximum weight stable sets. Nevertheless, there is a major drawback in this algorithm: it is very sensitive to the parameters of iterations such as the step size. In other words, the results often fall far short of the global optimum value in case that the parameters of iterations cannot be set properly. Since those parameters of iterations can only be chosen through trial and error methods, it is thus difficult to meet realtime requirements.
By addressing the problems mentioned above, to minimize the energy with satisfying the requirements of realtime multicast services in wireless multihop networks, we will construct a conflict graph based framework for designing joint NC subgraph and scheduling optimization algorithms. Here we would like to point out that this framework is mainly inspired from [10], but there are two important differences: it focuses on realtime services with a fixed and qualified data rate rather than nonrealtime services; the optimization target is to minimize the energy of transmissions rather than to maximize the network capacity. Afterwards, using Lagrangian relaxation, we will propose a joint optimization algorithm by sampling the maximum stable sets randomly in a conflict graph, which includes two steps: first, it randomly samples a certain number of maximum stable sets of the conflict graph at each iteration; then, it solves a scheduling optimization problem by searching for those random sets. This method transforms a NPhard scheduling subproblem into a plain linear programming problem so that it can be solved efficiently. More importantly, the accuracy can be adjusted not only through the parameters of iterations, but also by the number of random samplings. It thus can have good convergence and robustness.
The remainder of this paper is organized as follows. The network model with a directed hypergraph and its corresponding conflict graph is introduced in Section 2. A conflict graph based realtime multicast optimization framework is presented in Section 3. Using Lagrangian relaxation, a joint optimization algorithm is proposed in Section 4. The numerical simulation results are presented in Section 5. Finally, we conclude the paper in Section 6.
2. Network Model
Considering a wireless multihop network, we use a directed hypergraph (where denotes the set of nodes and is a collection of hyperarcs) to represent the model. We define as a hyperarc and as a set of neighbors of node . When node sends data, all nodes in are assumed to be within the receiving range. For any hyperarc , we have and . For each node, there are at most hyperarcs.
The scheduling problem studied here is a scheduling for all of the hyperarcs defined above. When scheduling multiple hyperarcs for transmissions, we must avoid the interference caused by those conflict nodes. The specific hyperarc conflict situation depends on the network interference model. We consider the following two commonly used protocol interference models: the primary interference model and the secondary interference model. It is assumed that each node can only receive data from one node every time in the primary interference model, while in the secondary interference model, besides the above constraints, it is also assumed that any receiving node can only receive data correctly when all other neighbors are in a dormant state. They are strictly defined as follows.
Interference Models. For any two simultaneous scheduling hyperarcs and , the necessary and sufficient conditions for no conflict are as follows:(1) and , ;(2)there is for the primary interference model; there is and , where is an empty set for the secondary interference model.
Note that in both the primary interference model and the secondary interference model parameters are defined symmetrically. Therefore, we can use an undirected graph to express the conflict among hyperarcs. In this paper, we use the method presented in [10] to formulate an undirected graph based conflict graph. It is defined as follows.
Conflict Graph . Given a directed hypergraph defined above, we can formulate an undirected graph , which represents conflicts among all hyperarcs under an interference model. The vertex set in is the collection of all the hyperarcs in . Each edge in represents the conflict between two connected vertices according to an interference model. That is, the vertices in represent the hyperarcs in ; each edge in represents a conflict between the two connected vertices. The two corresponding hyperarcs then cannot be scheduled at the same time.
Clearly, based on the above definitions and a specific interference model, we can easily formulate the corresponding conflict graph according to any hypergraph . For a conflict graph , we now define any subset that there is no edge connected between any two nodes as a stable set . A stable set can be indicated by a column vector of length , which is defined as
A maximal stable set is one that is not contained in any other stable set. A maximum stable set is a stable set of largest cardinality. The stability numbers of are the cardinality of the maximum stable set. The stable set polytope (denoted by CH_{SS}) is the convex hull of the incidence vectors of all stable sets of . For example, assuming that there are stable sets in and is the incidence vector of the stable set , then we have
Finally, we take an example to illustrate the notations. A directed hypergraph with five nodes is shown in Figure 1 and its corresponding conflict graph in Figure 2.
Note that both the primary and the secondary interference models lead to the same conflict graph. has a node for each hyperarc in ; that is, . Figure 2 shows that the stable set polytope for this example is the convex hull of the incidence vectors of the five stable sets , , and .
3. Optimization Framework
To integrate NC subgraph optimization problem with scheduling problem together in one framework, we need to define some essential NC variables. First, we define the rate of NC data packet injected to hyperarc by node as . It is referred to as the coded packet rate in this paper. Then, we use a collection , also known as a NC subgraph, to represent the coded packet rate injected to all the hyperarcs. For simplifying the optimization framework for supporting realtime multicast services requiring high reliabilities, we assume that the physical layer can ensure perfect reliable transmissions through appropriate power control, channel coding, modulation, and other methods when the coded packet rate on any hyperarc is not more than . This assumption is realistic for some emerging realtime services such as Tactile Internet, which requires an endtoend failure rate of 10^{−7} [11]. That is, when the transmission power is large enough within a certain range, this requires guaranteeing a very low error rate for each active link by suitable techniques. It thus leads to a reasonable assumption of perfect reliable links. Finally, we would like to point out that our framework can be extended for matching the realtime services with nonnegligible error rates, for example, a more realistic signaltointerferenceandnoise ratio (SINR) threshold based framework, which needs further studies in the future.
Then, we use the variable to indicate the transmission rate of the information flow transferred to the terminal when transmitting from to in case that the coded packet rate injected to the hyperarc is . Here represents the collection of all sink nodes in a multicast session. To simplify the notation, define the variable as follows: According to (3), the variable can be viewed as the transmission rate of the information flow transferred to the sink node when transmitting from to .
Now let us focus on the optimization target of minimizing the energy of realtime multicast. The objective function can be expressed as a function of the coded packet rate , which is defined similar to [3] asIn this function, represents the energy required when the coded packet transferred with the rate over the hyperarc . For realtime multicast services requiring low reliabilities, we can modify this function by introducing some factors for describing transmission costs under different link qualities. However, it is beyond the scope of this paper and left to future studies.
Note that the definition of the function in [3] is only suitable for a single communication link, while the definition used here is for the energy consumption of a hyperarc. Therefore, we assume that the value of depends on the link with the longest distance over a hyperarc . Let denote the distance between any two nodes ; then we have
Using the definitions above, we now can model the optimization problem of minimizing the energy of transmissions with a fixed multicast rate as a joint scheduling and NC subgraph optimization problem; that is,s.t.
Capacity constraints are
Flow constraints are
Scheduling constraints are
Note that the objective function in this optimization framework is a convex function and CH_{SS} is a convex set. Therefore, the optimization problem is a convex optimization problem. The three constraints are explained as follows.(1)Capacity constraints: when transmitting from to , the rate of the network information flow transferred to the sink node is not more than the sum of the rate of network coded packets injected to all the hyperarcs including the link .(2)Flow constraints : in a multicast session, the inputoutput relationship of the information flow rate at each node must meet the information flow balance constraints. Here the required fixed multicast rate is denoted by .(3)Scheduling constraints: in order to decompose the NC subgraph into a convex set of effective scheduling, the NC subgraph must be located in the stable set polytope CH_{SS} of a conflict graph.
Theoretical studies have shown that the optimal NC subgraph can be obtained by solving the optimization problems above. Then, the multicast transmission can meet all the constraints through RNC [2, 3]. Note that in this optimization framework we only consider the case using NC in a single multicast stream. The case of multiple multicast streams is thus similar to a single multicast stream. The case of using NC among multicast streams is left to future studies.
4. Optimization Algorithm
The global optimization in the stable set polytope CH_{SS} has been proved to be a NPhard problem [12]. Therefore, there is no efficient solution to solve it even with a centralized algorithm. We thus propose to construct a heuristic subgradient algorithm by using Lagrangian dual theory.
4.1. PrimalDual Subgradient Algorithm
Now let denote Lagrangian multipliers. Using Lagrangian relaxation, we can construct a dual function of the objective function of the original problem by moving the inequality constraint (7) into the objective function; that is,where and is defined as a flow constraint in (8). Here represents the corresponding Lagrangian function with , which is derived asIn (11), note that the last line can be decomposed into two parts that are coupled with Lagrangian multiplier . Substituting (11) into (10), we then can decompose this problem into two subproblems; that is, Then, the dual problem of the original problem becomesTo solve this dual problem, we use the projected subgradient algorithm with the iterative rulewhere is a suitable step size and is a subgradient at step , which is computed byIn (15), and are the solutions of subproblems 1 and 2, respectively, at step .
Note that the computation of each subgradient is only related to those local variables, that is, and . Therefore, all the computations can be carried out within a singlehop for each node. It indicates that a node based algorithm can be carried out in a distributed way without collecting the entire information of the network. That is, for updating Lagrangian multiplier by (14) and (15), each node only needs to communicate with its neighbor nodes rather than other nodes within two hops or more.
However, the intermediate solutions and might not be optimal. We thus use a primary recovery technique to enable nearoptimal results [13]. For simplifying the description, the optimal results at step for each subproblem are denoted byThe core idea of the primary recovery technique is to compute the convex set of all historical results as the current optimal solution of the main problem. In detail, the optimal solution of the main problem in step can be obtained through the following rules:whereFinally, the primaldual subgradient iterative optimization algorithm is summarized in Algorithm 1.

Larsson et al. have proved that the iterative algorithm can guarantee to converge to the global optimal solutions of the main problem if setting (where and ) and [13]. Other options also can be found in [3]. Considering the stopping criterion of the algorithm, it can be designed according to the gap between the optimal results of the main problem generated at each iteration and those at the last iteration. For example, the stopping criterion can be designed as , where is the accuracy requirement.
4.2. Solution of the Subproblems
As shown in Algorithm 1, the key of the algorithm lies in solving the two subproblems. First, let us focus on subproblem 1 in (12); that is,Note that is nonnegative. We thus can consider as an equivalent length of the corresponding communication link . Therefore, this subproblem can be treated as a classic shortest path searching problem: searching for the shortest path from the source node to each sink node , in which the length of the path is measured by . Obviously, this problem can be solved by the traditional asynchronous distributed BellmanFord algorithm [14]. Therefore the solution of subproblem 1 can be distributed.
Now let us focus on subproblem 2 in (12); that is,Then, we transform it into a maximum problem:Additionally, we define the weight of each hyperarc corresponding to each node asNote that both and are local variables which can be collected and computed within onehop of node . Therefore, can be obtained via a node based distributed algorithm. Because belongs to a stable set polytope CH_{SS} and can also be viewed as an equivalent weight of , this subproblem is a typical maximum weighted stable set problem in the scheduling optimization [15]. The difficulty of this subproblem lies in searching for the entire CH_{SS}, which results in a NPhard problem [12]. This means that, for both centralized algorithms and distributed algorithms, the bottleneck is the solution of the scheduling subproblem.
In order to avoid the NPhard complexity caused by searching for the entire CH_{SS}, a greedy algorithm named GWMIN was introduced in [10], in which the core idea of updating weights heuristically comes from [16]. The key steps of GWMIN are to update the weights in a heuristic way and search for a maximum stable set by certain rules, instead of searching for the entire CH_{SS}, at each iteration. However, the GWMIN has poor accuracy and robust performances. To overcome these drawbacks, this paper proposes a random sampling CH_{SS} based algorithm, given a conflict graph, sampling maximum stable sets randomly at each iteration, and solving subproblem 2 with random maximum stable sets. Now we use to represent the incidence vector of the maximum stable set obtained at the th random sampling. Then we can achieve the approximate solution of (21) by solving the following optimization problem:
Apparently, the optimization problem (23) is a linear programming problem so that an efficient polynomial algorithm such as the simplex method can be adopted. Therefore, using the maximum stable sets obtained by sampling CH_{SS} randomly, we can transfer the original NPhard problem into a general linear programming problem to obtain an approximate solution. Finally, given a conflict graph , the random sampling algorithm for each maximum stable set is presented in Algorithm 2.

As analyzed before, the computation of can be distributed. Furthermore, the random sampling algorithm for maximum stable sets also can be implemented in a distributed algorithm proposed in [17]. After collecting the necessary information of and maximal stable sets in a distributed way, problem (23) can be solved by a centralized algorithm. As a result, on the proposed solution for problem (12), only this part is centralized, and all of the others can be distributed.
5. Simulation Results
In order to analyze the performances of the proposed algorithm, we developed a simulation package in MATLAB. To solve problem (23), we used CVX, a package for specifying and solving convex programs [18]. For each simulation, we used a random network topology model in which all nodes were distributed within a rectangular area randomly and uniformly with a unit density. Any two nodes are viewed as reachable when the distance between them is less than a certain communication radius . We used multicast scenarios with one session and two sink nodes in simulations. We also used an auxiliary interference model to model a conflict graph. We selected the leftmost node as the source node and the rightmost two nodes as sink nodes for each multicast scenario within a random generated network topology.
5.1. Performance Comparisons
We simulate the following three algorithms: a full centralized algorithm indicating the global optimal results; the GWMIN algorithm presented in [10]; and the random maximum stable sets (RMSS) based algorithm proposed in this paper.
For a fair comparison, we used the same stopping criterion and iterative parameters for both GWMIN and RMSS. The step size is set as and the stopping criterion is set as . We simulated four different random scenarios with the numbers of nodes of 5, 10, 5, and 20, respectively. Note that the number of nodes in a conflict graph is the number of hyperarcs in a hypergraph, which is growing exponentially with the number of neighbors for each node. Accordingly, in order to reduce the number of neighbors for each node, is set to 1.8, in scenarios where the number of nodes is not more than 15; is set to 1.6 in scenarios where the number of nodes is 20. In addition, we also adopted different for simulations in scenarios with different numbers of nodes. Actually, the parameter has a great effect on the optimization performance. It will be discussed in the next section. For each case, we solved the optimization problem with different algorithms for 100 random topologies and then computed the average values of optimization results, as shown in Table 1.
The last row in Table 1 represents the normalized gain of the RMSS algorithm over the GWMIN algorithm. It indicates that RMSS outperforms GWMIN significantly. In most cases, this gain can be more than 20%. Further, we can see that the performance of RMSS is almost close to the optimal result when the number of nodes is 5, 10, and 15. The main reason is that GWMIN selects only one stable set with relative large weights at each iteration while RMSS selects stable subsets randomly at each iteration. Therefore, after a certain number of iterations, RMSS can better approximate the global optimal results. Comparing with GWMIN, the complexity of RMSS is mainly to select more maximal stable sets for the optimization. Fortunately, this is a linear programming problem, which can be solved efficiently.
Finally, note that the performance of RMSS is far from the global optimization result when the number of nodes is 20. This is because the number of hyperarcs in a conflict graph is growing exponentially with the network size. To achieve a good performance, therefore, the parameter needs to increase with the increase of the network size. The impacts of on the performance will be analyzed in the next subsection.
5.2. Effects of the Parameters
For studying the impacts of the parameter and the step size, we have chosen a multicast scenario with the number of nodes of 20 and random topologies. First, let us focus on the impacts of the parameter . For this case, the step size is still set to , the communication radius is set to 1.6, and the parameter in RMSS is set to 40, 50, and 80, respectively. The simulation results for 1000 iterations are shown in Figure 3.
From Figure 3, we can see that the performance of RMSS is always better than GWMIN. This is because the size of the maximum stable sets employed by RMSS is much larger than that employed by GWMIN, which results in the size of the search space in RMSS being also much larger than that in GWMIN. Figure 3 also shows that the performance of RMSS can be improved by increasing the value of the parameter . This is also because the size of the search space in RMSS increases with the increase of the parameter . Therefore, the larger the parameter is, the better the performance it can reach.
Additionally, Figure 3 shows that both RMSS and GWMIN have converged within 100 iterations. To compare their efficiencies with each other, we have recorded the running times of 1000 iterations for each simulation and then computed the average running time of one iteration for both RMSS and GWMIN, as concluded in Table 2. The simulations have been carried out under the following software and hardware environments: CVX (version 1.22), MATLAB (version 8.3.0), Windows 7 OS, 3.4 G CPU, and 8 G RAM.
From Table 2, we can see that the average running time of RMSS with is twice that of GWMIN due to larger search space. However, from Figure 3, we also find that by some improved methods the 80RMSS can converge within about 50 iterations, which is nearly half of that in GWMIN. It indicates that both GWMIN and 80RMSS can have similar efficiency for the convergence. Moreover, Table 2 shows that the average running time decreases with the decrease of the value of . Therefore, the parameter can provide a good way for the tradeoff between the performance and the efficiency.
Now let us focus on the influence of the step size on the performance of the algorithms. For this case, we modified the step size to and kept the other parameters similar to those in the above simulations. Then, we carried out simulations for GWMIN and RMSS with . The simulation results are shown in Figure 4.
From Figure 4, we can see that RMSS still can converge within 100 iterations, which is similar to that in Figure 3. However, Figure 4 shows that GWMIN cannot converge even within 1000 iterations. It indicates that GWMIN is much more sensitive to the step size than RMSS. Therefore, comparing with GWMIN, RMSS can adapt to the step size and achieve better robustness. It means that, to achieve fast and good convergence, GWMIN must choose the step size carefully in advance while RMSS can set it randomly within a reasonable range. Therefore, RMSS is suitable realtime services much more than GWMIN. Finally, as shown in Figure 4, the performances of GWMIN are better than RMSS before 400 iterations due to the randomness of initial values (i.e., the initial value of is set randomly in the simulation). From the trend of the convergent results, however, we can see that RMSS outperforms GWMIN significantly.
6. Conclusions
In wireless multihop networks, there is a challenge on minimizing the energy of transmissions with guaranteeing a fixed and qualified data rate for realtime multicast services. By addressing this issue, using a conflict graph model, we construct a joint network coding and scheduling optimization framework for modeling this optimization problem. Using Lagrangian decompositions, we also propose a primaldual subgradient algorithm to solve it. This method can transform a NPhard scheduling subproblem into a normal linear programming problem to obtain an approximate solution by sampling maximal stable sets randomly. Moreover, most calculations of the proposed algorithm can be carried out by a distributed algorithm, and only the solution of a linear programming problem needs to be done by a centralized algorithm. Simulation results show that the performance of this algorithm outperforms the existing algorithms significantly. The studies also show that the parameter provides a good way for reaching a tradeoff between the efficiency and the performance. Finally, compared with the existing algorithms, the proposed algorithm can converge by setting the step size of iterations randomly, which results in better robustness.
Since the parameter , the initial values, and the step size have significant impacts on performance of the proposed algorithm, how to set the most suitable parameters and update them online needs more research. Finally, we would like to point out that recent studies have shown that intersession NC can achieve better throughput performance than intrasession NC [19, 20]. Therefore, how to integrate intersession NC into the optimization framework is another worthy research direction in the future.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to thank the anonymous reviewers for valuable comments and suggestions for improving the quality of this paper. This work was supported in part by National Natural Science Foundation of China (61001068).