Abstract

For the high speed sensor networks applications such as Internet of Things, multimedia transmission, the realization of high-rate transmission under limited resources has become a problem to be solved. A high speed transmission and energy optimization model oriented to lifecycle maximization is proposed in this paper. Based on information-directed mechanism, the energy threshold set and the relay node distance selection will be done in the process of target tracing, as a result, retaining a balance between transmission rate and energy consumption. Meanwhile, multiagent coevolution is adopted to achieve the maximum of network lifecycle. Comparing with the relevant methods, indexes for network such as hops, throughput, and number of active nodes, standard deviation of remaining energy, and the network lifecycle are considered, and the simulated experiments show that the proposed method will promote the transmission rate effectively, prolong the network lifecycle, and improve network performance as a whole.

1. Introduction

The wireless sensor network (WSN) is a kind of multihop ad hoc network which is composed of a large number of tiny sensor nodes by means of wireless communications [1, 2]. The function of wireless sensor network is perception, acquisition, and processing of the information of monitored object in the coverage area, and then the information is sent to the observer. The network is widely used in military, intelligent transportation, environmental monitoring, health care, and so forth [36]. Due to the small size of sensor nodes, as a result, the energy of battery carried is limited. When sensor nodes are deployed in complex areas, how to optimize the network energy consumption and increase the network lifecycle is one of the key issues. Studies have focused on energy-efficient optimization of wireless sensor network, and some remarkable results such as minimizing the energy for routing and maximum lifecycle have been involved [710].

Nowadays facing real-time requirements of the new high speed sensor networks, such as wireless multimedia sensor networks and industrial real-time control, energy conservation cannot meet the needs of practical application completely. The high speed sensor network, as an emerging network, is to advance the traditional sensor network technology in the application of multimedia transmission, and so forth, with higher bandwidth and more computing and processing power, which has become a hotspot in recent years [11, 12]. Compared to the traditional networks, the wireless sensor network has a higher demand in aspects of throughput, real-time, computational complexity, power requirements, and so on. The most essential difference between them lies in throughput requirement which is the basis for high speed sensor network applications. Therefore how to meet the requirements of transmission performance while reducing energy consumption has become an urgent problem.

For wireless sensor networks, the issues of energy saving and lifecycle have been concerned. Literature [13] reveals the findings of lifecycle and predictive perception space in the route selection, in which the differences of node number and perception scale are used with linear programming equations to design the different scenarios to promote the lifecycle of networks by computing the normalization. But this method considers the single factor of lifecycle rather than the others. Literature [14] designs a cross-layer unified optimum convex model based on the integration of physics layer, MAC (Media Access Control) layer, and network layer, thus prolonging the network lifecycle. Literature [15] provides a dynamic tuning scheme for the transmission scope of sensor nodes against the emergence of energy hole, which adopts centralized algorithm and distributed algorithm to regulate the transmission scale, thus reducing the complexity of inquiry and promoting the lifecycle. Literature [16] predicts the next location of target-closed routing based on the information of target tracing, in which the nodes located adopt a method of inducing tracking to lead a distant node to move close to the center, thus completing local coverage positioning and tracking the energy consumption. All above studies of lifecycle for the sensor networks are based on the limited resources of sensor network, which do not consider the impact of transmission rate on the system, thus failing to meet with the application requirement of real-time high performance transmissions.

In this paper, some corresponding researches have been done under the background above; at the same time multiagent is used to optimize system. Literature [17] provides a simulation platform of multiagent oriented structure through formal modeling, on which the perception and rationalization can be made by decision-making and action prediction of agents, and the cooperation mechanism of each node can be simulated. Literature [18] analyzes the cooperation and competition relations of wireless sensor networks available for resource task and information processing, by utilizing the characters of multiagent such as distribution and being autonomous and self-organizing. The studies above only modeled and analyzed the network structure based on the multiagent properties, but not the systematic operation. Literature [19] focuses on the confliction and precision of target tracing sensors for the wireless sensor networks. Literature [20] designs a data convergence system with mobile multiagent for the wireless sensor networks, thus obtaining the effective convergence of time and energy to meet with users. The above two use the multiagent technology to realize target tracing and information convergence of networks but did not consider the transmission performance of data. Literature [21] proposed a multiagent energy balance strategy based on two different indicators in cases of wireless sensor networks hotspot being frequently inquired. Literature [22] introduced multiagent to realize the resource distribution in the course of task scheduling, on which the resource distribution and task were heavily correlative. Literature [23] introduces multiagent to locate the target and classify the task and uses the support vector machines to reduce communication load and optimized the system. The studies above analyze the transmission performance on resource distribution, energy consumption, and load, respectively, but did not consider the indicators comprehensively. All the results above are about the use of agent in wireless sensor network, they are all from the perspectives of the similarity of the network structures to formalize modeling or define unit individual, and the issue of system optimization by multiagent has not been considered.

