• <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.

    啦啦啦免费观看视频1| 国产熟女午夜一区二区三区| 在线天堂中文资源库| 一个人免费看片子| 久久人人爽人人片av| 夜夜骑夜夜射夜夜干| 成人影院久久| 国产精品亚洲av一区麻豆| 久久精品亚洲av国产电影网| 亚洲人成77777在线视频| 中文字幕人妻丝袜一区二区| 国产在线观看jvid| 天天躁夜夜躁狠狠躁躁| 欧美日韩黄片免| 一二三四在线观看免费中文在| 91精品国产国语对白视频| 色综合欧美亚洲国产小说| 视频区图区小说| 伊人亚洲综合成人网| 狠狠婷婷综合久久久久久88av| 午夜福利免费观看在线| 无遮挡黄片免费观看| 91字幕亚洲| 黑人巨大精品欧美一区二区mp4| 99国产精品一区二区三区| 亚洲色图综合在线观看| 伦理电影免费视频| 可以免费在线观看a视频的电影网站| 青草久久国产| a在线观看视频网站| 纵有疾风起免费观看全集完整版| 老司机影院毛片| 19禁男女啪啪无遮挡网站| 人人妻人人爽人人添夜夜欢视频| 国产91精品成人一区二区三区 | 男人爽女人下面视频在线观看| 亚洲五月婷婷丁香| 国产精品.久久久| www.999成人在线观看| 欧美黑人精品巨大| 在线观看人妻少妇| 人人妻,人人澡人人爽秒播| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品中文字幕在线视频| 侵犯人妻中文字幕一二三四区| 黄色视频在线播放观看不卡| 国产精品久久久久久精品电影小说| 国产在视频线精品| netflix在线观看网站| 中文字幕高清在线视频| 成人影院久久| 黑人欧美特级aaaaaa片| 丰满人妻熟妇乱又伦精品不卡| 国产精品 欧美亚洲| 久久久久国产一级毛片高清牌| 一级片'在线观看视频| 大码成人一级视频| 国产成人欧美在线观看 | 极品少妇高潮喷水抽搐| 丁香六月天网| 99久久人妻综合| 日韩欧美国产一区二区入口| 亚洲av成人不卡在线观看播放网 | 国产精品 国内视频| 亚洲国产精品一区二区三区在线| 午夜福利在线免费观看网站| 最近中文字幕2019免费版| 久久人妻福利社区极品人妻图片| 欧美日本中文国产一区发布| 日韩三级视频一区二区三区| 亚洲av美国av| 美女高潮喷水抽搐中文字幕| 91麻豆精品激情在线观看国产 | 国产主播在线观看一区二区| 亚洲成av片中文字幕在线观看| 午夜视频精品福利| 亚洲熟女精品中文字幕| 狠狠狠狠99中文字幕| 欧美黑人欧美精品刺激| 免费观看人在逋| 91成人精品电影| 国产日韩欧美视频二区| 国产精品久久久久久精品电影小说| 五月天丁香电影| 国产成+人综合+亚洲专区| 亚洲精品国产av成人精品| 精品国产国语对白av| 日日夜夜操网爽| 91大片在线观看| 黑人欧美特级aaaaaa片| 国产精品一区二区免费欧美 | 午夜免费成人在线视频| 自线自在国产av| 亚洲伊人色综图| 亚洲国产毛片av蜜桃av| 在线观看一区二区三区激情| 一区福利在线观看| 蜜桃在线观看..| 国产一卡二卡三卡精品| 久久中文看片网| 国产伦理片在线播放av一区| 亚洲男人天堂网一区| 国产一卡二卡三卡精品| 日韩欧美一区视频在线观看| 午夜免费观看性视频| 久久精品国产亚洲av香蕉五月 | 大片电影免费在线观看免费| 精品国产乱码久久久久久小说| 日韩有码中文字幕| 欧美日韩福利视频一区二区| 亚洲国产av新网站| 国产一区二区在线观看av| 日韩制服骚丝袜av| 午夜老司机福利片| 99热国产这里只有精品6| 麻豆av在线久日| 中亚洲国语对白在线视频| 在线观看免费日韩欧美大片| 香蕉国产在线看| 亚洲熟女精品中文字幕| 亚洲av男天堂| 777米奇影视久久| 国产亚洲午夜精品一区二区久久| 多毛熟女@视频| 国产av国产精品国产| 欧美日韩国产mv在线观看视频| 咕卡用的链子| 男女午夜视频在线观看| 日本91视频免费播放| tube8黄色片| 精品一品国产午夜福利视频| 人妻一区二区av| 久久国产精品大桥未久av| 天堂中文最新版在线下载| 99国产精品一区二区蜜桃av | 丝袜脚勾引网站| 国产精品国产三级国产专区5o| 午夜精品国产一区二区电影| 国产深夜福利视频在线观看| 少妇精品久久久久久久| 在线永久观看黄色视频| 老汉色av国产亚洲站长工具| 俄罗斯特黄特色一大片| 天天操日日干夜夜撸| 久久精品成人免费网站| 十八禁人妻一区二区| 18禁黄网站禁片午夜丰满| 777米奇影视久久| 午夜激情久久久久久久| 青青草视频在线视频观看| 精品少妇内射三级| 黑人操中国人逼视频| 久久久久网色| 久热爱精品视频在线9| 亚洲成人免费电影在线观看| 国产精品一二三区在线看| 嫩草影视91久久| av视频免费观看在线观看| 国产精品九九99| 99国产精品一区二区蜜桃av | 97精品久久久久久久久久精品| 韩国高清视频一区二区三区| 人人妻,人人澡人人爽秒播| 精品亚洲乱码少妇综合久久| 欧美av亚洲av综合av国产av| 1024香蕉在线观看| www.999成人在线观看| 精品福利永久在线观看| 久久九九热精品免费| 久久性视频一级片| 亚洲美女黄色视频免费看| 99久久国产精品久久久| 最新在线观看一区二区三区| 午夜激情久久久久久久| 搡老乐熟女国产| 国产99久久九九免费精品| 又黄又粗又硬又大视频| 亚洲成人免费电影在线观看| 自线自在国产av| 欧美在线黄色| 久久毛片免费看一区二区三区| 国产在线视频一区二区| 最近中文字幕2019免费版| 母亲3免费完整高清在线观看| 一级片免费观看大全| 后天国语完整版免费观看| 在线av久久热| 欧美另类亚洲清纯唯美| 制服人妻中文乱码| 精品卡一卡二卡四卡免费| 欧美性长视频在线观看| 亚洲中文字幕日韩| 精品人妻一区二区三区麻豆| 久久精品aⅴ一区二区三区四区| 男女床上黄色一级片免费看| 午夜福利在线免费观看网站| 久久免费观看电影| 国产成人av教育| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品一区三区| 天天躁日日躁夜夜躁夜夜| 亚洲欧美日韩另类电影网站| 黄色 视频免费看| 欧美av亚洲av综合av国产av| 中文字幕人妻丝袜一区二区| 国产精品.久久久| 狠狠狠狠99中文字幕| 青草久久国产| 精品久久蜜臀av无| 9色porny在线观看| 老司机影院成人| 男人舔女人的私密视频| 国产精品国产av在线观看| 国产一区二区三区在线臀色熟女 | 成在线人永久免费视频| 国产精品 国内视频| 国产一区二区在线观看av| 欧美久久黑人一区二区| 超碰成人久久| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久人人做人人爽| www.av在线官网国产| 纵有疾风起免费观看全集完整版| 亚洲美女黄色视频免费看| 欧美另类亚洲清纯唯美| 日日爽夜夜爽网站| 下体分泌物呈黄色| 久久精品国产a三级三级三级| 国产淫语在线视频| 老汉色av国产亚洲站长工具| 国产免费福利视频在线观看| 亚洲精品中文字幕一二三四区 | 黄频高清免费视频| 亚洲美女黄色视频免费看| 亚洲精品乱久久久久久| 久久久久久久精品精品| 黑人巨大精品欧美一区二区mp4| tube8黄色片| 国产精品二区激情视频| 满18在线观看网站| 日本vs欧美在线观看视频| 久久热在线av| 女性被躁到高潮视频| 午夜免费成人在线视频| kizo精华| 欧美xxⅹ黑人| 夜夜夜夜夜久久久久| 精品少妇久久久久久888优播| 狂野欧美激情性bbbbbb| 亚洲欧美色中文字幕在线| 女人久久www免费人成看片| 嫁个100分男人电影在线观看| 免费久久久久久久精品成人欧美视频| 91九色精品人成在线观看| 另类精品久久| 欧美av亚洲av综合av国产av| 大香蕉久久网| 十八禁人妻一区二区| 亚洲伊人久久精品综合| 精品少妇久久久久久888优播| 国产亚洲av高清不卡| 欧美在线黄色| 久久精品人人爽人人爽视色| 久久九九热精品免费| 777米奇影视久久| 精品国产国语对白av| 国产欧美日韩一区二区三区在线| 黄片小视频在线播放| 国产精品久久久久久人妻精品电影 | 99国产精品一区二区三区| 日韩中文字幕欧美一区二区| 亚洲国产看品久久| 日韩欧美免费精品| 日本一区二区免费在线视频| av天堂久久9| 老司机在亚洲福利影院| 天天躁狠狠躁夜夜躁狠狠躁| 免费女性裸体啪啪无遮挡网站| 中文欧美无线码| 欧美国产精品一级二级三级| 成年人午夜在线观看视频| 亚洲国产毛片av蜜桃av| 日韩视频一区二区在线观看| 中文字幕人妻熟女乱码| 国产高清videossex| 欧美国产精品va在线观看不卡| 三上悠亚av全集在线观看| 亚洲精品一二三| 男女下面插进去视频免费观看| 熟女少妇亚洲综合色aaa.| 老司机影院成人| av在线app专区| 国产精品自产拍在线观看55亚洲 | 亚洲专区国产一区二区| 亚洲欧美色中文字幕在线| 俄罗斯特黄特色一大片| 日本vs欧美在线观看视频| 日日摸夜夜添夜夜添小说| 国产欧美日韩一区二区三区在线| 欧美老熟妇乱子伦牲交| 国产伦人伦偷精品视频| 99热网站在线观看| 亚洲国产欧美网| 飞空精品影院首页| 久久久久精品人妻al黑| 在线观看免费日韩欧美大片| 精品一区在线观看国产| 色婷婷久久久亚洲欧美| 天天躁狠狠躁夜夜躁狠狠躁| 90打野战视频偷拍视频| 免费在线观看黄色视频的| 精品一区二区三区av网在线观看 | 777米奇影视久久| 欧美国产精品va在线观看不卡| 亚洲全国av大片| 亚洲精品美女久久av网站| 如日韩欧美国产精品一区二区三区| 国产一区二区 视频在线| 秋霞在线观看毛片| 久久毛片免费看一区二区三区| 丝袜人妻中文字幕| 免费一级毛片在线播放高清视频 | 夫妻午夜视频| 在线 av 中文字幕| 久久狼人影院| 欧美xxⅹ黑人| 亚洲国产欧美在线一区| 男男h啪啪无遮挡| 伊人久久大香线蕉亚洲五| 日韩三级视频一区二区三区| 如日韩欧美国产精品一区二区三区| 久久午夜综合久久蜜桃| 国产亚洲一区二区精品| 久久av网站| 国产精品成人在线| 国产精品一区二区在线不卡| 国产精品久久久久久精品古装| 深夜精品福利| 久久性视频一级片| 人人妻人人添人人爽欧美一区卜| 超碰97精品在线观看| 免费在线观看完整版高清| 国产精品成人在线| 精品一区在线观看国产| 久久久精品区二区三区| 在线天堂中文资源库| 欧美人与性动交α欧美精品济南到| 少妇的丰满在线观看| av不卡在线播放| 久热这里只有精品99| 1024视频免费在线观看| 国产精品一区二区在线不卡| 亚洲精品中文字幕在线视频| 亚洲中文av在线| 午夜老司机福利片| 美国免费a级毛片| 天天影视国产精品| 男女下面插进去视频免费观看| 91麻豆av在线| 老司机影院成人| 蜜桃在线观看..| 后天国语完整版免费观看| 黄色毛片三级朝国网站| 新久久久久国产一级毛片| 日本a在线网址| 一级片'在线观看视频| 777久久人妻少妇嫩草av网站| 久久久欧美国产精品| 久久精品亚洲av国产电影网| 人人澡人人妻人| 在线观看www视频免费| 99re6热这里在线精品视频| 美女主播在线视频| 精品少妇内射三级| 久久狼人影院| 在线观看免费日韩欧美大片| 少妇猛男粗大的猛烈进出视频| videos熟女内射| 国产伦人伦偷精品视频| 成人影院久久| 亚洲国产成人一精品久久久| 丝袜美腿诱惑在线| www.精华液| 亚洲中文字幕日韩| 国产亚洲欧美在线一区二区| 99热网站在线观看| 热99国产精品久久久久久7| 国产av一区二区精品久久| 国产精品偷伦视频观看了| 国产免费视频播放在线视频| 自线自在国产av| 天堂中文最新版在线下载| 午夜成年电影在线免费观看| 飞空精品影院首页| 视频区欧美日本亚洲| 中文字幕色久视频| 一级黄色大片毛片| 岛国毛片在线播放| 狂野欧美激情性xxxx| 久久国产精品影院| www.av在线官网国产| 十分钟在线观看高清视频www| 老司机影院成人| 蜜桃在线观看..| 国产国语露脸激情在线看| 岛国在线观看网站| 啦啦啦 在线观看视频| 大片电影免费在线观看免费| 婷婷成人精品国产| 国产伦人伦偷精品视频| 国产精品 欧美亚洲| 美女视频免费永久观看网站| 日韩视频一区二区在线观看| 国产成人啪精品午夜网站| 一区二区三区激情视频| 亚洲精品美女久久av网站| 亚洲免费av在线视频| 国产精品99久久99久久久不卡| 国产欧美日韩精品亚洲av| 亚洲成国产人片在线观看| 久久国产精品影院| 大香蕉久久网| 一级片'在线观看视频| 老汉色∧v一级毛片| 国产av精品麻豆| 久久久久精品国产欧美久久久 | 国产精品国产av在线观看| 黄色片一级片一级黄色片| 色老头精品视频在线观看| 精品少妇黑人巨大在线播放| 十分钟在线观看高清视频www| 一级黄色大片毛片| 大香蕉久久网| 美女扒开内裤让男人捅视频| 黄网站色视频无遮挡免费观看| 12—13女人毛片做爰片一| 俄罗斯特黄特色一大片| 国产亚洲午夜精品一区二区久久| 黄色 视频免费看| 日本a在线网址| 久久综合国产亚洲精品| 91成人精品电影| 这个男人来自地球电影免费观看| 黄色a级毛片大全视频| 亚洲精品成人av观看孕妇| 老汉色av国产亚洲站长工具| 午夜精品久久久久久毛片777| 午夜免费鲁丝| 日韩欧美国产一区二区入口| 最近最新免费中文字幕在线| 亚洲精品av麻豆狂野| 日本欧美视频一区| 国产一区二区三区av在线| 中文字幕高清在线视频| 丝袜喷水一区| 国产精品一区二区在线观看99| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美精品自产自拍| 亚洲国产日韩一区二区| 国产精品麻豆人妻色哟哟久久| 99久久人妻综合| 国产有黄有色有爽视频| 欧美少妇被猛烈插入视频| 免费在线观看影片大全网站| 高清黄色对白视频在线免费看| 国产成人av激情在线播放| 欧美在线黄色| 国产麻豆69| 亚洲欧美激情在线| 免费观看人在逋| 日日摸夜夜添夜夜添小说| 欧美黑人精品巨大| 免费不卡黄色视频| 精品少妇一区二区三区视频日本电影| 国产精品久久久人人做人人爽| 一区在线观看完整版| 下体分泌物呈黄色| 黑丝袜美女国产一区| 久久99热这里只频精品6学生| 国产淫语在线视频| 女人精品久久久久毛片| 美女视频免费永久观看网站| 一级片免费观看大全| 欧美人与性动交α欧美软件| 国产亚洲欧美在线一区二区| 午夜久久久在线观看| 亚洲全国av大片| 超色免费av| 亚洲国产精品成人久久小说| 老司机影院毛片| 黑人猛操日本美女一级片| 无遮挡黄片免费观看| 精品国产乱码久久久久久小说| 韩国高清视频一区二区三区| a级毛片在线看网站| 丰满少妇做爰视频| 中文字幕人妻丝袜一区二区| 一本色道久久久久久精品综合| 狠狠精品人妻久久久久久综合| 桃红色精品国产亚洲av| 欧美激情久久久久久爽电影 | 人妻人人澡人人爽人人| 人妻一区二区av| 法律面前人人平等表现在哪些方面 | 亚洲国产看品久久| 日本撒尿小便嘘嘘汇集6| 美女视频免费永久观看网站| 亚洲精品自拍成人| 免费日韩欧美在线观看| 亚洲精品久久成人aⅴ小说| 又大又爽又粗| 老司机深夜福利视频在线观看 | 日本av免费视频播放| 国产精品免费视频内射| 国产亚洲精品久久久久5区| 丝袜美足系列| 久久久水蜜桃国产精品网| 青青草视频在线视频观看| 午夜福利乱码中文字幕| 老司机午夜十八禁免费视频| 999精品在线视频| 亚洲国产欧美一区二区综合| 女性被躁到高潮视频| 亚洲人成77777在线视频| 脱女人内裤的视频| 国产91精品成人一区二区三区 | 99久久国产精品久久久| 12—13女人毛片做爰片一| 69精品国产乱码久久久| 麻豆乱淫一区二区| 精品福利永久在线观看| 老司机亚洲免费影院| 91av网站免费观看| 日韩免费高清中文字幕av| 免费高清在线观看日韩| 少妇猛男粗大的猛烈进出视频| 国产一区二区三区综合在线观看| 国产精品一区二区在线不卡| 美女午夜性视频免费| 亚洲成国产人片在线观看| 欧美+亚洲+日韩+国产| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 中文欧美无线码| 国产日韩欧美视频二区| 久久精品成人免费网站| 亚洲国产精品成人久久小说| 香蕉丝袜av| 午夜福利影视在线免费观看| 久久精品人人爽人人爽视色| 国产av国产精品国产| 欧美日本中文国产一区发布| 色视频在线一区二区三区| 亚洲av美国av| 黑人欧美特级aaaaaa片| kizo精华| 高清在线国产一区| 成年动漫av网址| 国产男人的电影天堂91| 日韩视频在线欧美| 极品人妻少妇av视频| 纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| 一边摸一边做爽爽视频免费| 最新的欧美精品一区二区| 一级黄色大片毛片| 麻豆av在线久日| 九色亚洲精品在线播放| 在线看a的网站| a级毛片黄视频| 国产精品 国内视频| 亚洲精品第二区| 免费高清在线观看视频在线观看| 777久久人妻少妇嫩草av网站| a级毛片在线看网站| 午夜激情av网站| 大型av网站在线播放| av免费在线观看网站| 高清黄色对白视频在线免费看| 日韩 亚洲 欧美在线| 男女边摸边吃奶| 亚洲欧美日韩另类电影网站| 国产成人影院久久av| 国产精品久久久久久精品古装| 成人手机av| 亚洲精品中文字幕在线视频| 欧美日韩一级在线毛片| 韩国精品一区二区三区| 久久久精品国产亚洲av高清涩受| av线在线观看网站| 亚洲av男天堂| 巨乳人妻的诱惑在线观看| 欧美日韩亚洲国产一区二区在线观看 | 久久久久久人人人人人| 国产精品成人在线| 秋霞在线观看毛片| 色精品久久人妻99蜜桃| 狂野欧美激情性bbbbbb| 桃红色精品国产亚洲av| 国产无遮挡羞羞视频在线观看| 国产熟女午夜一区二区三区| 12—13女人毛片做爰片一| 99热国产这里只有精品6| 久久精品国产a三级三级三级| 日本wwww免费看| 国产野战对白在线观看| 91麻豆av在线| avwww免费| 搡老乐熟女国产| 91精品三级在线观看| 国产1区2区3区精品| 午夜福利免费观看在线| 欧美人与性动交α欧美软件| 国产成人精品久久二区二区91| 国产精品免费视频内射|