<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Neurorobot.</journal-id>
<journal-title>Frontiers in Neurorobotics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Neurorobot.</abbrev-journal-title>
<issn pub-type="epub">1662-5218</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnbot.2022.840332</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>LEACH Protocol Optimization Based on Weighting Strategy and the Improved Ant Colony Algorithm</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Cheng</surname> <given-names>Xuezhen</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1460711/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Xu</surname> <given-names>Chuannuo</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Liu</surname> <given-names>Xiaoqing</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Li</surname> <given-names>Jiming</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1529881/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Zhang</surname> <given-names>Junming</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>College of Electrical Engineering and Automation, Shandong University of Science and Technology</institution>, <addr-line>Qingdao</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Shandong Senter Electronic Co</institution>, <addr-line>Zibo</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>College of Energy and Mining Engineering, Shandong University of Science and Technology</institution>, <addr-line>Qingdao</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Yimin Zhou, Shenzhen Institutes of Advanced Technology (CAS), China</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Avishek Nag, University College Dublin, Ireland; Alireza Mousavi, Brunel University London, United Kingdom</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Junming Zhang <email>truegis&#x00040;163.com</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>03</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>16</volume>
<elocation-id>840332</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>12</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>14</day>
<month>02</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2022 Cheng, Xu, Liu, Li and Zhang.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Cheng, Xu, Liu, Li and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license></permissions>
<abstract>
<p>This article aims to address problems in the current clustering process of low-energy adaptive clustering hierarchy (LEACH) in the wireless sensor networks, such as strong randomness and local optimum in the path optimization. This article proposes an optimal combined weighting (OCW) and improved ant colony optimization (IACO) algorithm for the LEACH protocol optimization. First, cluster head nodes are updated <italic>via</italic> a dynamic replacement mechanism of the whole network cluster head nodes to reduce the network energy consumption. In order to improve the quality of the selected cluster head nodes, this article proposes the OCW method to dynamically change the weight according to the importance of the cluster head node in different regions, in accordance with the three impact factors of the node residual energy, density, and distance between the node and the sink node in different regions. Second, the network is partitioned and the transmission path among the clusters can be optimized by the transfer probability in IACO with combined local and global pheromone update mechanism. The efficacy of the proposed LEACH protocol optimization method has been verified with MATLAB simulation experiments.</p></abstract>
<kwd-group>
<kwd>optimal combination weighting</kwd>
<kwd>improved ant colony optimization</kwd>
<kwd>path superiority</kwd>
<kwd>LEACH optimization</kwd>
<kwd>routing protocol</kwd>
</kwd-group>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content></contract-sponsor>
<counts>
<fig-count count="12"/>
<table-count count="5"/>
<equation-count count="31"/>
<ref-count count="29"/>
<page-count count="13"/>
<word-count count="6528"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>It is known that wireless sensor networks (WSNs) are composed of many spatially distributed sensor nodes with limited energy (Efe et al., <xref ref-type="bibr" rid="B4">2013</xref>), whereas the sensor nodes are usually powered by light batteries. Frequent charging or battery replacement of the sensor nodes would cause inconvenience to maintenance; hence, balancing the energy consumption of sensor nodes and prolonging the network lifetime are the two most important indicators to evaluate the performance of the WSNs (Tripathi et al., <xref ref-type="bibr" rid="B24">2013</xref>; Mukherjee et al., <xref ref-type="bibr" rid="B18">2018</xref>). The energy loss of the network can directly affect the performance and life of the network, which should be delicately dealt with to keep the low energy loss of the network in the communication process (Yan et al., <xref ref-type="bibr" rid="B26">2018</xref>; Mohar et al., <xref ref-type="bibr" rid="B17">2020</xref>; Lv et al., <xref ref-type="bibr" rid="B14">2021</xref>). Many routing protocols can be used in WSNs, where hierarchical routing protocols are the most widely adopted. The typical hierarchical routing protocols include low-energy adaptive clustering hierarchy (LEACH), power-efficient gathering in sensor information systems (PEGASIS), threshold-sensitive energy-efficient sensor network (TEEN) protocol, and hybrid energy-efficient distributed (HEED) clustering approach, which can gather the nodes into clusters to form a specific hierarchy (Galkin, <xref ref-type="bibr" rid="B6">2018</xref>). Particularly, the typical LEACH protocol was proposed by Heinzelman (<xref ref-type="bibr" rid="B7">2000</xref>) and adopted the &#x0201C;wheel&#x0201D; cycle mode for the first time, which is widely applied due to its low power consumption, node equality, and self-clustering adaptation. In such a protocol, each sensor node contains a clustering algorithm and a data transmission algorithm among nodes (Chen et al., <xref ref-type="bibr" rid="B1">2015</xref>), where the clustering algorithm can randomly change the cluster head nodes by comparing the size of the random number and threshold. After deployment, all sensor nodes self-organize to form different clusters. Generally, each cluster contains a cluster head node and multiple sensor member nodes. The nodes of the specific cluster can only hop one step inside the cluster to the cluster head. The cluster head node of each cluster will also hop one step to transmit the received information to the sink node to complete an iteration circle. The structure of the LEACH protocol is illustrated in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Network structure diagram of LEACH protocol.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0001.tif"/>
</fig>
<p>In the LEACH protocol, the energy of the sensor nodes is evenly distributed, and the network duration can be prolonged by balancing the energy consumption of the sensor nodes (Wang et al., <xref ref-type="bibr" rid="B25">2018</xref>). However, the random selection of the cluster head nodes in the protocol could lead to uneven distribution of the cluster head nodes, which is prone to a large or small number of member nodes in the cluster (Marappan and Rodrigues, <xref ref-type="bibr" rid="B16">2016</xref>). Each cluster head node and the base station adopt a single-hop transmission mode, which could easily cause premature death of the remote cluster head, poor expansibility, and uneven energy consumption phenomena in the network. The replacement mechanism of clusters using wheels will speed up the replacement frequency of the cluster head, which could result in more times of the information transmission in the network and shortened network life (Rahman et al., <xref ref-type="bibr" rid="B21">2013</xref>; Jameii and Maadani, <xref ref-type="bibr" rid="B10">2016</xref>). Then, the LEACH protocol has been improved with the fractional lion (FLION) algorithm to generate the optimal route (Sirdeshpande and Udupi, <xref ref-type="bibr" rid="B22">2017</xref>), whereas the fractional derivative is introduced to detect the neighbor solution, and the forward link algorithm is also used to select channels to improve the network survival duration. However, the distribution density of the nodes is not considered in the fitness function. Lalwani et al. (<xref ref-type="bibr" rid="B11">2018</xref>) used the harmony search algorithm (HSA) to determine the optimal routing, where the fitness function is designed considering the node density, energy, and distance factor, by choosing the smallest distance nodes for data transmission so as to achieve the energy consumption reduction of the nodes. However, when the cluster contains non-local nodes, the network performance will be affected. Ning et al. (<xref ref-type="bibr" rid="B20">2017</xref>), Ezhilarasi (<xref ref-type="bibr" rid="B5">2019</xref>) adopted an improved particle swarm optimization (PSO) to optimize the clusters of the WSN process. However, the selection rules of the cluster head node are relatively complex, and the cluster scale rapid update frequency consumes high energy due to fast convergence. At the same time, local optimum would easily occur. Hence, certain algorithms are proposed to tackle such problems, where the cuckoo search algorithm (CSA) (Huang and Hua, <xref ref-type="bibr" rid="B9">2020</xref>) and the fruit fly algorithm (FFA) (Dai et al., <xref ref-type="bibr" rid="B2">2020</xref>) are developed to optimize the cluster and routing protocol and to determine the transmission path according to the distance, energy, and trust value from the node to the base station such factors. However, the generalization capability of the surviving nodes during communication has not been considered (Maheshwari et al., <xref ref-type="bibr" rid="B15">2021</xref>).</p>
<p>In recent years, some improved heuristic intelligent algorithms have become research hotspots. For instance, the ant colony optimization (ACO) algorithm has been extensively studied and applied since it can assist to select the optimal path for the fused data transmission and the network life cycle extension (Ding, <xref ref-type="bibr" rid="B3">2020</xref>). Experts and scholars have explored the ACO algorithms to reduce the energy consumption of the nodes. An ACO multipath routing protocol has been proposed based on the angle factor and entropy to optimize the cluster head nodes with high reliability in the process of data transmission (Hou et al., <xref ref-type="bibr" rid="B8">2017</xref>; Liu and Li, <xref ref-type="bibr" rid="B13">2018</xref>). To tackle the energy consumption balance problem, Zou and Qian (<xref ref-type="bibr" rid="B29">2019</xref>) proposed improved ACO (IACO) with a sensor node transfer function and pheromone updating routing rules for optimal WSN routing but with little consideration of the impact factors. Nayyar and Singh (<xref ref-type="bibr" rid="B19">2020</xref>) proposed an energy-saving multipath routing protocol based on the ACO algorithm for dynamic WSNs, in which the optimal path for the adjacent nodes is designed, but the scalability is poor. Therefore, it is of great significance to investigate the method to improve the transmission performance of the WSNs and select the optimal data transmission path quickly and effectively.</p>
<p>This article aims to propose an IACO algorithm for the LEACH protocol optimization. The specific sections are as follows: the first section is the improvement of the clustering algorithm based on optimal combination weighting. Furthermore, the designed transmission path between clusters based on the IACO algorithm is proposed in section Improved Clustering Algorithm Based on OCW. Simulation experiments for the proposed LEACH protocol optimization verification are described in section Energy Consumption Analysis of the Sensor Node. The Conclusion is given in section Optimization of the Transmission Path Between Clusters Based on IACO. <xref ref-type="table" rid="T1">Table 1</xref> compares the energy-saving, cluster establishment time, uniform distribution of cluster head nodes, and path selection of the representative typical routing protocols.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>WSNs cluster routing protocol performance comparison.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Protocol</bold></th>
<th valign="top" align="left"><bold>Energy saving performance</bold></th>
<th valign="top" align="left"><bold>Cluster creation time</bold></th>
<th valign="top" align="left"><bold>Uniform distribution of cluster head nodes</bold></th>
<th valign="top" align="left"><bold>Path selection</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">LEACH</td>
<td valign="top" align="left">Poor</td>
<td valign="top" align="left">Quick</td>
<td valign="top" align="left">Difference</td>
<td valign="top" align="left">Single jump</td>
</tr>
<tr>
<td valign="top" align="left">TEEN</td>
<td valign="top" align="left">Good</td>
<td valign="top" align="left">Faster</td>
<td valign="top" align="left">Difference</td>
<td valign="top" align="left">Single jump</td>
</tr>
<tr>
<td valign="top" align="left">PEGASIS</td>
<td valign="top" align="left">Better</td>
<td valign="top" align="left">Slow</td>
<td valign="top" align="left">&#x02013;</td>
<td valign="top" align="left">More jumps</td>
</tr>
<tr>
<td valign="top" align="left">HEED</td>
<td valign="top" align="left">Better</td>
<td valign="top" align="left">Slower</td>
<td valign="top" align="left">Good</td>
<td valign="top" align="left">More jumps</td>
</tr>
<tr>
<td valign="top" align="left">EEUC</td>
<td valign="top" align="left">Good</td>
<td valign="top" align="left">Medium</td>
<td valign="top" align="left">Medium</td>
<td valign="top" align="left">More jumps</td>
</tr>
</tbody>
</table>
</table-wrap></sec>
<sec id="s2">
<title>Improved Clustering Algorithm Based on OCW</title>
<sec>
<title>The Update Mechanism of the Cluster Head Node</title>
<p>In the LEACH, cluster head nodes are replaced at each round, which would increase the network power consumption. Therefore, the cluster head nodes are updated within the existed clusters for the whole network clusters. The node whose energy is lower than the average energy <italic>E</italic><sub><italic>a</italic></sub> in the cluster loses the qualification to be the cluster head node, whereas the node with the highest remaining energy in the cluster is defined as the cluster head node of the next round, and the average energy <italic>E</italic><sub><italic>a</italic></sub> in the cluster is defined as,</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mo>&#x02211;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>i</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>E</italic><sub><italic>i</italic></sub> is the remaining energy of the surviving node <italic>i</italic> in the cluster and <italic>N</italic><sub><italic>alive</italic></sub> is the total number of the surviving nodes. If the energy of the current cluster head nodes is lower than this value, WSNs update the cluster size.</p>
<p>The average energy <italic>E</italic><sub><italic>all</italic></sub> of all the cluster head nodes in the whole network is calculated as,</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M2"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi></mml:mstyle><mml:mi>l</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mo>&#x02211;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>r</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>B</mml:mi></mml:mstyle></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>E</italic><sub><italic>r</italic></sub> and <italic>B</italic> are the residual energy and the total number of the cluster head nodes.</p></sec>
<sec>
<title>The Threshold of the Cluster Head Node Selection</title>
<sec>
<title>Improved Threshold for Cluster Head Node Selection</title>
<p>The improved LEACH protocol has the same rules as cluster head selection in LEACH protocol, while the selection threshold is defined as <italic>T</italic><sub><italic>imp</italic></sub>(<italic>n</italic>), written as,</p>
<disp-formula id="E3"><label>(3)</label><mml:math id="M3"><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>T</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>n</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>&#x003C9;</mml:mi><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:msub><mml:mfrac><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold"><mml:mn>0</mml:mn></mml:mstyle></mml:msub></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003C9;</mml:mi><mml:mstyle mathvariant="bold"><mml:mn>2</mml:mn></mml:mstyle></mml:msub><mml:mfrac><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>v</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>v</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x003C9;</mml:mi><mml:mstyle mathvariant="bold"><mml:mn>3</mml:mn></mml:mstyle></mml:msub><mml:mfrac><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mrow><mml:mi>max</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mrow><mml:mi>max</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02212;</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mrow><mml:mi>min</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;</mml:mtext><mml:mo>,</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle><mml:mo>&#x02208;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>G</mml:mi></mml:mstyle></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>0</mml:mn></mml:mstyle><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mo>,</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle><mml:mo>&#x02209;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>G</mml:mi></mml:mstyle></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:mrow></mml:math></disp-formula>
<disp-formula id="E4"><mml:math id="M4"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold-italic"><mml:mi>R</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>S</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mi>&#x003C0;</mml:mi><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi><mml:mi>p</mml:mi></mml:mstyle></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>i</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>S</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>|</mml:mo><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mstyle mathvariant="bold"><mml:mo>&#x02264;</mml:mo></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>R</mml:mi></mml:mstyle><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>S</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>&#x02208;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>S</mml:mi></mml:mstyle></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>E</italic><sub>0</sub> is the initial energy of the survived node, <italic>l</italic><sub><italic>alive</italic></sub> is the number of the existing neighboring nodes, <italic>d</italic><sub>max</sub> is the distance from node <italic>i</italic> to the sink node, <italic>d</italic><sub>max</sub> and <italic>d</italic><sub>min</sub> are the farthest and closest distance from the surviving node to the sink node, respectively, and {&#x003C9;<sub>1</sub>, &#x003C9;<sub>2</sub>, &#x003C9;<sub>3</sub>} are the weights of the node energy consumption, density, and distance, respectively.</p>
<p>In the improved threshold, the probability of being selected as the cluster head node is related to the energy, density, and distance from the sink node. As the remaining energy of the node decreases and the smaller the energy ratio is, the lower the probability of the node being selected as the cluster head node. If node <italic>i</italic> is the closest to the sink node, the distance factor is 1, otherwise, the distance factor is 0. Therefore, the range of the distance factor is between 0 and 1, and the closer the distance to the sink node, the closer the value to be 1.</p></sec>
<sec>
<title>Threshold Weight Selection</title>
<p>To ensure the rationality for the cluster head nodes selection, the optimal combined weighting (OCW) is used to adjust the three weights of the improved threshold, which combines the analytic hierarchy process (AHP) and the entropy method to ensure the objectivity of the threshold weight.</p>
<sec>
<title>Weight of AHP</title>
<p>The AHP method can decompose the problems into different levels with an evaluation index matrix to solve the maximum eigenvalues of the matrix and the corresponding eigenvector, so as to conduct a consistency test to obtain the weight of different evaluation indices (Yang et al., <xref ref-type="bibr" rid="B27">2017</xref>). However, if the AHP method is used alone, the weight cannot be reasonable due to the large subjective component and certain persuasiveness. Still, the AHP method can be applied to determine the threshold weight in WSNs, and the steps are described as follows:</p>
<list list-type="simple">
<list-item><p>(1) Construction of the structural model of the cluster head node selection. The hierarchical model of the threshold weight in Region 1 (<italic>d</italic><sub><italic>i</italic></sub> &#x02264; <italic>d</italic><sub>0</sub>) and Region 2 (<italic>d</italic><sub><italic>i</italic></sub> &#x0003E; <italic>d</italic><sub>0</sub>) are depicted in <xref ref-type="fig" rid="F2">Figures 2A,B</xref>.</p></list-item>
<list-item><p>(2) Determination of the evaluation index matrix. The residual energy, distribution density, and distance of the nodes are taken as evaluation indices, which are composed of the evaluation index matrix <italic>C</italic>, written as,
<disp-formula id="E5"><label>(4)</label><mml:math id="M5"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold-italic"><mml:mi>C</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable style="text-align:axis;" equalrows="false" columnlines="none none none none none none none none none" equalcolumns="false" class="array"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>12</mml:mn></mml:mstyle></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>13</mml:mn></mml:mstyle></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>21</mml:mn></mml:mstyle></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>23</mml:mn></mml:mstyle></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>31</mml:mn></mml:mstyle></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>32</mml:mn></mml:mstyle></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
</p>
<p>where <italic>c</italic><sub><italic>ij</italic></sub> is the measurement value of the evaluation index, and <inline-formula><mml:math id="M6"><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:math></inline-formula>. The value of the <italic>c</italic><sub><italic>ij</italic></sub> usually adopts Santy&#x00027;s 1&#x02013;9 scale method as listed in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<p>According to the importance of the impact factors in different regions, evaluation index matrices corresponding to influencing factors can be obtained from different hierarchical models, as shown in Eq. (5).
<disp-formula id="E6"><label>(5)</label><mml:math id="M7"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold-italic"><mml:mi>C</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable style="text-align:axis;" equalrows="false" columnlines="none none none none none none none none none" equalcolumns="false" class="array"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold"><mml:mn>2</mml:mn></mml:mstyle><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold"><mml:mn>5</mml:mn></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>2</mml:mn></mml:mstyle></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold"><mml:mn>3</mml:mn></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>5</mml:mn></mml:mstyle></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>3</mml:mn></mml:mstyle></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
</p></list-item>
<list-item><p>(3) The calculation of the AHP weight. According to Eq. (6), the maximum eigenvalue and eigenvector of the evaluation index matrix can be obtained, &#x003BB;<sub>max</sub> and <italic>W</italic>, respectively, where the corresponding eigenvector is <inline-formula><mml:math id="M8"><mml:mi>W</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>W</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mtext>&#x000A0;</mml:mtext><mml:msub><mml:mrow><mml:mi>W</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mtext>&#x000A0;</mml:mtext><mml:msub><mml:mrow><mml:mi>W</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>, and the obtained normalized weight is written in Eq. (7). The eigenvalue of Eq. (5) is 3.0037, and the weight of the energy, density, and distance are 0.581, 0.309, and 0.109 respectively.
<disp-formula id="E7"><label>(6)</label><mml:math id="M9"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold-italic"><mml:mi>C</mml:mi><mml:mi>W</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mo class="qopname">max</mml:mo></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E8"><label>(7)</label><mml:math id="M10"><mml:msub><mml:mi>&#x003C9;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>A</mml:mi><mml:mi>i</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:msub><mml:mi>W</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:msubsup><mml:mo>&#x02211;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold"><mml:mn>3</mml:mn></mml:mstyle></mml:msubsup><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:msub><mml:mi>W</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mstyle></mml:mrow></mml:mstyle></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;&#x000A0;</mml:mtext><mml:mo stretchy="false">(</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mn>1</mml:mn></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mn>2</mml:mn></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mn>3</mml:mn></mml:mstyle><mml:mo stretchy="false">)</mml:mo></mml:math></disp-formula>
</p></list-item>
<list-item><p>Consistency test. When the weights of the different impact factors are obtained, the final consistency test can be carried out to ensure the effectiveness of the obtained weights. The test expressions are expressed in Eqs. (8), (9), and (10). The consistency index of the weight values of the tested regions 1 and 2 is 0.0018, and the consistency ratio is 0.0032.
<disp-formula id="E9"><label>(8)</label><mml:math id="M11"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold-italic"><mml:mi>C</mml:mi><mml:mi>R</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>C</mml:mi><mml:mi>I</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>R</mml:mi><mml:mi>I</mml:mi></mml:mstyle></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E10"><label>(9)</label><mml:math id="M12"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="bold-italic"><mml:mi>C</mml:mi><mml:mi>I</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mo class="qopname">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle><mml:mo>-</mml:mo><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E11"><label>(10)</label><mml:math id="M13"><mml:mstyle mathvariant="bold-italic"><mml:mi>R</mml:mi><mml:mi>I</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mo>&#x003BB;</mml:mo><mml:mrow><mml:mi>max</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle><mml:mo>&#x02212;</mml:mo><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow></mml:mfrac></mml:math></disp-formula>
</p>
<p>where, <italic>RI</italic> is the average consistency indicator, and <italic>CR</italic> is the consistency ratio. If <italic>C</italic> meets the consistency test, <italic>CR</italic> &#x0003C; 0.1. <italic>CI</italic> is a consistency indicator and <italic>CI</italic> &#x0003D; 0 with a high degree of consistency <italic>C</italic>.</p></list-item>
</list>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Hierarchical model of threshold weight.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0002.tif"/>
</fig>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Evaluation index matrix value measurement table.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Scale values</bold></th>
<th valign="top" align="left"><bold><italic>Meaning</italic></bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">1</td>
<td valign="top" align="left">The two evaluation indicators are of the same importance</td>
</tr>
<tr>
<td valign="top" align="left">3</td>
<td valign="top" align="left">The previous evaluation index is more important</td>
</tr>
<tr>
<td valign="top" align="left">5</td>
<td valign="top" align="left">The previous evaluation index is important</td>
</tr>
<tr>
<td valign="top" align="left">7</td>
<td valign="top" align="left">The previous evaluation index is very important</td>
</tr>
<tr>
<td valign="top" align="left">9</td>
<td valign="top" align="left">The previous evaluation index is extremely important</td>
</tr>
<tr>
<td valign="top" align="left">2,4,6,8</td>
<td valign="top" align="left">The importance is the median value of the above adjacent indicators</td>
</tr>
<tr>
<td valign="top" align="left">Reciprocal</td>
<td valign="top" align="left">If the evaluation index i and j are measured as <italic>c</italic><sub><italic>ij</italic></sub>,but <inline-formula><mml:math id="M14"><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:math></inline-formula></td>
</tr>
</tbody>
</table>
</table-wrap></sec>
<sec>
<title>Weight of Entropy Method</title>
<p>The entropy method is an objective weighting method without considering the correlation among factors. Here, the entropy method is applied to determine the threshold weight in WSNs.</p>
<p>Assuming WSNs containing <italic>n</italic> (<italic>n</italic> = 1,2,&#x02026;, <italic>N</italic>) sensor nodes, the impact factor of the <italic>n</italic>th node is <italic>f</italic><sub><italic>nt</italic></sub> <italic>t</italic> (<italic>t</italic> = 1,2,&#x02026;, <italic>T</italic>), to affect the threshold of the cluster head node selection.</p>
<p>(1) Analysis of the impact factors. According to Eq. (11), the impact factors are determined as the energy <italic>f</italic><sub><italic>n</italic>1</sub>, distance <italic>f</italic><sub><italic>n</italic>2</sub>, and relative density <italic>f</italic><sub><italic>n</italic>3</sub> of the nth nodes,</p>
<disp-formula id="E12"><label>(11)</label><mml:math id="M15"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mn>0</mml:mn></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E13"><label>(12)</label><mml:math id="M16"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mo class="qopname">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mo class="qopname">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E14"><label>(13)</label><mml:math id="M17"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>i</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>i</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>(2) The entropy method for the weight of impact factors determination.</p>
<p>(a) The above three impact indicators are normalized as,</p>
<disp-formula id="E15"><label>(14)</label><mml:math id="M18"><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mstyle mathvariant="bold"><mml:mo>&#x02032;</mml:mo></mml:mstyle></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>&#x02212;</mml:mo><mml:mi>min</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mi>max</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x02212;</mml:mo><mml:mi>min</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula>
<p>where max(<italic>f</italic><sub><italic>nt</italic></sub>) and min(<italic>f</italic><sub><italic>nt</italic></sub>) are the maximum and minimum of the impact factors, respectively. The evaluation matrix <italic>R</italic> of the normalized impact factor is obtained as,</p>
<disp-formula id="E16"><label>(15)</label><mml:math id="M19"><mml:mstyle mathvariant="bold-italic"><mml:mi>R</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mn>11</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mn>12</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mn>1</mml:mn></mml:mstyle><mml:mn>3</mml:mn></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mn>21</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mn>22</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mn>23</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>&#x022EF;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x022EF;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x022EF;</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi><mml:mn>2</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi><mml:mn>3</mml:mn></mml:mstyle></mml:mrow><mml:mo>&#x02032;</mml:mo></mml:msubsup></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula>
<p>(b) The weight of the <italic>t</italic>th impact factor, <italic>E</italic><sub><italic>t</italic></sub> is expressed as,</p>
<disp-formula id="E17"><label>(16)</label><mml:math id="M20"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mo class="qopname">ln</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi></mml:mstyle></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi></mml:mstyle></mml:mrow></mml:munderover></mml:mstyle><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo class="qopname">ln</mml:mo><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E18"><label>(17)</label><mml:math id="M21"><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mstyle mathvariant="bold"><mml:mo>&#x02032;</mml:mo></mml:mstyle></mml:msubsup></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:msubsup><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi></mml:mstyle></mml:msubsup><mml:mrow><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>F</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mstyle mathvariant="bold"><mml:mo>&#x02032;</mml:mo></mml:mstyle></mml:msubsup></mml:mrow></mml:mstyle></mml:mrow></mml:mfrac></mml:math></disp-formula>
<p>where <italic>P</italic><sub><italic>nt</italic></sub> is the proportion of the <italic>t</italic> impact factor of the <italic>n</italic>th node under all indicators of the node.</p>
<p>(c) The entropy value (weight) of each impact factor is written as,</p>
<disp-formula id="E19"><label>(18)</label><mml:math id="M22"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn><mml:mo>-</mml:mo></mml:mstyle><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>T</mml:mi><mml:mo>-</mml:mo></mml:mstyle><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>T</mml:mi></mml:mstyle></mml:mrow></mml:msubsup></mml:mstyle><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>n</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold"><mml:mn>0</mml:mn></mml:mstyle><mml:mo>&#x0003C;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>S</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>&#x0003C;</mml:mo><mml:mn>1</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
</sec>
<sec>
<title>Optimal Combined Weighting</title>
<p>The OCW is a method where both the quantitative and qualitative analyses are used to reasonably allocate the weight, and how to allocate the weight of the AHP and entropy is the key of this method (Yang et al., <xref ref-type="bibr" rid="B27">2017</xref>). The weight calculation of the OCW method can be expressed as,</p>
<disp-formula id="E20"><label>(19)</label><mml:math id="M23"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>A</mml:mi><mml:mi>i</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mstyle mathvariant="bold"><mml:mo>&#x0002B;</mml:mo></mml:mstyle><mml:msub><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>2</mml:mn></mml:mstyle></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>&#x003C9;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003C9;<sub><italic>t</italic></sub>, &#x003C9;<sub><italic>Ai</italic></sub>, and &#x003C9;<sub><italic>St</italic></sub> are the weights of the optimal combination, AHP, and entropy method, respectively. &#x003BB;<sub>1</sub> and &#x003BB;<sub>2</sub> are the importance degrees of the AHP and entropy methods, respectively.</p></sec></sec></sec>
<sec>
<title>A Bunch of Rules</title>
<p>Once the cluster head nodes are determined, they are broadcast inside the WSNs so that the rest sensor nodes are invited to join specific clusters. If a sensor node receives multiple invitation messages within a period of time, it will determine the distance between the cluster head nodes and itself according to the strength of the received information, then make the decision to join the closer cluster and transmit the requested information to the specific cluster head node. The cluster head node receives the distributed information from each member node and decides which nodes can be joined. Then, the clusters can be established, and it enters the data transmission stage.</p></sec></sec>
<sec id="s3">
<title>Energy Consumption Analysis of the Sensor Node</title>
<p>The differences in energy consumption with different routing protocols can be used to evaluate the transceiver characteristics of the sensor nodes. Hence, a typical radio energy consumption model is usually adopted, including the energy consumption of the data transmitting circuit and the power amplifier circuit (Yang et al., <xref ref-type="bibr" rid="B27">2017</xref>), as demonstrated in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Wireless communication model of sensor nodes.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0003.tif"/>
</fig>
<p>There are mainly two energy consumption models for the WSNs, i.e., free-space models and multipath fading models, the adoption of which depends on the distance between the sender and the receiver (Li et al., <xref ref-type="bibr" rid="B12">2020</xref>) and written as,</p>
<disp-formula id="E22"><label>(20)</label><mml:math id="M25"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>T</mml:mi><mml:mi>X</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>T</mml:mi><mml:mi>X</mml:mi></mml:mstyle><mml:mo>&#x02212;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>e</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>T</mml:mi><mml:mi>X</mml:mi></mml:mstyle><mml:mo>&#x02212;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mo>&#x02212;</mml:mo><mml:msup><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mo>&#x02217;</mml:mo></mml:msup><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>e</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mtext>&#x000A0;</mml:mtext><mml:mo>+</mml:mo><mml:msup><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mo>&#x02217;</mml:mo></mml:msup><mml:msub><mml:mi>&#x003B5;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>&#x02217;</mml:mo><mml:msup><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold"><mml:mn>2</mml:mn></mml:mstyle></mml:msup><mml:mo>,</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>f</mml:mi></mml:mstyle><mml:mtext>&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mi>d</mml:mi><mml:mo>&#x02264;</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold"><mml:mn>0</mml:mn></mml:mstyle></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>&#x02212;</mml:mo><mml:msup><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mo>&#x02217;</mml:mo></mml:msup><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>e</mml:mi><mml:mi>l</mml:mi><mml:mi>e</mml:mi><mml:mi>c</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msup><mml:mstyle mathvariant="bold-italic"><mml:mi>l</mml:mi></mml:mstyle><mml:mo>&#x02217;</mml:mo></mml:msup><mml:msub><mml:mi>&#x003B5;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>m</mml:mi><mml:mi>p</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>&#x02217;</mml:mo><mml:msup><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold"><mml:mn>4</mml:mn></mml:mstyle></mml:msup><mml:mo>,</mml:mo><mml:mtext>&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>f</mml:mi></mml:mstyle><mml:mtext>&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mi>d</mml:mi><mml:mo>&#x0003E;</mml:mo><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold"><mml:mn>0</mml:mn></mml:mstyle></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>E</italic><sub><italic>TX</italic>&#x02212;<italic>elec</italic></sub>(<italic>l</italic>) is the loss of the transmitting circuit energy, <italic>E</italic><sub><italic>TX</italic>&#x02212;<italic>amp</italic></sub>(<italic>l, d</italic>) is the loss of the amplification circuit energy, and <italic>E</italic><sub><italic>elec</italic></sub> is the energy loss of the sender/receiver when performing data transmission. The calculation of the threshold of the data transmission distance is written as,</p>
<disp-formula id="E23"><label>(21)</label><mml:math id="M26"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>0</mml:mn></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B5;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>f</mml:mi><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003B5;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>m</mml:mi><mml:mi>p</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:msqrt></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003B5;<sub><italic>fs</italic></sub> and &#x003B5;<sub><italic>mp</italic></sub> represent the amplification energy consumption parameters of the free-space model and the multipath fading model, respectively. The type of the transmitter amplifier determines the values of the two parameters.</p></sec>
<sec id="s4">
<title>Optimization of the Transmission Path Between Clusters Based on IACO</title>
<sec>
<title>Flow Pattern Division</title>
<p>According to the energy consumption model, the energy consumption of the nodes takes the distance threshold <italic>d</italic><sub>0</sub> as the intermediate value, so the regions can be divided by the data transmission distance threshold <italic>d</italic><sub>0</sub> in the LEACH protocol optimization. The range (<italic>d</italic><sub>min</sub> &#x02264; <italic>d</italic><sub><italic>i</italic></sub> &#x02264; <italic>d</italic><sub>0</sub>) within the distance from the sink node <italic>d</italic><sub>0</sub> is area 1, and the range within the range of WSN nodes (<italic>d</italic><sub>0</sub> &#x02264; <italic>d</italic><sub><italic>i</italic></sub> &#x02264; <italic>d</italic><sub>max</sub>) is area 2. The partitioning diagram is illustrated in <xref ref-type="fig" rid="F4">Figure 4</xref>.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Schematic diagram of the zoning of wireless sensor network deployment node range.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0004.tif"/>
</fig>
<p>In WSNs, the node positions are usually fixed after deployment. When <italic>d</italic><sub><italic>i</italic></sub> &#x02264; <italic>d</italic><sub>0</sub>, the distance between nodes will not be regarded as the main impact factor, so the weight in the threshold will change. In contrast, when <italic>d</italic><sub>0</sub> &#x02264; <italic>d</italic><sub><italic>i</italic></sub>, the energy consumption of the node transmission and the distance between nodes increase to the fourth power. Hence, distance is the main impact factor. The weight in Eq. (19) will change with the importance of the impact factor.</p></sec>
<sec>
<title>Transmission Path Optimization Among Clusters Based on the IACO Algorithm</title>
<p>The ACO algorithm is a heuristic optimization method, which can be used to search the optimal path through individual efforts and group cooperation <italic>via</italic> accumulated pheromones on the path with positive feedback (Song and Yao, <xref ref-type="bibr" rid="B23">2017</xref>; Li et al., <xref ref-type="bibr" rid="B12">2020</xref>). In WSNs, the optimization of the transmission data path of the routing protocol has the same characteristics as the ACO algorithm to obtain the best foraging route for ants. The application of the ACO algorithm in WSNs still has some limitations, including low path efficiency and the appearance of local optimal solution. Here, the developed IACO algorithm is applied to WSNs to conduct path optimization from three aspects, improved transition probability, path superiority, and pheromone updating mechanism.</p>
<sec>
<title>Transfer Probability Improvement</title>
<p>The transition probability of the ACO algorithm considers only the path pheromone concentration and the distance between two nodes. The ant in the process of the optimal path search can also increase the probability of invalid path search and reduce the efficiency of path construction. Because the impact factors are not comprehensive, the flow of the data transmission of the nodes is larger, which could result in premature death affecting the network operation.</p>
<p>By using the ACO algorithm in WSN, the distance between two nodes is only considered in the next hop of the cluster node selection, which could cause more candidate nodes and a large amount of data transmission accompanied by more energy consumption. In this article, the transition probability is fully considered with the node energy for the ants&#x00027; node search ability enhancement, so as to speed up the convergence speed, and the developed transition probability function is formatted as,</p>
<disp-formula id="E24"><label>(22)</label><mml:math id="M27"><mml:msubsup><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle></mml:msubsup><mml:mo stretchy="false">(</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mfrac><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>&#x003C4;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>&#x003B1;</mml:mi></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>&#x003B7;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>&#x003B2;</mml:mi></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi><mml:mi>a</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>&#x003B3;</mml:mi></mml:msup></mml:mrow><mml:mrow><mml:mstyle displaystyle="true"><mml:msub><mml:mo>&#x02211;</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>s</mml:mi></mml:mstyle><mml:mo>&#x02208;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mi>e</mml:mi></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle></mml:msub></mml:mrow></mml:msub><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>&#x003C4;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>&#x003B1;</mml:mi></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mi>&#x003B7;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>&#x003B2;</mml:mi></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi><mml:mi>a</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mi>&#x003B3;</mml:mi></mml:msup></mml:mrow></mml:mstyle></mml:mrow></mml:mfrac><mml:mtext>&#x000A0;</mml:mtext><mml:mo>,</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold-italic"><mml:mi>s</mml:mi></mml:mstyle><mml:mo>&#x02208;</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mi>e</mml:mi></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle></mml:msub></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>0</mml:mn></mml:mstyle><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mo>,</mml:mo><mml:mtext>&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;&#x000A0;</mml:mtext><mml:mstyle mathvariant="bold-italic"><mml:mi>s</mml:mi></mml:mstyle><mml:menclose notation="updiagonalstrike"><mml:mo>&#x02208;</mml:mo></mml:menclose><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>w</mml:mi><mml:mi>e</mml:mi></mml:mstyle><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle></mml:msub></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:math></disp-formula>
<disp-formula id="E25"><label>(23)</label><mml:math id="M28"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi><mml:mi>a</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>a</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>allowed</italic><sub><italic>k</italic></sub> is the cluster head node set where <italic>k</italic> nodes have not been reached, <italic>E</italic><sub><italic>j</italic></sub> is the remaining energy of the next hop node, and <italic>E</italic><sub><italic>iave</italic></sub> is the average energy of the adjacent nodes of node <italic>i</italic>. &#x003C4;<sub><italic>ij</italic></sub> is the pheromone concentration of the path from node <italic>i</italic> to node <italic>j</italic>, and &#x003B7;<sub><italic>ij</italic></sub> is the heuristic function to be defined for the cluster head node</p>
<p>Based on the forwarding distance <italic>j</italic> node and the distance of the aggregation node, written as,</p>
<disp-formula id="E26"><label>(24)</label><mml:math id="M29"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003B7;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mstyle mathvariant="bold"><mml:mo>&#x0002B;</mml:mo></mml:mstyle><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>d</italic><sub><italic>ij</italic></sub> is the Euclidean distance between node <italic>i</italic> and the next hop node <italic>j</italic>, and <italic>d</italic><sub><italic>js</italic></sub> is the Euclidean distance between the next hop node <italic>j</italic> and the sink node.</p></sec>
<sec>
<title>Path Superiority Determination</title>
<p>When all ants are transferred to the sink node, each ant corresponds to a transmission path, so the path superiority degree can be used as the standard to measure the optimal path. The ant with a higher superiority degree is the optimal transmission path. In the previous path optimization process, the path with higher average energy is highly likely to be the best transmission path; however, the mean energy cannot represent the current node energy level. For instance, if certain nodes have high energy mean, but the energy difference between the actual nodes is substantially large, it will cause premature death of the nodes and failure of the transmission path. The path superiority is thus reflected from the lower hop count, higher mean energy, and uniform energy distribution. The higher the path superiority, the higher quality of the selected transmission path.</p>
<p>The variation coefficient is defined by the ratio between the SD and mean of the data, which is used to evaluate the difference degree of the data distribution. The smaller the coefficient of the variation, the more uniform the data distribution, and the smaller difference between the data (Zhu, <xref ref-type="bibr" rid="B28">2017</xref>). Therefore, the path variation coefficient <italic>P</italic><sub><italic>cv</italic></sub> is used to analyze the energy balance of the nodes.</p>
<disp-formula id="E27"><label>(25)</label><mml:math id="M30"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi><mml:mi>v</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>s</mml:mi><mml:mi>d</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>m</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msqrt><mml:mrow><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow></mml:mfrac><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mi>&#x003B4;</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>2</mml:mn></mml:mstyle></mml:mrow></mml:msup></mml:mrow></mml:msqrt></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>a</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>P</italic><sub><italic>sd</italic></sub> is the SD of the energy of all nodes in the path <italic>P</italic><sub><italic>m</italic></sub>, <italic>E</italic><sub><italic>ave</italic></sub> is the mean energy of all nodes in the path, and &#x003B4; is the number of nodes contained in the path.</p>
<p>Hence, the path superiority <italic>P</italic><sub><italic>s</italic></sub> can be written as,</p>
<disp-formula id="E28"><label>(26)</label><mml:math id="M31"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mo class="qopname">min</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>c</mml:mi><mml:mi>v</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mstyle mathvariant="bold"><mml:mo>&#x0002B;</mml:mo></mml:mstyle><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>J</mml:mi></mml:mstyle></mml:mrow></mml:mfrac></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>E</italic><sub>min</sub> is the minimum of the node energy, <italic>E</italic><sub><italic>con</italic></sub> is the sum of the node energy consumption, and <italic>J</italic> is the total hop number in the path.</p></sec>
<sec>
<title>Pheromone Updating</title>
<p>In the ACO algorithm, the pheromone concentration varies according to the length of the path, and ants plan the next route according to the pheromone concentration of different paths. Ants tend to choose paths with higher pheromone concentrations, and other paths are ignored. As a result, the path searching will fall into the local optimization, the ant will no longer search for new paths and the path searching process stops. Therefore, the IACO algorithm adds local and global pheromone updates to inter-cluster transmission path planning to avoid such local optimum problems. The developed pheromone updating mechanism is written as,</p>
<disp-formula id="E29"><label>(27)</label><mml:math id="M32"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003C4;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle><mml:mo>&#x0002B;</mml:mo><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle><mml:mo>-</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>&#x003C1;</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>&#x003C4;</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003C1;</mml:mi><mml:mo>&#x00394;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003C4;</mml:mi></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="E30"><label>(28)</label><mml:math id="M33"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mo>&#x00394;</mml:mo><mml:msub><mml:mrow><mml:mo>&#x003C4;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>N</mml:mi></mml:mstyle></mml:mrow></mml:munderover></mml:mstyle><mml:mo>&#x00394;</mml:mo><mml:msubsup><mml:mrow><mml:mo>&#x003C4;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle></mml:mrow></mml:msubsup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where &#x003C1; is the parameter regulating the pheromone volatilization speed, &#x00394;&#x003C4;<sub><italic>ij</italic></sub>(<italic>t</italic>) is the pheromone increment in the path, &#x00394;&#x003C4;<sub><italic>ij</italic></sub>(0) &#x0003D; 0 at the initial time, and <inline-formula><mml:math id="M34"><mml:mi>&#x00394;</mml:mi><mml:msubsup><mml:mrow><mml:mi>&#x003C4;</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula> is the pheromone concentration left by the <italic>k</italic>th ant in the path (i, j).</p>
<list list-type="bullet">
<list-item><p>Local pheromone update: If the ant node carries out data forwarding from the cluster head node <italic>i</italic> &#x02192; <italic>j</italic>, the pheromone concentrations of the corresponding paths should be locally updated,</p></list-item>
</list>
<disp-formula id="E31"><label>(29)</label><mml:math id="M35"><mml:mo>&#x00394;</mml:mo><mml:msubsup><mml:mi>&#x003C4;</mml:mi><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle></mml:msubsup><mml:mo stretchy="false">(</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mstyle mathvariant="bold-italic"><mml:mi>E</mml:mi></mml:mstyle><mml:mstyle mathvariant="bold-italic"><mml:mi>r</mml:mi></mml:mstyle></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi></mml:mstyle><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>d</mml:mi></mml:mstyle><mml:mo stretchy="false">(</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold-italic"><mml:mi>j</mml:mi></mml:mstyle><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula>
<p>where <italic>E</italic><sub><italic>r</italic></sub>(<italic>j</italic>) is the number of the member nodes in the next hop node.</p>
<list list-type="bullet">
<list-item><p>Global pheromone update: When all ants move to the sink node, each ant corresponds to a transmission path, and the pheromone concentration of the path is updated globally. Based on the path superiority, the global update rules is written as,</p></list-item>
</list>
<disp-formula id="E32"><label>(30)</label><mml:math id="M36"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mo>&#x00394;</mml:mo><mml:msubsup><mml:mrow><mml:mo>&#x003C4;</mml:mo></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>k</mml:mi></mml:mstyle></mml:mrow></mml:msubsup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mo stretchy="false">(</mml:mo></mml:mstyle><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>t</mml:mi></mml:mstyle></mml:mrow><mml:mstyle mathvariant="bold"><mml:mo stretchy="false">)</mml:mo></mml:mstyle></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mstyle mathvariant="bold"><mml:mn>1</mml:mn></mml:mstyle></mml:mrow><mml:mrow><mml:msub><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>P</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mstyle mathvariant="bold-italic"><mml:mi>s</mml:mi></mml:mstyle></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>The LEACH protocol is optimized with the improved clustering algorithm to cluster sensor nodes, and the influence factors of nodes in different network regions are different so as to replace the cluster head nodes dynamically. When the energy of the cluster head node reaches the limit value, the cluster is updated and the data are fused in the cluster head node. Although the IACO algorithm is used to find the optimal transmission path between clusters, it can effectively avoid the stagnation of the IACO algorithm in the local optimum. The specific process is depicted in <xref ref-type="fig" rid="F5">Figure 5</xref>.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Improved LEACH protocol flowchart.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0005.tif"/>
</fig></sec></sec></sec>
<sec id="s5">
<title>Experimental Verification and Analysis</title>
<sec>
<title>Experimental Environment Setting</title>
<p>The experiment verification is carried out in a MATLAB simulation environment, and 200 sensor nodes are deployed in a 200 &#x000D7; 200 m network area. The node deployment within the network is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, where the circle with sink node is the center, the area with the threshold distance <italic>R</italic> as the radius of the circle is area 1, and the rest is area 2. The parameters of the basic network and IACO algorithm in the experiments are listed in <xref ref-type="table" rid="T3">Tables 3</xref>, <xref ref-type="table" rid="T4">4</xref>.</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Sensor node deployment diagram.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0006.tif"/>
</fig>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>Basic parameter settings of the network.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Parameter</bold></th>
<th valign="top" align="left"><bold>Parameter value</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Network area</td>
<td valign="top" align="left">200 &#x000D7; 200 m</td>
</tr>
<tr>
<td valign="top" align="left">Number of nodes</td>
<td valign="top" align="left">200 individuals</td>
</tr>
<tr>
<td valign="top" align="left">The location of the sink node</td>
<td valign="top" align="left">(100, 100)</td>
</tr>
<tr>
<td valign="top" align="left">Node&#x00027;s initial energy <italic>E</italic><sub>0</sub></td>
<td valign="top" align="left">0.2J</td>
</tr>
<tr>
<td valign="top" align="left">Energy consumption at the receiving end <italic>E</italic><sub><italic>RX</italic></sub></td>
<td valign="top" align="left">50 nJ/bit&#x0002A;m<sup>2</sup></td>
</tr>
<tr>
<td valign="top" align="left">Energy consumption of the sender <italic>E</italic><sub><italic>TX</italic></sub></td>
<td valign="top" align="left">50 nJ/bit</td>
</tr>
<tr>
<td valign="top" align="left">Energy consumption of data fusion <italic>E</italic><sub><italic>DA</italic></sub></td>
<td valign="top" align="left">5 nJ/bit</td>
</tr>
<tr>
<td valign="top" align="left">Free space model energy consumption parameters &#x003B5;<sub><italic>fs</italic></sub></td>
<td valign="top" align="left">10 pJ/bit&#x0002A;m<sup>2</sup></td>
</tr>
<tr>
<td valign="top" align="left">Amplified energy consumption parameters of the multipath model &#x003B5;<sub><italic>mp</italic></sub></td>
<td valign="top" align="left">0.0013 pJ/bit&#x0002A;m<sup>2</sup></td>
</tr>
<tr>
<td valign="top" align="left">Control packet size</td>
<td valign="top" align="left">100 bits</td>
</tr>
<tr>
<td valign="top" align="left">Packet size</td>
<td valign="top" align="left">3,000 bits</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap position="float" id="T4">
<label>Table 4</label>
<caption><p>Basic parameter setting of IACO algorithm.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Parameter</bold></th>
<th valign="top" align="center"><bold>Parameter value</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Pheromone concentration weight factor &#x003B1;</td>
<td valign="top" align="center">1</td>
</tr>
<tr>
<td valign="top" align="left">Heuristic function weight factor &#x003B2;</td>
<td valign="top" align="center">5</td>
</tr>
<tr>
<td valign="top" align="left">Node energy consumption weight factor &#x003BB;</td>
<td valign="top" align="center">4</td>
</tr>
<tr>
<td valign="top" align="left">Pheromone Volatilization Coefficient &#x003C1;</td>
<td valign="top" align="center">0.1</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec>
<title>Performance Indicators</title>
<p>The life cycle of the network, the total energy consumption of the network, and the data received by the sink node are taken as indicators to evaluate the quality of the routing protocol. The specific analysis is as follows:</p>
<list list-type="order">
<list-item><p>Network life cycle: the duration from the normal operation of the WSN after the successful layout to the death of the last node. Three indicators were selected for evaluation, namely, the number of rounds in which the first node died (indicator 1), the number of rounds in which 10% of nodes died (indicator 2), and the number of rounds in which all nodes died (indicator 3). In most cases, the appearance of the dead nodes leads to the deterioration of the network detection quality, so the number of rounds where the first node dies is of particular importance.</p></list-item>
<list-item><p>Total network energy consumption: the total energy consumption of all nodes in the network during implementation. This index can reflect the balance degree of energy consumption. In the simulation, the size of the packet and control packet is set, and the energy consumptions of different nodes are calculated through the network operation.</p></list-item>
<list-item><p>Data received by the sink node: data received by the sink node after each round of the network operation.</p></list-item>
</list></sec>
<sec>
<title>Experimental Results and Analysis</title>
<p>The simulation results of different routing protocols are analyzed based on the above performance indicators to verify the effectiveness of the ILEACH protocol.</p>
<p><xref ref-type="fig" rid="F7">Figures 7</xref>, <xref ref-type="fig" rid="F8">8</xref> compare the number of surviving nodes and total network energy consumption of the LEACH protocol and the ILEACH protocol, respectively. To be specific, the impact of the ILEACH protocol on the network life cycle is analyzed <italic>via</italic> indicators 1 and 2, as shown in <xref ref-type="fig" rid="F8">Figure 8</xref>.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Comparison curve of the number of surviving nodes in the network.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0007.tif"/>
</fig>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>Comparison curve of total network energy consumption.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0008.tif"/>
</fig>
<p><xref ref-type="fig" rid="F9">Figure 9</xref> indicates that indicator 1 of the LEACH protocol and the ILEACH protocol appears in rounds 141 and 199, respectively, and indicator 2 appears when the network operates to rounds 596 and 986. When the LEACH protocol reaches &#x0007E;600 rounds, most nodes have no power, but the ILEACH protocol can be extended to &#x0007E;1,200 rounds. The network needs to be initialized in the early stage, which consumes energy quickly, and there are fewer surviving nodes in the later stage. However, ILEACH can effectively balance the energy consumption of nodes by considering the energy consumption of nodes in the process of clustering and data transmission, and the network energy consumption varies slowly in the later period. It is proven that the improved threshold and the replacement mechanism of the cluster head nodes can reduce the energy consumption of the nodes, and the ILEACH protocol can prolong the network life cycle much longer compared with the LEACH protocol.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p>Indicator analysis histogram.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0009.tif"/>
</fig>
<p><xref ref-type="fig" rid="F10">Figures 10</xref>, <xref ref-type="fig" rid="F11">11</xref> compare the number of the surviving nodes and total network energy consumption between the LEACH of ACO (ACO-LEACH) protocol and the LEACH of IACO (IACO-LEACH) protocol, where <xref ref-type="fig" rid="F11">Figure 11</xref> is the histogram of the index. <xref ref-type="fig" rid="F12">Figure 12</xref> illustrates that indicator 1 of the ACO-LEACH protocol and IACO-LEACH protocol occurs in rounds 199 and 329, respectively, and indicator 2 occurs in rounds 986 and 2,338. <xref ref-type="table" rid="T5">Table 5</xref> displays the number of the first node deaths, 10% node deaths, and all node deaths under different protocols. In conclusion, the developed IACO algorithm can gradually find the optimal transmission path for the WSN operation and save energy in the data transmission stage effectively.</p>
<fig id="F10" position="float">
<label>Figure 10</label>
<caption><p>Comparison curve of the number of surviving nodes in the network.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0010.tif"/>
</fig>
<fig id="F11" position="float">
<label>Figure 11</label>
<caption><p>Comparison curve of totalnetwork energy consumption.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0011.tif"/>
</fig>
<fig id="F12" position="float">
<label>Figure 12</label>
<caption><p>Indicator analysis histogram.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnbot-16-840332-g0012.tif"/>
</fig>
<table-wrap position="float" id="T5">
<label>Table 5</label>
<caption><p>Comparison table of the network life cycle of different protocols.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>Protocol</bold></th>
<th valign="top" align="left"><bold>Number of rounds in which the first node died</bold></th>
<th valign="top" align="left"><bold>Number of rounds in which 10% of the nodes died</bold></th>
<th valign="top" align="left"><bold>Number of rounds in which all nodes died</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">LEACH</td>
<td valign="top" align="left">141</td>
<td valign="top" align="left">596</td>
<td valign="top" align="left">659</td>
</tr>
<tr>
<td valign="top" align="left">ILEACH</td>
<td valign="top" align="left">199</td>
<td valign="top" align="left">986</td>
<td valign="top" align="left">1,296</td>
</tr>
<tr>
<td valign="top" align="left">ACO-ILEACH</td>
<td valign="top" align="left">329</td>
<td valign="top" align="left">2,338</td>
<td valign="top" align="left">3,487</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec></sec>
<sec sec-type="conclusions" id="s6">
<title>Conclusion</title>
<p>An algorithm based on OCW and IACO is proposed in the article to solve the problem of high energy consumption of the traditional LEACH protocol in WSNs. The ILEACH protocol adopts the cluster head node replacement mechanism to reduce the energy consumption considering the energy, density, and nodes distance for the threshold selection, which can effectively avoid the randomness of the clustering. Furthermore, the OCW is used to dynamically change the weight of the nodes according to the different impact factors of nodes in different regions. The developed IACO algorithm can optimize the transfer probability of the sensor nodes with the local update and global update strategies, which can prolong the life cycle of the network to a certain extent. The network environment is deployed through MATLAB simulation software to verify the feasibility of the ILEACH protocol.</p></sec>
<sec sec-type="data-availability" id="s7">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.</p></sec>
<sec id="s8">
<title>Author Contributions</title>
<p>XC and JZ proposed the idea in this study. CX and XL designed the experiment, performed the simulation experiments, analyzed the experiment results, and wrote the manuscript. JL corrected the manuscript. All authors contributed to the article and approved the submitted version.</p></sec>
<sec sec-type="funding-information" id="s9">
<title>Funding</title>
<p>The project (work) was supported by the National Natural Science Foundation of China Program (No. 62073198) and the Major Research Development Program of Shandong province of China (No. 2016GSF117009).</p></sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>XL is employed by Shandong Senter Electronic Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p></sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x00027;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p></sec>
</body>
<back>
<ack><p>We acknowledge the support of the College of Electrical Engineering and Automation and the College of Energy and Mining Engineering in Shandong University of Science and Technology and Shandong Senter Electronic Co., Ltd.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>S. X.</given-names></name> <name><surname>He</surname> <given-names>Z. S.</given-names></name> <name><surname>Xiong</surname> <given-names>H. L.</given-names></name> <name><surname>Liao</surname> <given-names>J. W.</given-names></name></person-group> (<year>2015</year>). <article-title>A wireless sensor signal reconstruction algorithm based on compressed sensing</article-title>. <source>Journal of computer science, Papers</source> <volume>38</volume>, <fpage>614</fpage>&#x02013;<lpage>624</lpage>. <pub-id pub-id-type="doi">10.3724/SP.J.1016.2015.00614</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dai</surname> <given-names>J. Y.</given-names></name> <name><surname>Deng</surname> <given-names>X. H.</given-names></name> <name><surname>Wang</surname> <given-names>B.</given-names></name> <name><surname>Wang</surname> <given-names>H. H.</given-names></name></person-group> (<year>2020</year>). <article-title>Clustering routing protocol for wireless sensor networks based on improved firefly optimized neural network</article-title>. <source>Journal of Beijing University of Posts and telecommunications, Papers</source> <volume>43</volume>, <fpage>131</fpage>&#x02013;<lpage>137</lpage>. <pub-id pub-id-type="doi">10.13190/j.jbupt.2019-161</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ding</surname> <given-names>H.</given-names></name></person-group> (<year>2020</year>). <article-title>Data fusion algorithm for wireless sensor networks based on ant colony optimization</article-title>. <source>Journal of Shenyang University of technology, Papers</source> <volume>42</volume>, <fpage>208</fpage>&#x02013;<lpage>212</lpage>. <pub-id pub-id-type="doi">10.7688/j.issn.1000-1646.2020.02.16</pub-id><pub-id pub-id-type="pmid">26861350</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Efe</surname> <given-names>K.</given-names></name> <name><surname>Ibrahim</surname> <given-names>K.</given-names></name> <name><surname>Cevdet</surname> <given-names>A.</given-names></name></person-group> (<year>2013</year>). <article-title>Active node determination for correlated data gathering in wireless sensor networks</article-title>. <source>Computer Networks</source> <volume>57</volume>, <fpage>1124</fpage>&#x02013;<lpage>1138</lpage>. <pub-id pub-id-type="doi">10.1016/j.comnet.2012.11.018</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ezhilarasi</surname> <given-names>M.</given-names></name> <name><surname>Krishnaveni</surname> <given-names>V.</given-names></name></person-group> (<year>2019</year>). <article-title>An evolutionary multipath energy-efficient routing protocol (emeer) for network lifetime enhancement in wireless sensor networks</article-title>. <source>Soft Computing A Fusion of Foundations Methodologies and Applications. Papers</source> <volume>23</volume>, <fpage>8367</fpage>&#x02013;<lpage>8377</lpage>. <pub-id pub-id-type="doi">10.1007/s00500-019-03928-1</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Galkin</surname> <given-names>P.</given-names></name></person-group> (<year>2018</year>). <article-title>Model of Reducing the Power Consumption for Node of Wireless Sensor Network in Embedded Control Systems. 2018 International Scientific-Practical Conference Problems of Infocommunications. Science and Technology (PIC SandT)</article-title>, pp. <fpage>252</fpage>-<lpage>256</lpage>. <pub-id pub-id-type="doi">10.1109/INFOCOMMST.2018.8631891</pub-id><pub-id pub-id-type="pmid">27295638</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Heinzelman</surname> <given-names>W. R.</given-names></name></person-group> (<year>2000</year>). <article-title>Energy-efficient communication protocol for wireless microsensor networks</article-title>. <source>Proc. IEEE Annu. Hawaii Int. Conf. Syst. Sci.</source> <volume>8</volume>:<fpage>8020</fpage>.</citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hou</surname> <given-names>M. T.</given-names></name> <name><surname>Zhao</surname> <given-names>Z. P.</given-names></name> <name><surname>Gao</surname> <given-names>M.</given-names></name> <name><surname>Zhang</surname> <given-names>N. N.</given-names></name></person-group> (<year>2017</year>). <article-title>Ant colony optimization multipath routing algorithm adopted angle factor</article-title>. <source>Computer Engineering and Applications. Papers</source> <volume>53</volume>, <fpage>107</fpage>&#x02013;<lpage>112</lpage>. <pub-id pub-id-type="doi">10.3778/j.issn.1002-8331.1604-0176</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>Y.</given-names></name> <name><surname>Hua</surname> <given-names>Y. Q.</given-names></name></person-group> (<year>2020</year>). <article-title>Routing optimization of wireless sensor networks with energy and path constraints</article-title>. <source>Journal of Xi&#x00027;an University of Electronic Science and technology, Papers</source> <volume>47</volume>, <fpage>113</fpage>&#x02013;<lpage>120</lpage>. <pub-id pub-id-type="doi">10.19665/j.issn1001-2400.2020.03.016</pub-id><pub-id pub-id-type="pmid">28368300</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jameii</surname> <given-names>S. M.</given-names></name> <name><surname>Maadani</surname> <given-names>M.</given-names></name></person-group> (<year>2016</year>). <article-title>Intelligent dynamic Connectivity Control Algorithm for cluster-based wireless sensor networks. 2016 11th International Conference for Internet Technology and Secured Transactions (ICITST)</article-title>, pp. <fpage>416</fpage>-<lpage>420</lpage>, <pub-id pub-id-type="doi">10.1109/ICITST.2016.7856744</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lalwani</surname> <given-names>P.</given-names></name> <name><surname>Das</surname> <given-names>S.</given-names></name> <name><surname>Banka</surname> <given-names>H.</given-names></name> <name><surname>Kumar</surname> <given-names>C.</given-names></name></person-group> (<year>2018</year>). <article-title>Crhs: clustering and routing in wireless sensor networks using harmony search algorithm</article-title>. <source>Neural Computing and Applications</source> <volume>30</volume>, <fpage>1</fpage>&#x02013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.1007/s00521-016-2662-4</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>X. Q.</given-names></name> <name><surname>Ma</surname> <given-names>R.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Hou</surname> <given-names>Y. Z.</given-names></name> <name><surname>Pei</surname> <given-names>Y. F.</given-names></name></person-group> (<year>2020</year>). <article-title>Improved design of Ant Colony Optimization and its application in Track Planning</article-title>. <source>Journal of Aeronautics, Papers</source> <volume>41</volume>, <fpage>213</fpage>&#x02013;<lpage>219</lpage>. <pub-id pub-id-type="doi">10.7527/S1000-6893.2020.24381</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>H. W.</given-names></name></person-group> (<year>2018</year>). <article-title>Non uniform clustering routing algorithm based on Ant Colony Optimization,&#x0201D; Journal of Huazhong University of science and Technology (Natural Science Edition), Papers 46(08)</article-title>, <fpage>50</fpage>&#x02013;<lpage>54</lpage>. <pub-id pub-id-type="doi">10.13245/j.hust.180810</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lv</surname> <given-names>L.</given-names></name> <name><surname>Peng</surname> <given-names>F.</given-names></name> <name><surname>Chen</surname> <given-names>W.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name></person-group> (<year>2021</year>). <article-title>Node Deployment of Wireless Sensor Networks Based on MOEA/P Algorithm. 2021 13th International Conference on Communication Software and Networks (ICCSN)</article-title>. pp. <fpage>173</fpage>-<lpage>178</lpage>. <pub-id pub-id-type="doi">10.1109/ICCSN52437.2021.9463630</pub-id><pub-id pub-id-type="pmid">27295638</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maheshwari</surname> <given-names>P.</given-names></name> <name><surname>Sharma</surname> <given-names>A. K.</given-names></name> <name><surname>Verma</surname> <given-names>K.</given-names></name></person-group> (<year>2021</year>). <article-title>Energy efficient cluster based routing protocol for wsn using butterfly optimization algorithm and ant colony optimization</article-title>. <source>Ad Hoc Networks</source> <volume>110</volume>, <fpage>102317</fpage>. <pub-id pub-id-type="doi">10.1016/j.adhoc.2020.102317</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marappan</surname> <given-names>P.</given-names></name> <name><surname>Rodrigues</surname> <given-names>P.</given-names></name></person-group> (<year>2016</year>). <article-title>An energy efficient routing protocol for correlated data using CLLEACH,nd in WSNs</article-title>. <source>Wireless Networks</source>. <volume>22</volume>(<issue>04</issue>): <fpage>1415</fpage>&#x02013;<lpage>1423</lpage>. <pub-id pub-id-type="doi">10.1007/s11276-015-1063-4</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mohar</surname> <given-names>S. S.</given-names></name> <name><surname>Goyal</surname> <given-names>S.</given-names></name> <name><surname>Kaur</surname> <given-names>R.</given-names></name></person-group> (<year>2020</year>). <article-title>Optimized sensor nodes deployment in wireless sensor network using bat algorithm</article-title>. <source>Wireless Personal Communications(</source>2), <fpage>1</fpage>&#x02013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1007/s11277-020-07823-z</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mukherjee</surname> <given-names>M.</given-names></name> <name><surname>Lu</surname> <given-names>Y. R.</given-names></name> <name><surname>Matam</surname> <given-names>R.</given-names></name> <name><surname>Choudhury</surname> <given-names>N.</given-names></name></person-group> (<year>2018</year>). <article-title>Energy trading in sleep scheduling for wireless rechargeable WSNs, 2018 10th International Conference on Communication Systems and Networks (COMSNETS)</article-title>. <source>pp</source>. <fpage>429</fpage>&#x02013;<lpage>431</lpage>. <pub-id pub-id-type="doi">10.1109/COMSNETS.2018.8328235</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nayyar</surname> <given-names>A.</given-names></name> <name><surname>Singh</surname> <given-names>R.</given-names></name></person-group> (<year>2020</year>). <article-title>Ieemarp- a novel energy efficient multipath routing protocol based on ant colony optimization (aco) for dynamic sensor networks</article-title>. <source>Multimedia Tools and Applications, Papers</source> <volume>79</volume>, <fpage>35221</fpage>&#x02013;<lpage>35252</lpage>. <pub-id pub-id-type="doi">10.1007/s11042-019-7627-z</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ning</surname> <given-names>W.</given-names></name> <name><surname>Yuan</surname> <given-names>Z.</given-names></name> <name><surname>Jinghao</surname> <given-names>L.</given-names></name></person-group> (<year>2017</year>). <article-title>A hierarchical clustering protocol for wireless sensor Networks based on improved particle swarm Optimization</article-title>. <source>Chinese Journal of Sensing Technology</source> <volume>30</volume>, <fpage>120</fpage>&#x02013;<lpage>125</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1004-1699.2017.01.022</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rahman</surname> <given-names>A. U.</given-names></name> <name><surname>Hasbullah</surname> <given-names>H.</given-names></name> <name><surname>Sama</surname> <given-names>N. U.</given-names></name></person-group> (<year>2013</year>). <article-title>Efficient energy utilization through optimum number of sensor node distribution in engineered corona-based (onsd-ec) wireless sensor network</article-title>. <source>Wireless Personal Communications</source> <volume>73</volume>, <fpage>1227</fpage>&#x02013;<lpage>1243</lpage>. <pub-id pub-id-type="doi">10.1007/s11277-013-1275-9</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sirdeshpande</surname> <given-names>N.</given-names></name> <name><surname>Udupi</surname> <given-names>V.</given-names></name></person-group> (<year>2017</year>). <article-title>Fractional lion optimization for cluster head-based routing protocol in wireless sensor network</article-title>. <source>Journal of the Franklin Institute</source> <volume>354</volume>, <fpage>4457</fpage>&#x02013;<lpage>4480</lpage>. <pub-id pub-id-type="doi">10.1016/j.jfranklin.2017.04.005</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>Y.</given-names></name> <name><surname>Yao</surname> <given-names>X.</given-names></name></person-group> (<year>2017</year>). <article-title>Design of Routing Protocol and Node Structure in Wireless Sensor Network Based on Improved Ant Colony Optimization Algorithm</article-title>. <source>2017 International Conference on Computer Network, Electronic and Automation (ICCNEA)</source> 2017, pp. <fpage>236</fpage>&#x02013;<lpage>240</lpage>. <pub-id pub-id-type="doi">10.1109/ICCNEA.2017.54</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tripathi</surname> <given-names>R. K.</given-names></name> <name><surname>Singh</surname> <given-names>Y. N.</given-names></name> <name><surname>Verma</surname> <given-names>N. K.</given-names></name></person-group> (<year>2013</year>). <article-title>Clustering algorithm for non-uniformly distributed nodes in wireless sensor network</article-title>. <source>Electronics Letters</source> <volume>49</volume>, <fpage>299</fpage>&#x02013;<lpage>300</lpage>. <pub-id pub-id-type="doi">10.1049/el.2012.3512</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>C.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Zhu</surname> <given-names>Y.</given-names></name> <name><surname>Han</surname> <given-names>Z.</given-names></name></person-group> (<year>2018</year>). <article-title>The Application of Data-Level Fusion Algorithm Based on Adaptive-Weighted and Support Degree in Intelligent Household Greenhouse</article-title>. In: Zhu Q., Na J., Wu X. (eds) Innovative Techniques and Applications of Modelling, Identification and Control. Lecture Notes in Electrical Engineering, vol 467. <publisher-loc>Springer, Singapore</publisher-loc>. <pub-id pub-id-type="doi">10.1007/978-981-10-7212-3_6</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yan</surname> <given-names>Z.</given-names></name> <name><surname>Goswami</surname> <given-names>P.</given-names></name> <name><surname>Mukherjee</surname> <given-names>A.</given-names></name> <name><surname>Yang</surname> <given-names>L.</given-names></name> <name><surname>Routray</surname> <given-names>S.</given-names></name> <name><surname>Palai</surname> <given-names>G.</given-names></name></person-group> (<year>2018</year>). <article-title>Low-energy pso-based node positioning in optical wireless sensor networks</article-title>. <source>Optik - International Journal for Light and Electron Optics. Papers</source> <volume>181</volume>, <fpage>378</fpage>&#x02013;<lpage>382</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijleo.2018.12.055</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>M.</given-names></name> <name><surname>Hu</surname> <given-names>W.</given-names></name> <name><surname>Wang</surname> <given-names>Y. S.</given-names></name></person-group> (<year>2017</year>). <article-title>Risk assessment of debris flow in the lower reaches of Zagunao River Based on Optimal Combination Weighting and fuzzy comprehensive evaluation method</article-title>. <source>Journal of natural science of Heilongjiang University, Papers</source> <volume>34</volume>, <fpage>271</fpage>&#x02013;<lpage>278</lpage>. <pub-id pub-id-type="doi">10.13482/j.issn1001-7011.2016.07.230</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>Z. C.</given-names></name></person-group> (<year>2017</year>). <article-title>Hybrid cross weighting method based on coefficient of variation G1 method</article-title>. <source>Statistics and decision making, Papers</source> <volume>12</volume>, <fpage>78</fpage>&#x02013;<lpage>81</lpage>. <pub-id pub-id-type="doi">10.13546/j.cnki.tjyjc.2017.12.018</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zou</surname> <given-names>Z. F.</given-names></name> <name><surname>Qian</surname> <given-names>Y.</given-names></name></person-group> (<year>2019</year>). <article-title>Wireless sensor network routing method based on improved ant colony algorithm</article-title>. <source>J Ambient Intell Human Comput</source> <volume>10</volume>, <fpage>991</fpage>&#x02013;<lpage>998</lpage>. <pub-id pub-id-type="doi">10.1007/s12652-018-0751-1</pub-id></citation></ref>
</ref-list> 
</back>
</article>