As to the problems above, this paper proposes a high speed transmission and energy optimization model oriented to lifecycle maximization (HS-EO), which tracks targets in the way that is information-directed, and selects the relay nodes and sets the energy threshold based on the wireless transmission performance, then gives a link selection based on multiagent coevolution algorithm, and finally realizes the maximum lifecycle of networks while high performance transmission is achieved.

The structure of this paper is organized as follows: Section 2 refers to the problem presentation; Section 3 presents the principles and schemes of HS-EO; Section 4 proposes a multiagent coevolution algorithm based on HS-EO; Section 5 carries on the correlative experimentations and analyses; Conclusions go last.

2. Problem Descriptions

HS-EO model is based on the information-directed mechanism and aiming at achieving high throughput and energy optimized transmission in the process of target tracing. Information-directed way can help to trace the target in wireless sensor network without knowing the location of the target node in advance; only some relevant information of the target node is needed. In this process, it is needed to obtain and update information with its neighbor nodes by information exchange, shown in Figure 1. denotes the start node, denotes the relay node, and denotes the target node. While selects the next hop, it needs to transmit information to the adjacent nodes, namely, , , and , and then it goes on close to . There are two conditions to be met: one is that the correlative information of target node can be perceived and obtained in wireless sensor networks; the other is that each node can have its own coordinate by GPS or anchors.

The route that is information-directed rests on routing strategies. HS-EO is designed to increase the transmission rate and reduce the energy consumption, in which the relay node selection of high performance transmission is shown in Figure 2.

For node , the radius which is able to achieve maximum throughput capacity within its communication radius can be obtained, and the nodes , , around it can provide high speed rate as well. In order to select high speed nodes while considering energy consumption, a multiagent coevolution algorithm is adopted to optimize the multiple objectives in this paper.

3. HS-EO Model Definition

In the process of target tracing, HS-EO is a high performance transmission model which is information-directed. For the problem of routing selection in data transmission, nodes are scheduled through energy threshold setting, and further the nodes with high throughput are selected.

3.1. Information-Directed

Since the information-directed way is suitable for the data transmission in wireless networks, it is often used to make the target location and tracking [24]. In the two-dimensional plane, the position information of the target node can be obtained by a series of measures which use sequence Bayesian filter on a set of nodes . Based on the prior information of the target node at the time t, the probability density function is used to obtain the estimates at time . Among them, the sequence Bayesian filter can measure and update through the following formula:

To quantify the contribution of each sensor node, the “Mutual Information (MI)” can be used, which can define the information contribution to the evaluation of relevant nodes. Using the random variable , , and joint probability density function , then the MI can be defined as follows:where shows the relationship between and , which explains how much information that transfers about . In sequential Bayesian filter, the contribution of measurement for node at the time can be expressed as follows:

Formula (3) shows how much information about location information of the target node is offered by , so reflects the variation of a priori information generated by node . Thus, if , this means node can continue to approach to the direction of target node.

3.2. Selection of Relay Nodes
3.2.1. Communication Distance Analysis

In multihop wireless sensor networks, the data are often forwarded by the adjacent nodes, so the relationship between throughput and transmission distance needs to be analyzed [25]. Assuming that the transmission power of the node is fixed, a half-duplex communication mode is adopted, and the location of the nodes in the region obeys the random Poisson distribution. When a node begins transmission, if the node coordinate is taken as the origin, during the transmission, the position of other nodes can be denoted as a set of location coordinates , which is defined as follows: where represents the time period of transmission, represents the nodes transmitted successfully in , and represents the coordinate origin.

Under the condition of , the probability of successful transmission from the origin to a node can be expressed as follows:

According to the definition of the transmission of wireless networks, the intensity of throughput can be defined as follows:where represents the effective traversal distance to the destination node in single hop.

