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

    Dwell Scheduling Algorithm for Digital Array Radar

    2018-04-16 06:49:46QunZhangDiMengYingLuoandYijunChen

    Qun Zhang, Di Meng,, Ying Luo and Yijun Chen

    (1.Institute of Information and Navigation, Air Force Engineering University, Xi’an 710077, China; 2.Collaborative Innovation Center of Information Sensing and Understanding, Xi’an 710077, China)

    With the improvement of radar digitization,digital array radar (DAR) has been widely researched in the radar industry[1]. Compared with traditional analog phased array radar, DAR not only has the advantages of high detection precision, complex target detection capability, and strong anti-jamming capability, but also has the characteristics of flexible signal processing mode. It can be used to search a number of airspaces, and track and image multiple targets simultaneously—resulting in effective time resource management of radar[2].

    Reasonable, flexible and effective scheduling strategies are the key to exploiting the high adaptive potential of DAR. Among existing studies, there are mainly two ways to implement resource management of DAR: template-based scheduling and adaptive scheduling. Adaptive scheduling methods flexibly adjust the resource scheduling strategy according to the working environment and task requirements, and are the most effective but also the most complex scheduling methods[3].

    So far, many scholars have studied the radar resource scheduling model. In Ref.[4], an adaptive algorithm of phased array radar based on time window is proposed to verify the rationality and validity of introducing time window in the scheduling process. A resource allocation model based on “Quality of Service” is proposed in Ref.[5] to optimize the target tracking performance of phased array radar systems. In Ref.[6], an adaptive scheduling algorithm based on information entropy is proposed to solve the optimal searching problem in target regions. In order to further improve the source utilization ratio of the radar system and give full play to the advantages of multifunction digital array radar, pulse interleaving technology is proposed. Its core idea is to schedule the transmission or reception pulse of other tasks in a single task pulse interval. Orman analyzed this method in Ref.[7] and proposed a heuristic algorithm to solve the problem of adaptive beam-dwell scheduling for phased array radar. In Ref.[8], an optimal algorithm based on pulse interleaving for radar task quadratic programming model is proposed to improve the scheduling success ratio of tasks with high priority. Aiming at the problem of beam-dwell scheduling for DAR, an algorithm based on analyzing scheduling interval is proposed in Ref.[9]. However, in almost all existing adaptive resource scheduling methods for DAR, the imaging mission is not taken into account. Joint scheduling imaging tasks with search, tracking and other tasks not only improve the recognition ability of the radar to the target and give full play to the advantages of multi-task cooperation, but also feedback the target characteristic information obtained by the image to the transmitter radar system.This enables the realization of dynamic adjustment of imaging tasks and improves the resource utilization ratio[10]. Traditional imaging algorithms require continuous observation of targets for a long time to obtain high-resolution images, but the alternation of different tasks inevitably leads to discontinuous synthetic aperture sampling of the imaging target azimuth.

    Based on the above issues, in this paper we propose a dwell scheduling algorithm for DAR, which takes the imaging task requirements into account in the radar resource scheduling model, by using the compressed sensing (CS) ISAR imaging method.

    1 Task Model and Constraints Based on Pulse Interleaving

    DAR is a kind of all-digital phased array radar, which can receive and transmit the beam in a digital way. It can be interleaved in the form as shown in Fig.1.txj,twj,trjdenote the radar beam transmitting duration, waiting duration of a dwell, and the radar beam receiving duration respectively forj(j=1,2).With the pulse interleaving technology in DAR, time resources of the waiting duration can transmit or receive other resident tasks,and furthermore, the receiving duration of dwell tasks can overlap each other at the same time.

    Fig.1 Form of pulse interleaving in DAR

    1.1 Sparse aperture imaging based on compressed sensing

    The traditional DAR needs to separate part of the continuous time resources to achieve imaging capabilities in the implementation of the target searching and tracking, which leads to the contradiction of the radar resource allocation and low radar efficiency. Under the framework of CS theory, continuous observation of the target image can be transformed into random sparse observation, and high-quality target ISAR image can be obtained under sparse aperture condition. This provides an effective technical support for incorporating imaging task requirements into the DAR resource scheduling model.

    In order to improve the adaptive ability of radar imaging, the method proposed in Ref.[10] can be used to recognize the characteristics of the targets after entering the stable tracking phase, and then calculate the demand of radar resources for target imaging based on the cognitive results. We estimate the distancei, speedi, headingi, target sizei, sparsity of echo in Doppler domainiand observation timeciof thei-th target. Assuming that the target is non-maneuvering, the full aperture imaging of the radar needs to transmitNi=Pfcipulses (Pfis the pulse repetition frequency), and the measurement dimensionMi(Mi

    Mi≥cilnNi

    (1)

    wherecis a constant associated with the recovery accuracy, usually taken as 0.5-2(we assumec=1 in this paper).

    1.2 Radar dwell model

    The DAR dwell model is defined as

    T={et,st,tx,tw,tr,ω,M,pri,Pt,P}

    (2)

    whereetdenotes the expected scheduling time;stdenotes the actual scheduling time;tx,tw,trdenote the transmitting duration, waiting duration and receiving duration respectively;ωdenotes the time window of the task;Mdenotes the number of pulse repetitions in the search and tracking tasks, and denotes the azimuth-oriented observation dimension in the imaging task; pri denotes the pulse repetition reputation interval;Ptdenotes the pulse transmitting power; andPdenotes the task priority[13].

    1.3 Scheduling constraints

    For tracking and imaging tasks, assuming the distance from theith target to the radar isi, the waiting duration for theith task dwell can be calculated from the predicted position information of the target, which is given by

    twi=2i/c

    (3)

    For the search tasks, it is generally impossible to obtain the arrival time of the echo in the absence of prior information. Therefore, in order to ensure that radar echo signals can be effectively received after the searching pulse is transmitted, once the transmitting duration ends, the antenna system must be in the receiving state until the maximum waiting time, i.e. the waiting duration of the search task dwell cannot be preempted.

    It should be noted that the arrival time of radar echo is affected by the range size of the target.In order to ensure the imaging quality, the received pulse should be broadened properly. Suppose that the distance from thei-t target to the radar isiand the range size of the target isyi, the width of the actual received pulse of theith imaging task can be expressed as

    twi=2(i+yi)/c

    (4)

    Although the overlap of multiple dwell tasks can improve the resource utilization of DAR, the energy consumption of the system is bound to increase as the radar is in the transmitting state for a long time. In order to avoid the damage of the transmitter due to the long working time, the restriction of energy constraint must be considered in the scheduling algorithm. The transient energy of the system attmoment can be expressed as

    (5)

    whereP(x) is the power parameter, andτis the look-back period, which is related to the heat transfer of the system. The energy constraint of the system can be defined as thatE(t) cannot exceed the maximum thresholdEmaxof instantaneous energy at any time, i.e.

    E(t)≤Emax

    (6)

    which meansE(t) cannot exceed the maximum thresholdEmaxof instantaneous energy at any time.

    In simulation, the energy consumption of radar beam and the variation of energy state in time Δtcan be estimated in prior by the parameters of antenna gain, transmission power, pulse width and pulse number, so as to reduce the complexity of algorithm.

    2 Dwell Scheduling Algorithm for DAR

    2.1 Design of task priority

    Without loss of generality, the target at a closer distance, higher speed and moving toward the radar is more threatening and requires a higher priority for tracking and imaging. If the tracking is divided into precision tracking and normal tracking, theith priority of the precision tracking task and the normal tracking task is defined as

    (7)

    (8)

    We divided the search task into high priority search (with priority of 3) and low priority search (with priority of 1). We assume only one verification beam is transmitted for confirmation after searching for a new target. Then the new target is added to the existing tracking task list to recognize the feature, calculate the tracking priority, and scheduled within the next scheduling interval.

    When a precise tracking task enters the stable tracking phase (assume radar enters the stable tracking phase after transmittinglibeam to track), the strategy of tracking and imaging are adopted at the same time in the next scheduling interval. The priority of theith imaging task that enters the stable tracking phase is defined as

    (9)

    Obviously, the priority of the imaging task is between 0 and 1, which is ensured to be lower than the priority of search and tracking tasks. The scheduling of an imaging task usually takes several scheduling intervals, in order to ensure that the pulses transmitted by the radar are not wasted during the imaging process, and a priority dynamic adjustment strategy is adopted for imaging tasks. Ifith imaging task is implemented in thekth scheduling interval, the priority of theith imaging task should be increased appropriately when allocating resources to the next scheduling interval.

    Pi,k+1=Pi,k+ΔP

    (10)

    where ΔPis the step value of priority.

    2.2 Adaptive adjustment of imaging accumulation time

    In order to make better use of the time resources of DAR to image as many precise tracking targets as possible, we propose the self-adaptive adjustment strategy of the imaging accumulation time of different imaging targets. After the end of each scheduling interval, the ISAR image of the target is reconstructed by using all of the previous observed sub-pulses up to the scheduling interval.

    Generally, as the imaging accumulation time increases, the quality of ISAR images increase accordingly and gradually come to a standstill. After reaching the standstill, the quality of ISAR image become worse if the imaging accumulation time continues to increase. Therefore, the mutual information can be taken as the similarity measure of adjacent reconstructed ISAR image. Mutual information represents the amount of information that two images contain about each other. For adjacent reconstructed imagesAandB, their mutual informationI(A,B) is expressed as

    (11)

    wherepi,pjdenote the gray probability distribution ofAandBrespectively, andpijdenotes the joint gray probability distribution. The larger the valueI(A,B) is, the higher the similarity of the two structure image is. After selecting an appropriate thresholdTαby reference to the desired resolution, if the mutual information of the target ISAR image obtained after the two adjacent scheduling intervals is less than the thresholdTα, continue to analyze the imaging task in the next scheduling interval.If the mutual information is more than the thresholdTα, it is considered that the imaging quality of target is up to the expected standard, and the imaging task is finished.

    2.3 Dwell scheduling algorithm

    In order to measure the performance of the radar resource scheduling, the scheduling success ratio (SSR), the hit value ratio (HVR), the time utilization ratio (TUR) and the energy utilization radio (EUR) are taken as the criterion. They can be expressed as

    SSR=N′/N

    (12)

    (13)

    (14)

    (15)

    whereNis the total number of tasks for requested scheduling,N′ is the number of scheduled tasks,Ttis total simulation time,Ptis the peak power of each transmitted pulse andPavis the average power delivered by the radar.

    In terms of the above performance index, assume that there areNdwell tasks in the scheduling interval [t0,te], an effective model of DAR based on pulse interleaving can be defined as

    (16)

    whereN′ is the number of search tasks andq1,q2,q3,q4are adjustment factors. The first constraint gives the scope of the actual execution time of each task; the second constraint indicates that there is no conflict between the transmitted pulses, i.e. the transmitting duration of the radar task cannot be preempted; the third constraint indicates that the search tasks cannot perform pulse interleaving; the fourth constraint indicates that the task dwell received pulse can be overlapped in the same time without colliding with the transmitted pulse; and the fifth constraint represents the energy constraint for task scheduling.

    The nonlinear programming problem in Eq.(16) is an N-P problem, and it is difficult to obtain the optimal solution. As a consequence, heuristic algorithms are usually used to obtain sub-optimal solutions[16]. With the self-adaptive adjustment strategy based on pulse interleaving, a heuristic method for solving the optimization problem is given. The effective algorithm of DAR based on pulse interleaving is summarized as follows.

    D= te-t0Δt ,

    Step1TakeNtasks for requested scheduling in the scheduling interval [t0,te], addKtasks with the latest scheduled start time less thant0to the delete list, discretize the system time; suppose that the length of each time slot is Δt, the number of time slot is introduce time pointertp=t0, initialize the time slot vector U={u1,u2,…uD}=0 and energy state vector E.

    Step2Let the remainingN-Ktasks join the application list according to the corresponding priority from high to low (tasks with the same priority are arranged according to the expected execution time), and leti=1.

    Step3Determine whether theith task can be scheduled at timetp.

    If the task meets the time and energy constraints in Eq.(16), it will be sent to the execution list and removed from the application list. The time slot vector U and time pointertpcan be rewritten as follows.

    Case 1: Search tasks

    uk=1,k∈( stiΔt , sti+txi+twi+triΔt )

    (17)

    tp=sti+txi+twi+tri

    (18)

    Case 2: Tracking tasks

    uk=2,k∈

    (19)

    tp=sti+txi

    (20)

    Case 3: Imaging tasks

    uk=3,k∈( stΔt , sti+txi+twi+triΔt )

    (21)

    tp=sti+txi

    (22)

    In order to distinguish different task types, the time slotsukare assigned 1, 2 and 3 respectively. Update the energy state vector E=E+ΔE(ΔE is the system energy consumption change), leti=i+1, and return to step 3. If the scheduling fails, adjust the actual execution time of the task in the time window, and lettp=tp+Δtp(Δtpis the minimum step size of the pointer slide).

    Step4Iftp

    Step5Ifi≤N-K, return to step 3, otherwise go to step 6.

    Step6The scheduling analysis of the scheduling interval ends. Use the sparse-aperture ISAR imaging method for target imaging by using all of the previously observed sub-pulses up to this scheduling interval. When the mutual information is more than the predefined threshold, this imaging task is no longer scheduled in the next scheduling interval. Otherwise, put the imaging task into the next scheduling interval for analysis.

    The flow chart of the scheduling process in one scheduling interval is as shown in Fig.2.

    Fig.2 Flow chart of the scheduling process

    3 Simulation and Analysis

    In the simulation, we choose three kinds of radar working mode: search, tracking and imaging. Typical parameters for various tasks are summarized in Tab. 1. The total simulation time was 6 s, the scheduling interval length was 50 ms, and the radar can provide an average power of 400 W. For search and tracking tasks, radar transmits narrowband signals, center frequencyfc=10 GHz, signal bandwidthB=10 MHz, and pulse repetition frequencyPf=1 000 Hz. For imaging tasks, radar transmits LFM signals, center frequencyfc=10 GHz, signal bandwidthB=300 MHz, pulse repetition frequencyPf=1 000 Hz, and priority sliding step ΔP=0.1.

    Tab.1 Parameters of radar tasks

    Traditional phased array radar algorithm (traditional algorithm)[15], the pulse interleaving scheduling algorithm(simple task algorithm) which does not consider the imaging tasks[9]and the adaptive dwell scheduling algorithm proposed in this paper (proposed algorithm) are compared in the simulation. Comparison curves of three different algorithms are shown in Tab.2-Tab.5.

    Tab.2 Scheduling success ratio %

    Tab.3 Hit value ratio %

    Tab.4 Time utilization ratio %

    Tab.5 Energy utilization ratio %

    From Tab.2, we see that when the number of tasks is 20, system resources are relatively abundant and the competition for resources among tasks is not obvious. All the three scheduling algorithms can successfully schedule all tasks. With the further increase of the tasks, the scheduling success radio of the traditional algorithm begins to decline significantly, while the other two algorithms based on the pulse interleaving can still successfully schedule all tasks. When the number of tasks increased to 60, the simple task algorithm cannot schedule more tasks, but the proposed algorithm can successfully schedule all tasks.This is due to the saturation of radar resources in the simple task algorithm.In the proposed algorithm, the imaging tasks have the lowest priority, so the remaining resources of the system can be fully utilized by the dynamic adjustment of the imaging accumulation time and the flexible sparse aperture allocation, without affecting the scheduling of the search and tracking tasks. Consequently, the task number of successful scheduling is improved.

    From Tab.3, we see that when the number of tasks is 20, the system resources are relatively abundant; but when the number of tasks increased to 40, the hit value ratio of traditional algorithm has fallen to a great extent. The precision tracking task can be imaged with the idle time of searching and tracking in proposed algorithm, so that it maintains a high hit value radio when the number of tasks reaches 80.

    From Tab.4 and Tab.5, we see that when the number of tasks is 20, the time utilization ratio and the energy utilization ratio plateau at 0.1 because of the resource bottleneck in traditional algorithm. Through the use of the pulse interleaving technology, the simple task algorithm can further utilize the radar system resources, the time utilization ratio can reach about 0.57 and energy utilization can reach about 0.5. On the basis of pulse interleaving, the idle time of search and tracking task is fully exploited in the proposed algorithm, so the time utilization ratio and energy utilization ratio can reach about 0.75 and 0.6 respectively.

    In order to validate that the DAR can achieve the target imaging while searching and tracking, the imaging results of three imaging tasks above are compared with the traditional full aperture ISAR imaging results, and the peak signal to noise ratio (PSNR) is used to measure the effect of the proposed algorithm. The PSNR is defined as

    (23)

    whereσ(i,j) denotes the sparse aperture imaging result in this paper,σ(i,j) denotes the traditional full aperture imaging result,mandndenote the number of rows and columns of the ISAR image matrix respectively. The larger the PSNR value, the better the imaging effect. The contrast between the traditional full aperture imaging and the sparse aperture imaging in this paper is shown in Tab.6. From Tab.6, it is shown that the proposed imaging method can be used to solve the problem of azimuth aperture sparsity and obtain satisfactory image quality.

    Tab.6Comparison of the image quality by the conventional full aperture and the sparse aperture

    4 Conclusion

    In this paper, a dwell scheduling algorithm has been proposed based on pulse interleaving for DAR. The time and energy constraints are set respectively, and the optimal resource scheduling algorithm is established. The simulation results show that the proposed algorithm can accomplish tracking and search task effectively as well as obtain high quality sparse aperture ISAR image. It not only makes use of the multi-task cooperative advantage of DAR, but also improves the resource utilization ratio of radar system. It should be noted that, in this paper, we assume that the target is non-maneuvering and the maneuverability of the target in the actual imaging process is ignored. For the target imaging of maneuvering flight, in addition to the sparsity of the azimuthal aperture, the optimal imaging time and the quality evaluation criteria of the imaging should be considered. Related content should be further studied.

    [1] Wang Feng, Fu Youguang. Performance comparison of digital phased array and conventional phased array radar[J]. Journal of CAEIT, 2012, 7(2): 148-151. (in Chinese)

    [2] Cheng Ting, He Zishu, Li Huiyong. An adaptive dwell scheduling algorithm for digital array radar[J]. Acta Electronica Sinica, 2009, 37(9): 2025-2029. (in Chinese)

    [3] Komorniczak W. Pietrasinski[J]. Selected problems of MFR resources management[C]∥The 3rd International Conference on Information Fusion, Paris, France, July 10-13, 2000.

    [4] Jang D S, Choi H L, Roh J E. A time-window-based task scheduling approach for multi-function phased array radars[C]∥IEEE International Conference on Control, Automation and Systems, Gyeonggi-do, South Korea, Oct. 26-29, 2011.

    [5] Lee C G, Kang P S, Shi C S, et al. Radar dwell scheduling considering physical characteristics of phased array antenna[C]∥Proceedings of the 24th IEEE International Real-Time Systems Symposium(RTSS’03), Cancun, Mexico, Dec 5, 2003.

    [6] Zhang Juan, Xia Zhongting. An adaptive scheduling searching method based on information entropy[J]. Modern Radar, 2015, 37(8): 33-36. (in Chinese)

    [7] Orman A J, Potts C N, Shahani A K. Scheduling for a multifunction phased array radar system[J]. European Journal of Operational Research,1996,90(1):13-25.

    [8] Zhao Yu, Li Jianxun, Cao Yinglan, et al, Adaptive scheduling algorithm based on quadratic programming for multifunction phased array radars[J]. Systems Engineering and Electronics, 2012,34: 698-703. (in Chinese)

    [9] Cheng Ting, He Zishu. Novel radar dwell scheduling algorithm based on pulse interleaving[J]. Journal of Systems Engineering and Electronics, 2009, 20(2): 247-253. (in Chinese)

    [10] Chen Yijun, Zhang Qun, Luo Ying, et al. Adaptive scheduling algorithm for radar based on sparse aperture ISAR imaging[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2013, 33(4): 171-176. (in Chinese)

    [11] Li Wenjing, Chen Hongwei. A kind of ISAR imaging algorithm based on compressed sensing[J]. Computer Simulation, 2015, 32(8): 10-13. (in Chinese)

    [12] Chen Y J, Zhang Q, Luo Y, et al. Measurement matrix optimization for ISAR sparse imaging based on genetic algorithm[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(12):1875-1879.

    [13] Chen J, Wang L, Zhang W, et al. Multifunction phased radar resource management via maximal pulse interleaving technique[J]. Arabian Journal for Science and Engineering, 2013, 38(11): 3081-3091.

    [14] Hou Q K, Liu Y, Fan L J, et al. Compressed sensing digital receiver and orthogonal reconstructing algorithm for wideband ISAR radar[J]. Science China Information Sciences, 2015, 58(2):1-10.

    [15] Wang Qiang, Xu Jungang, Wang Hongan, et al. A new priority table based real-time scheduling algorithm[J]. Acta Electronica Sinica, 2004, 32(2): 310-313. (in Chinese)

    [16] Focke R W, Villiers J P D, Inggs M R. Interval algebra-An effective means of scheduling surveillance radar networks[J]. Information Fusion, 2015, 23(3):81-98.

    国产极品天堂在线| 好男人视频免费观看在线| 深爱激情五月婷婷| 国产成年人精品一区二区| 亚洲av男天堂| 少妇被粗大猛烈的视频| eeuss影院久久| 免费观看a级毛片全部| 国产亚洲5aaaaa淫片| 欧美+亚洲+日韩+国产| 亚洲第一区二区三区不卡| 成人高潮视频无遮挡免费网站| 久久久a久久爽久久v久久| 国产精品久久久久久久久免| 麻豆av噜噜一区二区三区| 一进一出抽搐动态| 国产黄片视频在线免费观看| 久久久精品欧美日韩精品| АⅤ资源中文在线天堂| 99久久中文字幕三级久久日本| 久久精品久久久久久噜噜老黄 | 亚洲在久久综合| 亚洲av中文字字幕乱码综合| 一级毛片aaaaaa免费看小| 亚洲第一电影网av| 成年av动漫网址| 久久草成人影院| 美女 人体艺术 gogo| 最近手机中文字幕大全| 亚洲五月天丁香| 嫩草影院入口| 欧美xxxx黑人xx丫x性爽| 99国产精品一区二区蜜桃av| 日本av手机在线免费观看| 搡老妇女老女人老熟妇| 久久综合国产亚洲精品| 有码 亚洲区| 国产一区二区在线av高清观看| www.色视频.com| 色噜噜av男人的天堂激情| www日本黄色视频网| 亚洲成人中文字幕在线播放| 自拍偷自拍亚洲精品老妇| 欧美+亚洲+日韩+国产| 此物有八面人人有两片| 国产极品精品免费视频能看的| 国产高清不卡午夜福利| 精品熟女少妇av免费看| 夫妻性生交免费视频一级片| 九草在线视频观看| www.av在线官网国产| 久久鲁丝午夜福利片| 国产探花极品一区二区| 欧美高清成人免费视频www| 午夜老司机福利剧场| а√天堂www在线а√下载| 日韩欧美精品免费久久| 中文字幕精品亚洲无线码一区| 欧美极品一区二区三区四区| 男女那种视频在线观看| 久久久欧美国产精品| av天堂中文字幕网| 男女啪啪激烈高潮av片| avwww免费| 成年女人永久免费观看视频| 一级毛片久久久久久久久女| 国产亚洲精品久久久com| 成人二区视频| 欧美一级a爱片免费观看看| 欧美高清成人免费视频www| 九九在线视频观看精品| 国产精品一区二区三区四区久久| 非洲黑人性xxxx精品又粗又长| 免费一级毛片在线播放高清视频| 在线观看av片永久免费下载| 亚洲av成人av| 国产精品美女特级片免费视频播放器| 中国美白少妇内射xxxbb| 黄片无遮挡物在线观看| 卡戴珊不雅视频在线播放| 一卡2卡三卡四卡精品乱码亚洲| 成人高潮视频无遮挡免费网站| 免费观看a级毛片全部| 69人妻影院| 亚洲av不卡在线观看| 国内精品一区二区在线观看| 免费不卡的大黄色大毛片视频在线观看 | 一边亲一边摸免费视频| 国产黄色小视频在线观看| 国产精品蜜桃在线观看 | 九色成人免费人妻av| 禁无遮挡网站| 内地一区二区视频在线| 午夜福利在线观看免费完整高清在 | 婷婷色av中文字幕| 一个人看视频在线观看www免费| 美女大奶头视频| 丝袜美腿在线中文| 日本在线视频免费播放| 欧美精品一区二区大全| 免费大片18禁| 在线播放国产精品三级| 精品熟女少妇av免费看| 成人漫画全彩无遮挡| 中出人妻视频一区二区| 九九在线视频观看精品| 狂野欧美激情性xxxx在线观看| 国产亚洲欧美98| 天天一区二区日本电影三级| 国产精品一区二区三区四区久久| 麻豆av噜噜一区二区三区| 国内精品一区二区在线观看| 国产男人的电影天堂91| 美女 人体艺术 gogo| 成人永久免费在线观看视频| 男女下面进入的视频免费午夜| 亚洲高清免费不卡视频| 青春草视频在线免费观看| 日韩一本色道免费dvd| 深爱激情五月婷婷| 亚洲久久久久久中文字幕| 国产亚洲精品久久久久久毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲人成网站在线播放欧美日韩| 国产视频首页在线观看| 国产精品永久免费网站| 国内精品美女久久久久久| 亚洲国产欧美在线一区| 国产一级毛片在线| 免费观看在线日韩| 欧美日韩乱码在线| 别揉我奶头 嗯啊视频| 国产极品精品免费视频能看的| 久久精品综合一区二区三区| 国产精品蜜桃在线观看 | 亚洲av熟女| 天堂影院成人在线观看| 国产精品一区二区性色av| 嫩草影院精品99| 欧美成人免费av一区二区三区| 九九在线视频观看精品| 超碰av人人做人人爽久久| 精品久久久久久成人av| 天天躁夜夜躁狠狠久久av| 欧美日韩乱码在线| 亚洲欧美日韩高清在线视频| 久久国产乱子免费精品| 1024手机看黄色片| 黄色配什么色好看| 国产一区二区在线av高清观看| 免费看av在线观看网站| 国产极品精品免费视频能看的| or卡值多少钱| 丝袜喷水一区| 亚洲精品乱码久久久久久按摩| 午夜亚洲福利在线播放| 丝袜喷水一区| 久久99蜜桃精品久久| 毛片一级片免费看久久久久| 啦啦啦韩国在线观看视频| av在线观看视频网站免费| 69av精品久久久久久| 国产高清三级在线| 日韩欧美精品v在线| 好男人视频免费观看在线| 精品久久久噜噜| 国产av在哪里看| 欧美成人a在线观看| 精品久久久久久久久久免费视频| 日韩欧美精品v在线| 日韩高清综合在线| 欧美日韩国产亚洲二区| 长腿黑丝高跟| 精品久久久久久成人av| 一区二区三区高清视频在线| 超碰av人人做人人爽久久| 级片在线观看| 国产精品久久久久久亚洲av鲁大| 春色校园在线视频观看| 亚洲最大成人中文| 一级毛片电影观看 | 国产探花在线观看一区二区| 搡女人真爽免费视频火全软件| 久久久久网色| 变态另类成人亚洲欧美熟女| 色综合亚洲欧美另类图片| 99久久无色码亚洲精品果冻| 联通29元200g的流量卡| 色综合亚洲欧美另类图片| 亚洲人与动物交配视频| 在线观看免费视频日本深夜| 日本av手机在线免费观看| 亚洲人成网站高清观看| 欧美日韩精品成人综合77777| 99久久精品国产国产毛片| 国产伦理片在线播放av一区 | 久久草成人影院| 欧美三级亚洲精品| 欧美激情在线99| 中文字幕久久专区| 日本与韩国留学比较| 欧美成人一区二区免费高清观看| 亚洲精华国产精华液的使用体验 | 久久这里有精品视频免费| 亚洲自拍偷在线| 18禁裸乳无遮挡免费网站照片| 97超碰精品成人国产| 热99在线观看视频| 日本黄色片子视频| 日本免费一区二区三区高清不卡| 好男人视频免费观看在线| 少妇的逼水好多| 99久久精品热视频| 又爽又黄无遮挡网站| 欧美三级亚洲精品| 国产一区二区亚洲精品在线观看| 又爽又黄a免费视频| 亚洲av中文av极速乱| 六月丁香七月| kizo精华| 99在线人妻在线中文字幕| 如何舔出高潮| 人人妻人人澡欧美一区二区| 简卡轻食公司| 亚洲第一区二区三区不卡| 校园人妻丝袜中文字幕| 欧美高清成人免费视频www| 美女国产视频在线观看| 麻豆av噜噜一区二区三区| 日本黄色视频三级网站网址| 嫩草影院新地址| 欧美潮喷喷水| 免费不卡的大黄色大毛片视频在线观看 | 日韩强制内射视频| 亚洲三级黄色毛片| 久久久久久久久大av| 久久99热6这里只有精品| 成人特级黄色片久久久久久久| 午夜亚洲福利在线播放| 在线播放无遮挡| 97超视频在线观看视频| 别揉我奶头 嗯啊视频| 国产av在哪里看| 国内精品宾馆在线| 三级男女做爰猛烈吃奶摸视频| 最近视频中文字幕2019在线8| 亚洲国产精品国产精品| 日本与韩国留学比较| 精品欧美国产一区二区三| 大香蕉久久网| 丰满乱子伦码专区| 国产色婷婷99| 日韩欧美一区二区三区在线观看| 嘟嘟电影网在线观看| 精品人妻偷拍中文字幕| 国模一区二区三区四区视频| www.av在线官网国产| a级毛片a级免费在线| 午夜久久久久精精品| 97人妻精品一区二区三区麻豆| 麻豆成人午夜福利视频| av.在线天堂| 一区二区三区四区激情视频 | 毛片一级片免费看久久久久| 亚洲欧美日韩东京热| 少妇的逼好多水| 亚洲国产精品国产精品| 一区福利在线观看| 免费看日本二区| 在线观看av片永久免费下载| 麻豆成人午夜福利视频| 不卡一级毛片| 97热精品久久久久久| 国产探花极品一区二区| 国产精品女同一区二区软件| 精品久久久久久久人妻蜜臀av| 久久精品国产清高在天天线| 人妻久久中文字幕网| 神马国产精品三级电影在线观看| 国产精品美女特级片免费视频播放器| 成人毛片60女人毛片免费| 国产精品久久视频播放| 成年版毛片免费区| 丝袜美腿在线中文| www.av在线官网国产| 青春草亚洲视频在线观看| 精品久久久久久久久久免费视频| 日韩欧美三级三区| 中文字幕人妻熟人妻熟丝袜美| av国产免费在线观看| 久久久久性生活片| 日韩一本色道免费dvd| 给我免费播放毛片高清在线观看| 国产精品久久久久久久电影| 91久久精品国产一区二区成人| 亚洲成av人片在线播放无| 免费不卡的大黄色大毛片视频在线观看 | 欧美在线一区亚洲| 日韩,欧美,国产一区二区三区 | 欧美3d第一页| 最好的美女福利视频网| 日本撒尿小便嘘嘘汇集6| 两个人的视频大全免费| 精品欧美国产一区二区三| 在线a可以看的网站| 天美传媒精品一区二区| 精品少妇黑人巨大在线播放 | a级毛色黄片| 欧美一区二区精品小视频在线| 亚洲精品国产av成人精品| 亚洲av二区三区四区| 欧美日韩乱码在线| 亚洲在久久综合| 狂野欧美激情性xxxx在线观看| 成人综合一区亚洲| 国产男人的电影天堂91| 中文字幕久久专区| 深夜精品福利| 啦啦啦观看免费观看视频高清| 亚洲综合色惰| 国产综合懂色| 爱豆传媒免费全集在线观看| 乱码一卡2卡4卡精品| 亚洲欧洲国产日韩| 九九爱精品视频在线观看| 日本与韩国留学比较| 中文字幕精品亚洲无线码一区| 一进一出抽搐gif免费好疼| 亚洲四区av| 熟女电影av网| 毛片女人毛片| 1000部很黄的大片| 亚洲四区av| 国产精品久久久久久久电影| 春色校园在线视频观看| 国产极品精品免费视频能看的| 99riav亚洲国产免费| 成年版毛片免费区| 18禁在线无遮挡免费观看视频| 岛国毛片在线播放| 一个人观看的视频www高清免费观看| 国产探花极品一区二区| 亚洲av熟女| 亚洲欧美日韩东京热| 日本免费a在线| 日日摸夜夜添夜夜添av毛片| 一级毛片久久久久久久久女| 一个人看视频在线观看www免费| 欧美性感艳星| 亚洲在久久综合| 免费观看人在逋| 国产伦在线观看视频一区| 国产精品精品国产色婷婷| 久久久欧美国产精品| 成人特级av手机在线观看| 色播亚洲综合网| 日韩国内少妇激情av| 美女xxoo啪啪120秒动态图| 欧美日本视频| 午夜亚洲福利在线播放| 国内精品一区二区在线观看| 天天躁日日操中文字幕| 一个人看的www免费观看视频| 久久6这里有精品| 成人高潮视频无遮挡免费网站| 嫩草影院新地址| 精品少妇黑人巨大在线播放 | 精品日产1卡2卡| 亚洲精品日韩av片在线观看| 精品久久久久久成人av| 99热网站在线观看| 亚洲人成网站高清观看| 国产精品人妻久久久影院| 99热这里只有是精品50| 久久精品影院6| 女的被弄到高潮叫床怎么办| 在线观看av片永久免费下载| 人妻系列 视频| 夜夜看夜夜爽夜夜摸| 免费看美女性在线毛片视频| 国产午夜精品论理片| 亚洲av一区综合| 伊人久久精品亚洲午夜| 夜夜爽天天搞| 插逼视频在线观看| 国内精品宾馆在线| 日本免费一区二区三区高清不卡| 日韩一区二区三区影片| 精华霜和精华液先用哪个| 国产精品蜜桃在线观看 | 国产69精品久久久久777片| 欧美日韩在线观看h| 成人二区视频| 婷婷六月久久综合丁香| 男女那种视频在线观看| 亚洲欧洲日产国产| 91狼人影院| 一级毛片久久久久久久久女| 99久久九九国产精品国产免费| 天天躁日日操中文字幕| 成年女人看的毛片在线观看| 国产精品美女特级片免费视频播放器| 人人妻人人澡欧美一区二区| 成人亚洲精品av一区二区| 国内揄拍国产精品人妻在线| 亚洲美女搞黄在线观看| 亚洲av.av天堂| 国产在线精品亚洲第一网站| 亚洲内射少妇av| 午夜视频国产福利| 我要搜黄色片| 人人妻人人澡人人爽人人夜夜 | 国产精品人妻久久久影院| 国产探花在线观看一区二区| 我要看日韩黄色一级片| 国产成人一区二区在线| 女人十人毛片免费观看3o分钟| 久久久久久久久久久免费av| 欧美xxxx黑人xx丫x性爽| 亚洲成人av在线免费| 插阴视频在线观看视频| 高清午夜精品一区二区三区 | 哪个播放器可以免费观看大片| 欧美成人免费av一区二区三区| 草草在线视频免费看| 精品久久久久久久久久久久久| 亚洲18禁久久av| 日韩一区二区视频免费看| 国产视频首页在线观看| 欧美激情久久久久久爽电影| 一本久久中文字幕| 国产乱人视频| 极品教师在线视频| 国产精品一及| 亚洲av中文字字幕乱码综合| 搡老妇女老女人老熟妇| 欧美日韩综合久久久久久| 亚洲va在线va天堂va国产| 我的老师免费观看完整版| 狂野欧美白嫩少妇大欣赏| 日本与韩国留学比较| 午夜精品国产一区二区电影 | 晚上一个人看的免费电影| 日韩制服骚丝袜av| 久久久久久久亚洲中文字幕| 日韩强制内射视频| 少妇高潮的动态图| 午夜精品在线福利| 久久99精品国语久久久| 九九爱精品视频在线观看| 久久久久久久久久黄片| 午夜视频国产福利| 校园人妻丝袜中文字幕| 欧美日韩综合久久久久久| 麻豆成人午夜福利视频| 中文欧美无线码| 成人毛片60女人毛片免费| 三级毛片av免费| 91av网一区二区| 国内久久婷婷六月综合欲色啪| 深夜精品福利| 日日摸夜夜添夜夜添av毛片| 性插视频无遮挡在线免费观看| 中文亚洲av片在线观看爽| 午夜爱爱视频在线播放| 18禁黄网站禁片免费观看直播| 久久久久久久久久久丰满| 久久精品影院6| av在线亚洲专区| 嫩草影院入口| 国产高潮美女av| 性欧美人与动物交配| 亚洲欧美成人综合另类久久久 | 国产精品综合久久久久久久免费| 一级av片app| 日韩制服骚丝袜av| 国产精品野战在线观看| 69人妻影院| 全区人妻精品视频| 日本撒尿小便嘘嘘汇集6| 成人国产麻豆网| 午夜精品一区二区三区免费看| 青春草亚洲视频在线观看| 三级经典国产精品| 人妻少妇偷人精品九色| 亚洲精品成人久久久久久| 国内久久婷婷六月综合欲色啪| АⅤ资源中文在线天堂| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧美清纯卡通| 亚洲真实伦在线观看| 亚洲不卡免费看| 久久精品影院6| 熟妇人妻久久中文字幕3abv| 亚洲精品国产av成人精品| 波多野结衣高清无吗| 精品久久久久久成人av| 有码 亚洲区| 久久鲁丝午夜福利片| 久久九九热精品免费| 2021天堂中文幕一二区在线观| 色尼玛亚洲综合影院| 色哟哟·www| 亚洲国产精品合色在线| 日韩高清综合在线| 国产av在哪里看| 在线播放国产精品三级| a级毛片a级免费在线| 久久久久久久久久久免费av| 乱系列少妇在线播放| 五月伊人婷婷丁香| 啦啦啦观看免费观看视频高清| 一个人免费在线观看电影| 欧美日韩一区二区视频在线观看视频在线 | 小蜜桃在线观看免费完整版高清| 狠狠狠狠99中文字幕| av.在线天堂| 久久精品久久久久久久性| 我的女老师完整版在线观看| 中国美白少妇内射xxxbb| 日韩欧美精品免费久久| 国产亚洲5aaaaa淫片| 天堂影院成人在线观看| 听说在线观看完整版免费高清| 久久精品国产亚洲av涩爱 | 国产精品福利在线免费观看| 国产亚洲av片在线观看秒播厂 | 亚洲欧美日韩无卡精品| 国产大屁股一区二区在线视频| 我要看日韩黄色一级片| 高清午夜精品一区二区三区 | 日日撸夜夜添| 日日摸夜夜添夜夜添av毛片| 插阴视频在线观看视频| 亚洲成人精品中文字幕电影| 成人美女网站在线观看视频| 日本av手机在线免费观看| 不卡一级毛片| 人妻久久中文字幕网| 国产成人福利小说| 国语自产精品视频在线第100页| 在线免费十八禁| 国产不卡一卡二| 日韩一区二区视频免费看| 国产成人a∨麻豆精品| 日韩人妻高清精品专区| 亚洲第一电影网av| 欧美性猛交╳xxx乱大交人| 国产精品久久久久久久电影| 久久精品国产99精品国产亚洲性色| 亚洲自偷自拍三级| 人妻制服诱惑在线中文字幕| 插阴视频在线观看视频| 日韩欧美精品v在线| 精品久久久噜噜| 亚洲av熟女| 中文字幕av成人在线电影| 99久久人妻综合| 嫩草影院精品99| 中文亚洲av片在线观看爽| 性欧美人与动物交配| 国产成人精品婷婷| 国产精品.久久久| 国产黄色视频一区二区在线观看 | 午夜老司机福利剧场| 国产精品福利在线免费观看| 国产日韩欧美在线精品| 亚洲欧美日韩卡通动漫| 又粗又爽又猛毛片免费看| 国产在视频线在精品| 少妇裸体淫交视频免费看高清| 身体一侧抽搐| 中文字幕制服av| 久久久久久国产a免费观看| av天堂中文字幕网| 国产精品一二三区在线看| 亚洲av免费在线观看| 久久精品国产亚洲av涩爱 | av国产免费在线观看| 久久这里有精品视频免费| videossex国产| 亚洲av中文av极速乱| 久久久a久久爽久久v久久| 亚洲七黄色美女视频| 国内少妇人妻偷人精品xxx网站| 亚洲第一电影网av| 在线观看美女被高潮喷水网站| 1000部很黄的大片| 别揉我奶头 嗯啊视频| 美女xxoo啪啪120秒动态图| 精品一区二区三区人妻视频| 欧美极品一区二区三区四区| or卡值多少钱| 国产极品精品免费视频能看的| 校园春色视频在线观看| 国产一区二区三区在线臀色熟女| 国产极品精品免费视频能看的| 亚洲人成网站在线播| 嫩草影院精品99| 身体一侧抽搐| 别揉我奶头 嗯啊视频| 免费看美女性在线毛片视频| 欧美精品一区二区大全| 尤物成人国产欧美一区二区三区| 日韩成人av中文字幕在线观看| 午夜福利高清视频| 国产综合懂色| 日日摸夜夜添夜夜爱| 亚洲真实伦在线观看| 大又大粗又爽又黄少妇毛片口| 嫩草影院入口| 深夜a级毛片| 亚洲人成网站在线观看播放| 午夜福利视频1000在线观看| 日韩精品有码人妻一区|