• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A Layered Energy-Efficient Multi-Node Scheduling Mechanism for Large-Scale WSN

    2024-05-25 14:42:26XueZhaoShaojunTaoHongyingTangJiangWangandBaoqingLi
    Computers Materials&Continua 2024年4期

    Xue Zhao,Shaojun Tao,Hongying Tang,Jiang Wang and Baoqing Li

    Science and Technology on Microsystem Laboratory,Shanghai Institute of Microsystem and Information Technology,Chinese Academy of Sciences,Shanghai,200050,China

    ABSTRACT In recent years,target tracking has been considered one of the most important applications of wireless sensor network(WSN).Optimizing target tracking performance and prolonging network lifetime are two equally critical objectives in this scenario.The existing mechanisms still have weaknesses in balancing the two demands.The proposed heuristic multi-node collaborative scheduling mechanism (HMNCS) comprises cluster head (CH)election,pre-selection,and task set selection mechanisms,where the latter two kinds of selections form a two-layer selection mechanism.The CH election innovatively introduces the movement trend of the target and establishes a scoring mechanism to determine the optimal CH,which can delay the CH rotation and thus reduce energy consumption.The pre-selection mechanism adaptively filters out suitable nodes as the candidate task set to apply for tracking tasks,which can reduce the application consumption and the overhead of the following task set selection.Finally,the task node selection is mathematically transformed into an optimization problem and the genetic algorithm is adopted to form a final task set in the task set selection mechanism.Simulation results show that HMNCS outperforms other compared mechanisms in the tracking accuracy and the network lifetime.

    KEYWORDS Node scheduling;pre-selection;target tracking;WSN

    1 Introduction

    1.1 Background

    Wireless sensor network(WSN)is comprised of many miniature sensor nodes that are deployed in a monitored area.These nodes form a self-organizing network through multi-hop wireless communication.In military scenarios,WSN is often applied for target tracking.However,sensor nodes are mostly deployed in harsh field environments and difficult to replace batteries.Therefore,many researchers aim to extend the lifetime of WSN from different perspectives,such as regulating the data reporting process at the zone level through a general event-thresholding model in heterogeneous WSN[1],obtaining suitable routing protocols for WSN [2],finding efficient relay selection [3] and so on.In this paper,we mainly aim to extend the network lifetime from the perspective of scheduling nodes to remove redundancy in WSN for target tracking.The dilemma is that improving tracking accuracy requires more nodes to track the target and consumes more energy[4,5].Therefore,it is important to investigate how to balance tracking accuracy and energy consumption to prolong the network lifetime while obtaining favorable tracking performance[6].

    To save energy,nodes are put in a sleep state where only the sensing module is turned on,while the computation and communication modules are turned off until a target is sensed.After being awakened by the target,the nodes enter the active state and are organized into a cluster to track the target.How to select the appropriate task nodes from the awakened nodes to perform the tracking task,obtain better tracking accuracy,and extend the network lifetime is the research motivation of this paper.

    1.2 Contributions

    The main contributions of this paper are summarized below:

    Firstly,we propose a cluster head (CH) election mechanism,which innovatively introduces the changing trend of the distance between the target and the node to obtain the optimal CH.This reduces the frequency of CH rotation and thus reduces energy consumption.

    Secondly,we develop a pre-selection mechanism that pre-selects a portion of cluster members to apply for tracking the target.Some unsuitable nodes are prevented from transmitting data to CH,and some transmission consumption is reduced.

    Finally,a mathematical problem is formulated and solved to obtain the final task set from the candidate task set.The tracking accuracy,the total energy consumption,and the balance degree of the residual energy of cluster members are jointly considered.Simulation experiments are conducted and the results show that the proposed heuristic multi-node collaborative scheduling mechanism(HMNCS)has the most superior performance.

    The rest of this paper is organized as follows.Section 2 presents the related work.Section 3 introduces the system model.The proposed protocol is explained in Section 4.Section 5 includes simulation experiments and performance analysis.Finally,Section 6 concludes the paper.

    2 Related Work

    2.1 Existing Researches

    In WSN for target tracking,there are many researchers dedicated to research on node scheduling.In[3],Energy proficient clustering(EEPC)was proposed to reduce energy consumption.The placement and energy level are used to select the CH and the relay nodes are obtained by calculating particle fitness value with their velocity and location.In [7],researchers proposed an improved Unscented Kalman Filter (UKF) that considers the uncertainty of the network and unpredictability of target motion.They utilized a genetic adaptive multi-objective algorithm to obtain a better set of task nodes without knowing the number of sensors to be selected.In [8],a scheduling mechanism based on a Type-2 Fuzzy Logic System(T2FLS)and a software agent approach was developed.The mechanism integrates energy,the number of neighboring nodes,and available bandwidth.In [9],researchers applied target state dynamics to predict the target trajectory and used the particle swarm algorithm for network optimization.A coordinated tracking mechanism was designed to save overall network energy consumption.In [10],Zhu et al.proposed an environment-adaptive event-driven robust set affiliation fusion estimation mechanism to schedule the required number of nodes to track the target.Experiments demonstrated that it reduces energy consumption and computational complexity while ensuring localization accuracy.In [11],the sensor selection problem was transformed into a multiobjective optimization framework,and a new mutual information upper bound(MIUB)–based sensor selection mechanism was proposed.The authors studied several sensor selection strategies,such as how to obtain a task set when minimizing the number of selected sensors or minimizing the difference between the performance obtained when all sensors transmit data and the performance obtained when only task nodes transmit data.In[12],researchers proposed a GPB1-UIF(the Generalized Pseudo-Bayesian estimator of first order and the unscented information filter) algorithm with estimation feedback and an information quality metric framework to obtain the task set using a multi-objective optimization strategy.In[13],Fu et al.analyzed the target detection probability and residual energy of sensors,proposing a new energy-balanced sensor node scheduling algorithm to balance tracking quality and network lifetime.In[14],Alagha et al.considered factors such as residual energy,power cost,and data confidence,proposing a new data-driven active node selection mechanism that uses genetic and greedy algorithms to select the optimal task set.In [15],the scheduling problem was transformed into a partially observable Markov decision process(POMDP)optimal policy problem by analyzing the intrinsic relationship between tracking performance and energy consumption.A dynamic cluster membership scheduling (DCMS) algorithm was proposed to solve the trade-off problem.In[16],spatial and temporal correlations among nodes were utilized instead of just spatial correlation.With the support of reconstructed data from the previous moment,a constrained convex relaxation plus rounding method was applied to obtain the task set,and the proposed mechanism was experimentally shown to extend the network lifetime.In[17],an accuracy enhancement algorithm and an active node selection algorithm based on the robust interval fusion preference aggregation(IF&PA)method were developed to significantly extend the network lifetime while providing high-accuracy measurements.Measurement uncertainty can be reduced even with a small number of sensors.In[18],game theory was used to reduce redundancy,and only one active node was periodically selected in each region to transmit data.Extensive simulation results showed that the proposed mechanism effectively reduced energy waste and extended the network lifetime.Finally,in[19],a cooperative node selection mechanism considering the spatial correlation among nodes was proposed to dynamically select the optimal set of nodes for the tracking task,significantly reducing energy consumption without degrading tracking accuracy.In [20],the proposed mechanism selects the set of sensors based on the location of the target and the mobile sink,thereby extending the network lifespan and reducing complexity.In[21],the researchers proposed an innovative glowworm swarm optimization technique(SS-GSO)to appropriately divide the clusters and obtain the optimal CHs.The residual energy,total energy consumption,and other factors are considered to efficiently schedule the nodes.The simulation results show that SS-GSO performs better in energy consumption and lifetime compared with other techniques.In [22],the researchers developed a long-term dynamic resource allocation algorithm to obtain the optimal resource scheduling solution for accurate and consecutive tracking.The algorithm offloads time-sensitive tasks to sink to improve computing efficiency.

    2.2 Remaining Weaknesses

    It can be seen that researchers have explored different ways to eliminate redundancy from multiple perspectives in Section 2.1.However,certain limitations still exist within the current approaches.Firstly,there is a lack of comprehensive consideration for the rotation process.Since target tracking is an ongoing task,frequent rotations can lead to substantial energy wastage.Secondly,in the majority of existing mechanisms,all cluster members are required to transmit pertinent information to the CH for selection,without taking into account the energy consumption involved in this process.In reality,if certain criteria can be established for cluster members to halt the transmission of information from nodes that do not contribute to prolonging the network lifespan,it could result in significant energy savings.

    3 System Model

    3.1 Network Model

    Without losing generality,we consider an L ?L monitored area where N nodes are randomly deployed.The target enters the area from the boundary and moves within the area according to certain rules until leaving the network.To save energy and continuously monitor the target,before the node perceives the target,only the sensing module is turned on,and the computing and transmission modules are turned off,that is,the sensor node is in a dormant state.After sensing the target,the node is awakened and enters the wake-up state,which means all modules are turned on.After awakened nodes are self-organized into a cluster,CH collects and fuses the data about the target from task nodes,and then transmits it to the sink.Here the assumptions are as follows:

    1.Sink is in the center of the region and has unlimited energy.

    2.The transmission radius of the sensor is more than twice the sensing radius to ensure that all nodes in the cluster can communicate with each other.

    3.Nodes can obtain their positions through the global positioning system (GPS) or any other positioning algorithm.

    Fig.1 shows the schematic diagram of the network scene.

    Figure 1: The schematic diagram of the network scene

    3.2 Target Movement Model

    We denote the state of the target atkthmoment isXK=,k=1,2,3,···,wherexkandykare the position coordinates of the target in the horizontal and vertical directions,respectively,andandare the velocities of the target in the horizontal and vertical directions,respectively.

    The state transition equation of the target is as follows:

    3.3 Energy Consumption Model

    The most common energy consumption model is adopted,as shown below[13]:

    whereεfsandεampare the energy consumed by each bit of data through the power amplifier for the free space and multi-path wireless channel model,respectively,andd0is the threshold distance between the two models.Eelecis the energy consumed by each bit of data transmitted or received by the sensor.dis the distance between the sender and the receiver.

    4 The Proposed Protocol

    The proposed protocol HMNCS will be described in detail in this section.Section 4.1 describes the process of target tracking.Section 4.2 presents the CH selection mechanism.Section 4.3 includes the pre-selection mechanism.Finally,the task set selection mechanism is introduced in Section 4.4.

    4.1 The Overall Framework of HMNCS

    Throughout the journey of the target within the network,a set of strategies has been devised to enhance tracking accuracy and prolong the network lifespan.These strategies encompass the CH selection mechanism,the pre-selection mechanism,and the task set selection mechanism,collectively referred to as the heuristic multi-node collaborative scheduling mechanism(HMNCS),illustrated in Fig.2.

    Figure 2: The block diagram of HMNCS

    The target enters WSN from the area boundary,triggering the awakening of the node responsible for sensing the target.In the initial round,since no node possesses complete information about the cluster at the outset,the first awakened node assumes the role of CH,while the remaining nodes become cluster members.If an awakened node fails to receive the announcement message transmitted by other nodes,it broadcasts its announcement message,declaring itself as the temporary CH,and awaits join messages from other nodes.However,if the awakened node successfully receives the announcement message from other nodes,it sends a join message to the temporary CH,seeking admission into the cluster.The temporary CH responds by sending a confirmation message to the awakened nodes,verifying their inclusion in the cluster.Subsequently,the temporary CH employs the CH selection mechanism to determine the actual CH and broadcasts the information in the cluster.Once the cluster members are aware of the true CH,the temporary CH relinquishes its role and becomes a cluster member.

    After a cluster is formed,each node obtains the number of awakened nodes in its neighborhood by exchanging information between itself and neighbor nodes.Then it determines whether it applies to CH to be a task node according to the pre-selection mechanism,which will be introduced in Section 4.3.If it meets the pre-selection conditions,it will send an applying message including the residual energy and the observation of the target state to CH.CH derives the task set by applying the task set selection mechanism,which will be described in Section 4.4,and broadcasts the scheduling result in the cluster.Those selected task nodes keep awake to track the target,while those non-task nodes enter the dormant state to save energy.Task nodes keep tracking and send observations to CH.CH fuses the collected information,uses UKF to obtain the estimated target state,and then sends it to the sink.In the second and subsequent rounds,those nodes not aware of the target get out of the cluster,and those newly awakened nodes join the cluster.And as the target moves,CH and the task set may need to rotate.CH determines whether it needs rotating or not.If it does,CH obtains the new CH according to the CH selection mechanism and broadcasts it in the cluster.The newly designated CH assumes the responsibilities of the CH,while the former CH transitions into a cluster member or departs from the cluster entirely.CH also needs to determine whether the task set needs rotating or not.When the task set needs rotating,CH gets the new task set by using the task set selection mechanism and broadcasts it.Task nodes diligently track the target,while non-task nodes enter a dormant state to conserve energy.This process iterates until the target leaves the monitored area.

    In this paper,we assume that the target moves slowly,it is usually a long time from a task set selected to it is rotated.During the process,those nodes that meet the pre-selection conditions only send an apply message once.Hence the energy consumption of applying accounts for a small part of total consumption and is acceptable.

    4.2 The CH Election Mechanism

    The proposed CH election mechanism is shown in Algorithm 1.

    When CH meets the constraints in Eq.(10)or Eq.(11),it needs rotating.

    In Algorithm 1,the first and second parts of Eq.(8)are included to make the transmitting distance between cluster members and CH short to reduce energy consumption.The average residual energy of CLUSTER MEMBERs is used as a threshold to avoid nodes with low energy becoming CH.In Eq.(10),when the residual energy of CH is less than a certain value,CH needs to rotate to make sure the CH node has enough energy to complete its task in the next round.It should be noted that the energy consumption of a CH exceeds that of a cluster member.If nodes with low residual energy were to become CH,they would have a short lifespan,leading to an accelerated emergence of coverage holes,an increased frequency of CH rotation,and energy wastage[10].Moreover,the relative position between the target and the node is also taken into account.Nodes that are approaching the target are more likely to become CH,thus delaying CH rotation caused by failure to sense the target.

    4.3 The Pre-Selection Mechanism

    Considering redundancy,all nodes do not have to participate in tracking tasks.Instead,nodes transmit pertinent information to CH to express their interest in becoming task nodes,and CH selects a set of task nodes from these applicants.To minimize energy consumption during the application process and alleviate the computational load on the CH,a pre-selection mechanism is developed.This mechanism enables each node to decide whether or not to apply based on factors such as residual energy,estimation error,node density,and other relevant criteria.

    The requirements which nodes need to satisfy to apply are as follows:

    The nodes that meet the pre-selection requirements apply to CH to become a task node,and those nodes that do not meet the conditions turn into a dormant state directly until the next task set rotation or is awakened by the next target.

    The flow chart for determining whether to apply is shown in Fig.3.

    Figure 3: The flowchart of the pre-selection mechanism

    4.4 The Task Set Selection Mechanism

    The selection of the task node set takes into account both the estimation accuracy and the network lifespan,where the latter is influenced by the total energy consumption and the energy balance level.In this paper,UKF is employed to estimate the true state of the target.The trace of the error covariance matrix in UKF serves as an indicator of the average square error in state estimation [10].The total energy consumption is the sum of the energy consumption of each node in WSN,and the energy balance level is measured by the variance of the residual energy sequence of nodes in the cluster.Their expressions are as follows:

    whereεTSis the average square error obtained by using observations of task nodes to estimate the real target state through UKF,andPis the error covariance matrix in UKF.is the expected energy consumption ofithnode if it becomes a task node.and varTSare the expected total energy consumption of the network and the expected variance of the residual energy sequence of nodes in the cluster,respectively,after the task setTScompletes the tracking task in the current round.

    The optimization objective function for selecting the task set is as follows:

    Eq.(20) is an optimization problem.If the exhaustive algorithm is used to get the results,there will be++···+=2m-1 candidate task sets(assuming there aremnodes in the cluster),and the number increases withmat an exponential rate.Therefore,it is necessary to apply an algorithm with low complexity to solve this problem.

    In this paper,the genetic algorithm is adopted.CH is always in the task set,and task nodes are selected from those applying nodes.Hence the length of the chromosome is equal to the number of applying nodes,and each gene on the chromosome represents a node.Each node has two states:Selected and unselected,so binary encoding can be perfectly applied to this case.Each gene has two values,1 for selected and 0 for unselected.Each chromosome corresponds to a task set.Suppose that there are eight applying nodes whose numbers are 11,32,45,8,55,135,146,and 16,respectively.If the chromosome code is “01110100”,the second,third,and sixth nodes are selected,that is,the corresponding task set is 32,45,8,135,as shown in Fig.4.

    Figure 4: The comparison diagram of chromosomes and task sets

    According to the objective function,the fitness function of the candidate task set is set,as shown in Eq.(21).

    Due to significant differences in absolute values,it is necessary to normalize the average square estimation error,energy consumption,and variance.Normalization is carried out for all chromosomes in the population.The task set selection algorithm is shown in the Algorithm 2.

    5 Simulation Experiments and Performance Analysis

    5.1 Simulation Scenarios and Parameter Settings

    In this section,the performance of HMNCS,the dynamic chain-based collaboration mechanism(DCBC)[23],the greedy balance replace heuristic algorithm(GBRHA)[13]and the node scheduling algorithm for all-directional intrusion detection(NSAID)[24]are evaluated on MATLAB.In DCBC,the target tracking task is completed by forming a dynamic tracking chain around the target.GBRHA is an energy-balanced sensor node scheduling algorithm.By analyzing the target detection probability and the residual energy of the sensor nodes,multiple sensor nodes are scheduled for cooperative target tracking.In NSAID,the monitored region is divided into the outermost and inner regions,which can effectively save energy while ensuring the quality of all-around detection.

    In the simulation,the monitored area is set to 300 m ?300 m,and the sink is in the center of the area,with infinite energy.To evaluate the performance of HMNCS under different deployment densities,experiments with 400,600,and 800 nodes are carried out.The sensing radius of nodes is set to 40 m,and the transmission radius is set to twice the sensing radius to ensure that all nodes in the cluster can communicate with each other.The sampling interval is 1 s,and every 5 s is a round.The target enters the network from the random position of a boundary and moves according to the state transition Eq.(1),until it leaves the network.Fusion energy consumption means the energy consumed by fusing received data packets and its own data packets.The settings of simulation parameters are shown in Table 1.

    Table 1: The settings of simulation parameters

    The performance criteria used in this paper include the number of alive nodes,the total energy consumption,the tracking error,and the network lifespan.The lifespan can be usually measured in terms of FND (first node death),HND (half node death),and LND (last node death).However,in large-scale WSN,the death of the first node has little effect on the connectivity and coverage of the network,so HND and LND are adopted.The average value and variance of the tracking error sequence from deployment to the round when the last node runs out of energy are calculated to indicate the tracking performance.

    5.2 Performance Analysis

    The curves of the number of alive nodes and the total energy consumption of HMNCS,DCBC,GBRHA,and NSAID,when 400,600,and 800 nodes are deployed in the network,are plotted in Figs.5 and 6,respectively.In the long term,the number of alive nodes of HMNCS is higher than that of other compared mechanisms after the same rounds with different node deployment densities.In other words,the nodes of HMNCS die slower than those of other mechanisms,and the total energy consumption of HMNCS is the lowest after the same rounds.This is because HMNCS uses various methods to reduce energy consumption.Firstly,in the CH selection mechanism,not only the residual energy of nodes,the distance to other nodes,the distance to the sink,etc.,are considered,but also the relative movement trend between the target and each node is taken into account.The nodes that are closer and closer to the target are more likely to be CH,which reduces the rotation energy consumption.Further,using the event trigger mechanism avoids a part of applications and saves energy.Taking energy consumption and energy balance factors into account when the task set is finally obtained,all these operations result in lower energy consumption and longer network lifetime.As the deployment density of nodes increases,it is found that the advantage of HMNCS over the compared mechanisms increases,indicating that HMNCS is more suitable for scenarios with high node density.

    Figure 5: The number of alive nodes under different densities

    Figure 6: The total energy consumption under different densities

    However,in the early stage after deployment,the number of alive nodes of HMNCS is slightly lower than that of DCBC,and the gap decreases with the increase of node deployment density.In HMNCS,after the target appears,a cluster is formed in WSN.Then CH and a task set are selected,and only some nodes track the target,while the rest of the CLUSTER MEMBERs enter the dormant state.In DCBC,the nodes that sense the target form a tracking chain,and the chain head is selected.Data of all nodes is transmitted to the chain head from both ends of the tracking chain node-by-node and fused.Each node in the cluster transmits and fuses data,but the transmission distance is short in DCBC.In HMNCS,only the task set nodes in the cluster transmit data directly to CH,but the transmission distance is longer than that of DCBC.In DCBC,the energy consumption of each node is low,so nodes die slowly in the beginning.However,there are so many nodes participating in data transmission that the total energy consumption is large.In HMNCS,task nodes consume more energy,while non-task nodes hardly consume energy.As a result,the number of alive nodes in the network is slightly lower than that of DCBC in the early stage after deployment,but with the increase of rounds,the advantages of HMNCS are gradually emerging,and the number of alive nodes is the highest.

    Fig.7 shows the average tracking error and variance of tracking error from deployment to the round of last node death for 400,600,and 800 nodes.Fig.8 shows HND and LND under three deployment densities.The average value and variance of tracking error of HMNCS are lower than those of compared mechanisms,and HND and LND are the greatest,i.e.,HMNCS is the best in terms of tracking performance and network lifetime.

    Figure 7: The tracking error under different densities

    Figure 8: HND&LND under different densities

    The factors considered in HMNCS are relatively comprehensive.In the pre-selection mechanism,the tracking error of the previous moment is used to adjust the number of applying nodes at the current moment.Moreover,the average square estimation error is considered in the optimization objective function.The expected energy consumption and the residual energy of nodes are used to filter whether to apply or not.Energy consumption,and energy balance factors,which directly affect the network lifespan,are also included in the objective function.By setting appropriate weights for different factors,good performance can be achieved in both network lifespan and tracking performance.

    It is worth noting that DCBC cannot achieve high tracking accuracy in the later stage due to its fast node death rate.The early period of high tracking accuracy does not last for a long time,so tracking performance is not good in the long run.

    For all the mechanisms,the lifetime of WSN increases with the increase in node density.Among them,the WSN that uses HMNCS has the most extended lifespan,so it is more suitable for highdensity WSN.

    Based on the above simulation results,it can be concluded that the performance of HMNCS is superior to that of the compared mechanisms.In every step of HMNCS,such as the CH election mechanism,the pre-selection mechanism,and the optimization objective function of getting the final task set,the accuracy,energy consumption,and energy balance factors are considered in different ways.Therefore,HMNCS achieves a clever balance between the tracking accuracy and the network lifetime and obtains good results in both aspects of performance over the long term.

    6 Conclusion

    In this paper,an innovative two-layer node scheduling protocol called the Heuristic Multi-Node Collaborative Scheduling mechanism,is proposed to obtain optimal scheduling results during target tracking in WSN.HMNCS consists of three parts,which are cluster head election,pre-selection,and task set selection mechanism,respectively.By reducing the redundancy and making those necessary nodes track the target,this work balances the tracking performance and the network lifetime as much as possible.Simulation results confirm the superiority of the proposed mechanism.For HND,HMNCS is improved by 31.1%,83.9%,and 112.3% compared with DCBC,GBRHA and NSAID,respectively.For LND,HMNCS is improved by 34.3%,72.8%,and 89.0% compared with DCBC,GBRHA and NSAID,respectively.

    With the development of unmanned aerial vehicles(UAVs)in recent years,it is very common that wireless sensor network deployed on the ground cooperate with UAVs to fulfill tasks.UAVs can be regarded as movable sensor nodes,and all nodes form a three-dimensional heterogeneous wireless sensor network.In such circumstances,how to schedule nodes and plan the flight paths of UAVs efficiently to achieve reliable data collection and other requirements in the actual scenario is our future research work.

    Acknowledgement:The authors would like to express their gratitude to the members of the research group for their support.

    Funding Statement:This work is supported by the Project Program of Science and Technology on Micro-System Laboratory,No.6142804220101.The author who received the grant is identified by Xue Zhao.

    Author Contributions:The authors confirm contributions to the paper as follows: Study conception and design: Xue Zhao;simulation: Shaojun Tao and Hongying Tang;analysis,interpretation of results:Jiang Wang;draft manuscript preparation:Baoqing Li.All authors reviewed the results and approved the final version of the manuscript.

    Availability of Data and Materials:The authors confirm that the data and materials supporting the findings of this study are not applicable as there were no specific datasets or materials used.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    免费无遮挡裸体视频| 美女大奶头视频| 在线免费十八禁| 观看美女的网站| 一a级毛片在线观看| 欧美日韩一区二区视频在线观看视频在线 | 能在线免费观看的黄片| 深夜精品福利| 黄色配什么色好看| 亚洲第一电影网av| 麻豆成人午夜福利视频| 一个人免费在线观看电影| 亚州av有码| 蜜臀久久99精品久久宅男| 国产午夜精品论理片| 99久久成人亚洲精品观看| 99热全是精品| 六月丁香七月| 丰满乱子伦码专区| 哪里可以看免费的av片| 午夜影院日韩av| 国产高清不卡午夜福利| 黄色配什么色好看| 好男人在线观看高清免费视频| 日韩av在线大香蕉| 波多野结衣高清作品| 看非洲黑人一级黄片| 久久久午夜欧美精品| 色综合色国产| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品成人久久小说 | 婷婷精品国产亚洲av| 中文在线观看免费www的网站| 日本一二三区视频观看| 综合色丁香网| 国产高清视频在线观看网站| 热99在线观看视频| 人人妻人人澡欧美一区二区| av.在线天堂| .国产精品久久| 两个人视频免费观看高清| 丰满人妻一区二区三区视频av| 日本成人三级电影网站| 国产日本99.免费观看| 波多野结衣高清作品| 不卡一级毛片| 99精品在免费线老司机午夜| 乱系列少妇在线播放| 嫩草影院精品99| 少妇丰满av| 国产男靠女视频免费网站| 日本-黄色视频高清免费观看| 中国美白少妇内射xxxbb| 精品久久久噜噜| 日韩欧美国产在线观看| 国内精品宾馆在线| 国产一区二区在线av高清观看| 国产老妇女一区| 男女那种视频在线观看| 亚洲高清免费不卡视频| 久久久久久久久大av| 国产免费一级a男人的天堂| 内射极品少妇av片p| 国产综合懂色| 国产日本99.免费观看| 亚洲精品国产av成人精品 | 亚洲性夜色夜夜综合| 成人国产麻豆网| 12—13女人毛片做爰片一| 精品国内亚洲2022精品成人| 久久九九热精品免费| 日本a在线网址| 久久久久久大精品| 久久久久久久久久黄片| 少妇丰满av| 国产精品久久久久久亚洲av鲁大| 欧美一级a爱片免费观看看| 精品欧美国产一区二区三| 香蕉av资源在线| 亚洲乱码一区二区免费版| 日韩精品中文字幕看吧| a级一级毛片免费在线观看| 一卡2卡三卡四卡精品乱码亚洲| 最近2019中文字幕mv第一页| 一级毛片我不卡| 国产黄色小视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av熟女| av天堂在线播放| 成人一区二区视频在线观看| 成熟少妇高潮喷水视频| 亚洲aⅴ乱码一区二区在线播放| 一a级毛片在线观看| 又黄又爽又免费观看的视频| 亚洲在线自拍视频| 淫妇啪啪啪对白视频| 国产在线精品亚洲第一网站| 精品一区二区三区av网在线观看| 成人美女网站在线观看视频| 精品国内亚洲2022精品成人| 日韩欧美在线乱码| av专区在线播放| 一本精品99久久精品77| 亚洲欧美精品综合久久99| 91精品国产九色| 亚洲精品亚洲一区二区| 熟妇人妻久久中文字幕3abv| 99久久成人亚洲精品观看| 免费在线观看成人毛片| 国产亚洲精品久久久com| 国产成人a∨麻豆精品| 久久精品综合一区二区三区| 日韩三级伦理在线观看| 日韩三级伦理在线观看| 搡女人真爽免费视频火全软件 | 亚洲熟妇熟女久久| 国产精品,欧美在线| 丰满人妻一区二区三区视频av| 91麻豆精品激情在线观看国产| 免费av毛片视频| 久久精品国产自在天天线| 22中文网久久字幕| 亚洲人成网站在线播放欧美日韩| 三级经典国产精品| 熟妇人妻久久中文字幕3abv| 成年av动漫网址| 99久久精品国产国产毛片| 亚洲欧美清纯卡通| 最近在线观看免费完整版| 国产欧美日韩一区二区精品| 一进一出抽搐动态| 免费在线观看成人毛片| 美女高潮的动态| 成人二区视频| 免费黄网站久久成人精品| 一本精品99久久精品77| 大又大粗又爽又黄少妇毛片口| 免费一级毛片在线播放高清视频| 在线播放无遮挡| 内射极品少妇av片p| 色综合色国产| 人妻丰满熟妇av一区二区三区| 亚洲久久久久久中文字幕| 国产精品嫩草影院av在线观看| 日本a在线网址| 欧美一区二区亚洲| 亚洲国产精品sss在线观看| 亚洲一区二区三区色噜噜| 波野结衣二区三区在线| 免费一级毛片在线播放高清视频| 亚洲av不卡在线观看| 简卡轻食公司| 狂野欧美白嫩少妇大欣赏| 午夜激情欧美在线| 国产精品伦人一区二区| 桃色一区二区三区在线观看| 啦啦啦啦在线视频资源| 三级毛片av免费| 性色avwww在线观看| 久久精品影院6| 成人三级黄色视频| 亚洲成人久久性| 成人精品一区二区免费| 国产又黄又爽又无遮挡在线| 亚洲av电影不卡..在线观看| 久久精品国产亚洲网站| 能在线免费观看的黄片| 秋霞在线观看毛片| 亚洲国产精品成人久久小说 | 一进一出好大好爽视频| 少妇丰满av| 97人妻精品一区二区三区麻豆| 亚洲精品在线观看二区| 一个人看视频在线观看www免费| 在线播放国产精品三级| 国产91av在线免费观看| 色尼玛亚洲综合影院| 干丝袜人妻中文字幕| 国产日本99.免费观看| 能在线免费观看的黄片| 欧美日韩乱码在线| 麻豆国产97在线/欧美| 我的女老师完整版在线观看| 免费高清视频大片| 久久久精品欧美日韩精品| 观看美女的网站| 少妇猛男粗大的猛烈进出视频 | 狂野欧美激情性xxxx在线观看| 天堂√8在线中文| av在线天堂中文字幕| 69av精品久久久久久| 一进一出好大好爽视频| 神马国产精品三级电影在线观看| 欧美日本亚洲视频在线播放| 亚洲中文字幕一区二区三区有码在线看| 你懂的网址亚洲精品在线观看 | 国产精品人妻久久久影院| 欧美bdsm另类| 露出奶头的视频| 国产黄a三级三级三级人| 色综合站精品国产| 日韩成人av中文字幕在线观看 | 成人精品一区二区免费| 欧美绝顶高潮抽搐喷水| 在线播放无遮挡| 九九在线视频观看精品| 婷婷精品国产亚洲av在线| 国产精品久久久久久精品电影| 日韩欧美三级三区| 你懂的网址亚洲精品在线观看 | 国产三级中文精品| 国产一区二区三区av在线 | 18禁黄网站禁片免费观看直播| 久久精品国产亚洲av香蕉五月| 哪里可以看免费的av片| 一区二区三区四区激情视频 | 国产精品日韩av在线免费观看| 亚洲av第一区精品v没综合| 最近在线观看免费完整版| 国产真实乱freesex| 亚洲欧美清纯卡通| 午夜激情欧美在线| 又爽又黄无遮挡网站| 成人综合一区亚洲| 精品少妇黑人巨大在线播放 | 嫩草影院新地址| 插阴视频在线观看视频| 中文字幕精品亚洲无线码一区| 如何舔出高潮| 亚洲图色成人| 少妇的逼水好多| 小说图片视频综合网站| 亚洲精品一区av在线观看| 久久久久久伊人网av| 女生性感内裤真人,穿戴方法视频| 丝袜美腿在线中文| 国产麻豆成人av免费视频| 亚洲精品乱码久久久v下载方式| 日日摸夜夜添夜夜添小说| 亚洲av五月六月丁香网| 国产精品一区二区性色av| 一a级毛片在线观看| 91久久精品国产一区二区三区| 男人的好看免费观看在线视频| 精品免费久久久久久久清纯| 精品人妻熟女av久视频| 别揉我奶头~嗯~啊~动态视频| 麻豆精品久久久久久蜜桃| 欧美色视频一区免费| 国产精品一区www在线观看| 成人亚洲精品av一区二区| 色av中文字幕| 99热精品在线国产| 91狼人影院| 国产高清不卡午夜福利| 欧美区成人在线视频| 国产一区二区三区av在线 | 久久精品夜夜夜夜夜久久蜜豆| 噜噜噜噜噜久久久久久91| 日韩亚洲欧美综合| 久久人妻av系列| 99久久九九国产精品国产免费| 97碰自拍视频| 日韩中字成人| avwww免费| 免费电影在线观看免费观看| 久久久久久久午夜电影| 亚洲av第一区精品v没综合| 69av精品久久久久久| 久久九九热精品免费| 久久人人爽人人片av| 天堂影院成人在线观看| 搡老熟女国产l中国老女人| 在线观看免费视频日本深夜| 日本五十路高清| 免费高清视频大片| 嫩草影院入口| 成年免费大片在线观看| 国产在线精品亚洲第一网站| 两性午夜刺激爽爽歪歪视频在线观看| 嫩草影视91久久| 精品人妻视频免费看| 国产精品一及| 91久久精品国产一区二区三区| 亚洲va在线va天堂va国产| 亚洲av不卡在线观看| 亚洲精华国产精华液的使用体验 | 一级毛片aaaaaa免费看小| 三级国产精品欧美在线观看| 国产av不卡久久| 精品无人区乱码1区二区| 赤兔流量卡办理| 大型黄色视频在线免费观看| 哪里可以看免费的av片| 麻豆成人午夜福利视频| 精品久久久久久成人av| 成人一区二区视频在线观看| 久久久久国产精品人妻aⅴ院| 久久精品91蜜桃| 日韩制服骚丝袜av| 久久韩国三级中文字幕| 欧美中文日本在线观看视频| 国产精品不卡视频一区二区| 三级国产精品欧美在线观看| 国产伦精品一区二区三区视频9| 99久久无色码亚洲精品果冻| 亚洲av不卡在线观看| 亚洲国产精品合色在线| 欧美日韩精品成人综合77777| 国产黄色小视频在线观看| 女人被狂操c到高潮| 在线播放无遮挡| 97碰自拍视频| 女生性感内裤真人,穿戴方法视频| 久久国内精品自在自线图片| 亚洲天堂国产精品一区在线| 欧美一级a爱片免费观看看| 色综合色国产| 国内精品美女久久久久久| 久久精品国产亚洲av香蕉五月| 久久人人爽人人爽人人片va| 欧美xxxx黑人xx丫x性爽| 国产人妻一区二区三区在| 精品99又大又爽又粗少妇毛片| 亚洲av成人精品一区久久| 在线播放无遮挡| 无遮挡黄片免费观看| 变态另类丝袜制服| 亚洲欧美清纯卡通| 国产在视频线在精品| 亚洲中文字幕日韩| 97超视频在线观看视频| 日本一二三区视频观看| 一本精品99久久精品77| 一级黄色大片毛片| 91麻豆精品激情在线观看国产| 国产成人影院久久av| 九九爱精品视频在线观看| 亚洲av中文字字幕乱码综合| 真人做人爱边吃奶动态| 国产精品一区二区三区四区免费观看 | 久久久久国产精品人妻aⅴ院| 身体一侧抽搐| 香蕉av资源在线| 91久久精品国产一区二区成人| 国产片特级美女逼逼视频| 国产精品永久免费网站| 国语自产精品视频在线第100页| 午夜a级毛片| 亚洲精品影视一区二区三区av| 亚洲美女搞黄在线观看 | 一级毛片久久久久久久久女| 亚洲av不卡在线观看| 天天一区二区日本电影三级| 亚洲精品亚洲一区二区| 91久久精品国产一区二区三区| 男人和女人高潮做爰伦理| 成人一区二区视频在线观看| 12—13女人毛片做爰片一| 国产私拍福利视频在线观看| 插阴视频在线观看视频| 国产精品久久久久久久久免| 亚洲最大成人中文| 热99在线观看视频| 美女cb高潮喷水在线观看| 精品一区二区三区av网在线观看| 99久久久亚洲精品蜜臀av| 乱人视频在线观看| 少妇熟女欧美另类| 成人美女网站在线观看视频| av在线蜜桃| 99热只有精品国产| 禁无遮挡网站| av在线老鸭窝| 国产免费一级a男人的天堂| 床上黄色一级片| 亚洲五月天丁香| 晚上一个人看的免费电影| 国产高清不卡午夜福利| 五月玫瑰六月丁香| 一进一出抽搐动态| 久久久久久久久久成人| www日本黄色视频网| 成人无遮挡网站| 丝袜喷水一区| 精品熟女少妇av免费看| 久久精品综合一区二区三区| 国产69精品久久久久777片| 黄色配什么色好看| 亚洲av美国av| 亚洲av熟女| 成人鲁丝片一二三区免费| 精品久久久噜噜| 不卡一级毛片| 高清日韩中文字幕在线| 性插视频无遮挡在线免费观看| 成人二区视频| 精品少妇黑人巨大在线播放 | 日本a在线网址| 免费黄网站久久成人精品| 亚洲av二区三区四区| 免费不卡的大黄色大毛片视频在线观看 | 免费观看在线日韩| 成年女人永久免费观看视频| 女的被弄到高潮叫床怎么办| 亚洲最大成人中文| 国产一区二区三区av在线 | 亚洲丝袜综合中文字幕| 欧美一区二区精品小视频在线| 听说在线观看完整版免费高清| 老司机影院成人| 国产伦精品一区二区三区四那| 国内精品美女久久久久久| 变态另类丝袜制服| 久久99热6这里只有精品| 伊人久久精品亚洲午夜| 国产不卡一卡二| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区免费欧美| 国产爱豆传媒在线观看| 激情 狠狠 欧美| 色av中文字幕| 国内精品美女久久久久久| 久久久久久大精品| 亚洲欧美精品综合久久99| 亚洲中文字幕日韩| 男女边吃奶边做爰视频| 久久久国产成人精品二区| 亚洲经典国产精华液单| 免费高清视频大片| 三级毛片av免费| h日本视频在线播放| 亚洲国产欧洲综合997久久,| 一级黄片播放器| 一级黄色大片毛片| 女生性感内裤真人,穿戴方法视频| aaaaa片日本免费| 高清午夜精品一区二区三区 | 国产午夜精品久久久久久一区二区三区 | 97超视频在线观看视频| 在线免费观看不下载黄p国产| 成人高潮视频无遮挡免费网站| 国产私拍福利视频在线观看| 久久久国产成人免费| 亚洲性久久影院| www日本黄色视频网| 日韩高清综合在线| 午夜日韩欧美国产| 搡老岳熟女国产| 午夜精品在线福利| 欧美+日韩+精品| 久久久久久伊人网av| 狂野欧美白嫩少妇大欣赏| 少妇熟女aⅴ在线视频| 国产私拍福利视频在线观看| 最近的中文字幕免费完整| 内射极品少妇av片p| 亚洲熟妇中文字幕五十中出| 春色校园在线视频观看| 一卡2卡三卡四卡精品乱码亚洲| 人妻夜夜爽99麻豆av| 一个人免费在线观看电影| 国产 一区 欧美 日韩| 色在线成人网| 夜夜爽天天搞| 国产激情偷乱视频一区二区| 亚洲国产精品sss在线观看| 美女cb高潮喷水在线观看| 尾随美女入室| 国产亚洲91精品色在线| 欧美3d第一页| 寂寞人妻少妇视频99o| 嫩草影院精品99| 色视频www国产| 国产精品99久久久久久久久| 日韩欧美一区二区三区在线观看| 插阴视频在线观看视频| 在线国产一区二区在线| 国产高清激情床上av| 人妻少妇偷人精品九色| 欧美激情国产日韩精品一区| 高清日韩中文字幕在线| 色视频www国产| 在线天堂最新版资源| 中国国产av一级| 又黄又爽又刺激的免费视频.| 日本与韩国留学比较| 亚洲精品国产av成人精品 | 国产私拍福利视频在线观看| 麻豆av噜噜一区二区三区| 两个人的视频大全免费| 中国国产av一级| 精品一区二区三区av网在线观看| 亚洲精品日韩在线中文字幕 | 听说在线观看完整版免费高清| 午夜福利成人在线免费观看| 亚洲电影在线观看av| 听说在线观看完整版免费高清| 丝袜喷水一区| 日本五十路高清| 午夜老司机福利剧场| 麻豆一二三区av精品| 丰满的人妻完整版| 狂野欧美白嫩少妇大欣赏| 全区人妻精品视频| 天美传媒精品一区二区| 可以在线观看毛片的网站| 婷婷色综合大香蕉| 亚洲国产色片| 欧美人与善性xxx| 国产成年人精品一区二区| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看 | 九九热线精品视视频播放| 亚洲欧美精品综合久久99| 久久久久久九九精品二区国产| 少妇猛男粗大的猛烈进出视频 | 日本精品一区二区三区蜜桃| 亚洲欧美精品综合久久99| 免费看a级黄色片| 桃色一区二区三区在线观看| 亚洲真实伦在线观看| 97超级碰碰碰精品色视频在线观看| 男人舔奶头视频| 麻豆精品久久久久久蜜桃| 国产精品福利在线免费观看| 熟女人妻精品中文字幕| 在线观看免费视频日本深夜| 美女内射精品一级片tv| 大香蕉久久网| 精品福利观看| 精品乱码久久久久久99久播| 婷婷六月久久综合丁香| 伦精品一区二区三区| 亚洲av美国av| 直男gayav资源| 搡老岳熟女国产| 国产黄色视频一区二区在线观看 | 久久久成人免费电影| 午夜老司机福利剧场| 日本熟妇午夜| 久久综合国产亚洲精品| 国产亚洲精品综合一区在线观看| 亚洲熟妇中文字幕五十中出| 国产av一区在线观看免费| 日本色播在线视频| 国产亚洲欧美98| 中文字幕av成人在线电影| 亚洲欧美清纯卡通| 国产免费男女视频| 日韩欧美三级三区| 秋霞在线观看毛片| 日本黄大片高清| av在线播放精品| 国产乱人偷精品视频| 欧美又色又爽又黄视频| 国产淫片久久久久久久久| 九九在线视频观看精品| 日日摸夜夜添夜夜爱| 亚洲丝袜综合中文字幕| 国产蜜桃级精品一区二区三区| 99热这里只有精品一区| 国内少妇人妻偷人精品xxx网站| 在线观看午夜福利视频| 伊人久久精品亚洲午夜| 精品久久久久久久久av| 亚洲天堂国产精品一区在线| 亚洲自拍偷在线| 日韩 亚洲 欧美在线| 三级国产精品欧美在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲精品av在线| 国产亚洲91精品色在线| 久久久久精品国产欧美久久久| 久久国产乱子免费精品| 能在线免费观看的黄片| 国产爱豆传媒在线观看| 99久国产av精品国产电影| 欧美区成人在线视频| 91久久精品国产一区二区成人| 日韩精品中文字幕看吧| 精品免费久久久久久久清纯| 国产精品人妻久久久久久| a级毛片a级免费在线| 国产av麻豆久久久久久久| 亚洲内射少妇av| 人妻丰满熟妇av一区二区三区| 亚洲国产欧洲综合997久久,| 国产在线精品亚洲第一网站| 99热网站在线观看| 国产精品永久免费网站| 99久久精品热视频| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产自在天天线| 免费看av在线观看网站| 国产精品无大码| 精品日产1卡2卡| 少妇被粗大猛烈的视频| 99久国产av精品| 3wmmmm亚洲av在线观看| 99热这里只有精品一区| 一a级毛片在线观看| 内地一区二区视频在线| 一区福利在线观看| 亚洲美女搞黄在线观看 | 男插女下体视频免费在线播放| 欧美日韩一区二区视频在线观看视频在线 | 天堂√8在线中文| 麻豆一二三区av精品| 欧美性猛交╳xxx乱大交人| 久久精品国产亚洲av涩爱 | 在线观看美女被高潮喷水网站| 久久久久免费精品人妻一区二区| 偷拍熟女少妇极品色| 不卡视频在线观看欧美|