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

    An improved adaptive differential evolution algorithm for single unmanned aerial vehicle multitasking

    2022-01-05 09:41:04JianliSuHuaWang
    Defence Technology 2021年6期

    Jian-li Su, Hua Wang

    School of Astronautics, Beihang University, Beijing,100191, China

    Keywords:Unmanned aerial vehicle Multitasking Adaptive differential evolution Mutation factor Crossover factor

    ABSTRACT

    1. Introduction

    UAVs are increasingly widely used in the past decade, and its technology becomes more mature. Perform a single task can no longer meet the demand for UAV,particularly in the field of military[1]. The use of UAV has made warfare more efficient and has resulted casualties [2,3]. In consequence, research on multitasking of UAV has aroused worldwide attention.

    Although multi-UAVs cooperative multi-targets assignment is a wide-ranging topic, research on multitasking of single UAV is particularly important, such as the RQ-4, MQ-1 and MQ-9 are expensive to produce, and the high cost of their operations and maintenance, allowing large scale cooperation of these vehicles is not possible [4]. Single UAV multi-tasking assignment, like multi-UAVs cooperative multi-targets assignment, proved to be a NPhard problem that getting a global optimized solution [5]. Many intelligent algorithms have been put forward and improved by scholars to address the computing complexity, for example,mathematical algorithms, genetic algorithm (GA) [6-9], particle swarm optimization (PSO) [10-12], and ant colony optimization(ACO) [13-15]. GA, which is formed by the biological laws and natural genetic mechanisms search algorithm, works as a parallel algorithm. Professor J. Holland said GA can be mainly realized by the means of crossover, mutation, and selection operators in his proposed book in 1975 [16]. In the PSO algorithm, which is proposed by Kennedy and Eberhart [17], several particles working as potential solution of the search space are randomly generated in the first stage. Secondly, the iteration process can help to find the optimal solution.Thirdly,two“extreme values”are used to update the velocity and displacement of the particle,one of which is called the local optimal solution working for the particle itself as the optimal solution and another one is called the global optimal solution for the entire population. In Ref. [18], the spherical vectorbased SPSO algorithm was proposed to solve the UAV path planning problem and minimize the cost function. In Ref. [19], the proposed motion-encoded PSO algorithm is applied to the UAV's search for moving targets. In Ref. [20], Rauch-Tung-Striebel (RTS)smoother is used to improve the PSO algorithm,and the algorithm is applied to UAV path planning,and good results are obtained.The ACO imitates the process of ant foraging. At first, the ant colony randomly searches for food and finally gathers at the place where there is more food. The ant colony foraging will show a certain pattern. The algorithm that imitates this process is called ACO. In order to determine the trajectory of the UAV so that it can find the missing target in the shortest time,Sara Perez-Carabaza proposes a new method based on ACO[21].

    Those heuristic algorithms are put to be widely used in UAV multi-target attacking and path planning as digital technology and computer technology improvement, each one of which has its advantages and disadvantages. Although GA has good global searching ability, but its local searching turned out not very effective.During the early stage of PSO with the fast convergence speed,the low efficiency results in poor accuracy and divergence. ACO has many advantages, however, the complicated calculations cause long processing time and converge slowly,and easily fall into local optimal solution.

    In recent years, differential evolution (DE) [22] is widely concentrated and applied, especially for multi-targets assignment and path planning. Zhao and Huang present an improved DE algorithm,aiming at the problem that UAV cooperative multi-targets assignment in 3D environments which makes it difficult to model and calculate [23,24]. An approximate clustering dual-strategy DE algorithm is proposed to improve the efficiency of obtaining the optimal solution.An improved discrete DE algorithm is imported in the interval-valued intuition fuzzy decision making, aiming at multiple UAV task assignment [25]. Wang establishes a weapontarget assignment (WTA) model and proposes a self-adaptive DE algorithm to get better parameters of the WTA [26]. In addition,there are many studies about the UAV route planning based on DE algorithm. It is generally accepted that the collaborative path planning for UAV as inequality constrained optimization problems.A series of improved DE algorithms have been proposed to solve the UAV path planning problem, especially in 3D environments[27-29]. An improved DE algorithm adopts B-Spline curve is presented to obtain a 3D online optimum route[30].Zhang develops a DE algorithm with the exponential selection to get better planning route[31].Aiming at the UAVs dynamic path planning,an improved DE algorithm is adopted[32].Lei presents an effective DE algorithm with hybrid strategy to obtain an optimal path [33].

    Same as other heuristic algorithms,there are some drawbacks to DE algorithm, such as early convergence and stagnation. Hence,many scholars propose some improvements, especially in the combination of DE algorithm and adaptive mechanisms. Among them,JADE has been accepted by most researchers[34].In order to make the parameters of the objective function can be obtained quickly and accurately, Li proposes an enhanced adaptive DE algorithm [35]. In Ref. [36], a fitness-based adaptive DE algorithm is proposed, then the algorithm is optimized by utilizing three adaptive strategies. Tan proposes an improved adaptive DE algorithm with different strategies for mutation factor and crossover factor [37]. Shen presents an adaptive differential evolutionary strategy based on double mutation strategies [38]. Tao [39] introduces an adaptive differential evolutionary strategy using an improved elitist strategy. Meng proposes a novel parameter adaptive DE (PaDE) to solve the shortcomings of many DE variants[40,41].Based on the idea of greedy search,Miguel Leon proposes a new greedy adaptation method to dynamically adjust the mutation factor and crossover rate in DE [42]. M Leon combines RAM and JAPDE to propose RAM-JAPDE algorithm,which can simultaneously adjust the selection probability of control parameter pair and mutation strategy in DE [43]. In Ref. [44], it is pointed out that the performance of DE depends on the selection of mutation factor and crossover probability. As on this, a method of setting the control parameters of the differential evolution method is proposed.

    In this paper,aiming at the difficulties of single UAV multi-task execution in a three-dimensional environment,and answering how to make UAVs reasonably allocate the execution order of targets in UAV missions and improve efficiency under complex threedimensional terrain environments and multiple constraints, a new adaptive differential evolution algorithm is proposed. In the solution process, the three-dimensional voyage cost is used to express the allocation relationship between the UAV and the target,and the dynamic adaptive mutation and cross-differential evolution strategy are combined to adjust the balance of detection and development in algorithm search to avoid falling into local optimality.The accuracy of task execution is improved,and it is suitable for solving larger-scale target execution tasks. Firstly, a single UAV multitasking model and its related constraints are established in a three-dimensional environment. Secondly, an adaptive evolutionary algorithm based on population iteration times and individual fitness values to adjust the mutation factor and crossover factor is proposed to solve the problem of exploratory and exploitive balance existing in the current evolutionary process.Thirdly, through different scales of simulation experiments, the rationality of the IADE algorithm is proved.Finally,the comparative experiments with other algorithms show the superiority of the IADE algorithm.

    The rest of this article is organized as follows. Section 2 mainly explains the UAV multitasking model, various of constraints and corresponding mathematical formulations. Section 3 introduces the IADE algorithm. In Section 4, the simulation experiments are carried out and compared with other algorithms. Section 4 summarizes the whole article and proposes the work to be done in the future.

    2. UAV multi-targets assignment model

    2.1. Objective function

    UAV multitasking is a constrained optimization problem (COP).As usual,this kind of problem contains both equality and inequality constraints. According to the optimization theory, the mathematic description of COP is as follows:

    where x?D represents the decision variable of optimization problem,f(x)is the objective function,which depicts targets to be optimized, hi(x) and gj(x) respectively expresses the equality constraint and inequality constraint.

    From the description above,connecting with specific problem of UAV multitasking, assumption can be made that a single UAV attacks multiple targets or executes multi-tasks in different areas of the mission space.The key to solve this problem is to obtain attack and execution sequence, making the UAV have the shortest flight distance and minimum flight time that satisfy all constraints.

    The objective function of UAV multitasking model can be described as follows:

    where djis the flight distance from UAV to target j, wjis the objective weight of j, and its range is between 0 and 1, and the larger the value, the higher the priority, tjdenotes the UAV flight time to target,ckrepresents the constraint violation,α and β denote the scaling factors which are used to balance the polynomial of the objective function, keeping the voyage cost, flying time and constraint violation in the same order of magnitude. Djrepresents the decision variable,which determines the target to be performed by the UAV.As shown in Fig.1,the 3D environment and application scenarios are displayed.

    Fig.1. Application scenario in three-dimensional environment.

    2.2. Constraint relation

    (1) Maximum flying distance constraint Maximum flying distance constraint represents the performance and fuel consumption index of UAV, which is the UAV maximum range limit in the process of completing all missions.The key to this constraint is to ensure that the UAV can complete all targets, otherwise, it will be penalized. The constraint function is written as:

    where djis the estimated cost of voyage, D represents the UAV maximum range limit.C1violdenotes the constraint violation that it will be penalized if the UAV exceeds the maximum flying distance while performing the j-th task.

    (2) Airspeed constraint

    where V(j)minand V(j)maxrespectively represents the maximum speed and minimum speed of UAV; C2violrepresents that the UAV will be punished if it flies beyond the speed range of the specified minimum and maximum speed.

    (3) Maximum flying time constraint

    Maximum flying time constraint represents the total time elapsed for the UAV to complete all missions. The key to this constraint is to ensure that the UAV completes all missions within the specified flight time.

    where C3violrepresents that the UAV will be punished if it surpasses the maximum flying time.

    (4) Timing constraint between target points

    Timing constraint reflects the execution order among the target points.Priority to important target points must take precedence on execution, thus other target points will be performed in order according to their constraints.

    where Tj

    (5) Maximum number of tasks to execute constraint.

    The maximum number of missions performed by the UAV is represented by this constraint

    where Mj

    3. Improved adaptive differential evolution

    3.1. Standard differential evolution

    The DE algorithm first initializes the population,and then carry out mutation, crossover and selection operations to continuously iterate and update the population size. The steps of DE algorithm are as follows:

    (1) Initial population

    DE algorithm randomly generates NPindividuals within the dimension D of the search range of the solution according to

    (2) Mutation

    The DE algorithm often use some fixed mutation strategies to generate mutant individuals vito ensure the diversity of the population. The deterministic mutation strategies in DE algorithm are shown below:

    Where xbestdenotes the best individual,F is the scaling factor r1r2r3r4, r5represent distinct random integers on [1, NP] and are not equal to i.

    (3) Crossover

    The basic principle of crossover operation is to generate experimental vector after parameter hybrid of mutation vector and predetermined vector. The specific expression is as follows:

    where CR is the crossover probability; jrand represents a random integer on [1,D].

    (4) Selection

    The main idea of the DE algorithm is to select the next generation of individuals according to the law of greed,which means the excellent individuals are retained and the inferior ones are eliminated.For the optimization problem,the excellent individual has a smaller fitness value.

    3.2. JADE

    DE/current-to-pbest,as a new mutation strategy,forms the base of JADE which is an adaptive DE algorithm and is with optional archive.

    where randc is the Cauchy distribution,if Fi>1,let Fibe equal to1,or regenerate Fiif Fi≤0, initialize uFto 0.5 and updated by as follows:

    where c is the value within interval (0,1), SFdenotes the set of all excellent mutation factors.

    In JADE, CR is generated independently from a normal distribution which has a mean of uCRand a standard deviation of 0.1.

    where meanAdenotes the average value, and SCRis the set of all excellent CRi.

    3.3. PaDE

    In the whole evolution process,k groups are maintained where all individuals are classified in for evolution in PaDE algorithm.Then, uFcan be updated according to the following equation.

    Furthermore, the selection probability is updated according to the following equation.

    In order to adapt to dynamic changes for the population size, a new scheme for reducing the size of the parabolic population is proposed.

    Presenting a new mutation strategy with time stamp can help to initialize all the inferior individuals in the first stage, followed by the insertion into the external archive.

    3.4. RAM-JAPDE

    RAM-JAPDE is a combination of Rank-Based Adaptation of Mutation Strategies (RAM) and Joint Adaptation of Parameters in DE(JAPDE) to produce a new algorithm. The main process and calculation formulas of RAM-JAPDE algorithm are as follows:

    3.5. Improved mechanisms

    Among many adaptive DE algorithms, mutation and crossover operation are two important steps in DE family.The improvement is mainly aimed at F and CR,whose values are of great significance to whether the algorithm can avoid falling into the local optimum and find the global optimal solution.

    For F, it should be as large as possible in the pre-phase of evolution, and quickly approach to the global optimal solution,avoiding trapping in local optimum. The population is clustered around the global optimal solution in later period, so the F value should be as small as possible to improve the search accuracy and find the global optimal solution.

    CR determines the number of mutated individuals that can enter the selection operation.If the value of CR is large,the offspring will inherit less information from the parent, and the population renewal mainly depends on the process of mutation, which not only enriches the diversity of the population,but also improves the ability of global search.If the value of CR is small,the search of the offspring will mainly focus on the surrounding of the parent generation, which will narrow the optimization range, but will improve the convergence speed, which is beneficial for local optimization.

    Since the F and CR values of many DE variants are generated randomly, independent of the number of iterations and individual fitness values,its convergence is not as good as expected.For these problems,an improved adaptive DE algorithm is presented(IADE).In order to make the mutation factor adjust to the extent of evolution, a new dynamic adaptive mutation factor is proposed. The specific expression is as follows:

    where g represents gth generation,G denotes the total generations of evolution, Fmaxand Fminrespectively denotes respectively denotes the maximum and minimum values of the mutation factor.

    The crossover factor CR controls the degree to which each dimension of individual parameters participates in the crossover step, in order to make it change with the individual fitness values and ensure the stability of algorithm convergence, presenting a dynamic adaptive crossover factor. The specific expression is as follows:

    where CRuand CRlrespectively represents the minimum and maximum limits of CR, fiis the individual fitness value, fmaxrepresents the largest value, fbestdenotes the best value.

    In view of the above description, Fig. 2 is the algorithm flowchart of IADE.

    4. Experiments and analysis

    To verify the performance of the IADE algorithm solving how a single UAV attacks multiple targets or executes multi-tasks, simulation experiments have been carried out,and the IADE algorithm is simultaneously compared with other algorithms.

    4.1. Performance of the IADE algorithm

    To prove the validity of the ability of the algorithm to find the optimal solution without falling into the local optimal solution,two experiments are conducted with different target numbers: experiment 1 has a smaller target scale,while experiment 2 has a larger target scale. Before the experiment, we firstly established the 3D space topographic map, and set the initial position of UAV and target points, as well as some parameters of UAV itself.

    Table 1 shows some initializations of the UAV's own parameters,where Num is the UAV number,Position denotes the initial position of UAV, D represents the flight range of UAV, V describes the maximum and minimum flight speeds of UAV and N denotes the maximum number of missions performed by UAV per flight.

    Table 2 sets the number and location of target points, T is the number of targets.

    Parameters setting of the two groups of experiments are shown in Table 3. Where Target is target points, P is individuals in the population, G represents the number of evolutionary iterations,Num denotes the number of times the program is run.

    Fig. 3 shows the optimal simulation results of two groups of experimental UAV attacking targets in the three-dimensional environment. The yellow circle indicates the starting position of the UAV,the yellow stars indicate the location of the target points,the red line indicates the order in which the UAV executes the target points.

    Fig. 2. The algorithm flowchart of IADE.

    Table 1 Initialization of UAV parameters.

    Fig. 4 shows the convergence curves of the two groups of experiments. It can be seen that due to the small number of target points in Experiment 1, the algorithm gets the optimal solutionfrom the beginning till the end of evolution. Experiment 2 converges to the optimal solution at about generation 50, which indicates that the IADE algorithm has good convergence in the early stage of evolution.From the convergence curve of Experiment 2,it can be seen that the algorithm still converges half way through execution, which proves that the algorithm can effectively obtain the optimal solution under the joint action of the adaptive mutation factor and the adaptive crosssover factor, so as to ensure the balance between detection and development.

    Table 3 Set evolution parameters.

    Table 4 lists the statistical data of the results of the two groups of experiments. Where Avgtime denotes the average algorithm execution time, Bestcost represents the optimal total voyage cost,Avgcost represents the average total cost of the voyage, Optimal denotes the efficiency of obtaining the optimal solution, which is defined as a percentage.According to the results of thetwo groups of experiments,the IADE algorithm can solve the UAV multitasking problem well, which can get the optimal solution.

    By comparing the results of the two experiments, although the large-scale experiment 2 adds to the complexity of the UAV multitasking problem,the IADE algorithm can still find the optimal solution within limited time, and the optimal solution proportion reaches 92.5%.The scale of Experiment 1 is relatively small,which indicates that the algorithm can quickly get the optimal solution,and the running time is short, and the efficiency of the optimal solution reaches 97.5%.Although the tasks in experiment 2 is twice as large as that in experiment 1, the average running time of the optimal solution obtained by the two groups of experiments is at the same order of magnitude,which indicates that it is suitable for solving multitasking problems of different scales with high efficiency, and the Avgcost is close to the Bestcost.

    Table 5 shows the optimal allocation results of the two groups of experiments.

    4.2. Comparison with other algorithms

    To further verify the validity and rationality of IADE algorithm,comparison of this algorithm, DMDE algorithm and standard DE algorithm is also made. The comparative experiments are carried out in the exact same simulation environment. Table 6 is the experimental results. The simulation results and convergence curves of each algorithm are respectively shown in Fig.5,Fig.6 and Fig. 7.

    Through the comparative experiments, it can be seen that the IADE algorithm has the shortest execution time, the fastestconvergence speed, the lowest average flight cost and the highest efficiency in obtaining the optimal solution. Although the DEDM algorithm[17]performs well in multi-UAV cooperative multi-target assignment, it performs poorly in dealing with single-UAV multitasking,especially when the scale of the task is relatively large.TheDEDM algorithm has the longest running time and the highest average flight cost, and the efficiency of obtaining the optimal solution is only 30 %. Although standard DE algorithm is better than DMDE algorithm, its execution time is relatively longer compared with IADE algorithm. The efficiency of obtaining the optimal solution is 72.5 %, which is 20 % lower than IADE algorithm.

    Table 2 Initialization of target coordinates.

    Fig. 3. Simulation results of Experiment 1 and Experiment 2.

    Fig. 4. Convergence curves of Experiment 1 and Experiment 2.

    Table 4 Statistics of the results.

    Table 6 Statistics of the results.

    Table 5 Assignment results.

    Fig. 5. Results and convergence curves of IADE algorithm.

    Fig. 6. Results and convergence curves of DMDE algorithm.

    Fig. 7. Results and convergence curves of DE algorithm.

    5. Conclusions

    In this paper,a novel algorithm is proposed to solve the problem of difficult modeling and complex calculation of UAV during multitask execution. The main contributions are as follows:

    (1) A single UAV multitasking model and its related constraints are established in a three-dimensional environment.

    (2) An adaptive evolutionary algorithm based on population iteration times and individual fitness values to adjust the mutation factor and crossover factor is proposed to solve the problem of exploratory and exploitive balance existing in the current evolutionary process.

    (3) Through different scales simulation experiments, the rationality of the IADE algorithm is proved.

    (4) The comparative experiments with other algorithms show the superiority of the IADE algorithm.

    The future work is to apply the IADE algorithm to UAV path planning.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    This work is not funded and supported by any organization or institute.

    在线观看舔阴道视频| 此物有八面人人有两片| 宅男免费午夜| 两性夫妻黄色片| 免费看日本二区| 免费在线观看日本一区| 久久香蕉精品热| 午夜福利欧美成人| 精品久久久久久久末码| 国产精品 欧美亚洲| 日韩国内少妇激情av| 日本与韩国留学比较| 亚洲av免费在线观看| 午夜福利在线观看吧| 亚洲第一欧美日韩一区二区三区| 亚洲中文日韩欧美视频| 欧美+亚洲+日韩+国产| 久久精品影院6| 啪啪无遮挡十八禁网站| 色尼玛亚洲综合影院| 国内精品一区二区在线观看| 亚洲国产精品sss在线观看| or卡值多少钱| 欧美3d第一页| 久久久久国产精品人妻aⅴ院| 午夜影院日韩av| 90打野战视频偷拍视频| 国产97色在线日韩免费| 国产亚洲精品av在线| 999久久久国产精品视频| www.999成人在线观看| 高清毛片免费观看视频网站| 免费观看精品视频网站| 99久久精品一区二区三区| 国产激情欧美一区二区| 成人18禁在线播放| 精品一区二区三区av网在线观看| 亚洲av成人精品一区久久| 夜夜躁狠狠躁天天躁| 色老头精品视频在线观看| 国产精品国产高清国产av| 亚洲真实伦在线观看| 久久久久久久午夜电影| 两人在一起打扑克的视频| 亚洲av片天天在线观看| 国产精品影院久久| 51午夜福利影视在线观看| 日本黄色片子视频| 一个人看视频在线观看www免费 | 国产精品九九99| 久久性视频一级片| 日日摸夜夜添夜夜添小说| 每晚都被弄得嗷嗷叫到高潮| 露出奶头的视频| 在线观看免费午夜福利视频| 99国产精品99久久久久| 18禁观看日本| 叶爱在线成人免费视频播放| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av| 亚洲片人在线观看| 青草久久国产| 一区二区三区国产精品乱码| 2021天堂中文幕一二区在线观| 精品久久久久久久毛片微露脸| 三级毛片av免费| 免费一级毛片在线播放高清视频| ponron亚洲| 国产精品九九99| 99久久国产精品久久久| netflix在线观看网站| 欧美又色又爽又黄视频| 狠狠狠狠99中文字幕| 99久国产av精品| www日本在线高清视频| 99re在线观看精品视频| 18禁裸乳无遮挡免费网站照片| 欧美大码av| 99视频精品全部免费 在线 | 99国产精品一区二区三区| 巨乳人妻的诱惑在线观看| 搡老妇女老女人老熟妇| 国内久久婷婷六月综合欲色啪| av福利片在线观看| 日韩三级视频一区二区三区| 国产69精品久久久久777片 | 麻豆国产av国片精品| 精品国内亚洲2022精品成人| 2021天堂中文幕一二区在线观| 搡老妇女老女人老熟妇| 最近最新中文字幕大全电影3| 亚洲国产高清在线一区二区三| 99热只有精品国产| 老司机午夜十八禁免费视频| 成人高潮视频无遮挡免费网站| 99re在线观看精品视频| 欧美乱妇无乱码| 免费观看精品视频网站| 少妇人妻一区二区三区视频| 久久久久国产精品人妻aⅴ院| 精品人妻1区二区| 成在线人永久免费视频| 熟女少妇亚洲综合色aaa.| h日本视频在线播放| 国产91精品成人一区二区三区| 欧美成狂野欧美在线观看| 国产精品电影一区二区三区| 午夜精品久久久久久毛片777| 观看免费一级毛片| 97碰自拍视频| 亚洲精品中文字幕一二三四区| 国产精品久久电影中文字幕| 又紧又爽又黄一区二区| 久久香蕉国产精品| 色av中文字幕| 两个人的视频大全免费| 久久久精品欧美日韩精品| 首页视频小说图片口味搜索| 日本成人三级电影网站| 国产精品九九99| 欧美丝袜亚洲另类 | 欧美在线一区亚洲| 99热这里只有是精品50| 99国产极品粉嫩在线观看| 制服丝袜大香蕉在线| 一卡2卡三卡四卡精品乱码亚洲| 老司机在亚洲福利影院| 久久香蕉国产精品| 老司机福利观看| 一二三四在线观看免费中文在| 色吧在线观看| xxxwww97欧美| 久久香蕉精品热| 亚洲国产欧美人成| 亚洲国产中文字幕在线视频| 亚洲国产高清在线一区二区三| 亚洲精品粉嫩美女一区| 熟女少妇亚洲综合色aaa.| 一区二区三区激情视频| 免费大片18禁| 国产在线精品亚洲第一网站| 最好的美女福利视频网| 国产成人影院久久av| 一级毛片精品| 男人舔女人下体高潮全视频| 熟女人妻精品中文字幕| 熟妇人妻久久中文字幕3abv| 亚洲在线观看片| 久久久久国内视频| 日韩高清综合在线| 法律面前人人平等表现在哪些方面| 国产毛片a区久久久久| 男女午夜视频在线观看| 国产精品爽爽va在线观看网站| 久久天躁狠狠躁夜夜2o2o| 可以在线观看的亚洲视频| 国产精品电影一区二区三区| 亚洲成人免费电影在线观看| 一级a爱片免费观看的视频| 床上黄色一级片| 午夜免费激情av| 色av中文字幕| 亚洲午夜精品一区,二区,三区| 好男人电影高清在线观看| 国产黄a三级三级三级人| 老司机午夜福利在线观看视频| 欧美成人免费av一区二区三区| 国产亚洲精品一区二区www| 亚洲国产精品sss在线观看| 日韩欧美免费精品| e午夜精品久久久久久久| 99久久精品一区二区三区| 国产爱豆传媒在线观看| 亚洲av成人av| 亚洲五月天丁香| 在线播放国产精品三级| 亚洲欧美精品综合一区二区三区| 一级a爱片免费观看的视频| 亚洲欧美一区二区三区黑人| 色av中文字幕| 亚洲欧美精品综合一区二区三区| 1024香蕉在线观看| 在线观看午夜福利视频| 亚洲国产中文字幕在线视频| 全区人妻精品视频| 国产av不卡久久| 国产成人aa在线观看| 韩国av一区二区三区四区| 日本与韩国留学比较| 国产午夜精品久久久久久| 久久性视频一级片| 国产久久久一区二区三区| 国产成人av激情在线播放| 白带黄色成豆腐渣| 欧美又色又爽又黄视频| 99在线人妻在线中文字幕| 欧美激情在线99| 夜夜夜夜夜久久久久| 天堂网av新在线| 亚洲专区国产一区二区| 国产三级中文精品| 欧美av亚洲av综合av国产av| 久久久精品欧美日韩精品| tocl精华| 久久午夜综合久久蜜桃| 成人永久免费在线观看视频| 免费在线观看影片大全网站| 午夜a级毛片| 综合色av麻豆| 精品一区二区三区视频在线 | 丝袜人妻中文字幕| www国产在线视频色| 波多野结衣巨乳人妻| 757午夜福利合集在线观看| 男人舔女人的私密视频| 日本精品一区二区三区蜜桃| 好男人在线观看高清免费视频| 国产探花在线观看一区二区| 日本黄色视频三级网站网址| 欧美色欧美亚洲另类二区| 熟妇人妻久久中文字幕3abv| 美女免费视频网站| 叶爱在线成人免费视频播放| 女人被狂操c到高潮| 亚洲最大成人中文| 2021天堂中文幕一二区在线观| 久久久久久大精品| 午夜精品久久久久久毛片777| 精品一区二区三区av网在线观看| 亚洲色图av天堂| 男女下面进入的视频免费午夜| 欧美黑人巨大hd| 日韩欧美国产在线观看| 一二三四在线观看免费中文在| 亚洲片人在线观看| 麻豆av在线久日| 色综合欧美亚洲国产小说| 精品久久久久久久人妻蜜臀av| 国产一区二区三区视频了| 88av欧美| 毛片女人毛片| 99久久成人亚洲精品观看| 老司机深夜福利视频在线观看| 日本五十路高清| 日韩欧美免费精品| www国产在线视频色| 亚洲自偷自拍图片 自拍| 午夜福利高清视频| 久久久精品欧美日韩精品| 成年免费大片在线观看| 脱女人内裤的视频| 老熟妇仑乱视频hdxx| 久久中文字幕一级| 欧美色欧美亚洲另类二区| 亚洲国产日韩欧美精品在线观看 | 久久久久免费精品人妻一区二区| 日本黄大片高清| 草草在线视频免费看| 亚洲黑人精品在线| 观看免费一级毛片| 神马国产精品三级电影在线观看| 欧美日韩乱码在线| 两人在一起打扑克的视频| 99国产精品一区二区三区| 欧美日本亚洲视频在线播放| 噜噜噜噜噜久久久久久91| 精品一区二区三区四区五区乱码| 欧美乱色亚洲激情| av欧美777| 欧美成狂野欧美在线观看| 日韩大尺度精品在线看网址| 精品国产亚洲在线| 午夜福利免费观看在线| 久久精品91无色码中文字幕| 国产又黄又爽又无遮挡在线| 国产精品乱码一区二三区的特点| 全区人妻精品视频| 久久久久亚洲av毛片大全| 少妇丰满av| 亚洲精品粉嫩美女一区| 男女之事视频高清在线观看| av女优亚洲男人天堂 | 国产97色在线日韩免费| 欧美乱色亚洲激情| 在线永久观看黄色视频| 美女cb高潮喷水在线观看 | 亚洲国产看品久久| 嫩草影院入口| 久久伊人香网站| 香蕉久久夜色| 日本免费一区二区三区高清不卡| 国产精品九九99| 男人的好看免费观看在线视频| 国产97色在线日韩免费| 夜夜躁狠狠躁天天躁| 午夜两性在线视频| 欧美在线黄色| 女生性感内裤真人,穿戴方法视频| 51午夜福利影视在线观看| 99久久综合精品五月天人人| 国产伦一二天堂av在线观看| 欧美日韩黄片免| 美女高潮喷水抽搐中文字幕| 国模一区二区三区四区视频 | 人人妻人人澡欧美一区二区| 久久久国产精品麻豆| 亚洲国产色片| 中文字幕最新亚洲高清| 亚洲精品一卡2卡三卡4卡5卡| av欧美777| 色播亚洲综合网| 精品久久久久久成人av| 国产精品98久久久久久宅男小说| 精品国产三级普通话版| 18禁美女被吸乳视频| 无人区码免费观看不卡| 国产欧美日韩精品一区二区| 又紧又爽又黄一区二区| 美女黄网站色视频| 国产精品 国内视频| 精品久久久久久久久久久久久| 国产伦一二天堂av在线观看| 丁香欧美五月| 午夜免费激情av| 国产av一区在线观看免费| 不卡av一区二区三区| www.自偷自拍.com| 精品久久久久久久毛片微露脸| 欧美成人免费av一区二区三区| 久久精品综合一区二区三区| 九九热线精品视视频播放| 可以在线观看的亚洲视频| 亚洲av片天天在线观看| 久久人人精品亚洲av| 成人午夜高清在线视频| 99久国产av精品| 老熟妇乱子伦视频在线观看| 18禁国产床啪视频网站| 午夜激情福利司机影院| 欧美绝顶高潮抽搐喷水| 国产不卡一卡二| 桃红色精品国产亚洲av| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 国产黄色小视频在线观看| 国产单亲对白刺激| 久久国产精品影院| 国产黄a三级三级三级人| 成年人黄色毛片网站| 国产97色在线日韩免费| 中文字幕久久专区| 又黄又粗又硬又大视频| 色噜噜av男人的天堂激情| а√天堂www在线а√下载| 国产人伦9x9x在线观看| 国产精品98久久久久久宅男小说| 欧美av亚洲av综合av国产av| or卡值多少钱| 国产日本99.免费观看| 日韩国内少妇激情av| 搡老岳熟女国产| 国产高清有码在线观看视频| 最新在线观看一区二区三区| 国内精品一区二区在线观看| 免费av毛片视频| 这个男人来自地球电影免费观看| 亚洲人成网站高清观看| 久久人人精品亚洲av| 国产亚洲精品久久久com| 少妇熟女aⅴ在线视频| 老汉色av国产亚洲站长工具| 精品99又大又爽又粗少妇毛片 | 婷婷精品国产亚洲av| 亚洲自拍偷在线| 真人一进一出gif抽搐免费| 亚洲国产欧美网| 99re在线观看精品视频| 久久九九热精品免费| 搡老熟女国产l中国老女人| 久久久成人免费电影| 久99久视频精品免费| 国产精品98久久久久久宅男小说| 国产亚洲欧美在线一区二区| 99久久精品国产亚洲精品| 欧美日韩综合久久久久久 | 欧美性猛交黑人性爽| 色尼玛亚洲综合影院| av天堂在线播放| 国产成人影院久久av| 国产精品亚洲美女久久久| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 欧美一区二区精品小视频在线| 天天躁日日操中文字幕| 男女那种视频在线观看| 日韩欧美国产一区二区入口| 99riav亚洲国产免费| 中文资源天堂在线| 99热只有精品国产| 99久久国产精品久久久| 日韩国内少妇激情av| 国产三级黄色录像| 亚洲第一电影网av| 欧美日本亚洲视频在线播放| 国产欧美日韩一区二区精品| 国产精品乱码一区二三区的特点| 欧美日韩中文字幕国产精品一区二区三区| 在线播放国产精品三级| 亚洲自拍偷在线| 亚洲欧洲精品一区二区精品久久久| av天堂在线播放| 成人性生交大片免费视频hd| 麻豆成人av在线观看| 欧美成人性av电影在线观看| 欧美最黄视频在线播放免费| 巨乳人妻的诱惑在线观看| 可以在线观看的亚洲视频| 一级黄色大片毛片| 精品免费久久久久久久清纯| 99精品在免费线老司机午夜| 精品国产三级普通话版| 97超级碰碰碰精品色视频在线观看| 在线观看舔阴道视频| 哪里可以看免费的av片| 成人18禁在线播放| 好男人电影高清在线观看| 亚洲五月天丁香| 免费观看精品视频网站| 久久这里只有精品中国| 国产成人一区二区三区免费视频网站| 亚洲午夜精品一区,二区,三区| 19禁男女啪啪无遮挡网站| 国产精品永久免费网站| 久久国产乱子伦精品免费另类| 母亲3免费完整高清在线观看| 亚洲无线观看免费| 波多野结衣高清作品| 日韩欧美在线二视频| 人人妻人人澡欧美一区二区| 久久午夜亚洲精品久久| 国产精品乱码一区二三区的特点| 欧美日韩中文字幕国产精品一区二区三区| 国产人伦9x9x在线观看| 99re在线观看精品视频| 老司机午夜十八禁免费视频| 午夜福利在线观看免费完整高清在 | 亚洲 欧美 日韩 在线 免费| 搡老妇女老女人老熟妇| 一级毛片女人18水好多| 黄色女人牲交| 精品久久蜜臀av无| 丁香欧美五月| 女人被狂操c到高潮| 国产精品影院久久| 99精品久久久久人妻精品| 九九在线视频观看精品| 51午夜福利影视在线观看| 日韩欧美在线乱码| 黄频高清免费视频| 一本一本综合久久| 久久国产精品影院| 51午夜福利影视在线观看| 天堂av国产一区二区熟女人妻| 宅男免费午夜| 欧美日韩亚洲国产一区二区在线观看| 伊人久久大香线蕉亚洲五| 给我免费播放毛片高清在线观看| 欧美日韩瑟瑟在线播放| 亚洲狠狠婷婷综合久久图片| 国产综合懂色| 日韩欧美精品v在线| 国产淫片久久久久久久久 | 精品国产三级普通话版| 三级国产精品片| 男的添女的下面高潮视频| 亚洲国产最新在线播放| 欧美人与善性xxx| 色哟哟·www| 亚洲精品一区蜜桃| av在线观看视频网站免费| 中文乱码字字幕精品一区二区三区 | 午夜激情欧美在线| 欧美日韩综合久久久久久| 狂野欧美激情性xxxx在线观看| 国产 一区 欧美 日韩| 国产午夜精品久久久久久一区二区三区| 日本一二三区视频观看| 成年版毛片免费区| 乱人视频在线观看| 中文字幕久久专区| 久久精品夜夜夜夜夜久久蜜豆| 麻豆精品久久久久久蜜桃| 国产精品一区二区性色av| 女人十人毛片免费观看3o分钟| 97超碰精品成人国产| 在线免费观看的www视频| 干丝袜人妻中文字幕| 一区二区三区高清视频在线| 男女视频在线观看网站免费| 欧美+日韩+精品| 色综合站精品国产| av在线老鸭窝| 亚洲综合色惰| 国产色爽女视频免费观看| 黑人高潮一二区| 国产免费又黄又爽又色| 九九久久精品国产亚洲av麻豆| av线在线观看网站| eeuss影院久久| 久热久热在线精品观看| 亚洲精品亚洲一区二区| 国内精品宾馆在线| av福利片在线观看| av视频在线观看入口| 伦精品一区二区三区| 中文字幕熟女人妻在线| 免费人成在线观看视频色| 日韩欧美三级三区| 在线播放国产精品三级| 国产精品无大码| 热99re8久久精品国产| 亚洲精品aⅴ在线观看| 午夜免费激情av| av福利片在线观看| 日日摸夜夜添夜夜爱| 国产亚洲精品av在线| 欧美性猛交黑人性爽| 联通29元200g的流量卡| 99久久中文字幕三级久久日本| 国产v大片淫在线免费观看| 久久欧美精品欧美久久欧美| 水蜜桃什么品种好| 人妻夜夜爽99麻豆av| 大香蕉97超碰在线| 国产一级毛片七仙女欲春2| 午夜老司机福利剧场| 老司机影院毛片| 91久久精品国产一区二区成人| 亚洲在久久综合| 91在线精品国自产拍蜜月| 黄色配什么色好看| 欧美成人a在线观看| 岛国毛片在线播放| 久久6这里有精品| 国产毛片a区久久久久| 日本av手机在线免费观看| 国产成人一区二区在线| 免费观看人在逋| 亚洲性久久影院| 亚洲av.av天堂| 乱码一卡2卡4卡精品| 亚洲av免费在线观看| 免费人成在线观看视频色| 国产精品,欧美在线| 国产精品国产三级国产专区5o | 久久精品综合一区二区三区| 欧美丝袜亚洲另类| 人人妻人人澡人人爽人人夜夜 | 欧美日本亚洲视频在线播放| 人妻系列 视频| 午夜福利高清视频| 在线免费十八禁| 午夜福利在线观看免费完整高清在| 18禁在线无遮挡免费观看视频| 网址你懂的国产日韩在线| 一区二区三区免费毛片| 三级经典国产精品| 免费观看的影片在线观看| 国产精品一区二区性色av| 性插视频无遮挡在线免费观看| 少妇高潮的动态图| 99久久无色码亚洲精品果冻| 久久久久久久久久久免费av| 嫩草影院新地址| 天天一区二区日本电影三级| 精品久久久久久久末码| 91狼人影院| 国产亚洲最大av| 人人妻人人看人人澡| 久久久精品欧美日韩精品| 久久鲁丝午夜福利片| 欧美zozozo另类| 国产精品1区2区在线观看.| 亚洲精品国产av成人精品| 国产女主播在线喷水免费视频网站 | 3wmmmm亚洲av在线观看| 内地一区二区视频在线| 国产精品,欧美在线| 久久久精品大字幕| 国产亚洲午夜精品一区二区久久 | 99久久中文字幕三级久久日本| 国产 一区 欧美 日韩| 色吧在线观看| 99热这里只有是精品50| 国产男人的电影天堂91| 国产免费又黄又爽又色| 日本猛色少妇xxxxx猛交久久| 日韩在线高清观看一区二区三区| 天美传媒精品一区二区| 两个人的视频大全免费| 18禁在线无遮挡免费观看视频| 欧美三级亚洲精品| 国产精品久久久久久精品电影小说 | 国产精品一区二区三区四区免费观看| 免费看美女性在线毛片视频| 白带黄色成豆腐渣| 禁无遮挡网站| 午夜福利高清视频| 国产v大片淫在线免费观看| www.色视频.com| 国产色婷婷99| 国产视频首页在线观看| 看免费成人av毛片| 免费观看的影片在线观看| 婷婷色麻豆天堂久久 |