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

    Joint Optimization of Resource Allocation and Trajectory Based on User Trajectory for UAV-Assisted Backscatter Communication System

    2024-03-11 06:28:42PeizhongXieJunjieJiangTingLiYinLu
    China Communications 2024年2期

    Peizhong Xie ,Junjie Jiang ,Ting Li ,Yin Lu

    1 School of Communication and Information Engineering,Nanjing University of Posts and Telecommunications,Nanjing 210003,China

    2 Key Lab of Broadband Wireless Communication and Sensor Network Technology(Nanjing University of Posts and Telecommunications),Ministry of Education,Nanjing 210003,China

    Abstract: The Backscatter communication has gained widespread attention from academia and industry in recent years.In this paper,A method of resource allocation and trajectory optimization is proposed for UAV-assisted backscatter communication based on user trajectory.This paper will establish an optimization problem of jointly optimizing the UAV trajectories,UAV transmission power and BD scheduling based on the large-scale channel state signals estimated in advance of the known user trajectories,taking into account the constraints of BD data and working energy consumption,to maximize the energy efficiency of the system.The problem is a non-convex optimization problem in fractional form,and there is nonlinear coupling between optimization variables.An iterative algorithm is proposed based on Dinkelbach algorithm,block coordinate descent method and continuous convex optimization technology.First,the objective function is converted into a non-fractional programming problem based on Dinkelbach method,and then the block coordinate descent method is used to decompose the original complex problem into three independent sub-problems.Finally,the successive convex approximation method is used to solve the trajectory optimization sub-problem.The simulation results show that the proposed scheme and algorithm have obvious energy efficiency gains compared with the comparison scheme.

    Keywords: energy efficiency;joint optimization;UAV-assisted backscatter communication;user trajectory

    I.INTRODUCTION

    With the appearance and development of the Internet of Things and its large-scale commercial use,the energy problem of node batteries has become a major factor restricting its development,and the backscatter communication technology came into being[1-3].The backscatter device does not need to use any complex and high-power active RF wave components.It can not only reflect the incident RF wave by adjusting the impedance of the device to achieve information transmission,but also collect energy from the incident RF wave[4,5].This will reduce the power consumption of the node circuit and extend its service time,thus providing an energy-saving solution for the communication of the Internet of Things devices[6].

    Backscattering communication systems can be divided into the following categories: single-static backscattering communication systems,double-static backscattering communication systems and environmental backscattering communication systems [7].The single static backscatter communication system is that the carrier wave transmitter and signal receiver are set on the same equipment,the received backscatter signal has experienced two signal fading,and is also affected by the double near-far problem.The carrier transmitter and signal receiver of the dual-static backscatter communication system are separate.Selecting the appropriate location of the carrier transmitter can greatly improve the performance of the dual-static backscatter communication system.However,the double-static backscatter carrier transmitter is bulky and expensive to deploy,and will cover huge economic costs in the actual production environment.The carrier transmitter used in the environmental backscatter communication system is the environmental RF source with other functions and utilization value,including a TV tower,cellular base station,Wi-Fi signal,etc.Since the modulated environmental RF signal is a unpredictable dynamic signal,it will directly interfere with the signal receiver,which greatly limits the performance of the environmental backscatter communication system[8-10].

    Unmanned aerial vehicle (UAV) assisted wireless communication has gained widespread attention from academia and industry in recent years benefiting from its wide coverage,high mobility,low cost and high probability of line-of-sight links[11-13],Lei Lang et al.[14]maximized the minimum average rate of users by optimizing the flight path and bandwidth allocation of UAVs in the downlink wireless transmission system of UAV-assisted mobile user communication.In[15],the author put forward a UAV broadcast coverage strategy that can maximize energy efficiency.

    Considering the disadvantages of the above backscatter communication system,the UAV is introduced into the backscatter communication network.In [16],the author studied UAV assisted environment backscatter system,analyzed the outage probability of downlink backscatter communication and derived the asymptotic outage probabilities of infinite RF signal power and infinite Nakagami fading parameter.In [17],the author studied UAV assisted NOMA backscatter communication system,optimized the scheduling of energy collection time,the gateway amplification factor,and the location of the UAV to maximize the sum rate of system.In[18],the author proposed a backscatter communication network with a UAV equipping a carrier wave transmitter and a signal receiver at the same time.The UAV collects the BDs data and transmits it back to the base station.The throughput of the backscatter communication system is maximized by optimizing the trajectory,time distribution and backscatter coefficient of the UAV.Authors in [19] and [20] considered that UAVs are equipped with a carrier wave transmitter and a signal receiver at the same time,and a non-orthogonal multiple access communication scheme is selected to transmit information.In [19],authors maximized the number of bits successfully decoded in the uplink and minimized the flight time by optimizing the flight altitude of UAVs.In [20],authors optimized the hover position,reflection coefficient,and transmission power of UAV to maximize the total rate of nodes.Fan et al.[21]considered that multiple UAVs are used for time division multiple access communication,and the average node rate is maximized by optimizing UAV trajectory,reflection coefficient,and transmission power.

    The UAV-assisted backscatter communication network has shown better performance improvement in communication quality,but UAV battery energy consumption cannot be ignored.Therefore,issues related to energy consumption have been studied.In [22],authors considered that the scenario of UAV-assisted backscatter and HTT hybrid network model.UAV provides energy for the hybrid network.By optimizing the transmission power and trajectory of the unmanned aerial vehicle,the backscatter reflection coefficient of the Internet of Things nodes and the transmission power of active transmission,the total energy efficiency of all Internet of Things nodes is improved.However,the energy consumption considered in this paper is the BDs and active transmission,The flight energy consumption of UAV is ignored which is the largest power consumption.In [23],authors considered that The UAV flew from the base station to the hover position to collect data,and then flew back to the base station to transmit data.The energy efficiency was maximized by determining the optimal hover position.its energy consumption is UAV flight energy consumption and hover energy consumption.Since UAV flight power and hover power are fixed values during constant speed flight and hover,it is sufficient to study flight time and hover time.In [24],authors used the energy consumption model of UAV dynamic flight speed to calculate energy consumption accurately,which considers the impact of trajectory and speed on flight energy consumption.This paper studies the backscatter communication system that UAV is with a signal receiver and the position of the carrier transmitter is fixed,authors improved the system energy efficiency by optimizing the trajectory of UAV,node scheduling and transmission power of the carrier transmitter.

    At present,the research on UAV-assisted backscatter communication systems is based on the fixed position of BDs,while BDs may be mobile in reality.For mobile BDs,it is difficult to obtain instantaneous channel state information dynamically,but when the BD dynamic path is predictable,large-scale channel state information can be obtained in advance,and it can solve the problem effectively.

    In this paper,considering the backscatter communication system model of UAV equipped with a carrier wave transmitter and a signal receiver,and the moving path of ground backscatter equipment can be predicted,the trajectory and resource allocation of UAV are researched to improve the energy efficiency of the system.Firstly,the UAV transmits radio frequency carrier wave to the moving ground backscatter equipment,and then the BDs use part of the received sine signal for power supply and part of the sine signal for modulation to send to the UAV,and finally the UAV receives the signal to complete data acquisition.Setting the minimum energy consumption requirements of the ground reflecting equipment and the minimum data collection limit,the energy efficiency of the system is maximized by jointly optimizing the trajectory of the UAV,the launch power of the UAV and the mobile BD scheduling.When designing the optimization problem solving algorithm,the original fractional optimization problem is first converted into an equivalent integer problem based on Dinkelbach algorithm,and then the block coordinate descent method is used to decompose into sub-problems,and finally the successive convex approximation algorithm is used to solve the sub-problems.

    The main contributions are as follows:

    1) Since the backscatter equipment is mobile,a UAV-assisted backscatter communication system model consisting of a UAV that equips a carrier wave transmitter and a receiver and multiple mobile BDs is adopted.An optimization objective function is constructed to maximize the total energy efficiency of the system and meet the minimum energy collection constraints,minimum transmission data volume constraints,and maximum transmission power constraints.The problem is a non-convex,nonlinear,multivariable coupled optimization problem,and it is difficult to obtain the analytical solution directly.

    2) This paper will propose the DBS algorithm(Dinkelbach-Block Coordinate Descent-Successive Convex Approximations)to solve the problem.Firstly,the objective function is converted into a nonfractional programming problem based on Dinkelbach method,and then the original problem is converted into three sub-optimization problems to solve the multivariable coupling problem using the Block Coordinate Descent method;Finally,the non-convex optimization problem is converted into a convex optimization problem by using the Successive Convex Approximations algorithm.This paper will introduce the specific algorithm steps.

    3)This paper will build the simulation environment to verify the feasibility of the scheme.The simulation results show that the proposed scheme has obvious energy efficiency gain compared with the comparison scheme,and the proposed algorithm can effectively solve the problem and has good convergence.

    II.SYSTEM MODEL

    As shown in Figure 1,This UAV-assisted backscatter communication system includes K mobile BDs and a UAV with a carrier wave transmitter and a signal receiver,and K BDs are represented byKall ∈{1,2,...,N.}

    Figure 1.The UAV-assisted backscatter communication system.

    To facilitate analysis,a two-dimensional coordinate system is introduced into the system model.Setting the task time toT.In timet ∈[0,T],because the dynamic trajectory of BD is known,the k-th BD coordinate is defined aswk(t)=[xk(t),yk(t)],wherexk(t)andyk(t)respectively represent the abscissa and ordinate of the k-th BD at time t.The plane coordinates of the UAV are defined asq(t)∈[x(t),y(t)],wherexk(t) andyk(t) respectively represent the abscissa and ordinate of the UAV at timet,and the fixed flight height of the UAV is set asH.To simplify the analysis,the flight time is divided intoNtime slots that are small enough,and the duration of each time slot is represented byNall ∈{1,2,...,N}.The dynamic trajectory of BD can be approximately expressed as,n ∈{1,2,...,N}.The trajectory of the UAV can be approximately expressed as,n ∈{1,2,...,N},and the maximum speed of the UAV is expressed asVmax,so there are constraints on the position and speed of the UAV.

    where ||·||represents the two-norm calculation symbol,δ=T/Nrepresents the length of each time slot,qIrepresents the initial position,andqFrepresents the final position.

    Unlike the traditional terrestrial cellular network communication,the UAV can flexibly change its position,so the channel from the backscatter equipment BD to the UAV is dominated by the line-of-sight link.In the n-th time slot,the distance between the UAV and the k-th BD can be expressed as

    In this paper,the time division multiple access(TDMA)dynamic transmission strategy is adopted to provide a communication scheme for the UAV auxiliary backscatter communication system,that is,the moving BD sends data to the UAV in turn,and the UAV flies according to the designed trajectory,and receives data from the BDs during the flight.Define a binary variableak[n]∈{0,1}to indicate whether the k-th BD sends data to the UAV in the n-th time slot.Ifak[n]=1,the k-th BD will send data;Otherwiseak[n]=0.At most one BD establishes communication with UAV at the same time,so there is the following constraint.

    BDs can change the matching impedance of the antenna to modulate the data to the incident wave,and make the modulated signal to transmit to the UAV in the way of antenna reflection.The intensity of the reflected signal is related to the power reflection coefficient of BD and the transmitted power of the UAV.To simplify the model,whenak[n]=1,BD is in the state of total reflection,that is all energy can be used for reflected communication when the UAV establishes communication with the k-th node;Whenak[n]=0,BDs are in the state of energy collection.The carrier transmitter has power constraints.

    wherePmaxrepresents maximum transmitting power of UAV.

    In the process of energy collection,the backscatter device BD must have some energy loss due to internal circuit reasons.The wireless energy collection efficiency of the k-th BD is expressed byηk ∈(0,1).Therefore,the energy collected by the k-th BD in the whole mission cycle can be expressed as

    where Ω0represents the path loss when the reference distance is one meter.

    According to Shannon’s theory,the amount of data transmitted by the k-th BD in the whole task cycle can be expressed as

    whereβ0represents the channel power gain when the reference distance is one meter,andσ2represents the additive Gaussian white noise at the UAV receiver.

    In the UAV-assisted backscatter communication system,the total energy consumption of the system mainly includes carrier transmitter energy consumption and UAV flight energy consumption.The flight propulsion energy consumption of rotorcraft is determined by the flight speed and acceleration of UAV.In this paper,because the variable speed flight time of UAV only accounts for a small part of the total flight time of UAV,the additional energy consumption caused by UAV acceleration is ignored.The flight propulsion energy consumption power of UAV[14]is expressed as

    Therefore,the energy efficiency of the UAV-assisted backscatter communication system can be expressed as

    From the perspective of energy consumption,the optimization problem of maximizing the energy efficiency of UAV assisted backscatter communication system is established.The goal is to maximize the energy efficiency of the system through BD scheduling A,carrier transmission power P,UAV trajectory Q and UAV speed V,while ensuring the minimum communication throughput and working energy consumption of BD.The optimization problem(10)is as follows

    III.PROPOSED ALGORITHM

    3.1 Conversion of Equivalent Optimization Problems

    The objective function of the original optimization problem has uncertain parameters and is difficult to deal with.Because Dinkelbach method is widely used to deal with nonlinear fractional optimization problems,this paper applies this method to deal with the energy efficiency objective function.With the aid of auxiliary variablesη*,the original objective function can be equivalently converted into

    whereη*is the optimal energy efficiency.

    For any given non-negative valueη ≥0,iff(η*)≤0,η ≥η*,it means that the current value needs to be reduced to make it tend to the optimal energy efficiency under the condition of ensuring(11).Similarly,iff(η*)≥0,η ≤η*,it means that under the condition of ensuring(11),it is necessary to increase the current value to make it tend to the optimal energy efficiency.Onlyη=η*then can we obtain the optimal energy efficiency.Therefore,the original optimization problem can be transformed into the following equivalent parameter optimization problem(12).

    It is still a non-convex optimization problem because of the coupling relationship between optimization variables.For such problem,the block coordinate descent method is used to decompose it into three sub-optimization problems: scheduling optimization,power optimization,and trajectory optimization.

    3.2 Scheduling Optimization

    For the scheduling optimization subproblem,given the feasible carrier transmission power,flight path,speed,and auxiliary variables,the original problem is simplified to solve the scheduling of UAV and backscatter equipment.Since the uplink access variable of BD is a binary variable and needs to be processed continuously,the equivalent parameter problem of energy efficiency optimization(13)can be expressed as

    Obviously,the optimization problem is a linear programming problem,which can be solved by existing optimization tools(such as CVX).

    3.3 Power Optimization

    For the power optimization subproblem,given the feasible scheduling variables,flight path,speed,and auxiliary variables,the original problem is simplified to solve the carrier transmission power of the optimal UAV,and the energy efficiency optimization equivalent parameter problem(14)can be expressed as

    The objective function is a convex function,and the constraint conditionsC8 andC9 are convex functions,so the optimization problem is a convex optimization problem,which is solved by using convex optimization tools.

    3.4 Trajectory Optimization

    For the trajectory optimization subproblem,given feasible scheduling variables,launch power,and auxiliary variables,the original problem is simplified to solve the optimal UAV trajectory Q and the flight speed V varying with the trajectory.The equivalent parameter problem of energy efficiency maximization can be expressed as (15).The optimization problem has non-convexity and is difficult to solve because of nonconvex objective function and non-convex constraint termsC10,C11.The non-convex optimization problem needs to be converted into a convex optimization problem.It can be seen from formula(8)that the first term and the third term ofPUAV[n] are convex functions ofV[n],but the second term is non-convex.The non-convexity is solved by introducing relaxation variables{y[n]}.The process is as (16) while this formula can be equivalently transformed as(17).The optimization problem can be further obtained by substituting formula(18)into the optimization problem.

    However,the optimization problem is still nonconvex because there are non-convex terms in the objective function,and the constraintsC10,C11 andC12 are non-convex.These four non-convex terms can be solved by continuous convex optimization technology.Specifically,the objective function and the constraint conditionC10 are convex functions ofH2+||wk[n]-q[n]||2.Because the first order Taylor expansion function of the convex function is the lower bound of the original function,the first order Taylor expansion obtains the lower bound of the original function.

    The non-convex portion of the optimization problem is locally linearized using a first-order Taylor series expansion.The equation is considered valid when the expansion point within the local region corresponds to a local optimum.

    For the non-convex constraint condition (18),the expression on the right side of the equal sign is a joint convex function of variablesy[n]andV[n].The lower bound of the expression on the right side of the constraint conditionC12 can be obtained by using the continuous convex optimization technique.

    where{y{l}[n]}is obtained after thel-1 iteration,and{y{l}[n]}is a constant.The lower bound of inequality(20)is a linear function of variablesy[n]andV[n].The equation is satisfied wheny{l}[n]andV {l}[n]are local optimums.For constraint conditionC6,its lower bound function is obtained by Taylor first-order expansion

    Replace the logarithmic function of the objective function (18) and the constraint conditionC10 with the corresponding lower bound according to formula (19);Replace the right part of the constraint conditionC12 with its lower bound according to formula(20);According to formula(21),replace the left function of constraint conditionC11 with its lower bound,and then problem(18)is optimized to(22)by continuous convex optimization technology.

    The global optimal solution of problem (22) is the lower bound of the global optimal solution of problem(15).First,the feasible region of problem (22) is a subset of the feasible region of problem(15),because once the constraintsC14,C15 andC16 are satisfied,the constraintsC10,C11 andC12 must be satisfied.Secondly,because problem (22) uses the continuous convex optimization technique to take the lower bound of the objective function of problem(15),the optimal energy efficiency value of problem (22) must not exceed the value of problem(15).

    The proposed DBS algorithm is summarized and listed in Algorithm 1.The number of iterations for the DBS algorithm is defined asKite,according to literature [25],when the convergence threshold is?,The computational complexity of the Dinkelbach algorithm isO(1/?2log(3K2N3)).Set the convergence threshold of the interior point method toζ,the complexity of each complex iteration of steps 2,3,and 4 isO(log(1/ζ)).In summary,the computational complexity of DBS isO(Kite·1/?2log(3K2N3)·3log(1/ζ)).

    IV.SIMULATION RESULTS

    The UAV-assisted backscatter communication systems work at the frequency of 900MHz.It is assumed that there are four users moving along the known track in the square area,as shown in Figure 2.The starting coordinates of user 1,user 2,user 3,and user 4 are(50,50),(50,450),(450,450),(450,50),and (450,50)respectively moving in a uniform straight line towards(50,450),(450,450) and (450,50),and starting at the same time.The maximum power of the carrier transmitter is set asPmax=6W.Each node must meet the lower energy limit with=1*1014J.The flight altitude of the UAV is set asH=50m,its maximum flight speed is set asVmax=20m/s.The other parameters as follow:β0=0.1,η=0.5,σ2=-144dBm.This paper divides the flight time into 100 time slots.The parameters of rotorcraft are set according to the parameters in reference [24].Table 1 lists the values of various parameters involved in the flight dynamics model of rotorcraft.

    Table 1.Parameters of UAV.

    Figure 2.The user trajectory.

    The simulation first considers the trajectory optimization of UAV,where Figure 3 and Figure 4 show the flight path optimization results of UAV under Case 1 and Case 2,taking into account two situations of UAV at different starting points and ends and at different mission times.

    Figure 3.UAV optimal trajectory in Case 1.

    Figure 4.UAV optimal trajectory in Case 2.

    Case 1.The starting point and ending point of the UAV are set as(50,50).

    Case 2.The starting point and ending point of the UAV are set to(250,250). At the same time,adjust the known trajectory of user 1 to a uniform linear motion from(50,450)to(50,50).

    It can be seen from Figure 3 and Figure 4 that although the starting and ending positions of the UAV are different,and the user trajectory is slightly adjusted,the optimal trajectory of the UAV in both cases is as close as possible to the user trajectory,and the distance between the UAV and the user is shortened,which indicates that the optimized trajectory is closer to the user,and better large-scale channel status is obtained,and larger data transmission volume is obtained.In addition,when the task time is changed,the optimal flight path is also different.With the increase of the task time,the complexity of the optimal trajectory of the UAV also increases,and the optimal trajectory is also closer to the user’s trajectory.This shows that the longer the mission time,the larger the flight path optimization space,and the better the trajectory optimization effect.

    Figure 5 shows the time resources allocated by the UAV to four users in Case 1.Figure 6 shows the reflected power distribution of four users when the task time is set to 80 s.At the same time,it can be seen from the optimal trajectory of the UAV in Figure 3 that the user’s scheduling depends on the distance between the UAV and the user,and the UAV will establish communication transmission with the user closest to the current time.Observe Figure 3 and Figure 5.With the increase of task time,the trajectories close to User 1 and User 4 become longer and the scheduling time also becomes longer,indicating that under the condition of meeting the data volume constraints of User 2 and User 3,increasing the scheduling time of User 1 and User 4 can achieve better energy efficiency.

    Figure 5.Time resource allocation in Case 1.

    Figure 6.User scheduling distribution in Case 1.

    Figure 7 shows the time resources allocated by the UAV to four users in Case 2.Figure 8 shows the reflected power distribution of four users when the task time is set to 80 s.Observe Figure 4 and Figure 8.With the increase of task time,the trajectories close to user 2,user 3 and user 4 become longer and the scheduling time also becomes longer,indicating that under the condition of meeting the data volume constraint of user 1,increasing the scheduling time of other users can achieve better energy efficiency.

    Figure 7.Time resource allocation in Case 2.

    Figure 8.User scheduling distribution in Case 2.

    To verify the convergence of the algorithm,Figure 9 shows the relationship between different task times and convergence.It can be observed that when the task time is 60s,the energy efficiency will converge after 6 iterations;When the task duration is 80 s,it will converge after 10 iterations;When the task time is 100s,it will converge after 12 iterations.The experimental results show that the convergence decreases slightly with the increase of task duration.At the same time,it was found that under the initial trajectory setting,the longer the task time,the lower the energy efficiency.This is due to the relatively unreasonable initial trajectory,the longer the mission time increases the flight energy consumption,but does not improve the quality of data transmission;however,when the algorithm converges,the longer the task time is,the higher the energy efficiency is.This is because the optimal trajectory improves the communication transmission quality on the one hand,and optimizes the UAV energy consumption on the other hand.

    Figure 9.Effect of different task times on convergence.

    Finally,the performance of the methods used in this paper is evaluated and two sets of comparison schemes are proposed.Scheme 1 is based on the design presented in[24],which is different from the setup in reference[24]that employs multiple carrier transmitters.Considering the cost implications in practical applications,this paper only utilizes a single carrier transmitter,aligning with the assumptions made in this study.A high-power carrier transmitter is set at(250,250)to achieve optimal energy efficiency by optimizing BD scheduling,UAV transmission power and UAV trajectory.Scheme 2 represents a scenario without trajectory optimization for the unmanned aerial vehicle(UAV),assuming stationary hovering of the UAV.This approach is consistent with the common practice observed in UAV-assisted communication systems.It can be seen from Figure 10 that the energy efficiency of the proposed scheme in this paper is the best.The energy efficiency of the second scheme decreases with the increase of time.This is because the hovering flight of the UAV is a state of high energy consumption,and the growth rate of energy consumption over time exceeds the rate of data volume.

    Figure 10.Comparison of energy efficiency of different schemes.

    V.CONCLUSION

    This paper researches UAV-assisted backscatter communication systems,which is composed of multiple mobile BDs and a UAV with carrier transmitter and receiver.Firstly,by jointly optimizing BD scheduling,UAV launch power and UAV trajectory,the system model of energy efficiency optimization is established under the constraints of BD throughput and working energy consumption,and then a DBS algorithm based on Dinkelbach algorithm,block coordinate descent algorithm and continuous convex optimization technology is proposed.Because the original problem is a nonlinear fractional optimization problem,Dinkelbach algorithm transforms the original problem into an equivalent integral optimization problem;The block coordinate descent method decomposes the equivalent integer optimization problem into three subproblems and solves three coupling problems of solving variables;The successive convex approximation technology completes the approximate solution of nonconvex subproblems of trajectory optimization.Finally,the system is simulated,and the results show that the energy efficiency of the proposed scheme is significantly higher than that of the original scheme,and the proposed algorithm has better convergence.Under the research background of this paper,the choice of communication scheme using NOMA is also a problem worth studying.

    一级a爱视频在线免费观看| tube8黄色片| 视频在线观看一区二区三区| 一级片免费观看大全| 久久精品国产亚洲av高清一级| 亚洲,欧美精品.| 999久久久国产精品视频| 国产精品熟女久久久久浪| 亚洲欧洲精品一区二区精品久久久 | 国产精品秋霞免费鲁丝片| 亚洲中文av在线| 国产成人精品一,二区| 99热全是精品| 久久99热这里只频精品6学生| 国产精品三级大全| av又黄又爽大尺度在线免费看| 18在线观看网站| 亚洲,欧美,日韩| 99久久人妻综合| 国产片特级美女逼逼视频| 中文字幕最新亚洲高清| 又大又黄又爽视频免费| 国产精品av久久久久免费| 久久这里有精品视频免费| 国产无遮挡羞羞视频在线观看| 色婷婷久久久亚洲欧美| 久久综合国产亚洲精品| xxx大片免费视频| 国产男女内射视频| 欧美少妇被猛烈插入视频| 日韩av免费高清视频| 伦理电影大哥的女人| 成人午夜精彩视频在线观看| 国产亚洲av片在线观看秒播厂| 少妇精品久久久久久久| 美女午夜性视频免费| 午夜免费观看性视频| 国产xxxxx性猛交| 飞空精品影院首页| 中文字幕制服av| 97在线人人人人妻| 亚洲第一av免费看| 日韩精品免费视频一区二区三区| 国产亚洲一区二区精品| 可以免费在线观看a视频的电影网站 | 超碰成人久久| 美女国产高潮福利片在线看| 黄色视频在线播放观看不卡| 国产成人精品婷婷| 中文字幕人妻熟女乱码| 欧美日韩综合久久久久久| 人妻 亚洲 视频| 国产精品麻豆人妻色哟哟久久| 国产精品香港三级国产av潘金莲 | 亚洲熟女精品中文字幕| 高清黄色对白视频在线免费看| 人成视频在线观看免费观看| 亚洲精品aⅴ在线观看| 久久久久久久久久久免费av| 久久人妻熟女aⅴ| 天天躁夜夜躁狠狠久久av| 亚洲人成电影观看| 亚洲熟女精品中文字幕| 精品久久久久久电影网| 26uuu在线亚洲综合色| 国产欧美亚洲国产| 另类亚洲欧美激情| 黄色一级大片看看| 波野结衣二区三区在线| 中文字幕人妻丝袜制服| 亚洲视频免费观看视频| 国产日韩欧美亚洲二区| 精品国产超薄肉色丝袜足j| 五月天丁香电影| 精品少妇黑人巨大在线播放| 狂野欧美激情性bbbbbb| 中文字幕精品免费在线观看视频| 国产成人a∨麻豆精品| 日本爱情动作片www.在线观看| 色播在线永久视频| 丝袜美足系列| 午夜免费鲁丝| 各种免费的搞黄视频| 男人操女人黄网站| 日韩av在线免费看完整版不卡| 国产熟女欧美一区二区| 在线天堂中文资源库| 亚洲男人天堂网一区| 亚洲久久久国产精品| 成人免费观看视频高清| 色网站视频免费| 亚洲精品美女久久久久99蜜臀 | 伊人亚洲综合成人网| 汤姆久久久久久久影院中文字幕| 最新中文字幕久久久久| av免费在线看不卡| 亚洲精品久久久久久婷婷小说| 男的添女的下面高潮视频| 欧美精品一区二区免费开放| 哪个播放器可以免费观看大片| 欧美人与善性xxx| 国产一区二区激情短视频 | 热re99久久精品国产66热6| 欧美精品亚洲一区二区| 在线观看www视频免费| 久久精品国产亚洲av天美| av在线老鸭窝| 久久这里有精品视频免费| 丝袜喷水一区| 男人舔女人的私密视频| 精品国产超薄肉色丝袜足j| 蜜桃国产av成人99| 国产精品久久久久久精品电影小说| 侵犯人妻中文字幕一二三四区| 97精品久久久久久久久久精品| 在线观看www视频免费| 777久久人妻少妇嫩草av网站| 久久99热这里只频精品6学生| 高清不卡的av网站| 一级片'在线观看视频| 亚洲国产av影院在线观看| 亚洲熟女精品中文字幕| 大香蕉久久成人网| 最近最新中文字幕免费大全7| 欧美精品av麻豆av| 99久久综合免费| 一区二区三区乱码不卡18| 国产片内射在线| 欧美bdsm另类| 欧美中文综合在线视频| 毛片一级片免费看久久久久| 国产精品免费视频内射| 一本久久精品| 黄色怎么调成土黄色| 久久久久久久精品精品| 黄片播放在线免费| 黄色怎么调成土黄色| xxxhd国产人妻xxx| 亚洲精品av麻豆狂野| 91精品国产国语对白视频| 嫩草影院入口| 一区二区三区乱码不卡18| 99久国产av精品国产电影| 999精品在线视频| 欧美成人精品欧美一级黄| 九色亚洲精品在线播放| 日本欧美国产在线视频| 天天躁日日躁夜夜躁夜夜| 在线观看美女被高潮喷水网站| 欧美精品av麻豆av| 国产精品三级大全| 亚洲内射少妇av| 啦啦啦在线观看免费高清www| 日韩大片免费观看网站| 成人国产av品久久久| 夫妻午夜视频| 精品午夜福利在线看| 成年女人毛片免费观看观看9 | 极品少妇高潮喷水抽搐| 老女人水多毛片| 在现免费观看毛片| 2021少妇久久久久久久久久久| 婷婷色综合www| 一级爰片在线观看| 午夜免费观看性视频| 亚洲人成77777在线视频| 制服丝袜香蕉在线| 寂寞人妻少妇视频99o| 这个男人来自地球电影免费观看 | 久久人妻熟女aⅴ| 青草久久国产| 免费在线观看黄色视频的| 亚洲精品日韩在线中文字幕| 天天躁夜夜躁狠狠久久av| 国产极品天堂在线| 日韩一卡2卡3卡4卡2021年| 97在线视频观看| 男男h啪啪无遮挡| 亚洲欧美一区二区三区国产| 2021少妇久久久久久久久久久| 在现免费观看毛片| 男人爽女人下面视频在线观看| 国产日韩欧美视频二区| 丝袜脚勾引网站| 如何舔出高潮| 国产伦理片在线播放av一区| 成人毛片a级毛片在线播放| 一本色道久久久久久精品综合| 毛片一级片免费看久久久久| 热re99久久精品国产66热6| 午夜福利视频在线观看免费| 亚洲精品一区蜜桃| 男人舔女人的私密视频| 免费在线观看视频国产中文字幕亚洲 | 国产视频首页在线观看| 日韩不卡一区二区三区视频在线| 久久99蜜桃精品久久| 一级片免费观看大全| 制服丝袜香蕉在线| 制服丝袜香蕉在线| 久久久久久久久久久免费av| 男人添女人高潮全过程视频| 亚洲欧美精品综合一区二区三区 | 丝袜美足系列| 国产日韩欧美视频二区| 成年女人在线观看亚洲视频| 国产午夜精品一二区理论片| 黄片小视频在线播放| 久久久久国产网址| 日韩熟女老妇一区二区性免费视频| 色网站视频免费| 国产福利在线免费观看视频| 香蕉国产在线看| 观看av在线不卡| 极品人妻少妇av视频| 国产精品 欧美亚洲| av在线播放精品| 免费看av在线观看网站| 国产在线视频一区二区| 男的添女的下面高潮视频| 波多野结衣一区麻豆| 精品国产一区二区久久| 中国三级夫妇交换| 黑人猛操日本美女一级片| 在线观看国产h片| 少妇猛男粗大的猛烈进出视频| 国产成人av激情在线播放| 色网站视频免费| 久久精品国产亚洲av高清一级| 一区在线观看完整版| 欧美少妇被猛烈插入视频| 国产免费一区二区三区四区乱码| 丝瓜视频免费看黄片| 亚洲,欧美,日韩| 欧美激情高清一区二区三区 | 精品国产乱码久久久久久男人| 日韩伦理黄色片| 日日爽夜夜爽网站| 蜜桃国产av成人99| 国产又爽黄色视频| 久久狼人影院| 久久久精品免费免费高清| 日本黄色日本黄色录像| 久久99一区二区三区| 国产极品粉嫩免费观看在线| 黑人欧美特级aaaaaa片| 97精品久久久久久久久久精品| 亚洲第一av免费看| 丝袜喷水一区| 少妇的逼水好多| 国产高清不卡午夜福利| 汤姆久久久久久久影院中文字幕| 亚洲综合色惰| 国产日韩一区二区三区精品不卡| 侵犯人妻中文字幕一二三四区| 成年动漫av网址| 亚洲精品aⅴ在线观看| 777米奇影视久久| 亚洲婷婷狠狠爱综合网| 人妻系列 视频| 99久久综合免费| 三级国产精品片| 亚洲国产毛片av蜜桃av| 熟女av电影| 最新的欧美精品一区二区| 国产视频首页在线观看| 亚洲精品成人av观看孕妇| 秋霞伦理黄片| 寂寞人妻少妇视频99o| 国产精品秋霞免费鲁丝片| 一本色道久久久久久精品综合| 一区二区三区四区激情视频| 欧美成人午夜精品| 国产精品不卡视频一区二区| 国产精品偷伦视频观看了| 久久午夜综合久久蜜桃| 一级毛片我不卡| 超碰97精品在线观看| 一级爰片在线观看| 国产一区二区激情短视频 | 中文字幕另类日韩欧美亚洲嫩草| 观看av在线不卡| 男女国产视频网站| 啦啦啦中文免费视频观看日本| 中文乱码字字幕精品一区二区三区| 久久这里有精品视频免费| 日韩人妻精品一区2区三区| 成年女人毛片免费观看观看9 | 十分钟在线观看高清视频www| 精品国产国语对白av| 亚洲欧美精品综合一区二区三区 | 国产又色又爽无遮挡免| 久久精品人人爽人人爽视色| 亚洲国产精品一区二区三区在线| av在线app专区| 亚洲国产精品一区三区| 最新中文字幕久久久久| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩精品成人综合77777| xxx大片免费视频| 久久99热这里只频精品6学生| av女优亚洲男人天堂| 男女无遮挡免费网站观看| av网站免费在线观看视频| 桃花免费在线播放| 国产高清国产精品国产三级| 精品国产超薄肉色丝袜足j| 大香蕉久久网| 亚洲国产av新网站| 中文字幕av电影在线播放| 三上悠亚av全集在线观看| 国产探花极品一区二区| 亚洲国产av影院在线观看| 男男h啪啪无遮挡| 久久久久久人妻| 在线天堂中文资源库| 亚洲av男天堂| av电影中文网址| 欧美黄色片欧美黄色片| 高清黄色对白视频在线免费看| 18在线观看网站| 在线观看免费高清a一片| 国产在线一区二区三区精| 母亲3免费完整高清在线观看 | 人体艺术视频欧美日本| 自拍欧美九色日韩亚洲蝌蚪91| 日本av手机在线免费观看| 在线亚洲精品国产二区图片欧美| 精品人妻熟女毛片av久久网站| 99热国产这里只有精品6| 午夜福利影视在线免费观看| 久久精品熟女亚洲av麻豆精品| 丝袜美腿诱惑在线| 一区二区三区精品91| 一本—道久久a久久精品蜜桃钙片| 中国国产av一级| 亚洲国产欧美网| 一区二区三区激情视频| 久久av网站| 久久免费观看电影| 大码成人一级视频| videossex国产| 男女边吃奶边做爰视频| 看免费av毛片| 亚洲三区欧美一区| 日韩,欧美,国产一区二区三区| 波多野结衣一区麻豆| 成人18禁高潮啪啪吃奶动态图| 超碰97精品在线观看| 国产午夜精品一二区理论片| 欧美精品一区二区大全| 男女高潮啪啪啪动态图| www日本在线高清视频| 久久久久久久国产电影| 免费高清在线观看日韩| 美女国产视频在线观看| 欧美国产精品va在线观看不卡| 18禁动态无遮挡网站| 精品亚洲成a人片在线观看| 日韩欧美一区视频在线观看| 日韩在线高清观看一区二区三区| 亚洲av欧美aⅴ国产| 日韩av免费高清视频| www日本在线高清视频| 久久久精品免费免费高清| 婷婷色麻豆天堂久久| 成人影院久久| 国产片特级美女逼逼视频| 国产精品麻豆人妻色哟哟久久| 国产成人精品福利久久| 99热国产这里只有精品6| 国产成人a∨麻豆精品| 在线亚洲精品国产二区图片欧美| 精品一区二区三卡| 大陆偷拍与自拍| 中国三级夫妇交换| 亚洲国产最新在线播放| 亚洲图色成人| 日韩熟女老妇一区二区性免费视频| 99香蕉大伊视频| 国产一区二区 视频在线| 中文字幕av电影在线播放| 久久久久视频综合| 18禁裸乳无遮挡动漫免费视频| 亚洲av福利一区| 夫妻午夜视频| 久久这里只有精品19| av免费观看日本| 亚洲国产av新网站| xxx大片免费视频| 成人18禁高潮啪啪吃奶动态图| 午夜91福利影院| 热re99久久精品国产66热6| 国产精品国产三级国产专区5o| 免费在线观看视频国产中文字幕亚洲 | 超碰成人久久| 少妇人妻精品综合一区二区| 欧美精品高潮呻吟av久久| 90打野战视频偷拍视频| 亚洲美女视频黄频| 大片免费播放器 马上看| 国产亚洲精品第一综合不卡| 最近中文字幕2019免费版| 美女国产视频在线观看| 欧美 亚洲 国产 日韩一| 久久国产精品男人的天堂亚洲| 电影成人av| 日韩成人av中文字幕在线观看| 亚洲美女搞黄在线观看| 熟女电影av网| 国产片内射在线| 丝袜在线中文字幕| 国产精品av久久久久免费| 少妇的丰满在线观看| 国产成人aa在线观看| 伦精品一区二区三区| 欧美日韩精品成人综合77777| 9191精品国产免费久久| 91成人精品电影| 亚洲成人av在线免费| 自拍欧美九色日韩亚洲蝌蚪91| 天堂俺去俺来也www色官网| 国产探花极品一区二区| 国产97色在线日韩免费| 曰老女人黄片| 亚洲国产日韩一区二区| 欧美精品人与动牲交sv欧美| 欧美日韩亚洲国产一区二区在线观看 | 青春草视频在线免费观看| 丝袜人妻中文字幕| 中文字幕精品免费在线观看视频| 午夜福利一区二区在线看| 国产精品一国产av| 国产精品久久久久久av不卡| 新久久久久国产一级毛片| 国产免费福利视频在线观看| 国产精品免费大片| 最近2019中文字幕mv第一页| 国产成人午夜福利电影在线观看| av在线观看视频网站免费| 国产综合精华液| 久久久久久人人人人人| 国产白丝娇喘喷水9色精品| 日日撸夜夜添| www日本在线高清视频| 久热久热在线精品观看| 成年动漫av网址| 亚洲av电影在线进入| 又大又黄又爽视频免费| 久久久久久久精品精品| 国产 精品1| 精品人妻一区二区三区麻豆| 亚洲欧美成人综合另类久久久| 免费高清在线观看视频在线观看| 欧美日韩国产mv在线观看视频| 伦精品一区二区三区| 日韩av免费高清视频| 成年av动漫网址| 成人毛片60女人毛片免费| 美女xxoo啪啪120秒动态图| 国产深夜福利视频在线观看| 欧美xxⅹ黑人| 中国国产av一级| 亚洲av在线观看美女高潮| 精品人妻偷拍中文字幕| 日本av手机在线免费观看| 97在线人人人人妻| a级毛片在线看网站| 久久青草综合色| 涩涩av久久男人的天堂| 日韩一本色道免费dvd| 在线观看免费高清a一片| 久久婷婷青草| 99香蕉大伊视频| a级毛片黄视频| 日韩,欧美,国产一区二区三区| 2018国产大陆天天弄谢| av线在线观看网站| 超色免费av| 中文乱码字字幕精品一区二区三区| 成人亚洲精品一区在线观看| 久久人人爽av亚洲精品天堂| 久久久欧美国产精品| av又黄又爽大尺度在线免费看| 丁香六月天网| 国产精品一国产av| 水蜜桃什么品种好| 大陆偷拍与自拍| 两个人看的免费小视频| 久久狼人影院| 亚洲美女视频黄频| 麻豆精品久久久久久蜜桃| 最近中文字幕2019免费版| 美女xxoo啪啪120秒动态图| 国产又爽黄色视频| 中文欧美无线码| 免费日韩欧美在线观看| 国产一区二区激情短视频 | 国产极品粉嫩免费观看在线| 9热在线视频观看99| 国产成人精品在线电影| av一本久久久久| 亚洲第一av免费看| 日本wwww免费看| 欧美亚洲日本最大视频资源| 男的添女的下面高潮视频| 一级毛片我不卡| 亚洲av电影在线进入| 欧美日韩视频高清一区二区三区二| 国产一区有黄有色的免费视频| 哪个播放器可以免费观看大片| 少妇熟女欧美另类| 国产免费福利视频在线观看| 一区二区三区乱码不卡18| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 大码成人一级视频| 久久久久久人妻| 观看av在线不卡| 不卡av一区二区三区| 老司机影院毛片| 国产成人av激情在线播放| 亚洲国产精品999| 欧美精品av麻豆av| 日本91视频免费播放| 哪个播放器可以免费观看大片| 久久久久精品性色| 日本免费在线观看一区| 热re99久久国产66热| 亚洲国产日韩一区二区| 校园人妻丝袜中文字幕| 久热这里只有精品99| 少妇人妻久久综合中文| 建设人人有责人人尽责人人享有的| 久久免费观看电影| 中文乱码字字幕精品一区二区三区| 最近中文字幕2019免费版| 在线观看免费高清a一片| 国产精品二区激情视频| 欧美国产精品一级二级三级| 日韩电影二区| av片东京热男人的天堂| 亚洲国产欧美日韩在线播放| 这个男人来自地球电影免费观看 | 精品少妇久久久久久888优播| 国产成人午夜福利电影在线观看| 久久久久视频综合| 丝瓜视频免费看黄片| 日本av免费视频播放| 狂野欧美激情性bbbbbb| 亚洲欧美精品自产自拍| 王馨瑶露胸无遮挡在线观看| 久久久久精品人妻al黑| 色视频在线一区二区三区| 午夜免费男女啪啪视频观看| 日韩av免费高清视频| 亚洲欧美色中文字幕在线| 极品人妻少妇av视频| 久久久久久人妻| 观看av在线不卡| 另类精品久久| 日韩av在线免费看完整版不卡| 亚洲欧美一区二区三区国产| 永久免费av网站大全| 又黄又粗又硬又大视频| 久久鲁丝午夜福利片| 777米奇影视久久| 国产深夜福利视频在线观看| 日本wwww免费看| 国产无遮挡羞羞视频在线观看| 国产极品天堂在线| 国产精品久久久久久精品电影小说| 日本免费在线观看一区| 老司机亚洲免费影院| 亚洲第一区二区三区不卡| 国产精品熟女久久久久浪| 久久人妻熟女aⅴ| 777久久人妻少妇嫩草av网站| 纵有疾风起免费观看全集完整版| 成年人午夜在线观看视频| 叶爱在线成人免费视频播放| 好男人视频免费观看在线| 国产精品av久久久久免费| 久久 成人 亚洲| av有码第一页| 99久久综合免费| 丝袜美腿诱惑在线| www日本在线高清视频| 2022亚洲国产成人精品| 国产一级毛片在线| 少妇人妻久久综合中文| 青春草视频在线免费观看| 日本猛色少妇xxxxx猛交久久| 美女中出高潮动态图| 午夜激情久久久久久久| 日韩电影二区| www.熟女人妻精品国产| 精品人妻在线不人妻| 丰满乱子伦码专区| 亚洲精品一二三| 国产精品女同一区二区软件| 性高湖久久久久久久久免费观看| 啦啦啦视频在线资源免费观看| 女人久久www免费人成看片| 男男h啪啪无遮挡| 叶爱在线成人免费视频播放| 成年人午夜在线观看视频| 日韩人妻精品一区2区三区| 色哟哟·www| 日日啪夜夜爽| 一本色道久久久久久精品综合| 新久久久久国产一级毛片| 久久久久视频综合| 成年女人在线观看亚洲视频| 多毛熟女@视频| 伦精品一区二区三区|