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

    Targetdetection for low angle radar based on multi-frequency order-statistics

    2015-12-23 10:08:53

    National Laboratory of Radar Signal Processing,Xidian University,Xi’an 710071,China

    Targetdetection for low angle radar based on multi-frequency order-statistics

    Yunhe Cao?,Shenghua Wang,Yu Wang,and Shenghua Zhou

    National Laboratory of Radar Signal Processing,Xidian University,Xi’an 710071,China

    Forradartargets fl ying atlow altitude,multiple pathways produce fade or enhancement relative to the level that would be expected in a free-space environment.In this paper,a new detection method based on a wide-ranging multi-frequency radar for low angle targets is proposed.Sequentialtransmitting multiple pulses with differentfrequencies are fi rstapplied to decorrelate the coherence ofthe directand re fl ected echoes.After receiving allechoes, the multi-frequency samples are arranged in a sortdescending according to the amplitude.Some high amplitude echoes in the same range cell are accumulated to improve the signal-to-noise ratio and the optimalnumber ofhigh amplitude echoes is analyzed and given by experiments.Finally,simulation results are presented to verify the effectiveness ofthe method.

    multipath,signal detection,order statistic,multifrequency,low angle.

    1.Introduction

    The problem of detecting and tracking low angle targets in the presence of multipaths has attracted lots of interest in the radarcommunity formany years[1-10,14-20].In lowangle targetdetection[3,5-7,9],radars receive two ormore coherentechoes via multipaths(typically four propagation paths),resulting in fade or enhancementof the targetecho amplitude.Severe signal cancellation arises and the target cannotbe detected effectively when the directand re fl ected components are nearly outof phase.

    To alleviate multipath effects,the multi-frequency technique is used in some literature[3,7].The technique is effective to decorrelate the coherence of the direct and re fl ected echoes.The carrier frequency of radar changes within a large bandwidth from pulse to pulse.This kind of multi-frequency radar provides some characters of robustness to the signal fades and makes a proper use of multipath propagation to improve the target detection capability[8].Considering fourpropagation paths(one direct path and three re fl ected paths),ifthe directpath and the refl ected paths are in phase atsome frequency point,the echo amplitude in a multipath environmentshould be four times greater than that in the free space(only direct path)when the amplitude of the ground re fl ection coef fi cient is 1.0. Consequently,the signal-to-noise ratio(SNR)of the received signal is expected to get a gain of 12 dB.If the carrier frequency changes uniformly,the differentialphase between the directpath and the specularre fl ection path can be regarded as a uniformly distributed random variable and the gain of the received amplitude can range from 0 to 4. Consequently,the SNR of the return is bene fi cialin statistics for the detection of targets in multipaths[3].

    The noncoherent accumulation method is common in radarsignalprocessing.However,during the low angle targetdetection,because of the severe multipath cancellation atsome frequencies,extremely low SNR echoes should be excluded from the accumulation.A new detection method, which is based on the order statistics and noncoherentaccumulation,is described in this paper.In the method,only some test cells with strong energy are employed in the detection procedure and the optimal number of test cells is given accurately.With regard to the detection performance, ourmethod achievesa higherdetection probability than the M-out-of-N detectorproposed by Wilson and Carlson[3].

    This paper is organized as follows.In Section 2,a traditionalmultipath modeland the probability density function(PDF)ofthe targetecho forfourpropagation paths are given.In Section 3,the proposed multi-frequency orderstatistics method is described in detail.The expression of optimal parameter K is also introduced in this Section. Performance comparison of this method with the M-outof-N detector and the coherent accumulation method are presented in Section 4.Finally concluding remarks are given.

    2.Multipath model

    A simple fl at-earth multipath modelwith perfectre fl ection is used as shown in Fig.1.The radarreceives the echo signals from a low elevation targetand four components are taken into account.The fi rstcomponentgoes straightto the target and returns straight to the antenna,whereas the second componentgoes straight to the target and returns to the antenna via the re fl ecting plane.The third component travels to the target via the re fl ecting plane and returns straight to the antenna,whereas the last component travels to the target via the re fl ecting plane and so does the return[10].

    Fig.1 Plat-earth modelin multipath

    Suppose thatthe radar cross section(RCS)of the target is the same for the four propagation paths and the direct echo amplitude is normalized.The diffuse re fl ection is neglected for simplicity.Then,the voltage multipath factor of the received signalis given by

    whereα=2πfδ/c,which is converted from the distance differenceδbetween the directand the specularre fl ections, denotes the phase difference,ρis the complex re fl ection coef fi cient,f and c are the carrier frequency and the speed of light,respectively.

    Assuming thatthe frequency f is hopped uniformly over the radar work bandwidthΔf.It is expected to enable the echo amplitude G to traverse a cycle within the bandwidth Δf(and soαundergoes a phase shiftofπ),i.e.

    whereδcan be estimated as follows:

    where Rois the distance between the radar and the target, R1is the distance between the radar and the targetimage, haand htare the heightofthe radarantenna and the height of the target,respectively,and D is the ground range ofthe target.

    For a low angle target,ha,ht<<D.Using Taylor expansion,(3)can be written approximately as

    In this case,D≈R0.Equation(4)is often written[13]as

    For example,assuming that ha=20 m,ht=50 m, Ro=10 km,the work frequency bandwidthΔf is expected to be 0.75 GHz.In this case,phase differenceαcan be treated as a uniformly distributed random variable over [0,π].Therefore,W=α/2πis uniformly distributed over [0,1/2],and the PDF of W can be expressed as

    The phase of re fl ection coef fi cientρis approximately constant in a certain frequency range for horizontal polarization[11].Hence,we consider thatρis a complex constantnumberfordifferentcarrier frequencies(in the whole radar work bandwidth).Then we have

    After substituting(6)and(8)into the following PDF formula,we obtain the PDF of G

    3.Multi-frequency order-statistics method

    In this section,the expression of multi-frequency orderstatistics false alarm probability function is deduced and a detection threshold is given.Subsequently,the detection probability of the multi-frequency order-statistics method in the multipath environment is presented.Furethermore, the optimal number of high amplitude echoes is analyzedby experiments.Severalassumptions are made in the algorithm as follows:i)the noise is temporally white and circularly symmetric zero-mean complex Gaussian random variable;ii)the noise in different cells is mutually statistically independent;iii)all the multipath delayed signals from the target lie within the range cell under consideration;iv)the sample data in varied frequency pulses are mutually statistically independent since the frequency bandwidth is wide enough.

    3.1 Statisticalanalysis of test cells

    Assuming thata frequency agile radar has received N differentfrequency echoes Xt(t=1,2,...,N),the testcells are chosen from the same range cellof the N echoes.The sample data obtained through the linear-law detector in the testcells can be denoted as Yt=|Xt|.

    In the hypothesis H0,there is only noise in the testcell. Thus,the PDF of Yt,which is the Rayleigh PDF,can be written as

    In the multipath environment,the effects of multipath fl uctuation and target fl uctuation need to be taken into deduction of the PDF of Yt.Here,we consider the Swerling III target model for the following analysis,and the amplitude of targetecho A'is

    where A represents the amplitude of the target fl uctuation. The PDF of A for Swerling IIItargetis given by

    where f(a,g)is the joint PDF of a and g.The PDF of A'can be obtained by the derivative of FA'(a')with respect to a':

    Since the target fl uctuation A and the target multipath echoes amplitude G are independent,the PDF of A'can be rewritten as

    Substituting(9)and(12)into(15),and noticing that g∈is a positive number,then

    In the hypothesis H1,the samples in test cells contain the target signal.The PDF of the output signal of linearlaw detector can be expressed as

    where I0(·)is the modi fi ed Bessel function of the zero order.

    By(16)and(17),we can get the PDF of Ytin the hypothesis H1

    3.2 Detection probability of the proposed method

    Due to the intense variation of the return power of the different frequency echoes in the multipath environment,it is not optimal to accumulate all the received frequency echoes.It is feasible to choose only K strong frequency echoes to accumulate in the detection procedure.

    Let Y(t)(t=1,2,...,N)be a descending order,i.e.

    In this way,only K test cells with strong energy are involved in the non-coherent integration.De fi ne the teststatistic as

    According to the detection theory,the targetis detected when the test statistic Z exceeds the threshold value ZT, where the decision hypotheses are

    Since Y1,...,YNare independent identically distributed(IDD),we obtain from[12]that the joint PDF of the fi rst K large values after sorting is

    where y(t)is the order statistics,satisfying y(1)≥y(2)≥...≥y(K),fY(y(t))is the individual PDF of Y, and FY(y(K))denotes the CDF of Y(FY(y(K))=By using the iterated integraltransformation,we can getthe PDF of Z

    The function fZ(z)is complicated in calculation,so the numerical integration is generally needed.The function fY(·)in the hypothesis H1and H0have been given in (10)and(18),respectively.

    Given that the probability of false alarm Pfais known previously,the voltage threshold ZTcan be obtained by the following equation

    where fZ(z|H0)represents the PDF of Z in the hypothesis H0,which can be obtained by substituting(10)into(23).

    By using the PDF of Z in the hypothesis H1and the voltage threshold ZT,the detection probability Pdis

    where fZ(z|H1)is the PDF of Z in the hypothesis H1, which can be obtained by substituting(18)into(23).

    3.3 The optimum value ofK

    In what follows,the important problem is that how to judiciously choose the value of K and substantially enhance the multipath detection performance.For speci fi ed values of N,Pfaand Pd,the optimal value of K is expected to be picked outaccording to the minimum SNR required.

    Here,a point-by-pointgraphicalprocedure has been carried out by using the Monte-Carlo simulations based on 108independenttrials fora Swerling IIItarget.The simulation results of minimum SNR required are shown in Fig.2 to Fig.5,for N=5,11,21,and 31,respectively,and for the various indicated combinations of Pfaand Pdat two differentvalues.

    Fig.2 Minimum SNR required per pulse as a function of K for N=5

    Fig.3 Minimum SNR required per pulse as a function of K for N=11

    Fig.4 Minimum SNR required per pulse as a function of K for N=21

    Fig.5 Minimum SNR required per pulse as a function of K for N=31

    It can be seen from Fig.2 to Fig.5 that the value of K willaffectthe performance of the proposed method and there is an optimalvalue of K fora fi xed N attwo differentvalues.They also show that the traditional noncoherent accumulation method using allfrequency echoes(K=N) is notthe bestchoice in the multipath environment.

    For a particular N with re fl ection coef fi cient amplitudespeci fi ed,the optimal value of K can be reasonably considered to be independent of Pfaand Pdwhen 10?6≤Pfa≤10?4and 0.6≤Pd≤0.9.We can see that the optimal value of K seems to shift slowly with increasing N and be approximately proportional to√N(yùn).Thus, the region of optimal K may be expressed as

    where Z+denotes the positive integer.Fig.6 plots K versus the number of N.It can be seen from Fig.6 that the optimal value of K lies in the region(26).Additionally, the value ofˉρhas slightin fl uence on the choice of K.Generally speaking,the larger the value ofˉρ,the smaller the optimal K.That is to say,the value ofˉρis close to 0.6, we can selecta big integer satisfying(26),whereas a small integer is chosen when the value ofis close to 0.9.

    Fig.6 Parameter K versus N

    4.Performance comparison

    In this section,we perform Monte Carlo simulations based on 106independenttrials to realize the performance comparisons with the following speci fi cations:the number of carrier frequencies N=11,the number of reference cells L=64,false-alarm probability Pfa=10?4.The target’s RCS fl uctuation is assumed to obey a Swerling III distribution.

    Fig.7 depicts the variations of detection probability as a function of SNR per echo for the proposed detector and the M-out-of-N detector attwo differentˉρvalues.

    Fig.7 Detection probability oftwo methods atdifferentvalues

    According to the reference[3],the optimum of M for the M-out-of-N detector is 2,while Fig.3 shows thatthe optimum of K for the proposed method is differentat two differentˉρvalues.Itis obvious thatthe performance of theproposed detectoris betterthan thatofthe M-out-of-N detectorregardlessofthe value ofThe required SNRofthe proposed method is decreased by up to 2 dB for Pd=0.9 compared with that of the M-out-of-N detector for two cases plotted0.9.The value ofhas little in fl uence on the gain of SNR,thatis,the performance of the proposed algorithm is stable and reliable in multipath detection.

    Now we compare the proposed algorithm with the multiple pulse coherent accumulation method in a low angle circumstance.For the sake of fairness,the transmit power is the same for the two methods.The proposed method sequentially transmits 11 pulses with different frequencies and then four(i.e.K=4)large amplitude echoes are noncoherently integrated as test statistic.The coherent accumulation method transmits 11 pulses with the same frequency and then all the pulses are coherently integrated as test statistics.Fig.8 depicts the detection performance of the proposed method and the coherent accumulation method.The fi gure includes two group curves where differentre fl ection coef fi cientamplitudevalues are used.It can be seen from the fi gure thatthe proposed method has a good performance when asked for a high detection probability in the multipath environment.For example,the value ofrequested SNRmay reduce from 20 dBto 0 dBatdetection probability 0.9 when re fl ection coef fi cient amplitudeis equavalentto 0.9.When the re fl ection coef fi cien tis 0.6,the value of requested SNR may reduce from 13 dB to 2 dB atdetection probability 0.9.

    Fig.8 Effects of different SNR values on detection probability

    5.Conclusions

    In order to improve targetdetection performance at a low angle multipath environment,this paper proposes an order statistics detection method based on the multi-frequency technique.Only K large amplitude frequency samples rather than all are involved into the noncoherent integration,achieving a high detection probability at the same SNR.Moreover,a simple statistic model is given to obtain the optimalvalue of K in the multipath environment. Moreover,the proposed algorithm is robust to the multipath re fl ection coef fi cient.The simulation results demonstrate thatthe proposed method is feasible and effective.

    [1]J.Kim,H.J.Yang,N.Kwak.Low-angle tracking of two objects in a three-dimensional beamspace domain.IET Radar Sonar Navigation,2012,6(1):9-20.

    [2]Y.G.Hou,Y.Y.Shen,Z.Z.Zhang.Elevation estimation for low-angle targetbased on re fl ection paths suppression.Journal ofSystems Engineering and Electronics,2008,19(1):71-75.

    [3]S.L.Wilson,B.D.Carlson.Radar detection in multipath.IET Radar Sonar Navigation,1999,146(1):45-52.

    [4]A.Zaimbashi,M.Derakhtian.A.Sheikhi.GLRT-based CFAR detection in passive bistatic radar.IEEE Trans.on Aerospace and Electronic Systems,2013,49(1):134-157.

    [5]J.Ding,H.Chen,B.You.Fluctuating targetdetection in lowgrazing angle with multi-input multiple-output radar.Proc. of the International Conference on Signal Processing,2012: 1714-1719.

    [6]H.T.Hayvaci,A.D.Maio,D.Erricolo.Improved detection probability of a radar targetin the presence of multipath with prior knowledge of the environment.IET Radar Sonar Navigation,2013,7(1):6-46.

    [7]Z.Tarique,W.Q.Malik,D.J.Edwards.Bandwidth requirements for accurate detection of direct path in multipath environment.Electronics Letters,2006,42(2):100-102.

    [8]K.B.Yu.Recursive super-resolution algorithm for lowelevation targetangle tracking in multipath.IET Radar Sonar Navigation,1994,141(4):223-228.

    [9]P.E.Berry,D.P.Finch.Low-angle target detection with interference in a multipath environment.Proc.of the International Conference on Radar,2013:440-445.

    [10]W.D.Blair,M.Brandt-Pearce.Statistics of monopulse measurements for tracking targets in the presence of sea-surface induced multipath.Proc.of the IEEE Aerospace Conference, 1998:49-58.

    [11]M.I.Skolnik.Introduction to radar System.3rd ed.New York: McGraw-Hill Companies,2006.

    [12]H.A.David,H.N.Nagaraja.Order statistics.3rd ed.New Jersey:Wiley,2003.

    [13]D.K.Barton.Low-angle radar tracking.Proceedings of the IEEE,1974,62(6):687-704.

    [14]S.Sen,A.Nehorai.Adaptive OFDM radar for target detection in multipath scenarios.IEEE Trans.on Signal Processing, 2011,59(1):78-90.

    [15]C.S.Chae,J.T.Johnson.A study of sea surface rangeresolved Doppler spectra using numerically simulated lowgrazing-angle backscatter data.IEEE Trans.on Geoscience and Remote Sensing,2013,51(6):3452-3460.

    [16]D.Poullin,M.Flecheux,M.Klein.Elevation angle estimation for low-altitude targets using DVB(SFN broadcasters).IEEE Aerospace and Electronic Systems Magazine,2012,27(1): 27-35.

    [17]J.He,G.Liao,Z.Yang.Moving targettrack detection in radar low angle environment.Proc.of the IET International Radar Conference,2013:1-4.

    [18]I.Walterscheid,A.R.Brenner,J.Klare.Radar imaging with very low grazing angles in a bistatic forward-looking con fi guration.Proc.ofthe IEEE InternationalGeoscience and Remote Sensing Symposium,2012:327-330.

    [19]A.Hizal,S.Koc.Mitigating the multipath effects oflow angle monopulse tracking by even difference pattern.Proc.of the European Radar Conference,2012:290-293.

    [20]Y.D.Zhang,J.J.Zhang,M.G.Amin.Maneuvering target altitude tracking in over-the-horizon radars exploiting multipath Dopplersignatures.Proc.of the IEEE International Conference on Acoustics,Speech and Signal Processing,2013: 4056-4060.

    Biographies

    Yunhe Cao was born in 1978.He received his Ph.D.degree from Xidian University in 2006.He is now an associate professor at National Laboratory of Radar Signal Processing,Xidian University.His research interests include radar array signalprocessing,wideband signalprocessing and MIMO radar signalprocessing.

    E-mail:cyh xidian@163.com

    Shenghua Wang was born in 1979.She received her B.S.degree and M.S.degree from Xidian University in 2001 and 2004,respectively.She is now a Ph.D.candidate ofthe National Laboratory of Radar Signal Processing,Xidian University.Her research interests include radar array signal processing and low angle tracking technique.

    E-mail:wshh 2011@163.com

    Yu Wang was born in 1991.He received his B.S. degree from Xidian University in 2013.He is now a Ph.D.candidate ofthe National Laboratory of Radar Signal Processing,Xidian University.His research interests include radar signal processing,parameter estimation and wideband array signalprocessing.

    E-mail:xdwangyu@163.com

    Shenghua Zhou was born in 1981.He received his Ph.D.degree from Xidian University in 2011.He is now an associate professor at National Laboratory of Radar Signal Processing,Xidian University.His research interests include multi-station radar signal processing,parameter estimation and target detection.

    E-mail:shzhou@mail.xidian.edu.cn

    10.1109/JSEE.2015.00032

    Manuscriptreceived July 04,2014.

    *Corresponding author.

    This work was supported by the National Natural Science Foundation of China(61372136;61372134;61172137),the Fundamental Research Funds for the Central Universities(K5051202005),and the China Scholarship Council(CSC).

    九色成人免费人妻av| 精品人妻视频免费看| 99久国产av精品| 欧美一区二区精品小视频在线| 最近在线观看免费完整版| 91久久精品电影网| 少妇的逼水好多| 两人在一起打扑克的视频| 男人的好看免费观看在线视频| 免费看a级黄色片| 成人欧美大片| 午夜免费激情av| 嫩草影院精品99| 深夜精品福利| 精华霜和精华液先用哪个| 在线免费观看的www视频| 午夜影院日韩av| 久久99热这里只有精品18| 精品久久久久久久久亚洲 | 国产精品乱码一区二三区的特点| 亚洲精品在线观看二区| 97超视频在线观看视频| 亚州av有码| 露出奶头的视频| 午夜亚洲福利在线播放| 久久久国产成人免费| 麻豆av噜噜一区二区三区| 久久久久精品国产欧美久久久| 91在线观看av| 动漫黄色视频在线观看| 中文字幕人成人乱码亚洲影| 一二三四社区在线视频社区8| 亚洲av第一区精品v没综合| 两个人视频免费观看高清| 91麻豆精品激情在线观看国产| 国产男靠女视频免费网站| 麻豆久久精品国产亚洲av| 亚洲乱码一区二区免费版| 久久这里只有精品中国| 精品久久久久久久久av| 99riav亚洲国产免费| 一本一本综合久久| 欧美中文日本在线观看视频| 免费在线观看日本一区| 国产精品98久久久久久宅男小说| 欧美日韩综合久久久久久 | 蜜桃久久精品国产亚洲av| 婷婷丁香在线五月| 国产69精品久久久久777片| 成人亚洲精品av一区二区| 午夜精品一区二区三区免费看| 国产伦一二天堂av在线观看| 国内精品一区二区在线观看| 亚洲成人中文字幕在线播放| 国产精品女同一区二区软件 | 亚洲av成人av| 亚洲欧美清纯卡通| 成熟少妇高潮喷水视频| 在线播放无遮挡| 一个人看的www免费观看视频| 亚洲国产日韩欧美精品在线观看| 免费人成在线观看视频色| 非洲黑人性xxxx精品又粗又长| 亚洲av免费在线观看| 国产男靠女视频免费网站| 中文字幕av成人在线电影| 国产精品1区2区在线观看.| 国产亚洲精品综合一区在线观看| 免费在线观看亚洲国产| 老司机午夜福利在线观看视频| 国产精品亚洲av一区麻豆| 国产淫片久久久久久久久 | 国产一区二区激情短视频| 久久久久久久午夜电影| 五月玫瑰六月丁香| 亚州av有码| 一级黄片播放器| 色综合亚洲欧美另类图片| 婷婷色综合大香蕉| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美日韩高清在线视频| 久久99热6这里只有精品| or卡值多少钱| 欧美午夜高清在线| 亚洲av五月六月丁香网| 啦啦啦观看免费观看视频高清| 国产一区二区三区在线臀色熟女| 99国产极品粉嫩在线观看| 国产精品亚洲美女久久久| 欧美成人免费av一区二区三区| 国产日本99.免费观看| 丰满的人妻完整版| av天堂在线播放| 中文字幕av成人在线电影| 免费高清视频大片| 亚洲内射少妇av| 国产一区二区在线av高清观看| 18禁在线播放成人免费| 国产黄a三级三级三级人| 久久精品影院6| 欧美一区二区亚洲| 国产成人aa在线观看| 中国美女看黄片| 欧美在线一区亚洲| 男女之事视频高清在线观看| 99在线视频只有这里精品首页| 十八禁国产超污无遮挡网站| 国产高潮美女av| 午夜a级毛片| 18禁在线播放成人免费| 在线看三级毛片| 性欧美人与动物交配| 国产成人aa在线观看| 久久久久九九精品影院| 18美女黄网站色大片免费观看| 欧美黄色片欧美黄色片| 国产乱人视频| 内射极品少妇av片p| 少妇的逼水好多| 高清在线国产一区| 国产色婷婷99| 老司机福利观看| 成人性生交大片免费视频hd| 久久久久久大精品| 国产精品一区二区三区四区久久| 成人av在线播放网站| 啦啦啦观看免费观看视频高清| 中文亚洲av片在线观看爽| 国产精品综合久久久久久久免费| 日韩欧美国产一区二区入口| 少妇熟女aⅴ在线视频| 精品国产亚洲在线| 欧美高清性xxxxhd video| 国产真实乱freesex| www.www免费av| 嫩草影院入口| 国产精品久久久久久人妻精品电影| 亚洲av熟女| 97热精品久久久久久| 亚洲国产精品成人综合色| 国产一区二区三区在线臀色熟女| 国产伦人伦偷精品视频| 麻豆成人午夜福利视频| 校园春色视频在线观看| 亚洲第一电影网av| 真人做人爱边吃奶动态| 成人一区二区视频在线观看| 别揉我奶头 嗯啊视频| 两人在一起打扑克的视频| 一本精品99久久精品77| 久久国产精品人妻蜜桃| 非洲黑人性xxxx精品又粗又长| 中出人妻视频一区二区| 禁无遮挡网站| 久久草成人影院| 天堂av国产一区二区熟女人妻| av天堂中文字幕网| 午夜福利在线在线| 亚洲18禁久久av| 丝袜美腿在线中文| 在线免费观看不下载黄p国产 | 一级毛片久久久久久久久女| 欧美激情在线99| 一进一出抽搐动态| 一个人看视频在线观看www免费| 国产主播在线观看一区二区| 18禁黄网站禁片午夜丰满| 国产精品亚洲美女久久久| 热99在线观看视频| 我要搜黄色片| 桃红色精品国产亚洲av| 免费观看的影片在线观看| 老女人水多毛片| 成人毛片a级毛片在线播放| 俄罗斯特黄特色一大片| 内地一区二区视频在线| 国产高清视频在线播放一区| 我要搜黄色片| 成人一区二区视频在线观看| 18禁黄网站禁片免费观看直播| 最后的刺客免费高清国语| 乱码一卡2卡4卡精品| 亚洲av成人不卡在线观看播放网| 国内毛片毛片毛片毛片毛片| 又黄又爽又刺激的免费视频.| av中文乱码字幕在线| 欧美在线一区亚洲| 日日摸夜夜添夜夜添小说| av在线观看视频网站免费| 久久午夜福利片| 真人做人爱边吃奶动态| 国产高清有码在线观看视频| 天堂影院成人在线观看| 成人av一区二区三区在线看| 搡老妇女老女人老熟妇| 国产在线男女| 日韩欧美三级三区| xxxwww97欧美| 免费无遮挡裸体视频| 久久人人爽人人爽人人片va | 一区二区三区激情视频| 淫秽高清视频在线观看| 久久伊人香网站| 成人国产一区最新在线观看| 男插女下体视频免费在线播放| 少妇的逼水好多| 黄色女人牲交| 欧美绝顶高潮抽搐喷水| 老司机午夜福利在线观看视频| 直男gayav资源| 国产aⅴ精品一区二区三区波| 在线观看舔阴道视频| 精品一区二区三区视频在线观看免费| 麻豆国产97在线/欧美| 国产亚洲av嫩草精品影院| 老女人水多毛片| 亚洲国产日韩欧美精品在线观看| 国产一级毛片七仙女欲春2| 99久久久亚洲精品蜜臀av| 久久久成人免费电影| 欧美一区二区国产精品久久精品| 99riav亚洲国产免费| 精品一区二区三区av网在线观看| 国产乱人视频| 日本精品一区二区三区蜜桃| 99久久无色码亚洲精品果冻| 成年人黄色毛片网站| 国产精品伦人一区二区| 亚洲色图av天堂| 丰满的人妻完整版| 美女xxoo啪啪120秒动态图 | 很黄的视频免费| 欧美+日韩+精品| 国产免费av片在线观看野外av| 久久午夜亚洲精品久久| 小蜜桃在线观看免费完整版高清| 午夜福利欧美成人| 欧美成狂野欧美在线观看| 国产精品人妻久久久久久| av在线观看视频网站免费| 国产欧美日韩精品亚洲av| 国产精品自产拍在线观看55亚洲| 欧美日韩乱码在线| 非洲黑人性xxxx精品又粗又长| 在线国产一区二区在线| 男人舔奶头视频| 丰满乱子伦码专区| 免费在线观看成人毛片| 91在线观看av| 亚洲欧美日韩无卡精品| av视频在线观看入口| 色综合欧美亚洲国产小说| 人妻制服诱惑在线中文字幕| 在线观看一区二区三区| 日本 欧美在线| 可以在线观看的亚洲视频| 亚洲色图av天堂| 最近视频中文字幕2019在线8| 日韩欧美国产在线观看| 久久国产乱子免费精品| 狂野欧美白嫩少妇大欣赏| 成年版毛片免费区| 国产高清有码在线观看视频| 看黄色毛片网站| 97超视频在线观看视频| 精品久久久久久久久久久久久| 国产午夜精品久久久久久一区二区三区 | 99视频精品全部免费 在线| 在线观看美女被高潮喷水网站 | 国产欧美日韩一区二区三| 国产69精品久久久久777片| 免费在线观看影片大全网站| 国产成人影院久久av| 又黄又爽又免费观看的视频| 日韩欧美免费精品| 亚洲成av人片在线播放无| 国产精品免费一区二区三区在线| 亚洲人成网站在线播| 国产精品影院久久| 最近在线观看免费完整版| 亚洲专区国产一区二区| 亚洲人成网站在线播| 韩国av一区二区三区四区| АⅤ资源中文在线天堂| 在线播放国产精品三级| 精品国产三级普通话版| 91麻豆av在线| 一级a爱片免费观看的视频| av女优亚洲男人天堂| 国产亚洲精品久久久久久毛片| 久久6这里有精品| 午夜免费男女啪啪视频观看 | 男人舔女人下体高潮全视频| 婷婷丁香在线五月| 美女免费视频网站| 男人狂女人下面高潮的视频| 久久欧美精品欧美久久欧美| 国产欧美日韩一区二区精品| 国产午夜福利久久久久久| 免费人成视频x8x8入口观看| 一个人看视频在线观看www免费| 色综合站精品国产| 日日夜夜操网爽| 99热只有精品国产| 日本一本二区三区精品| 国内揄拍国产精品人妻在线| 久久精品国产99精品国产亚洲性色| 欧美+亚洲+日韩+国产| 毛片女人毛片| 永久网站在线| 男人和女人高潮做爰伦理| 婷婷精品国产亚洲av| 日本免费a在线| 午夜视频国产福利| 国产精品三级大全| 成人一区二区视频在线观看| 国产一区二区在线av高清观看| 黄色丝袜av网址大全| 少妇被粗大猛烈的视频| 免费在线观看成人毛片| 嫩草影视91久久| 亚洲av熟女| 国产熟女xx| 久久精品综合一区二区三区| av专区在线播放| 97超视频在线观看视频| 99久久成人亚洲精品观看| 最新中文字幕久久久久| 免费黄网站久久成人精品 | av在线老鸭窝| 在线免费观看的www视频| 国产av一区在线观看免费| 精品人妻1区二区| 给我免费播放毛片高清在线观看| 亚洲国产高清在线一区二区三| 69av精品久久久久久| 一本综合久久免费| 一个人免费在线观看电影| 国产精品人妻久久久久久| 国产av不卡久久| 国产成人a区在线观看| 国产免费av片在线观看野外av| 老司机福利观看| 97人妻精品一区二区三区麻豆| 成人鲁丝片一二三区免费| 一区二区三区高清视频在线| 国产人妻一区二区三区在| 亚洲成人免费电影在线观看| 美女cb高潮喷水在线观看| 看十八女毛片水多多多| 久久亚洲真实| 一级黄色大片毛片| 免费人成视频x8x8入口观看| 国产乱人伦免费视频| 69av精品久久久久久| 又黄又爽又免费观看的视频| 久久久久久久精品吃奶| .国产精品久久| 91久久精品国产一区二区成人| 久久久久国产精品人妻aⅴ院| 深爱激情五月婷婷| 欧美激情国产日韩精品一区| 99国产综合亚洲精品| 9191精品国产免费久久| 亚洲av一区综合| 久久久久久久久久成人| 嫩草影院新地址| 特级一级黄色大片| 日日干狠狠操夜夜爽| 成人高潮视频无遮挡免费网站| 亚洲av第一区精品v没综合| 亚洲va日本ⅴa欧美va伊人久久| 99国产综合亚洲精品| 好看av亚洲va欧美ⅴa在| 18禁在线播放成人免费| 人妻制服诱惑在线中文字幕| 日韩av在线大香蕉| 国内精品一区二区在线观看| 精品久久久久久久末码| 国产精品亚洲一级av第二区| 午夜福利视频1000在线观看| 全区人妻精品视频| 麻豆国产av国片精品| 1000部很黄的大片| 99热6这里只有精品| 日本五十路高清| 欧美黄色淫秽网站| 久久国产乱子免费精品| 欧美极品一区二区三区四区| 国产高清激情床上av| 国产精品久久电影中文字幕| 亚洲综合色惰| 夜夜躁狠狠躁天天躁| 少妇的逼好多水| 夜夜躁狠狠躁天天躁| 少妇的逼好多水| 午夜视频国产福利| 制服丝袜大香蕉在线| 在线十欧美十亚洲十日本专区| 51午夜福利影视在线观看| 国产伦在线观看视频一区| 观看美女的网站| 久久久久久久久久成人| 毛片一级片免费看久久久久 | 日本免费一区二区三区高清不卡| 美女 人体艺术 gogo| 2021天堂中文幕一二区在线观| 欧美中文日本在线观看视频| 久久这里只有精品中国| 精品久久久久久久久久免费视频| 久久人人爽人人爽人人片va | 日本精品一区二区三区蜜桃| 国产精品野战在线观看| av在线天堂中文字幕| 99国产精品一区二区三区| 色播亚洲综合网| 国产黄片美女视频| 免费观看的影片在线观看| 悠悠久久av| 久久午夜亚洲精品久久| 波多野结衣高清无吗| 日本撒尿小便嘘嘘汇集6| 成人三级黄色视频| 日本在线视频免费播放| www.999成人在线观看| 午夜视频国产福利| 超碰av人人做人人爽久久| 亚洲国产精品合色在线| 日韩欧美精品免费久久 | 99热这里只有是精品在线观看 | 天堂网av新在线| 久久久久久国产a免费观看| 在线观看舔阴道视频| 最近中文字幕高清免费大全6 | 大型黄色视频在线免费观看| 欧美3d第一页| 国产伦精品一区二区三区视频9| 直男gayav资源| 精品午夜福利在线看| 精品久久久久久久人妻蜜臀av| 偷拍熟女少妇极品色| 国产精品一区二区三区四区免费观看 | 久久人人爽人人爽人人片va | 黄片小视频在线播放| 国产av麻豆久久久久久久| 舔av片在线| 黄色一级大片看看| 青草久久国产| 97碰自拍视频| 欧美另类亚洲清纯唯美| 我的女老师完整版在线观看| 成人精品一区二区免费| 国产高清有码在线观看视频| 亚洲 欧美 日韩 在线 免费| 精品无人区乱码1区二区| 婷婷色综合大香蕉| 男女那种视频在线观看| 中文字幕熟女人妻在线| 亚洲午夜理论影院| 嫩草影院新地址| 日本成人三级电影网站| 午夜精品久久久久久毛片777| 最好的美女福利视频网| 欧美高清性xxxxhd video| 亚洲av免费在线观看| 中文字幕av在线有码专区| 亚洲av第一区精品v没综合| 欧美日韩综合久久久久久 | 91麻豆av在线| 久久久久久久久久黄片| 在线观看舔阴道视频| 久久亚洲真实| 欧美一区二区精品小视频在线| 欧美中文日本在线观看视频| 久久中文看片网| 中文字幕高清在线视频| 午夜福利免费观看在线| 免费人成在线观看视频色| 国产精品一区二区三区四区久久| 国产在线精品亚洲第一网站| 国产高清视频在线观看网站| 亚洲一区二区三区不卡视频| 丝袜美腿在线中文| 色哟哟哟哟哟哟| 最新在线观看一区二区三区| 国产私拍福利视频在线观看| 99久久99久久久精品蜜桃| 亚洲黑人精品在线| 成人毛片a级毛片在线播放| 国产精品伦人一区二区| 亚州av有码| 亚洲成av人片在线播放无| 一个人观看的视频www高清免费观看| 精品国产亚洲在线| 老司机午夜十八禁免费视频| 国产蜜桃级精品一区二区三区| 可以在线观看毛片的网站| 亚洲av第一区精品v没综合| 国产精品久久久久久久电影| xxxwww97欧美| 久久精品国产清高在天天线| 国产亚洲精品综合一区在线观看| 久久天躁狠狠躁夜夜2o2o| 91久久精品电影网| 午夜福利成人在线免费观看| 国产精品一及| 中文字幕高清在线视频| 嫩草影院入口| 欧美+亚洲+日韩+国产| 国产精品久久电影中文字幕| 国产精品三级大全| 亚洲色图av天堂| 精品熟女少妇八av免费久了| 男女下面进入的视频免费午夜| 免费在线观看日本一区| 两个人视频免费观看高清| 日本黄大片高清| 亚洲一区高清亚洲精品| 久久久久久久久久黄片| 特级一级黄色大片| 男人和女人高潮做爰伦理| 深夜a级毛片| 国产精品一区二区免费欧美| 制服丝袜大香蕉在线| 成年女人看的毛片在线观看| 国产成人欧美在线观看| 波多野结衣高清作品| 两性午夜刺激爽爽歪歪视频在线观看| 欧美成狂野欧美在线观看| 国产真实伦视频高清在线观看 | 三级国产精品欧美在线观看| 女人被狂操c到高潮| 91麻豆精品激情在线观看国产| 中文字幕高清在线视频| 欧美xxxx黑人xx丫x性爽| 国产乱人视频| www.熟女人妻精品国产| 久久久久久大精品| 在线观看免费视频日本深夜| 成年免费大片在线观看| 日本精品一区二区三区蜜桃| 天天一区二区日本电影三级| av女优亚洲男人天堂| 男女下面进入的视频免费午夜| av专区在线播放| 久久亚洲真实| 亚洲精品在线观看二区| 国产又黄又爽又无遮挡在线| 深爱激情五月婷婷| 欧美日韩黄片免| 国产精品久久久久久久久免 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜久久久久精精品| 欧美性猛交╳xxx乱大交人| 午夜两性在线视频| 一级av片app| АⅤ资源中文在线天堂| 国产乱人伦免费视频| 非洲黑人性xxxx精品又粗又长| 香蕉av资源在线| 好看av亚洲va欧美ⅴa在| 给我免费播放毛片高清在线观看| 精华霜和精华液先用哪个| 成人欧美大片| 18美女黄网站色大片免费观看| 最好的美女福利视频网| 天堂网av新在线| 99国产精品一区二区蜜桃av| 日韩免费av在线播放| 直男gayav资源| 99国产综合亚洲精品| 午夜福利在线观看免费完整高清在 | 精品一区二区三区人妻视频| 中文字幕免费在线视频6| 成年女人永久免费观看视频| 在线a可以看的网站| 久久久久性生活片| 五月玫瑰六月丁香| 99久久九九国产精品国产免费| 国产成人aa在线观看| 久久久久久久久大av| 内射极品少妇av片p| 精品熟女少妇八av免费久了| 黄色女人牲交| 亚洲,欧美,日韩| 成人av在线播放网站| 久久久成人免费电影| 亚洲,欧美,日韩| 成人av在线播放网站| www.色视频.com| 熟妇人妻久久中文字幕3abv| 好男人在线观看高清免费视频| 热99在线观看视频| 亚洲一区二区三区不卡视频| 亚洲国产色片| 男女做爰动态图高潮gif福利片| 国产熟女xx| 亚洲国产色片| 他把我摸到了高潮在线观看| 内射极品少妇av片p| 国产成人aa在线观看| 免费电影在线观看免费观看| 欧美精品啪啪一区二区三区| 午夜激情福利司机影院| 日本精品一区二区三区蜜桃| 看十八女毛片水多多多| 日本黄色片子视频| 88av欧美| 亚洲av二区三区四区| 女人被狂操c到高潮| 国产一区二区激情短视频| 在线天堂最新版资源| 精品午夜福利视频在线观看一区| 露出奶头的视频| 我的老师免费观看完整版|