There are two factors that need to be considered in selecting the next relay node: one is distance, and the other is throughput. Therefore, the throughput intensity needs to be calculated under the condition of distance , as shown in formula (7). Among them, represents the minimum SINR (Signal to Interference Plus Noise Ratio) value required to receive the data.

According to the derivation for formula (7), the maximum throughput and the corresponding transmission distance can be obtained, shown as follows:

Therefore, the nodes with high throughput can be further obtained.

3.2.2. Analysis of Network Energy Consumption

In wireless sensor networks, energy consumption is the key problem of data transmissions. In order to make the network energy consumption distribution more reasonable and then extend the effective work time of networks, the threshold of energy consumption of sensor nodes needs to be set, which is the foundation of node scheduling. On the basis, network lifecycle optimization can be achieved under the multiconstrained conditions.

(1) Scheduling of Node Energy Consumption. In order to avoid shortening network lifecycle due to premature exhaustion of energy of nodes, the node is not allowed to be the relay node in two cases, so two thresholds are defined accordingly.

Firstly, in all the neighbor nodes of , a node can not complete the transfer task due to the low energy, shown as the following formula: where represents the transmission energy consumption of data flow, and represents the neighbor nodes of , which are defined as follows: where represents the set of all nodes in the area, and represents the communication radius between nodes.

Secondly, in order to avoid inappropriate relay nodes chosen in transmission and then resulting excessive unnecessary energy consumption, the limitation of transmission energy consumption needs to be set. When node selects node for data forwarding, there is a difference between energy consumption of transmission between and and the expectation of that between and other nodes. We scale the difference by the standard deviation of energy consumption between and others and obtain the normalized difference. The probability that the normalized difference is greater than is greater than , as shown in formula (12), which means transmission cost is too high when is taken as the relay node, and the reason for that is often due to transmission distance or interference.where denotes the threshold set according to the scene; represents the probability that transmission energy consumption exceeds the setting threshold; represents transmission energy consumption between nodes and ; represents the transmission energy consumption between node and the other neighbor node ; indicates the expectation of the transmission energy.

As mentioned above, the neighbor nodes which meet the two conditions simultaneously may be chosen to participate in the data transmission, and the purpose of these conditions is to avoid the unexpected scheduling of nodes and shorten the lifecycle of network.

(2) Network Lifecycle. Lifecycle of wireless sensor network is defined as the survival time of the first failure node caused by energy depletion, so the network lifecycle should be analyzed based on each individual. Assume that all nodes have the same initial energy and communication radius, and a fixed transmission power is used. When the distance between nodes is , the transmission energy consumption with data can be expressed as below: where represents circuit power consumption of data transceiver and represents circuit power consumption for signal amplification during sending and receiving data.

The total amount of data forwarding between nodes can be expressed aswhere represents the data flow between the current node and its neighbor node .

Through the definition of node data flow and energy consumption, we can get the description of the survival time of node , as shown in

The lifecycle of wireless sensor network is relative to the survival time of the shortest one, so the corresponding network lifecycle of data flow is defined as

Therefore, the maximum network lifecycle can be defined as

(3) Set of Alternative Relay Nodes. According to the above definition, we can get a set of alternative relay nodes for the transmission scheme, which is defined as follows:

Condition ① requires that the relay nodes should satisfy the two thresholds of energy; condition ② requires that the relay nodes can be aware of the information change of the target node and then track the target. Among them, represents the set of nodes with high throughput which meets the energy threshold and is information-directed; represents the distance between the current node and neighbor node ; represents the area where the node around has high throughput. For convenience, assume that the area only contains no more than four nodes, which are closest to , and the number of nodes depends on the density of neighbor nodes. Shown as Figure 3, there are six nodes; namely, , , , , , and meet with the constraints, but in the valid area of , only , , , , four nodes, at most are selected as elements of the alternative relay nodes set, and then we can choose the proper nodes for routing.

4. Model Optimization Algorithms

According to the transmission requirements of HS-EO, a multiagent coevolution algorithm based on self-adjusting search space is proposed. Multiagent system is an important branch in the field of distributed artificial intelligence, which is a set composed of multiple agents. The aim of that is working together to solve the problem which can not be done by single agent in the way of collaborating with each other and environments. So an important issue of the system is to use the right collaborative strategy to complete the assigned task in the process of their interaction. Multiagent system takes each agent as a unit, and each individual has its own place; through the comparison with the neighbor agents, the specific behavior is carried out to achieve the population evolution, so as to obtain the process of finding optimal solutions.

