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

    Transmitter-Receiver Path Selection for Cell Range Extension Using Multi-Hop D2D

    2021-12-11 13:30:50FarahAkifKiranSultanAqdasMalikIjazQureshiandSabaMahmood
    Computers Materials&Continua 2021年8期

    Farah Akif,Kiran Sultan,Aqdas N.Malik,Ijaz M.Qureshi and Saba Mahmood

    1Department of Electrical Engineering,International Islamic University,Islamabad,44000,Pakistan

    2Department of CIT,FAS,King Abdulaziz University,Jeddah,21551,Saudi Arabia

    3Department of Electrical Engineering,Air University,Islamabad,44000,Pakistan

    4Department of Computer Science,Bahria University,Islamabad,44000,Pakistan

    Abstract:Conventional approach of dealing with more users per coverage area in cellular networks implies densifying the amount of(Access Point)AP which will eventually result in a larger carbon footprint.In this paper, we propose a base station off-loading and cell range extension (CRE) scheme based on multi-hop device-to-device(MHD2D)path selection between transmitter and receiver node.The paper also provides derivations of upper and lower bounds for energy efficiency,capacity,and transmit power.The proposed path selection scheme is inspired by the foraging behavior of honey bees.We present the algorithm as a modified variant of the artificial bee colony algorithm(MVABC).The proposed optimization problem is modeled as a minimization problem where we optimize the Energy Efficiency (EE).The proposed path selection MVABC is compared with the Genetic Algorithm (GA) and also with classical artificial bee colony (ABC) through simulations and statistical analysis.The student’s t-test, p-value, and standard error of means (SEM)clearly show that MVABC based path selection out-performs the GA and classical ABC schemes.MVABC based approach is 66%more efficient when compared with classic ABC and about 62%efficient when compared with GA based scheme.

    Keywords: Multi-hop; ultra-dense; D2D; path-selection; ABC; GA

    1 Introduction

    Conventional cellular networks provide coverage to a user by having a direct connection established between a user and the base station.Although the scheme provides communication with less delay, the recent trend of increased wireless data users and devices makes it difficult for the conventional scheme to fulfill the QoS (Quality of Service) requirements [1].This apparent increase in the data requirement or increase in the number of users reduces network reliability and efficiency when it is not able to cater to new users or provide dependable services to the edge users.Hence, the effect is similar to having a cell with reduced coverage.Increasing the number of BS in a cell may solve this issue but at the expense of increased BS and equipment installation charges along with an overall increase in power consumption.D2D communication promises to provide an alternative to this problem by allowing mobile devices to communicate directly with each other using low power in quite possibly the same band being used by some other cellular device [2].This concept is based on space division multiple access.

    From 3GPP (third generation global partnership project) LTE (long term evolution) release 12 and onwards, it is quite apparent that a MHD2D network is set to become a necessary part of 5G networks.Standards for D2D communication along with two-hop communication are already part of the latest release [3].Moreover, the cell range of base stations can be extended by utilizing MHD2D relay nodes.To ensure that the cell range extension (CRE) offered by D2D links meet the required QoS demand, the energy efficiency (EE) metric is commonly used [4] which is generally defined as:

    1.1 Related Work

    MHD2D is an emerging technology capable of providing service to user-driven demand.Catering for the similar user demand, the authors in [5] developed an algorithm for long-range D2D links which is done through multi-hop communication.However, the study lacks any insight regarding optimal path selection.In [6] the authors developed an algorithm for power control and routing in MHD2D where Dijkstra shortest path algorithm is used for routing.A drone-based D2D aided multi-hop network for emergency services is proposed in [7] where stochastic geometry is used for link coverage probability to minimize delay, however, the path selected for emergency message transmission is not considered in this paper.

    A standard Dijkstra shortest path-based routing algorithm is developed in [8] for D2D for multi-hop communication in the Internet of Things (IoT) application.D2D based solution for offloading traffic in a cellular network is proposed in [9], the proposed framework offers a multihop routing scheme for different types of data transmissions.In [10] a MHD2D energy-efficient routing protocol is proposed to perform resource allocation while improving throughput and packet delivery ratio.MHD2D routing for the millimeter-wave network is proposed in [11] where received signal strength (RSS) is measured for relay nodes and an algorithm is designed to select a route based on the probability of route selection as generated by their RSS values.

    Interference aware D2D routing is discussed in [12] where an adaptive algorithm is proposed which selects a route based on node location and amount of interference at next hop.Optimum path selection for 5G vehicular networks based on mmWave D2D links is discussed in [13] where the optimum path is found based on vehicular information rather than link quality.The authors in [14] investigate multi-hop data transfer between D2D transmitter and receiver by scheduling cellular resource blocks for optimizing transmission power which is done by maximizing the number of relay-based paths.In [15] a greedy algorithm is proposed for MHD2D relay selection.The algorithm works on a distributed local search.

    The above-presented literature review on path selection, routing, or relay selection in MHD2D network clearly shows a lack of heuristic approach towards problem-solving, this gap is largely addressed in the proposed research.Literature provides evidence of using an ABC (artificial bee colony) based algorithm for optimization.ABC based scheduling is proposed in [16] where modified ABC is implemented for reducing travel distance and waiting time for a manufacturing workshop.The authors in [17] utilized ABC for the vehicle-to-vehicle emergency information dissemination by using the algorithm minimum route selection from source to destination.In [18]the ABC algorithm is implemented for intrabody nanonetwork where it is used for optimal selection of biosensors along the transmission path.The low energy reliable routing protocol problem is presented as an NP-hard problem in [19] which is solved using an ABC algorithm.

    The literature survey supports the argument that heuristic approaches can be applied for optimum path selection, at the same time application for same in CRE for MHD2D remains unaddressed.

    1.2 Contribution

    The major contributions of the research are:

    ? An MHD2D link-based cell range extension scheme is proposed.

    ? Single antenna transmit power and capacity upper and lower bounds calculations for the network without interference are derived.

    ? Upper and lower bounds of single antenna energy efficiency for a system without interference are calculated are derived.

    ? Transmit power and capacity upper and lower bound calculations for single antenna network with interference are derived.

    ? Energy efficiency upper and lower bound calculations for a network with interference derived.

    ? MVABC proposed for source-destination path selection in a MHD2D based link for CRE.

    ? Comparison between MVABC and genetic algorithm (GA) based path selection for MHD2D link between the source and destination node through simulation over a large data set of 30 runs of 200 iterations each is presented.

    ? Comparison between MVABC and classic ABC based path selection for MHD2D link between the source and destination node through simulation over a large data set of 30 runs of 200 iterations each is presented.

    ? Statistical analysis of simulation results using student’s t-test, p-value, and SEM (standard error of means).

    1.3 Paper Organization

    Section 2 describes the system model and assumptions.In Section 3 we provide detailed derivations for transmitting power, capacity, and energy efficiency for two cases i.e., case one provides derivations without interference while case two provides derivations while considering the effects of interference in a single antenna network.Section 4 gives the mathematical formulation of the problem.In Section 5 we describe MVABC based path selection for MHD2D, which is followed by simulation and results with discussion in Section 6, and the paper is concluded in Section 7.

    2 System Model and Assumptions

    We consider a single input single output system without interference.The mobile users are distributed around the Macro Base Station (MBS) as Homogeneous Poisson Point Process (HPPP)and downlink communication from MBS to a far-off destination node is considered.A ‘far-off’node is a node that lies just outside the cell range of MBS.Since the destination node D is located outside the coverage area of MBS, MHD2D links are used for CRE.It is also assumed that nodes independently decide about their willingness to participate in relaying between source and destination.The decision is based on the incentives offered to mobile relay nodes from the service providers.Fig.1 shows that MBS has several options in terms of finding a path to node D via MHD2D users.

    Figure 1:Proposed D2D multi-hop source-to-destination communication

    For the sake of simplicity, we assume that MBS is already operating at optimized power.Hence, we have to find upper and lower bound for transmitting power of relaying devices and capacity along the path from MBS to node D to optimize an energy-efficient path as described by Eq.(1).It is also assumed that the dissipation in transceiver hardware and baseband processing is negligible.We consider a scenario (refer to Fig.1) where a source node S wants to communicate with node D.Node S sends a connection request (Con.Req.) to MBS.Since MBS cannot directly locate node D in its cell range so it forwards a route request to the users at its cell boundary who then forwards the request to their neighbors and ultimately the request reaches the destination node.In this centralized model of route discovery, the MBS discovers the multi-hop route between node S and node D.In a decentralized case node S will directly broadcast a route request to its neighbors which will then further broadcast it to their neighbor until node D is reached.Both centralized and decentralized modes have their advantages and disadvantages, however, the comparison of the two is beyond the scope of the current research.

    3 Energy Efficiency and Transmit Power Limits

    We will consider two cases wherein the first case we find the transmit power and EE limits for the channel without interference while in the second case we derive the limits for the channel with interference.

    3.1 Single Antenna Transmit Power Without Interference

    The lower bounds of transmit powerPtlbis required to ensure that a required threshold signal to noise ratioγthcan be ensured during transmission.As:

    wherePrthin Eq.(2) is the received power threshold which is the minimum acceptable received power that can ensure a signal to noise ratio (SNR) at least equal to the thresholdγth, soPrthcan be used for derivingPtlbi.e., the lower bound of transmit power to ensureγth, which is given by Eq.(4):

    whereNis the noise power,dis the distance between nodes,dois the reference distance which is taken equal to 1 meter for indoor communication and 100 m to 1 Km for outdoor communication,Lis a constant which is the path loss computed atdoandαis the pathloss exponent.

    For upper bound ofPt(transmit power) we will receive maximum powerPrmaxand the received SNR will beγmaxwhile staying withinPtlbandPtub(transmit power lower and upper bound respectively), so the maximum receiver power is given by Eq.(5):

    As it is obvious that received power increases with a decrease in distance between communicating nodes, hence, we denote the minimum distance asdmin.So, the upper bound of transmit power can be derived as in Eq.(6):

    So, transmit power upper bound is given as in Eq.(7):

    3.2 Single Antenna Capacity and EE Without Interference

    Here we are going to derive the expressions for upper and lower bounds of capacity and EE in a system without interference.Since the capacity is directly proportional to SNR so for the lower bound of capacity the transmit power will also be the lower bound and hence the SNR lower bound ofγthwill be used for the expression of lower bound of capacity as given by Eq.(8):

    Similarly, the upper bound for capacity can be derived usingγmaxas in Eq.(9):

    Now, the upper bound of EE can be written as in Eq.(10):

    Similarly, the lower bound for EE can be written as in Eq.(11):

    3.3 Single Antenna Transmit Power with Interference

    In this case, we will use the Signal to Noise and Interference (SINR) ratio represented byγ′which can be written as Eq.(12):

    whereP′is the interference power received from neighboring nodes.After substituting values of received power as in the case of single antenna transmit power without interference, we can write the expression for lower bounds of transmit power for nodei() and upper bound of transmit power for nodei() as given in Eqs.(13) and (14) respectively:

    3.4 Single Antenna Capacity and EE with Interference

    The expressions for lower and upper bounds of capacity for a network with interference can be derived similarly and is given by Eqs.(15) and (16) respectively:

    Energy efficiency upper and lower bounds can be expressed in terms of the ratio of capacity and transmit power for a network with interference and can be expressed in Eqs.(17) and(18) respectively:

    4 Problem Formulation

    LetV= {v0,v1,...,vn} is the set of nodes,E= {e0,e1,...,en} is the set of edges whereE?{|i,j=0,1,2,...,n}.Path selection problem for maximizing EE can mathematically be represented as:

    Such that:

    Eq.(19) states the energy efficiency for a particular pathpwritten as (EEp) is maximization problem or it’s negative (?EEp) is minimization problem.Eqs.(20) and (21) enforce the capacity constraint on the selected path for the case without interference and the case with interference respectively.Similarly, Eqs.(23)-(25) impose the constraints of transmit power without and with interference on SNR and SINR lower and upper bounds for path selection.

    5 MVABC for MHD2D Path Selection

    The ABC algorithm is a nature-inspired algorithm where the natural behavior of honey bee for rich nectar source exploration is used for real-world problem-solving.The idea was first proposed by Karaboga in 2005 [20] and ever since then the idea has been refined and also it has found numerous problem-solving applications [21-27].Each food source corresponds to the solution for the optimization problem.The ABC algorithm consists of three types of bees based on their roles in the exploration process.

    Scout:All potential food sources are explored by the scout bee.

    Onlooker:Evaluate the nectar quality and quantity as learned from scouts traversing.

    Employed Bees:Responsible for nectar collection from sources identified by the onlooker bees.If a food source is exhausted, the corresponding employed bee becomes a scout and looks for a new food source.

    In MVABC, we take inspiration from the ABC algorithm and propose an algorithm for the implementation of MVABC for multi-hop path selection.Fig.2 shows a flow chart of the pseudocode for the MVABC algorithm for path selection.The proposed algorithms are presented in the Algorithms 1-4.

    Algorithm 1:Initialization 1.For i=Np 2.For j=1 to D 3.M BS generates xi solutions 4.Where xi represents the nodes from SN 5.End 6. failurei=0 7.End

    When node S sends a connection request to MBS, MBS broadcasts to users/devices at the cell edge for node discovery.Through MHD2D relays, node D is discovered and a reply is generated by each intermediate node (which previously received MBS request for node discovery)and sent back to MBS.Based on the received replies, MBS generates an initial population withSNsolutions representing the nodes which sent replies.Each of these D2D relay (food source positions) represent a possible solutionXi=(i=1,2,...,SN)which is aDdimensional vector representing the number of optimization parameters that are taken as two to represent capacity and transmit power.In the algorithmfailureirepresents the abandonment limit of a solutionxiand is equal to the number of unimprovedxi.

    The following sequence takes place during the initialization step:

    ? Fix swarm sizeSNand number of cycles

    ? Determine the number of employed bees (EB), onlooker bees (OB), and food sources, i.e.,

    For path selection problem this will equal the number of nodes.

    ? MBS generates random variables that will serve as the initial population within the domain of decision variables.

    ? MBS then calculates the fitness of the population.

    ? Generate initial trial vector for the population.

    In this phase, MBS sends a route discovery message onNpwhereNprepresents all the possible nodes between S and D, each message sent from MBS behaves similarly to an EB.This route discovery message performs a neighborhood search similar to the EB.Since the node may already have forwarded node discovery request and received a reply, it will compare previous interaction with new interaction, similar to comparing the nectar amount of previous food source in EB memory with the nectar amount of the new food source.The old food location in memory is replaced by a new location if the new location has better nectar quality.The updated location is given by Eq.(27):

    where ?is randomly chosen between [?1,1] andj∈{1,2,...,D}.

    Figure 2:Flow chart of MVABC algorithm

    Algorithm 2:MBS initiates employed bee phase 1.For i=Np 2.Node i receives route request from MBS 3.Node i generates new solution using (27)4.Node i compares fitness using (29)5.If fit(xi+1)>fit(xi)6.newsol=currentsol 7.failurei=0 8.Else 9.newsol=oldsol 10.failurei=failurei+1 11.End

    To generate a new solution, first, we randomly choose a variable to change for the first food source i.e., first neighbor node, then we randomly choose a partner solutionXpand generate a new solution using Eq.(27).The solution is checked for obeying upper and lower bounds and after that fitness of the new solution is calculated using the objective function of Eq.(28):

    If the objective function is within the bound specified in Eqs.(10) and (11) for the case without interference or satisfies Eqs.(17) and (18) for the case with interference, then, the fitness of the population is calculated using Eq.(29):

    A greedy selection is performed between old and new solutions and the fittest of the two solutions is retained.The process continues until all the nodes who sent route request-reply (food sources) have been exhausted.

    images/BZ_621_613_1560_716_1622.pngAlgorithm 3:MBS Initiates Onlooker Bee Phase 1.For i=Np 2.MBS calculates pi using (31)3.End 4. trial=0, i=1 5.Repeat 6.If pi<λ 7.trial=trial+1 8.For j=1toD 9.MBS produces Xjnew using (27)10.End 11.If newpi

    Once the search process is complete and all possible paths have been visited by EB, i.e., MBS has sent route request messages on all possible paths, the D2D relays share their neighborhood information with MBS, this is similar to nectar amount being shared with OB through dance.The probability of a food source is calculated by OB (processing within MBS) using Eq.(30):

    wherefitiis the fitness of solutionifor the EE ratio (nectar amount) it represents, where a comparison is made between old and new EE ratio values and the new one replaces the old position (node) if the new position (node) has a better fitness.

    If MBS does not receive a route request reply from one of the nodes which previously sent a node discovery reply, after a wait ofttime intervals, that particular node will be abandoned.In classical ABC algorithm, in the onlooker bee phase, once the probabilities are calculated for all the possible fitness values, a random number is generated.The randomly generated number is compared with the probability ofith food source (pi).If the probability is less than the random number, then a new solution is generated using Eq.(27), otherwise, the old solution is maintained.In MVABC, we replace this randomly generated number with a number (λ) for making intelligent decision whether new solution should be generated or not.

    Initially the generation ofλdepends on food source generated using lower bounds ofptandC, for these values we find the objective function, fitness and probability using Eqs.(28)-(30)respectively.We label this particular value of probability for lower bounds of transmit power and capacity asλ.As shown in algorithm 3, ifpi<λ, i.e., probability of ith food source is less than a predefined threshold for acceptable ranges ofptandC, then a new solution is generated using Eq.(27), otherwise the old solution is retained and probability of next food source is compared.In the simulation we have chosen three difference values ofλnamely 0.2, 0.8 and 1 for the sake of demonstration of results for its different values.For simulationλmay be assigned any value within in acceptable bounds ofptandC.This value can be fixed to generated most optimum results.

    Algorithm 4:MBS Initiates Scout Bee Phase 1.If trial>maxt rial 2.MBS generates new solution using (31)

    The SB phase is entered when any solution having a trial limit greater than a maximum exceeds.When MBS cannot improve an established path any further through route discovers, then,the previous best solution is stored in the memory and a new solution is generated randomly between upper and lower bounds using Eq.(31):

    whereXminandXmaxare lower and upper bound of population respectively.

    The objective function and fitness of the new solution is calculated using Eqs.(28) and(29) respectively and it replaces the old solution even if it has poor fitness compared to the old solution.

    6 Simulation and Results

    The comparison and performance analysis of the proposed algorithm was carried out in Matlab R2018b.The simulation parameters are presented in Tab.1.

    6.1 Performance Metrices

    A statistical approach similar to one utilized in [28] was used to evaluate the results.We will briefly describe the performance metrices in this sub-section and their details are given in the sub-section below along with discussion on results.

    Table 1:Simulation parameters

    Best Cost(Mean abs(-EE)):As described in Section 1, in this research we consider maximizing the EE or minimizing the (-EE) and is presented for proposed algorithm MVABC along with classic ABC and GA.

    Standard Deviation(SD):Describes the variation in mean data values.

    Standard Error of Means(SEM):Measures deviation of sample’s mean from population mean.

    6.2 Results

    The proposed algorithm was compared with classic ABC and GA based path selection.Both algorithms were tested using large data set with 30 and 60 runs with 200 iterations per run.The population size and cycle limit were kept the same for all the simulated cases of MVABC, GA and classic ABC.

    For GA based algorithm we kept the crossover percentage at 0.7, mutation percentage at 0.3,and mutation rate at 0.1.A population size of 100 was used for all three algorithms.We calculated the mean best cost, SD and the SEM for each run.

    We compared GA, classic ABC and MVABC based path selection through simulation over a large data set with randomly generated seed values.The mean best, standard deviation, and SEM are given in Tab.2 for 30 runs of each algorithm.

    To show whether the experimental results obtained from MVABC, classic ABC and GA for path selection are significantly different, we perform Student’s t-test hypothesis [28].The t-test statistically verifies if the difference in the mean values is due to the variation in sampling or does it present evidence that MVABC provides better convergence to the solution.

    The t-value for the hypothesis is calculated using Eq.(32):

    where grand standard deviation,SX1X2is given by Eq.(33):

    Table 2:Statistical results from 30 runs of MVABC, classic ABC and GA

    Table 2(continued).

    whereis mean of data for MVABC,is mean of data for classic ABC for obtaining results as in Tab.3 and in subsequent calculations it represents mean of GA where results are presented in Tab.4.SX1is the standard deviation of data for MVABC, andSX2is the standard deviation of data for classic ABC and represent the same for GA in subsequent calculations as given in Tab.4.

    Table 3:Summary of statistical T-test results for MVABC and classic ABC

    The t-test will test the null hypothesis to check whether the mean of MVABC, classic ABC and GA are equalvs.an alternate hypothesis that the mean of MVABC, classic ABC and GA are not equal and the MVABC test results hold more significance.The result of the t-test performed on the data set labeled as Group1 (MVABC) and Group 2 (classic ABC) is given in Tab.3.

    Table 4:Summary of statistical T-test results for MVABC and GA

    whereNis the total number of runs of each algorithm,Tis the t-test value (who’s absolute is compared to the critical value), d.o.f (degree of freedom) refers to maximum logically independent values that can vary with freedom in a sample of data and is given byd.o.f=2n?2, wherenis the total number of iterations.Its value is important because it helps in determining the critical value against which a hypothesis is tested.The test shows that the means of MVABC and classic ABC are significantly different atp<0.05 significance level.Wherepis the probability of having a value less than the significance level indicates that the null hypothesis is rejected and the alternate hypothesis is accepted, in other words, we can say that the mean path cost calculated by MVABC is less than that using classic ABC.This implies that the nodes selected between source to destination using MVABC represented a more energy efficient path.

    Similar t-tests were performed between MVABC and GA.The statistical results are given in Tab.4.As in the case of comparison with classic ABC, MVABC compared to GA is again capable of providing a more energy efficient path.

    This can further be realized by examining the error bar graph in Fig.3.MVABC was tested for different values ofλ.Since the error bars are not overlapping (Fig.3) hence it can be seen that there is a significant difference in the results generated by MVABC and those by classic ABC and GA, so, the null hypothesis is rejected and the alternative hypothesis is accepted.Moreover,the difference in means calculated using the MVABC, classic ABC and GA can also be observed for various values ofλ, the results can be seen in Fig.4.The results represent better convergence of MVABC and better overall best cost value compared to classical ABC.In Fig.5a comparison is drawn between different values ofλapplied in MVABC with classic ABC.Various values ofλfor MVABC clearly offer least cost solution compared to the classic ABC.

    Figure 3:SEM comparison between MVABC, Classic ABC and GA

    Figure 4:Comparison of means between MVABC, Classic ABC and GA

    Figure 5:Comparison of different values of λ in MVABC with classic ABC

    7 Conclusion

    This paper proposed an MHD2D link for CRE or offloading of the base station.An MVABC algorithm was proposed for path selection between source and destination for an MHD2D link.Upper and lower bounds of transmit power, capacity, and energy efficiency for the network with and without interference were derived.The proposed MVABC approach for MHD2D path selection was compared with classical ABC and GA-based MHD2D path selection for a large data set with the same population size and same upper and lower bounds for population generation.The simulated results signified the improved performance of MVABC over classic ABC and GA based approach.The simulated results were further verified using statistical techniques namely the student’s t-test, p-value, and SEM.The detailed statistical analysis confirmed that MVABC outperformed GA based MHD2D based path selection.

    As future work, we will consider multi-antenna system with both centralized and decentralized approaches for path finding.We also plan on improving MVABC further in this respect and develop a generalized algorithm for single and multi-antenna systems.

    Funding Statement:The authors received no specific funding for this study.

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

    国产亚洲一区二区精品| 纵有疾风起免费观看全集完整版| 热99re8久久精品国产| 人人澡人人妻人| 欧美精品一区二区大全| 久久免费观看电影| 两性夫妻黄色片| 99热网站在线观看| 麻豆国产av国片精品| 久久午夜亚洲精品久久| 久久性视频一级片| 久久香蕉激情| 一区二区av电影网| 777米奇影视久久| 久久人妻av系列| 亚洲精品av麻豆狂野| 亚洲成人手机| 精品国产国语对白av| 午夜激情av网站| 韩国精品一区二区三区| 老司机影院毛片| 最黄视频免费看| 在线播放国产精品三级| 欧美成人午夜精品| 免费在线观看完整版高清| 亚洲精品自拍成人| 国产片内射在线| 岛国在线观看网站| 考比视频在线观看| 久久精品亚洲av国产电影网| 黄色怎么调成土黄色| 亚洲精品在线美女| 美女高潮到喷水免费观看| 亚洲九九香蕉| 91国产中文字幕| 法律面前人人平等表现在哪些方面| 两个人看的免费小视频| 久久久久久久久免费视频了| 多毛熟女@视频| 中文字幕色久视频| 国产精品久久电影中文字幕 | 一夜夜www| 成人手机av| 国产单亲对白刺激| 精品一区二区三区av网在线观看 | 亚洲成人免费av在线播放| 久久av网站| 久久精品91无色码中文字幕| 亚洲avbb在线观看| 亚洲欧美激情在线| 亚洲精品美女久久av网站| 波多野结衣一区麻豆| 午夜免费鲁丝| 韩国精品一区二区三区| 久久国产精品大桥未久av| 久久久精品国产亚洲av高清涩受| 老司机深夜福利视频在线观看| 色在线成人网| 久久精品国产亚洲av高清一级| 狂野欧美激情性xxxx| 一本大道久久a久久精品| 亚洲成人免费av在线播放| 欧美日韩亚洲综合一区二区三区_| 性少妇av在线| 97人妻天天添夜夜摸| 久热这里只有精品99| 亚洲人成77777在线视频| 国产免费av片在线观看野外av| 丁香六月欧美| 日日摸夜夜添夜夜添小说| 亚洲午夜精品一区,二区,三区| 久久久久久久国产电影| 99re在线观看精品视频| 精品人妻在线不人妻| 国产精品香港三级国产av潘金莲| 一本一本久久a久久精品综合妖精| 亚洲av电影在线进入| 汤姆久久久久久久影院中文字幕| 人人妻人人澡人人看| 9色porny在线观看| 欧美人与性动交α欧美软件| 不卡一级毛片| av在线播放免费不卡| 亚洲精品粉嫩美女一区| 亚洲国产欧美一区二区综合| 桃花免费在线播放| 国产av国产精品国产| 狠狠精品人妻久久久久久综合| 欧美日韩中文字幕国产精品一区二区三区 | 丝袜人妻中文字幕| 久久 成人 亚洲| 丝瓜视频免费看黄片| 国产精品自产拍在线观看55亚洲 | 亚洲va日本ⅴa欧美va伊人久久| 在线观看一区二区三区激情| 十八禁人妻一区二区| 久久热在线av| 免费在线观看影片大全网站| 免费观看av网站的网址| av天堂久久9| 飞空精品影院首页| 国产精品av久久久久免费| 久久久久国产一级毛片高清牌| 曰老女人黄片| 在线天堂中文资源库| 三级毛片av免费| av电影中文网址| avwww免费| 久久精品91无色码中文字幕| 少妇猛男粗大的猛烈进出视频| 日韩欧美国产一区二区入口| 啪啪无遮挡十八禁网站| 9色porny在线观看| 精品亚洲成a人片在线观看| 精品免费久久久久久久清纯 | 久久国产精品影院| 日韩熟女老妇一区二区性免费视频| aaaaa片日本免费| av不卡在线播放| 欧美黑人精品巨大| 欧美成狂野欧美在线观看| 淫妇啪啪啪对白视频| 久久久久久久大尺度免费视频| 精品久久久久久久毛片微露脸| 中文字幕色久视频| 亚洲精品国产区一区二| 国产亚洲精品久久久久5区| 午夜久久久在线观看| 国产日韩欧美视频二区| 黑丝袜美女国产一区| av又黄又爽大尺度在线免费看| 亚洲精品一卡2卡三卡4卡5卡| 欧美黄色片欧美黄色片| 亚洲精品乱久久久久久| 精品国产超薄肉色丝袜足j| 国产精品.久久久| 十分钟在线观看高清视频www| 下体分泌物呈黄色| 日本黄色日本黄色录像| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品影院| 99精国产麻豆久久婷婷| 搡老乐熟女国产| 久久久久久久大尺度免费视频| 日韩欧美免费精品| 我的亚洲天堂| 啦啦啦免费观看视频1| 午夜福利在线观看吧| 天堂8中文在线网| av线在线观看网站| 在线看a的网站| 国产日韩一区二区三区精品不卡| 国产一区有黄有色的免费视频| 一边摸一边抽搐一进一小说 | 亚洲 国产 在线| 狠狠婷婷综合久久久久久88av| 国产亚洲午夜精品一区二区久久| 久久人人97超碰香蕉20202| 狠狠婷婷综合久久久久久88av| 亚洲精品国产色婷婷电影| 亚洲视频免费观看视频| 99riav亚洲国产免费| 午夜激情久久久久久久| h视频一区二区三区| 国产av一区二区精品久久| 国产精品九九99| 国产精品成人在线| 正在播放国产对白刺激| 老熟妇乱子伦视频在线观看| 丁香六月天网| 露出奶头的视频| 国产欧美日韩综合在线一区二区| 黑人欧美特级aaaaaa片| 国产精品99久久99久久久不卡| 欧美精品高潮呻吟av久久| 日韩欧美免费精品| 正在播放国产对白刺激| 国产成人系列免费观看| 嫁个100分男人电影在线观看| 男男h啪啪无遮挡| 国产精品香港三级国产av潘金莲| 亚洲一卡2卡3卡4卡5卡精品中文| av免费在线观看网站| 欧美日本中文国产一区发布| e午夜精品久久久久久久| 精品少妇内射三级| 久久国产精品影院| 国产在线观看jvid| 午夜福利免费观看在线| 视频在线观看一区二区三区| 国产单亲对白刺激| 日韩免费高清中文字幕av| 新久久久久国产一级毛片| 国产一区二区三区在线臀色熟女 | 老熟妇仑乱视频hdxx| 99在线人妻在线中文字幕 | 日本wwww免费看| 亚洲成人手机| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品在线电影| 最近最新中文字幕大全电影3 | 午夜福利视频在线观看免费| 免费在线观看黄色视频的| 美女国产高潮福利片在线看| 美女高潮到喷水免费观看| 国产精品久久电影中文字幕 | 国产不卡一卡二| 麻豆国产av国片精品| 免费日韩欧美在线观看| 国产成人精品在线电影| 91老司机精品| 午夜福利一区二区在线看| 18禁黄网站禁片午夜丰满| 国产精品一区二区精品视频观看| 叶爱在线成人免费视频播放| 免费一级毛片在线播放高清视频 | 亚洲精品av麻豆狂野| av在线播放免费不卡| 精品国产乱码久久久久久小说| 欧美日韩av久久| 国产精品.久久久| www.熟女人妻精品国产| 99re在线观看精品视频| 国产97色在线日韩免费| 天天躁夜夜躁狠狠躁躁| 成人手机av| 在线 av 中文字幕| 香蕉丝袜av| 我的亚洲天堂| 1024香蕉在线观看| 国产亚洲午夜精品一区二区久久| 日韩成人在线观看一区二区三区| 91国产中文字幕| 五月天丁香电影| 国产在线免费精品| 丰满迷人的少妇在线观看| 国产亚洲精品第一综合不卡| 国产精品美女特级片免费视频播放器 | 久久国产精品大桥未久av| 久久午夜综合久久蜜桃| 人妻 亚洲 视频| 极品人妻少妇av视频| 午夜两性在线视频| 国产精品 国内视频| 在线观看舔阴道视频| 人人妻人人爽人人添夜夜欢视频| 夫妻午夜视频| 国产真人三级小视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 激情视频va一区二区三区| 亚洲av日韩在线播放| 国精品久久久久久国模美| 久久久精品免费免费高清| 两人在一起打扑克的视频| 国产一区二区三区视频了| 欧美精品人与动牲交sv欧美| 超碰成人久久| 国产有黄有色有爽视频| 国产成人精品久久二区二区免费| 亚洲熟妇熟女久久| 天堂8中文在线网| 亚洲精品久久成人aⅴ小说| 久久这里只有精品19| 丰满饥渴人妻一区二区三| 欧美在线一区亚洲| 欧美 亚洲 国产 日韩一| 国产精品.久久久| 两个人看的免费小视频| 国产一卡二卡三卡精品| 国产片内射在线| 天天影视国产精品| 国产区一区二久久| 99国产极品粉嫩在线观看| 亚洲熟女精品中文字幕| 丰满饥渴人妻一区二区三| 亚洲av欧美aⅴ国产| 午夜福利影视在线免费观看| 亚洲性夜色夜夜综合| 99精品在免费线老司机午夜| 久久久久精品人妻al黑| 我的亚洲天堂| 国产精品二区激情视频| 亚洲欧洲日产国产| 美女午夜性视频免费| 国产一区二区在线观看av| 另类精品久久| 久久性视频一级片| 国产亚洲一区二区精品| 一边摸一边抽搐一进一小说 | 日韩欧美一区视频在线观看| 熟女少妇亚洲综合色aaa.| 男女午夜视频在线观看| 99热网站在线观看| 极品少妇高潮喷水抽搐| 菩萨蛮人人尽说江南好唐韦庄| 超碰97精品在线观看| a级毛片在线看网站| 国产精品久久久久成人av| 久久久水蜜桃国产精品网| 久久影院123| 一区二区日韩欧美中文字幕| 午夜成年电影在线免费观看| 一区二区日韩欧美中文字幕| 大片电影免费在线观看免费| 免费少妇av软件| 夜夜夜夜夜久久久久| 国产一卡二卡三卡精品| 午夜老司机福利片| 高清欧美精品videossex| 欧美av亚洲av综合av国产av| 狠狠婷婷综合久久久久久88av| 少妇的丰满在线观看| 亚洲欧美精品综合一区二区三区| 一本一本久久a久久精品综合妖精| 日本wwww免费看| 激情在线观看视频在线高清 | 黄色视频不卡| 中文字幕色久视频| 一级黄色大片毛片| 亚洲av日韩精品久久久久久密| 亚洲第一欧美日韩一区二区三区 | 精品少妇黑人巨大在线播放| 成人国产av品久久久| 国产成人影院久久av| 777久久人妻少妇嫩草av网站| 狠狠狠狠99中文字幕| 久久精品国产综合久久久| av有码第一页| 精品卡一卡二卡四卡免费| 国产亚洲欧美在线一区二区| 国产1区2区3区精品| 老汉色av国产亚洲站长工具| 怎么达到女性高潮| 欧美激情极品国产一区二区三区| 久久国产亚洲av麻豆专区| 成年人免费黄色播放视频| 人人妻人人添人人爽欧美一区卜| 日本黄色日本黄色录像| 黄色毛片三级朝国网站| 90打野战视频偷拍视频| 91九色精品人成在线观看| 女性生殖器流出的白浆| 国产成人精品久久二区二区免费| 黑丝袜美女国产一区| 色播在线永久视频| 俄罗斯特黄特色一大片| 成人永久免费在线观看视频 | 侵犯人妻中文字幕一二三四区| 啪啪无遮挡十八禁网站| 久热爱精品视频在线9| 交换朋友夫妻互换小说| 国产精品98久久久久久宅男小说| av超薄肉色丝袜交足视频| 精品福利观看| 一边摸一边抽搐一进一小说 | 精品乱码久久久久久99久播| 天堂8中文在线网| 99精品欧美一区二区三区四区| 久久热在线av| 99在线人妻在线中文字幕 | 91字幕亚洲| 国产极品粉嫩免费观看在线| 自线自在国产av| 99热国产这里只有精品6| 少妇粗大呻吟视频| 人人妻人人澡人人看| 自线自在国产av| 国产成人一区二区三区免费视频网站| 免费看a级黄色片| 视频区图区小说| 欧美+亚洲+日韩+国产| 一区二区av电影网| 无遮挡黄片免费观看| e午夜精品久久久久久久| 女人高潮潮喷娇喘18禁视频| 大码成人一级视频| 人妻 亚洲 视频| 黄色丝袜av网址大全| 国产欧美日韩一区二区三| 91大片在线观看| 久久中文看片网| 9色porny在线观看| 国产在线免费精品| 天天躁日日躁夜夜躁夜夜| 精品国产乱码久久久久久男人| 又黄又粗又硬又大视频| 国产成人精品久久二区二区91| 一边摸一边抽搐一进一出视频| 啪啪无遮挡十八禁网站| 亚洲三区欧美一区| 欧美日韩一级在线毛片| 人人妻,人人澡人人爽秒播| 免费不卡黄色视频| 日韩 欧美 亚洲 中文字幕| 久久久国产一区二区| 国产又色又爽无遮挡免费看| 男女高潮啪啪啪动态图| 高清视频免费观看一区二区| 午夜福利视频精品| 精品国产一区二区三区久久久樱花| 亚洲欧美一区二区三区久久| 久久精品国产亚洲av高清一级| 手机成人av网站| 男女边摸边吃奶| 国产精品亚洲一级av第二区| 欧美成狂野欧美在线观看| 12—13女人毛片做爰片一| 精品国产乱码久久久久久小说| 嫩草影视91久久| 最新在线观看一区二区三区| 精品一区二区三卡| av天堂久久9| 99国产精品免费福利视频| 啦啦啦在线免费观看视频4| 91成年电影在线观看| 日韩免费高清中文字幕av| 99精国产麻豆久久婷婷| 久久国产亚洲av麻豆专区| 午夜福利一区二区在线看| 久久亚洲真实| 看免费av毛片| 国产高清激情床上av| 免费不卡黄色视频| 蜜桃国产av成人99| 久久亚洲真实| 高清欧美精品videossex| 国产精品.久久久| 91老司机精品| 两个人免费观看高清视频| 成人国产av品久久久| 男女免费视频国产| 精品卡一卡二卡四卡免费| 亚洲成av片中文字幕在线观看| 国产精品 欧美亚洲| 日韩成人在线观看一区二区三区| 99久久国产精品久久久| 纯流量卡能插随身wifi吗| 精品国产一区二区久久| 精品国产一区二区三区四区第35| 国产日韩一区二区三区精品不卡| 国产亚洲欧美在线一区二区| 亚洲va日本ⅴa欧美va伊人久久| 日韩人妻精品一区2区三区| 亚洲熟女毛片儿| 最黄视频免费看| 国产亚洲精品第一综合不卡| 多毛熟女@视频| 高清视频免费观看一区二区| 黑人操中国人逼视频| 高清黄色对白视频在线免费看| 变态另类成人亚洲欧美熟女 | 国产野战对白在线观看| 多毛熟女@视频| 99国产精品免费福利视频| 色综合婷婷激情| 男女之事视频高清在线观看| 日韩一区二区三区影片| 性少妇av在线| 国产成人系列免费观看| 制服人妻中文乱码| 女性生殖器流出的白浆| www日本在线高清视频| 又大又爽又粗| 天天躁夜夜躁狠狠躁躁| 成人三级做爰电影| 亚洲男人天堂网一区| 真人做人爱边吃奶动态| 狠狠狠狠99中文字幕| 亚洲熟女毛片儿| 黄色片一级片一级黄色片| 亚洲国产毛片av蜜桃av| 久久影院123| 久久精品人人爽人人爽视色| 曰老女人黄片| 久9热在线精品视频| 99国产极品粉嫩在线观看| 欧美日韩亚洲国产一区二区在线观看 | 丝袜喷水一区| 91大片在线观看| 国产免费现黄频在线看| 曰老女人黄片| 中文字幕人妻熟女乱码| 久久国产精品男人的天堂亚洲| 久久人妻av系列| 露出奶头的视频| 国产麻豆69| 日韩视频在线欧美| 50天的宝宝边吃奶边哭怎么回事| 国产在线免费精品| 久久 成人 亚洲| 成人免费观看视频高清| 日日夜夜操网爽| 蜜桃在线观看..| 久久久精品94久久精品| 啦啦啦免费观看视频1| 国产成人欧美在线观看 | 一个人免费看片子| 夜夜夜夜夜久久久久| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色 | 窝窝影院91人妻| 精品少妇黑人巨大在线播放| 人人澡人人妻人| 精品少妇黑人巨大在线播放| 成人18禁高潮啪啪吃奶动态图| 蜜桃在线观看..| 丰满人妻熟妇乱又伦精品不卡| 国产又色又爽无遮挡免费看| 91老司机精品| 免费观看av网站的网址| 在线观看www视频免费| 精品国产乱码久久久久久小说| 美女午夜性视频免费| 亚洲欧美日韩高清在线视频 | 女同久久另类99精品国产91| 精品一区二区三区四区五区乱码| 少妇被粗大的猛进出69影院| 亚洲国产欧美网| 午夜福利在线观看吧| 日韩制服丝袜自拍偷拍| 在线播放国产精品三级| 国产一区二区三区在线臀色熟女 | 91麻豆av在线| 国产欧美日韩精品亚洲av| 中文字幕高清在线视频| 极品教师在线免费播放| 欧美性长视频在线观看| 麻豆国产av国片精品| 久久国产精品男人的天堂亚洲| 黄色怎么调成土黄色| 男女免费视频国产| 纯流量卡能插随身wifi吗| 国产精品免费大片| 亚洲av日韩精品久久久久久密| 水蜜桃什么品种好| 亚洲人成伊人成综合网2020| 国产免费av片在线观看野外av| 午夜日韩欧美国产| 十八禁人妻一区二区| 国产亚洲精品第一综合不卡| 黄片播放在线免费| 精品免费久久久久久久清纯 | 五月开心婷婷网| 亚洲va日本ⅴa欧美va伊人久久| 男人操女人黄网站| 精品视频人人做人人爽| 男男h啪啪无遮挡| 视频在线观看一区二区三区| videos熟女内射| 久久久久久久久免费视频了| 狠狠狠狠99中文字幕| 午夜视频精品福利| 午夜久久久在线观看| 91成年电影在线观看| 国产日韩欧美视频二区| 成人国语在线视频| 天堂动漫精品| 手机成人av网站| 日本wwww免费看| 宅男免费午夜| 欧美人与性动交α欧美软件| 久久久国产成人免费| 在线观看免费视频网站a站| 精品高清国产在线一区| 又大又爽又粗| 国产在线观看jvid| 国产成人免费无遮挡视频| 亚洲精品自拍成人| 一个人免费看片子| 亚洲色图综合在线观看| 男女边摸边吃奶| 在线观看人妻少妇| 亚洲 欧美一区二区三区| 母亲3免费完整高清在线观看| 黑人巨大精品欧美一区二区蜜桃| 成人18禁高潮啪啪吃奶动态图| 久久人妻熟女aⅴ| 亚洲国产成人一精品久久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲九九香蕉| 老司机靠b影院| 国产成人啪精品午夜网站| 国产男女超爽视频在线观看| 丰满饥渴人妻一区二区三| 午夜激情久久久久久久| 人人妻人人澡人人爽人人夜夜| 一区二区三区激情视频| 亚洲第一欧美日韩一区二区三区 | 91老司机精品| 欧美精品一区二区大全| 人人妻,人人澡人人爽秒播| 国产精品一区二区免费欧美| 天天操日日干夜夜撸| av天堂在线播放| xxxhd国产人妻xxx| 麻豆av在线久日| 国产一区有黄有色的免费视频| 啦啦啦在线免费观看视频4| 久久精品成人免费网站| 少妇精品久久久久久久| 亚洲 欧美一区二区三区| 亚洲av国产av综合av卡| 亚洲第一av免费看| 亚洲av片天天在线观看| 久久九九热精品免费| 精品视频人人做人人爽| 欧美乱码精品一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久欧美国产精品| 老司机深夜福利视频在线观看| 日本一区二区免费在线视频| 久久中文看片网| 手机成人av网站| 欧美黄色淫秽网站| 国产精品.久久久|