HS-EO model is able to perceive all the design requirements to obtain information to track the target node and elect a collection of high throughput relay nodes under the condition of no abnormality of energy. On this basis, the transmission scheme of network life cycle maximization can be achieved. As a result, the optimization of transmission performance and overall energy consumption can be obtained. In fact, for wireless sensor networks, if without considering the meaningless consumption, energy consumption of nodes will occur only in transmission. Thus, the problem of network lifecycle maximization can be converted into a problem of lifecycle maximization for transmission path, defined as where represents the data flow through node , and all the nodes in the set of alternative relay nodes are composed of alternative paths in the solution space.

4.1. Definition of Agents

Agent represents a candidate solution for routing, which can be defined as the form of a multivariate vector, as shown in where represents the lifecycle of node , , , and    represents the number of relay nodes. The energy possessed by agent can be expressed as follows:

4.2. Behaviors of Agents

The key to the multiagent system is whether it is able to develop a suitable coordination strategy. The implementation of strategy is capable of generating a series of behaviors, and each agent will enhance its own energy through behaviors interacting with external. To this end, three kinds of behaviors are defined, namely, variation, competitive, and self-organized criticality behavior. First of all, agent can enhance its own energy by using knowledge, and variation is a kind of behavior using knowledge of itself; secondly, to obtain better solutions while computational cost is spent as little as possible, we can take the amount of calculation as limited environmental resources of agent, making each agent obtain more resources to increase its own energy by using competition. In addition, the interaction between agents often has nonlinear characteristics; in order to produce emergent phenomenon in the multiagent system, self-organized criticality behavior should be defined to improve the relationship between agents.

(1) Variation Behavior. The change of individual can promote the evolution of population, for the realization of collaborative evolution of multiagent, and variation behavior of agent is needed. The variation behavior will produce a new agent, is used to replace the existing , and the specific transformation is shown as follows:where denotes a random number allocated in interval of to node ; can get a random value, and the value belongs to .

(2) Competitive Behavior. Agent usually needs to compete with its neighbor nodes. Shown as Figure 4, each agent is taken as a unit, each agent will compare the energy with the surrounding neighbor nodes, and if there is no energy of neighbor nodes higher than its own, the nodes may continue to exist; otherwise, it will be replaced by other agents, and the new agent replacing the position is the individual derived by the agent with the maximum energy in the region.

Assuming that the current agent is , the agent with maximum energy is expressed as . Energy comparison among agents is needed. If , then agent does not change and continue to exist. Otherwise, will be replaced by offspring which is produced by where represents a random value that belongs to , is a preset value used to represent the degree of difference between vector and , and usually the value of is small, so that most of the information of will be retained.

(3) Self-Organized Criticality Behavior. Self-organized behavior refers to the formation of the state without outside intervention, produced entirely by the relationship of the occurrence of the various subsystems within the system. Criticality, as a kind of state of the system, denotes that small local changes will also affect the global situation, and all subsystem behaviors in the critical state are influencing and associated. The occurrence of the state of self-organized criticality is in accordance with the power-law distribution [26], and the definition of the behavior aims to control the individual variation and selection in population, so as to avoid stagnation and make the system update continuously.

Self-organized criticality can be described as follows: represents the agent with the maximum energy for tth generation; represents power-law distribution of tth generation, of which the corresponding numerical interval is [0.2, 0.8]. The significance of this interval is that the population is in standstill when the optimal agent has no change in a number of iterations.

The steps are as follows.

Step 1. The Agent is evolved by using variation and competitive behavior.

Step 2. If the value of has no change in consecutive times, then go to Step 3, or Step 1 is performed.

Step 3. For each agent , if , then a new individual will be generated to replace agent using formula (23).

4.3. Self-Adjustment of Search Space

With the number of nodes in the area increasing, leading to the expansion of the solution space, the optimal solution will be solved in larger search space, and the search speed will be reduced. In order to improve the efficiency of the algorithm, the search space needs to be converged properly. Firstly, determine the set of alternate relay nodes, in descending order according to the lifecycle. Then the nodes mentioned before will be selected, generally, the selected probability of which is relative high. According to the operation condition, the value of is extended dynamically, eventually extending to the whole set of alternative nodes. The step of self-adjustment for the search space is as follows.

Step 1. Initialize ; commonly the value of is small.

Step 2. Iteration is performed for the agent algorithm.

Step 3. If the optimal result is not updated in consecutive times, then go to Step 4, or Step 2 is performed.

Step 4. Let ; then go to Step 2.

Searching for the whole solution space directly is avoided by using the method. As a result, the optimal solution will be found quickly with high probability when the solution space is relatively large, and usually the probability of achieving the optimum can meet the requirements when .

4.4. Solution of Routing with Optimal Lifecycle

The specific steps of the multiagent coevolution algorithm are shown below.

Step 1. Initialize agent ; let and . Initialize .

Step 2. Competitive behavior is performed for each agent , the result of which is .

Step 3. If , variation behavior will be performed for each , and the corresponding result is .

Step 4. will be selected from to update .

Step 5. If , let ; then update , or let , , and go to Step 8.

Step 6. Self-organized criticality behavior is performed for each , and the result is .

Step 7. is obtained from ; then update .

Step 8. If , let , or let , .

Step 9. will be output and the algorithm ends when the end conditions are satisfied, or let , and go to Step 2.

represents the optimal agent update to the tth generation; represents the optimal agent in the tth generation; represents the dynamic adjustment parameter; represents the times that the optimal agent keeps the same; represents the probability of performing variation behavior.

5. Simulations and Analysis

In order to verify the validity of the HS-EO model, a series of simulation experiments are performed on the OPNET14.5 platform. The Shortest Path First (SPF) is a classical routing method, which is simple and performs with high efficiency, and is often used in the network transmission performance testing; Minimum Energy Consumption (MEC) is the basic method of energy consumption analysis for wireless sensor network. Therefore, SPF and MEC are adopted to compare with HS-EO proposed in this paper. In addition, a variety of scenarios are set up for simulation, and indexes of hops, throughput, number of active nodes, standard deviation of residual energy, and network lifecycle are used here. The experimental parameters are set as shown in Table 1.

5.1. Hops

The hop is a basic index of routing, which affects the performance of network transmission and energy consumption and so on. 100, 200, and 300 nodes are deployed in the testing area, respectively; 30 times tests for each condition are performed. As shown in Figure 5, the value of each selected node (except for the first node) represents the average of the recent five tests. In three cases, with the number of sensor nodes in unit area increasing, hops of all methods are increased. Among them, the hops of SPF are the least. The reason of that is that the shortest path is considered, so under the same condition SPF will select path with less hops; MEC has more hops, which is due to the energy consumption considered in wireless transmission, and more hops are performed to reduce the transmission cost between point and point. In contrast, HS-EO is relatively complex, by formula (8), the node with appropriate distance is selected to communicate in order to improve the transmission performance, and so the test results show that it is between the other two methods. In particular, when the density of sensor nodes is relatively large, HS-EO has some fluctuations on hops; one of the reasons is that the interference of the node number raise makes change, and on the other hand, it is the result of considering multiple nodes scheduling to achieve the energy load balancing.

5.2. Throughput

Throughput is an important index to affect the performance of wireless sensor networks. The throughput is the effective data traffic (unit: Kbits/sec) in the case of no packet loss. Test traffic for the link is the saturation flow, and the measured values can be effective to show the performance of each method; the results of the experiment are shown as Figure 6.

When the density of nodes is low, the throughput of the three methods is not large because the selected path is less. With the number of nodes increasing in the area, throughput of each method grows in different degree. Under the condition of same power in wireless communication, the transmission rate of long distance nodes is lower than the nodes with short distance. For SPF, although the hops reduce, the transmission distance increases accordingly, so that the test results should be lower than MEC. Compared with the previous two methods, the multiagent coevolution is adopted in HS-EO in process of tracing in order to achieve the multiobjective optimization; that is to say, to obtain load balancing between improving throughput and energy consumption of the node, nodes within the solution space are optimized for throughput, so the transmission rate is improved obviously. From the experimental data we can know that the throughput of HS-EO is consistent with its design objectives. At the same time, three kinds of methods decline with different degree after peaks in the interval (200, 250); this is because when the node density is too large, the interference between nodes becomes the main influencing factor of transmission, so the throughput will decrease slightly, which is related to the characteristics of wireless transmission.

5.3. Number of Activate Nodes

Active nodes refer to the nodes that are used in process of data transmission in the network lifecycle, and the number of active nodes is directly related to the performance of scheduling. If few parts of the network nodes are used with high frequency, then the energy of nodes will be depleted quickly, so as to make network enter a period in which energy fades fast. So how to make nodes schedule efficiently and uniformly is very important for improving network lifecycle and transmission performance, shown as Figure 7.

The number of nodes activated in SPF is only relevant to the number of relay nodes in routing process, so the number is less. Because of considering the problem of saving energy consumption for transmission end to end, MEC is capable of selecting relay node dynamically according to the energy of neighbor nodes, so there will be more nodes involved in the transmission than SPF. For HS-EO, with the increase of node density, as well as the energy threshold setting, more nodes will be chosen in scheduling due to multiagent coevolution participating in routing, so the experimental result of HS-EO is higher than that of the other two methods.

5.4. Standard Deviation of Residual Energy

Standard deviation of residual energy refers to residual energy standard deviation of other nodes when the energy of some node participating in routing is depleted, which is used to evaluate the energy consumption of the algorithm and reflects the distribution of the residual energy of each node. If the residual energies of all nodes are close to each other, namely, the standard deviation is smaller, this indicates that the network energy consumption distribution is uniform; otherwise, the energy distribution is not uniform and may lead to network lifecycle shortening. Standard deviation of residual energy can be defined as where represents the total number of sensor nodes in network, represents the current residual energy of node , represents the initial energy value of nodes, and represents the average residual energy of nodes.

To test the distribution of energy consumption of nodes, 30 times tests are performed in the case of 200 nodes, as shown in Figure 8. SPF fails to schedule the network energy consumption effectively, so the test result is relatively high and large fluctuations appear. MEC takes into account the optimization of transmission energy consumption of relay nodes and adjusts according to the energy consumption of neighbor nodes, so the value of is more uniform than SPF. Among the three methods, standard deviation of residual energy in HS-EO is the lowest, which reflects that the method can make load sharing of energy consumption more uniform. The reasons are as follows: firstly, threshold scheduling can avoid the extremes of energy consumption of nodes, thereby reducing the possibility of standard deviation being too large; secondly, because the purpose of multiagent coevolution is to select the node which has better energy in the alternative set of nodes and achieve the optimization of lifecycle for link, there is no significant difference of standard deviation among the residual energies of nodes; namely, the value of that is not high; thirdly, the experimental results of number of active nodes show that the number of relay nodes involved in HS-EO is more than the other two methods; therefore the sample space is large and the measurement data are more accurate. Thereby, in the aspects of energy scheduling, HS-EO is superior to the other two methods.

5.5. Network Lifecycle

Network lifecycle is an important index to measure the effectiveness of nodes. We test 50 times in the case of 200 nodes, and the experimental results are shown in Figure 9. For SPF, because the relay nodes are relatively small, the distribution of energy is not uniform and fails to optimize the communication cost as well, so the lifecycle of network is relatively short. MEC adopts the strategy with the low energy consumption; at the same time, the experimental results of standard deviation of residual energy show that the distribution of energy consumption is more uniform than SPF, and therefore the network lifecycle is relatively long. For HS-EO, although high throughput transmission will bring the energy loss partly, multiple nodes can be scheduled with energy thresholds used in multiagent coevolution, so as to balance the energy consumption and prolong the network lifecycle timely.

6. Conclusions

In wireless sensor network, the research of target tracing always focuses on object recognition, positioning, tracing, and other issues, whereas there are no further studies on routing methods. In view of the situation, according to the current needs of energy saving and high speed transmission for wireless sensor network, a high speed transmission and energy optimization model HS-EO oriented to lifecycle maximization is proposed in the paper. Meanwhile multiagent coevolution is adopted to implement optimization for maximum lifecycle. Through the comparison with SPF and MEC, in the case of five different indexes, the results show that HS-EO can improve the transmission rate effectively, prolong the network lifecycle, and enhance the performance of the network as a whole. In further studies, we will combine network function virtualization to do more related researches.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (61403109) and the Scientific Research Fund of Heilongjiang Provincial Education Department (12541169).