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

    Detection accuracy of target accelerations based on vortex electromagnetic wave in keyhole space

    2024-02-29 09:17:04KaiGuo郭凱ShuangLei雷爽YiLei雷藝HongPingZhou周紅平andZhongYiGuo郭忠義
    Chinese Physics B 2024年2期
    關(guān)鍵詞:忠義

    Kai Guo(郭凱), Shuang Lei(雷爽), Yi Lei(雷藝), Hong-Ping Zhou(周紅平), and Zhong-Yi Guo(郭忠義)

    School of Computer and Information,Hefei University of Technology,Hefei 230009,China

    Keywords: vortex electromagnetic waves,detect accelerations,keyhole space,spectrum spread

    1.Introduction

    Electromagnetic wave can carry not only spin angular momentum (SAM) but also orbital angular momentum(OAM),[1,2]which are related to circular polarization and helical phase structure.Electromagnetic wave carrying OAM is called as vortex wave[3]and shows huge advantages over conventional electromagnetic wave in practical application owing to its unique features.The orthogonality between different OAM modes provides extra degrees of freedom in addition to conventional physical components of electromagnetic wave, increasing the capacity of communication channel[4,5]and improving the resolution of imaging.[6–8]In addition,the coupling between OAM of electromagnetic wave and mechanical momentum has been utilized to manipulate particles.[9]In recent years, the rotational Doppler effect of vortex electromagnetic wave has offered a simple method to obtain threedimensional information of a moving object, such as velocity and acceleration.[10–12]So far, there have been extensive researches on vortex electromagnetic wave, including generation of high-quality OAM modes,[13,14]efficient detection method,[15]and its propagation characteristic through inhomogeneous medium.[16]

    Practically,environment has great influence on the transmission of electromagnetic vortex beams, causing crosstalk among OAM modes.In 2005,the concept of OAM spectrum was proposed in optics.[17]Subsequently,the OAM spectrum characteristics of vortex beams were extensively studied over a broadband electromagnetic spectrum in different transmission environments.On the one hand, the OAM spectrum spreading of optical beams has been carefully studied after passing through atmosphere[18,19]and seawater turbulences.[20,21]On the other hand,in the radio frequency range,the transmission of vortex electromagnetic wave is more obviously affected by the scattering of obstacles.[22–24]The scattering of high-order Bessel vortex electromagnetic waves on dielectric spheres was studied.[22]Fortunately,the OAM spectrum maintains well after the vortex electromagnetic wave has transmitted through conductor sphere and cone.[23]Very recently, it was reported that the off-axis vortex waves illuminate electrically large conductors of different shapes.[24]

    These results indicate that the OAM spectrum spreading due to inhomogeneous transmission environment will inevitably bring negative influence on the application of vortex beam.However,as a unique feature of vortex electromagnetic wave,the rotating Doppler effect and its application in detecting a rotating target have been studied mostly within the ideal line-of-sight scenarios until now.[25–30]Especially, the OAM spectrum spreading caused by non-line-of-sight scenario and non-cooperative target simultaneously is difficult to predict directly.[12]Therefore,it is urgent in practice to study the detection of rotating targets in non-line-of-sight scenarios.A“keyhole space”model for a typical non-line-of-sight scenario has been investigated in communication system with vortex electromagnetic wave.[31–33]For example,for indoor users to detect outdoor targets,the opacity of the wall to the signal will cause energy to leak out only through windows, doors, and other keyholes.Even in outdoor scenarios with abundant scattering environment,the keyhole effect can occur when the ring of the scattering object is smaller than the distance between the transmitter and the receiver.[34]However,the performance of detecting rotating targets with vortex electromagnetic wave still remains to be studied in the keyhole space model.

    In this work, we use the vortex electromagnetic wave in the keyhole space to detect a rotating target.The rest of this paper is organized as follows.In Section 2 the detection principle and analysis method are described.In Section 3 are studied the OAM spectrum spreading of vortex electromagnetic wave in several typical keyhole space situations.In Section 4,the detection errors of both longitudinal and angular accelerations of the target are calculated and analyzed.Finally, some conclusions are drawn from the present study in Section 5.

    2.Detection principle and analysis method

    Figure 1 schematically shows the detection of the rotating target by using the model of vortex electromagnetic wave in the keyhole space.The vortex electromagnetic wave is generated by an eight-arm Archimedes spiral antenna (ASA)[35]with a radius of 150 mm,propagates along thezaxis,and encounters the keyhole space that is symmetric about thezaxis 450 mm away.Considering the common diffraction aperture,a circular keyhole with radiusRis studied and its surrounding material is set as a perfect electric conductor(PEC).[33]After passing through the keyhole space,the vortex electromagnetic wave with spread OAM spectrum reaches the target 900 mm away, which rotates perpendicular to thezaxis with angular velocity?0+g′tand moves along thezaxis with longitudinal velocityv0+gt.?0andv0are the initial angular velocity and longitudinal velocity, respectively,g′andgare the angular acceleration and longitudinal acceleration,respectively.In the variable acceleration case,the longitudinal variable acceleration of the target isg(t)=g0+pt,and the angular variable acceleration isg′(t)=g′0+p′t,withpandp′being longitudinal variable acceleration and angular variable acceleration,respectively.Without loss of generality,the target can be treated as a scattering point in the far field.

    Fig.1.Detection model of longitudinally and angularly accelerated target with vortex electromagnetic wave in keyhole space.

    In a line-of-sight scenario,the signal received by antenna can be expressed as[12]

    wherefcis the carrier frequency,kis the wave number,ris the distance between the antenna and the target,s(t)is the relative distance difference,?(t)is the relative azimuth difference,σis the scattering coefficient,?0is the initial azimuth,andais the antenna radius.

    Based on Eq.(1)and time-frequency analysis method.we use a commonly used time-frequency analysis method to analyze time-varying non-stationary signals.It has an exponential kernel function, which can effectively suppress the interference term and maintain a high time-frequency resolution.[12]The slope of the time-frequency signal could be expressed as

    where the first part(Eq.(2a))and second part(Eq.(2b))correspond to uniform acceleration and variable acceleration, respectively, ?f(t) is the Doppler frequency shift,gandg′are the longitudinal and angular uniform acceleration,respectively,pandp′are the longitudinal variable acceleration and angular variable acceleration, respectively, and the value ofmcould be estimated by fitting method.When transmitting signals of different OAM modes with topological charges ofl1andl2simultaneously, the angular uniform acceleration(AUA),longitudinal uniform acceleration(LUA),the angular variable acceleration(AVA)and longitudinal variable acceleration(LVA)of target can be expressed as

    However, in keyhole space, OAM beam is blocked and does not fully illuminate the object, which results in the spread of OAM spectrum.Therefore, the transmitting signal when reaching the target can be expressed as

    whereAlis the weight of the OAM mode with the topological chargel.

    Without considering the attenuation of the transmission path, the frequency shift information carried by the moving target in the echo signal is only time-dependent.When the phase is differentiated with respect to time, all the phase effects independent of time are eliminated.Therefore,the slope of the time-frequency signal for the uniform acceleration can be expressed as

    and the slope of the time-frequency signal for the variable acceleration be expressed as

    We obtain the slopem′that can be expressed by using data fitting that is affected by different values ofm.The AUA,LUA,AVA,and LVA of target can be expressed as

    Both the keyhole space and the attitude of the target will lead to the impurities of the OAM modes and increased errors in detection of angular and longitudinal accelerations.We will study these issues in detail in future work.

    3.Results and discussion

    3.1.OAM spectra spreading

    Without loss of generality,the OAM spectrum spreading of vortex electromagnetic wave with the topological charge ofl=1 andl=2 atf=2.45 GHz are analyzed under different keyhole conditions by numerical simulation.In an actual detection,we only need to predict the error range of the acceleration and do not need to obtain a specific spectrum spread.However,in order to analyze the source of the errors,we measure the OAM spectrum spreading of vortex electromagnetic wave by using a uniform circular antenna (UCA) array with 256 sampling points.[15]

    Firstly, we investigate the keyhole spaces with different radii ofR=0.5λ,R=0.7λ,R=λ,R=1.5λ.Besides, two other cases are studied: one is the keyhole radius equal to that of the transmitting antenna ASA(rof ASA),and the other is the free space.The thickness of the keyhole is 5 mm[33]and the radius of circular observation surface is 85 mm, which is equal to the radius of the radiating area of the transmitting antenna ASA.

    Figures 2(a) and 2(b) show the mode purities of OAM beams with topological charge ofl=1 andl=2,respectively,after passing the keyhole spaces with different radii.The purities of main modes become higher when the radius of keyhole increases, because the influence of diffraction at the keyhole on the spiral phase construction of OAM beams becomes less.However,figure 2(b)shows that the OAM spectra spreading is more obvious for the case withl=2.It is because the energy of the OAM beams mainly concentrates into a ‘donut’, most of energy locates at a ring with radius of 90 mm and 170 mm for the case withl=1 andl=2,respectively, at the keyhole space.Therefore,the destruction of spiral phase is more serious withl=2 for the same radius of keyhole.

    Fig.2.OAM spectra of vortex electromagnetic waves with topological charge of l=1(a)and l=2(b),exiting from Archimedes spiral antenna and after passing a keyhole space with different radii.

    Secondly,the thickness of the keyhole is investigated,and the value ofwis set to 5 mm, 30 mm, 60 mm, 100 mm, and 200 mm in simulation.Without loss of generality, the radius of keyhole is set asλand the radius of circular observation surface is 85 mm.

    Figures 3(a) and 3(b) show the mode purities of OAM beam withl=1 andl=2,respectively,after passing the keyhole space with different thicknesses.The OAM beam withl= 1 can pass through the keyhole completely due to the small size of its energy‘donut’.Therefore,the OAM spectrum spreading is not obvious withl=1 as the thickness increases.For the OAM beam withl=2,the range of energy can exceed the size of keyhole,fortunately,the divergence of OAM beam will be limited during propagation in the keyhole space,resulting in a focusing effect.Therefore,the purity of main mode is improved more obviously withl=2.

    Fig.3.OAM spectra of vortex electromagnetic wave with l=1(a)and l=2(b)after passing keyhole with different thicknesses.

    Thirdly,it is considered that the attitude of the target may exert an influence on OAM mode in the keyhole space.As shown in Fig.1,the target is moving along the UCA array.The radius of keyhole is set asλand its thickness is 5 mm,and different radii of target trajectory are investigated:R′=85 mm,105 mm,150 mm,250 mm,and 400 mm.

    Figures 4(a) and 4(b) show the mode purities of OAM beam withl=1 andl=2,respectively,in the case of different radii of the target trajectories.It can be seen that the purity of the main mode withl=1 first increases and then decreases gradually, while the purity of main mode withl=2 is improved obviously.Because most of energy is located at a ring with radius of 130 mm and 300 mm for the case withl=1 andl=2, respectively, at the observation surface.If the target rotates at the ring where the OAM beam has the strongest energy,the purity of the main mode is the best.

    Finally, the centers of the target at different off-axis distances are investigated.As shown in Fig.1, its centers are located at(0,d, 900 mm), wheredindicates the off-axis distance.In simulation, the value ofdis set to 0, 100 mm,200 mm, 300 mm, and 400 mm.The thickness of the keyhole is 5 mm and the radius of circular observation surface is 85 mm.

    Fig.4.OAM spectra of vortex electromagnetic wave with l=1(a)and l=2(b)after passing keyholes with different radii of target trajectory.

    Fig.5.OAM spectra of vortex electromagnetic wave with l=1(a)and l=2(b)after passing through keyhole when target is off-axis.

    Figures 5(a) and 5(b) show the mode purities of OAM beam withl=1 andl=2,respectively,when the target is offaxis.In the cases withl=1 andl=2,the OAM spectra spread both obviously,consequently,purities of their adjacent modes increase.On the one hand,the purity of main mode decreases gradually withl=1 as thedincreases from 0 mm to 400 mm,because the trajectory of the target moves away from the ring of the OAM beam gradually.On the other hand, for the case withl=2, the target does not rotate at the ring whend=0 due to the larger range of the hollow region and moves then far away from the ring whendexceeds 200 mm.Therefore,the purity of main model=2 is improved whendincreases from 0 to 200 mm and then decreases whendincreases from 200 mm to 400 mm.

    3.2.Measurement errors and analysis

    In this subsection, the detection errors of a target that moves along the observation surface with a longitudinal acceleration and an angular acceleration in four situations are investigated by using composite OAM beam withl=±1 andl=±2.We consider the OAM modes with topological charge ofl=±1 andl=±2 for the following reasons.Firstly,the OAM spectrum is dependent on the geometrical relation among the radiation region of vortex electromagnetic wave,the keyhole and the target.For higher-order OAM modes,the variation trend of OAM spectrum will not change.Secondly,it can be easier to generate lower-order vortex electromagnetic waves with high quality at the radio frequency.Meanwhile,the using of OAM modes with opposite topological charges makes it easy to analyze the OAM spectrum and investigate the detection errors.For the uniform acceleration case, the LUA and AUA are set asg=15 m/s2andg′=100πrad/s2,respectively, with initial velocity ofv0=0 and?0=0.For the variable acceleration case, the LVA and AVA arep=30 andp′=300(m/s3),respectively,withg0=0 andg′0=0.

    Firstly,figures 6(a)and 6(b)show the detection errors for the cases that the different radii of keyhole withl=±1 andl=±2,respectively.As shown in Fig.6(a),the detection errors of the AUA and the AVA with different radii of keyhole keep the same as those without keyhole existing, for the energy of the OAM beam withl=1 andl=-1 illuminate the rotating target completely.The detection errors of AUA are all below 5%, while the detection errors of AVA keep below 15%withl=±1.Meanwhile,all the detection errors of LUA and LVA remain below 1%.Figure 6(b) shows the detection errors for thel=2 andl=-2.It is known that the higher mode will bring larger rotational Doppler shift,resulting in the lower detection errors of AUA and AVA.As a result, the detection errors of angular accelerations can be reduced by using higher modes when the radius of keyhole is larger.It can be seen that the errors of the AUA decrease from 100.4%to 1.3%and the errors of AVA fall from 80.4% to 3.78% gradually.It can be well explained by the fact that when the aperture is small,most of the energy of the higher-order mode is reflected back,and only a small potion of radiation irradiates the target,thus causing purity of mode to diverge into the value in the casel=0.Then the purity of main mode is increased obviously as the aperture size turns larger.When the aperture size is infinite, it is equivalent to the transmission of vortex electromagnetic wave in free space, and the higher the order of OAM mode used for detection,the smaller the errors of angular acceleration is.[12]Therefore,for angular accelerations,we should select the maximum OAM mode so that the detection errors can be reduced through using the aperture.However,the errors of LUA and LVA do not change linearly.According to Eq.(3), the detection of LUA and LVA are affected by not only the slopem, but also the detection values of AUA and AVA.In fact,we can emit a spherical wave withl=0 to eliminate the influence of obstacles on the vortex electromagnetic wave spectrum spreading, which are discussed in the following cases.

    Secondly, figures 6(c) and 6(d) show the detection errors for the different thickness of keyhole withl=±1 andl=±2, respectively.It can be found that the distribution of OAM spectrum is similar to those of different aperture sizes in Figs.2 and 3.Therefore, the changing trends of errors are the same in two cases.Due to the small range of energy withl=±1,the purities of main modes are almost unaffected when passing through keyholes with different thickness.As shown in Fig.6(c),the detection errors of AUA are below 5%,those of AVA keep below 15%,and those of LUA and LVA are below 1.5%.Figure 6(d)shows the plots of the detection errors forl=2 andl=-2, the errors of AUA are decrease from 89.325%to 1.25%and the errors of AVA fall from 72.48%to 5.45%apparently as the thickness of keyhole increases.It can be seen in Fig.3(b)that the larger aperture thickness limits the divergence of the higher-order modes and improves the purity of the main mode withl=2.Therefore, the errors of angular accelerations decrease gradually.Interestingly, it is noted that the detection errors increase a little whenw=200 mm compared with whenw=100 mm.It is because the higher modes can bring larger Doppler frequency shifts,which indicates whether the frequency shift fitting exceeds the display range in the time–frequency diagrams should be considered when the purities of higher modes increase to a certain value.It can be reduced by increasing the sampling time or adding a deblurring function before signal processing.

    Thirdly, figures 6(e) and 6(f) show the detection errors for the cases that the different radii of target trajectories withl=±1 andl=±2, respectively.For detecting the target accelerations, the detection errors can be minimized only when the target is illuminated by the main radiation region of the vortex electromagnetic wave.It can be seen in Fig.6(e) that the errors of AUA increase from 4%to 89.4%and the errors of AVA are raised from 13.17%to 106.49%withl=±1.Meanwhile, the errors of LUA and LVA remain above 20% whenR′=400 mm.It can be explained that when the target moves outside the main radiation region of OAM beam,the purity of main mode spreads significantly to higher-order modes.As shown in Fig.4(a), the purity ofl=1 is lowest and the purities of the high-order modes such asl=3 increase obviously whenR′=400 mm.According to Eq.(5), the higher-order modes have a strong influence on the slope estimation,which causes the estimated value ofm1to increase.Therefore, the detection errors increase apparently whenR′=400 mm withl=±1.Figure 6(f)shows the plots of the errors for the case ofl=±2.It can be seen that the errors of AUA decrease to 1.625%and those of AVA fall to 7.2%due to the fact that the purities of main modes increase as the value ofR′rises.Therefore,for different radii of target trajectory,we should comprehensively consider the radiation region of the vortex electromagnetic wave and the size of the aperture, then choose the appropriate OAM mode, so that the target can be illuminated by the energy of the vortex electromagnetic wave as much as possible,which can reduce the diffusion of main OAM purity.

    Fig.6.Errors of uniform and variable accelerations,caused by difference among keyhole radii with l=±1(a)and l±2(b),difference among keyhole thickness with l=±1(c)and l=±2(d),difference among radii of target trajectories with l=±1(e)and l=±2(f),and difference among offset positions of target with l=±1(g)and l=±2(h).

    Finally, figures 6(g) and 6(h) show the detection errors for the cases that the different off-axis positions of the target withl=±1 andl=±2, respectively.It is shown that the OAM spreading is obvious withl=1 for the case that position of the target is off-axis due to only part of the trajectory is illuminated by the main radiation region of the OAM beam.As a result, figure 6(g)shows that the detection errors of AUA increase from 3.65%to 78.4%and those of AVA rise from 13.03%to 83.13%apparently withl=±1.Meanwhile,the errors of LUA and LVA reach 5% asdrises.Figure 6(h)shows the detection errors for the cases ofl=±2.It is worth noting that the detection errors of AUA and AVA decrease to 3.9% and 6.23%, respectively, even though the OAM spectrum spreading is also obvious whend=200 mm, 300 mm,and 400 mm.According to Eq.(5),l=3 andl=1 may have opposite effects on the slope estimation withl=2.Therefore, errors in the detection of angular accelerations may be rather lower withl=±2.To demonstrate it, figure 7 shows the time–frequency diagrams ofd=400 mm withl=±2.As shown before, there are many interference items in Figs.7(a)and 7(b), but the errors of accelerations are rather small in Fig.6(h).

    Fig.7.Time–frequency diagrams of echo signal of (a) uniformly and(b) variably accelerated target with l =±2 when d =400 mm is detected.

    In order to consider the practicability of the results, we change the value of signal-to-noise ratio (SNR) to 15 dB.Without loss of generality, the cases of different radii of keyhole are studied, and the detection errors are calculated.As shown in Fig.8, the changing trends of detection errors are stable at an SNR of 15 dB in comparison with the results in Figs.6(a)and 6(b).

    Fig.8.Errors of uniform and variable accelerations at signal-to-noise ratio of 15 dB,caused by difference among radii of keyhole with l=±1(a)and l=±2(b).

    According to the above analysis, the discrete time–frequency signal leads to the detection errors of the target accelerations,which are caused by the spread of the OAM spectrum when OAM beam does not completely illuminate the rotating object.In fact,the detection ability of the OAM beam is determined by the relative position among the main radiation region of the vortex electromagnetic wave,the keyhole and rotating object.When the selected OAM mode is not blocked by the keyhole and the rotating object moves in its main radiation region,the detection errors of the target accelerations are minimized.Therefore,by pre-measuring the amplitudes of OAM beams with different modes transmitted by the antenna and by comparing with the aperture size and the attitude of the target, then the reasonable selection of the OAM modes, which are less affected by keyhole and illuminate the trajectory of the detected target as much as possible,can effectively reduce the error of angular acceleration,because the error is directly affected by the purities of main OAM modes.For example,when the size of the aperture or the trajectory of the target is large, we need to use higher-order modes to obtain higher OAM purities and thus smaller angular acceleration errors are obtained.Instead,smaller OAM modes should be used to ensure the accuracy of the detection.For the longitudinal acceleration of the target, we can emit vortex electromagnetic wave withl=0 to detect target and eliminate the influence of OAM spectrum, from which the obtained conclusion is obviously different from the conclusion that the detection accuracy by using high-order OAM modes is higher than that of loworder OAM modes when the target accelerations are detected in free space.[12]

    In particular,we consider a special case that the acceleration is set to 0.According to Eq.(2),wheng=0 andg′=0,the Doppler shift formula can be expressed as

    The longitudinal velocity and the angular velocity are set asv0=1.5 m/s and?0=40πrad/s,respectively.

    For example, the OAM spectral spreading ofR=0.5λwithl=1 in Fig.2(a)is chosen.When the fast Fourier transform(FFT)is performed on the echo signal with respect tot,the frequency spectrum is shown in Fig.8.It is found that the angular velocity and longitudinal velocity can be calculated by transmitting a single OAM mode.The angular velocity can be accurately calculated by using the frequency shift difference (20 Hz) between adjacent modes, and the longitudinal velocity can be calculated by frequency shift(24.5 Hz)of the component ofl=0 in simulation.[36]

    Fig.9.Simulation frequency spectrum of the echo when v0 =1.5 m/s and ?0=40π rad/s.

    Finally, it is necessary to point out that spectral components may be missing in the experiment.[37]To reduce the errors, the minimum frequency shift difference should be used to estimate the angular velocity.

    4.Conclusions

    In conclusion,this work studies the errors in detection of accelerations in keyhole space.Based on the effect of OAM spectrum spreading which is simulated and analyzed, the errors of accelerations of uniformly and variably accelerated targets in keyhole space under different composite OAM beams are calculated and compared by the multi-mode method.The results show that the errors can be reduced effectively by using a combination of appropriate modes in different cases.The change in keyhole and attitude of the target will embody in the relativeness between radiation region of vortex electromagnetic wave and the attitude of receiving surface.

    Acknowledgement

    Project supported by the National Natural Science Foundation of China(Grant Nos.11804073 and 61775050).

    猜你喜歡
    忠義
    玉米轉(zhuǎn)錄因子ZmbHLH91對非生物逆境脅迫的應(yīng)答
    傳話游戲
    數(shù)據(jù)挖掘技術(shù)在內(nèi)河航道維護(hù)管理中的應(yīng)用研究
    正史忠義類傳的新變化
    ——論《宋史·忠義傳》的纂修體例、入傳標(biāo)準(zhǔn)及文學(xué)色彩
    這條忠義街
    “忠義”立新人
    抱樸立新 德行忠義——江陰市立新中心小學(xué)“新忠義的力量”少先隊(duì)活動掠影
    全國人大代表米忠義
    朋友的情誼
    Enhancing Structural Stability and Pervaporation Performance of Composite Membranes by Coating Gelatin onto Hydrophilically Modified Support Layer*
    亚洲第一区二区三区不卡| 免费久久久久久久精品成人欧美视频| 91精品国产国语对白视频| 亚洲精品成人av观看孕妇| 人人妻,人人澡人人爽秒播 | 99精国产麻豆久久婷婷| 国产一区二区在线观看av| 久久久久久久精品精品| 免费黄色在线免费观看| 一区福利在线观看| 精品国产露脸久久av麻豆| 成人黄色视频免费在线看| 久久久久精品国产欧美久久久 | 日韩大片免费观看网站| av国产久精品久网站免费入址| 国产一级毛片在线| 国产毛片在线视频| 免费高清在线观看视频在线观看| 男女国产视频网站| 亚洲欧美精品综合一区二区三区| 久久狼人影院| 成人毛片60女人毛片免费| 最新在线观看一区二区三区 | 亚洲欧美中文字幕日韩二区| 久久国产精品男人的天堂亚洲| 成人亚洲欧美一区二区av| 日韩伦理黄色片| 免费观看性生交大片5| 妹子高潮喷水视频| 男女床上黄色一级片免费看| 久久精品亚洲av国产电影网| 久久久久精品国产欧美久久久 | 伦理电影免费视频| 一二三四在线观看免费中文在| 日韩视频在线欧美| 午夜福利网站1000一区二区三区| 国产男人的电影天堂91| av不卡在线播放| 女人高潮潮喷娇喘18禁视频| 99精国产麻豆久久婷婷| 一本色道久久久久久精品综合| 精品国产国语对白av| 亚洲专区中文字幕在线 | 又黄又粗又硬又大视频| 无遮挡黄片免费观看| 欧美最新免费一区二区三区| av片东京热男人的天堂| 人妻人人澡人人爽人人| 丝袜喷水一区| 久久人妻熟女aⅴ| 七月丁香在线播放| 欧美乱码精品一区二区三区| 另类精品久久| 狠狠精品人妻久久久久久综合| 在线观看人妻少妇| 国产成人精品在线电影| 精品久久久久久电影网| 日韩欧美精品免费久久| 国产又爽黄色视频| 欧美精品一区二区大全| 免费av中文字幕在线| 亚洲一区中文字幕在线| 秋霞伦理黄片| 这个男人来自地球电影免费观看 | 国产成人一区二区在线| 国产精品蜜桃在线观看| 国产精品香港三级国产av潘金莲 | 亚洲成国产人片在线观看| 我要看黄色一级片免费的| 在线观看人妻少妇| 人成视频在线观看免费观看| svipshipincom国产片| 国产 精品1| 秋霞伦理黄片| 成人手机av| 黄片播放在线免费| 晚上一个人看的免费电影| a级片在线免费高清观看视频| 又大又爽又粗| 免费观看性生交大片5| 中文字幕制服av| 亚洲国产精品999| 亚洲国产精品成人久久小说| 啦啦啦在线观看免费高清www| 高清av免费在线| 9色porny在线观看| 97人妻天天添夜夜摸| 精品久久久精品久久久| 高清在线视频一区二区三区| 成人免费观看视频高清| 国产一区二区三区av在线| 久久久国产一区二区| 亚洲成人手机| 日日啪夜夜爽| 国产熟女欧美一区二区| 日韩制服骚丝袜av| 午夜免费鲁丝| 欧美 日韩 精品 国产| 免费看不卡的av| 国产免费福利视频在线观看| 久久久久视频综合| 色婷婷久久久亚洲欧美| 久久婷婷青草| 精品福利永久在线观看| 国产一区亚洲一区在线观看| 女性被躁到高潮视频| 欧美av亚洲av综合av国产av | 国产1区2区3区精品| 黑人猛操日本美女一级片| 高清黄色对白视频在线免费看| 国产精品久久久av美女十八| 一边摸一边抽搐一进一出视频| 久久精品久久精品一区二区三区| 人人澡人人妻人| 免费av中文字幕在线| 一级毛片电影观看| 久久ye,这里只有精品| 黄片播放在线免费| 亚洲色图综合在线观看| 少妇猛男粗大的猛烈进出视频| 男男h啪啪无遮挡| 无限看片的www在线观看| 国语对白做爰xxxⅹ性视频网站| 欧美日韩精品网址| 亚洲精华国产精华液的使用体验| 在线观看免费日韩欧美大片| 久久97久久精品| 日韩制服丝袜自拍偷拍| 国产伦人伦偷精品视频| √禁漫天堂资源中文www| 久热爱精品视频在线9| 久久天堂一区二区三区四区| 黄片无遮挡物在线观看| 精品少妇内射三级| 国产日韩欧美在线精品| 男女边摸边吃奶| 在线观看国产h片| 国产一级毛片在线| 国产男人的电影天堂91| 国产 精品1| 亚洲视频免费观看视频| 别揉我奶头~嗯~啊~动态视频 | 多毛熟女@视频| 久久久久精品久久久久真实原创| 青春草视频在线免费观看| 精品国产一区二区三区四区第35| 国产乱人偷精品视频| 欧美 日韩 精品 国产| 最近的中文字幕免费完整| 精品国产国语对白av| 极品人妻少妇av视频| 在线观看免费午夜福利视频| 一本大道久久a久久精品| 国产又爽黄色视频| 电影成人av| 欧美精品av麻豆av| 女人被躁到高潮嗷嗷叫费观| 国产在线视频一区二区| 操美女的视频在线观看| 在现免费观看毛片| 久久久精品免费免费高清| 只有这里有精品99| 久久韩国三级中文字幕| 国产日韩欧美在线精品| 日韩大码丰满熟妇| 国精品久久久久久国模美| 丝瓜视频免费看黄片| 丝袜美腿诱惑在线| svipshipincom国产片| 亚洲三区欧美一区| 90打野战视频偷拍视频| 啦啦啦在线观看免费高清www| 97精品久久久久久久久久精品| 久久av网站| 国产极品粉嫩免费观看在线| 热99国产精品久久久久久7| 在线观看免费视频网站a站| 日韩大码丰满熟妇| av卡一久久| 日本91视频免费播放| 美女大奶头黄色视频| 国产高清国产精品国产三级| 欧美日韩亚洲高清精品| 国产免费福利视频在线观看| 国产一区二区 视频在线| 性色av一级| 国产一区二区激情短视频 | 亚洲av欧美aⅴ国产| 男女高潮啪啪啪动态图| 男人舔女人的私密视频| 永久免费av网站大全| 久久久久人妻精品一区果冻| 亚洲美女搞黄在线观看| 午夜福利网站1000一区二区三区| 日韩成人av中文字幕在线观看| 精品亚洲乱码少妇综合久久| 亚洲av国产av综合av卡| 丰满少妇做爰视频| 亚洲精华国产精华液的使用体验| 99热国产这里只有精品6| 高清欧美精品videossex| 欧美黑人精品巨大| av福利片在线| 亚洲国产精品成人久久小说| 国产在线免费精品| 一边亲一边摸免费视频| 欧美成人精品欧美一级黄| 成人亚洲欧美一区二区av| 一边亲一边摸免费视频| 熟妇人妻不卡中文字幕| 母亲3免费完整高清在线观看| 我的亚洲天堂| 天天躁狠狠躁夜夜躁狠狠躁| 黄片无遮挡物在线观看| www日本在线高清视频| 久久亚洲国产成人精品v| 菩萨蛮人人尽说江南好唐韦庄| 不卡av一区二区三区| kizo精华| 最近最新中文字幕免费大全7| 少妇 在线观看| 午夜福利免费观看在线| 国产亚洲最大av| 999久久久国产精品视频| 亚洲欧美精品自产自拍| 啦啦啦中文免费视频观看日本| 欧美av亚洲av综合av国产av | 日日啪夜夜爽| 中文欧美无线码| 亚洲国产最新在线播放| 亚洲成av片中文字幕在线观看| 最新在线观看一区二区三区 | 嫩草影视91久久| 国产成人精品无人区| 免费观看性生交大片5| 天天影视国产精品| 日本午夜av视频| 热99国产精品久久久久久7| 99国产精品免费福利视频| 最新在线观看一区二区三区 | av在线app专区| 1024视频免费在线观看| 一本大道久久a久久精品| 欧美黑人欧美精品刺激| www.精华液| 中文字幕亚洲精品专区| 制服丝袜香蕉在线| 国产1区2区3区精品| 久久精品aⅴ一区二区三区四区| 又大又黄又爽视频免费| 一区二区三区激情视频| 亚洲欧美精品自产自拍| 热re99久久国产66热| 国产av一区二区精品久久| 久久久久精品久久久久真实原创| 亚洲熟女毛片儿| 国产av精品麻豆| 国产亚洲最大av| 精品亚洲成国产av| 美女主播在线视频| 免费黄色在线免费观看| 狠狠婷婷综合久久久久久88av| 欧美久久黑人一区二区| 国产xxxxx性猛交| 国产黄频视频在线观看| 在线 av 中文字幕| 国产精品久久久人人做人人爽| 老鸭窝网址在线观看| 久久97久久精品| 国产欧美日韩一区二区三区在线| 51午夜福利影视在线观看| 欧美日韩视频高清一区二区三区二| 免费黄频网站在线观看国产| 99香蕉大伊视频| 欧美国产精品va在线观看不卡| 九草在线视频观看| 丰满迷人的少妇在线观看| 精品亚洲成国产av| 视频区图区小说| 亚洲精品成人av观看孕妇| 国产精品av久久久久免费| 久久久久久久精品精品| 日韩精品免费视频一区二区三区| 日本vs欧美在线观看视频| 国产一级毛片在线| 黄片无遮挡物在线观看| 日韩一区二区三区影片| 久久久久久久久久久久大奶| 久久久久网色| 免费观看av网站的网址| 精品国产国语对白av| 十八禁网站网址无遮挡| 久久久久久久久久久久大奶| 老司机深夜福利视频在线观看 | 午夜福利,免费看| 美女扒开内裤让男人捅视频| 国产成人系列免费观看| 性少妇av在线| 亚洲图色成人| 又大又黄又爽视频免费| 超碰97精品在线观看| 伊人亚洲综合成人网| 欧美日韩国产mv在线观看视频| 亚洲精品美女久久久久99蜜臀 | 日韩一区二区三区影片| 亚洲精品日韩在线中文字幕| 久热爱精品视频在线9| 最近最新中文字幕大全免费视频 | 日日啪夜夜爽| 免费在线观看黄色视频的| 国产精品一二三区在线看| 国产成人精品无人区| 午夜福利影视在线免费观看| 国产熟女欧美一区二区| 日韩免费高清中文字幕av| 19禁男女啪啪无遮挡网站| 久久ye,这里只有精品| 精品一品国产午夜福利视频| 中文字幕色久视频| 一个人免费看片子| av片东京热男人的天堂| 人人妻,人人澡人人爽秒播 | 亚洲欧洲日产国产| 男女无遮挡免费网站观看| 国产成人欧美在线观看 | 国产av码专区亚洲av| 久久韩国三级中文字幕| videosex国产| 中文字幕高清在线视频| 女人久久www免费人成看片| 校园人妻丝袜中文字幕| 亚洲精华国产精华液的使用体验| 久久久久人妻精品一区果冻| 黄片小视频在线播放| 国产亚洲午夜精品一区二区久久| 国产精品香港三级国产av潘金莲 | 国产精品嫩草影院av在线观看| 色网站视频免费| 少妇猛男粗大的猛烈进出视频| 深夜精品福利| 青青草视频在线视频观看| 久久婷婷青草| 国产亚洲午夜精品一区二区久久| 国产av一区二区精品久久| av又黄又爽大尺度在线免费看| 国产亚洲av片在线观看秒播厂| 国产国语露脸激情在线看| 老汉色∧v一级毛片| 午夜日本视频在线| 精品少妇一区二区三区视频日本电影 | 香蕉丝袜av| 亚洲婷婷狠狠爱综合网| 国产精品av久久久久免费| 国产成人a∨麻豆精品| 黄色一级大片看看| 青青草视频在线视频观看| 欧美激情 高清一区二区三区| 精品一区在线观看国产| 国产又色又爽无遮挡免| 中国国产av一级| 亚洲精品一二三| 丰满乱子伦码专区| 1024香蕉在线观看| 成人18禁高潮啪啪吃奶动态图| 宅男免费午夜| 天美传媒精品一区二区| 91精品三级在线观看| 丰满迷人的少妇在线观看| 亚洲第一青青草原| 午夜免费鲁丝| 成人免费观看视频高清| av国产久精品久网站免费入址| 天天添夜夜摸| 我要看黄色一级片免费的| 国产精品av久久久久免费| 天天躁日日躁夜夜躁夜夜| 久久国产精品大桥未久av| 王馨瑶露胸无遮挡在线观看| 肉色欧美久久久久久久蜜桃| 曰老女人黄片| 日韩av在线免费看完整版不卡| 在线免费观看不下载黄p国产| 日韩av在线免费看完整版不卡| 一二三四在线观看免费中文在| 亚洲激情五月婷婷啪啪| 在线观看www视频免费| 一区二区三区精品91| 韩国av在线不卡| 大香蕉久久成人网| 欧美亚洲日本最大视频资源| 少妇被粗大猛烈的视频| 亚洲图色成人| 国产成人免费观看mmmm| 精品一区二区三卡| 久久久精品94久久精品| 久久av网站| 久久久久久久精品精品| 人人妻人人澡人人爽人人夜夜| 国产免费一区二区三区四区乱码| 男女之事视频高清在线观看 | 丰满饥渴人妻一区二区三| 亚洲精品,欧美精品| 亚洲成人一二三区av| 国产精品一国产av| 纯流量卡能插随身wifi吗| 看非洲黑人一级黄片| 亚洲专区中文字幕在线 | 国产精品香港三级国产av潘金莲 | 青春草国产在线视频| 欧美日韩成人在线一区二区| 18禁动态无遮挡网站| 午夜老司机福利片| 天天影视国产精品| 国产亚洲最大av| 人妻一区二区av| 日本爱情动作片www.在线观看| 五月开心婷婷网| 亚洲成人手机| 大话2 男鬼变身卡| 亚洲欧美日韩另类电影网站| 欧美日韩视频精品一区| 黑人巨大精品欧美一区二区蜜桃| 亚洲免费av在线视频| tube8黄色片| 国产精品免费大片| 亚洲精品一二三| 天美传媒精品一区二区| 午夜福利在线免费观看网站| 又黄又粗又硬又大视频| 亚洲av欧美aⅴ国产| 国产亚洲av高清不卡| 在线观看免费日韩欧美大片| 国产爽快片一区二区三区| 精品国产一区二区久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久久精品94久久精品| 日韩,欧美,国产一区二区三区| 高清不卡的av网站| 国产av精品麻豆| 五月天丁香电影| a级毛片在线看网站| 久久久精品区二区三区| 久久久精品国产亚洲av高清涩受| 一个人免费看片子| 伦理电影免费视频| 亚洲在久久综合| 美女福利国产在线| 免费女性裸体啪啪无遮挡网站| 嫩草影院入口| 成年人免费黄色播放视频| 大片电影免费在线观看免费| 另类亚洲欧美激情| 最近中文字幕高清免费大全6| 色94色欧美一区二区| 夜夜骑夜夜射夜夜干| 永久免费av网站大全| 日韩一区二区三区影片| 亚洲av成人精品一二三区| 91成人精品电影| √禁漫天堂资源中文www| 久久av网站| 中文字幕高清在线视频| 成人毛片60女人毛片免费| 麻豆乱淫一区二区| 操出白浆在线播放| 色精品久久人妻99蜜桃| 热re99久久精品国产66热6| 国产欧美日韩综合在线一区二区| 亚洲第一av免费看| 国产一区二区 视频在线| 久久久国产精品麻豆| 亚洲自偷自拍图片 自拍| 男的添女的下面高潮视频| 亚洲精品在线美女| e午夜精品久久久久久久| 欧美黑人精品巨大| 亚洲精品久久午夜乱码| 国产高清不卡午夜福利| 久久久久精品性色| 日本av免费视频播放| 国产乱人偷精品视频| 日本猛色少妇xxxxx猛交久久| 亚洲,欧美,日韩| 91精品伊人久久大香线蕉| 亚洲av日韩在线播放| 天堂俺去俺来也www色官网| 亚洲av成人精品一二三区| 精品国产一区二区三区四区第35| 亚洲情色 制服丝袜| 国产精品.久久久| 色婷婷久久久亚洲欧美| 日本av免费视频播放| 极品少妇高潮喷水抽搐| 久久久久精品国产欧美久久久 | 嫩草影视91久久| 一级片'在线观看视频| 高清视频免费观看一区二区| 精品国产国语对白av| a级毛片在线看网站| 欧美日韩亚洲高清精品| 中文精品一卡2卡3卡4更新| 久久亚洲国产成人精品v| 亚洲婷婷狠狠爱综合网| 一级,二级,三级黄色视频| 欧美亚洲日本最大视频资源| 成年女人毛片免费观看观看9 | 中文字幕另类日韩欧美亚洲嫩草| 最近的中文字幕免费完整| 久久久久精品久久久久真实原创| 日韩精品免费视频一区二区三区| 亚洲色图综合在线观看| 男女下面插进去视频免费观看| 久久久精品区二区三区| 最近最新中文字幕免费大全7| kizo精华| 女人爽到高潮嗷嗷叫在线视频| 色精品久久人妻99蜜桃| 老熟女久久久| 伦理电影免费视频| 老司机影院毛片| 国产 精品1| 久久精品久久久久久噜噜老黄| 久久狼人影院| 国产成人免费无遮挡视频| 亚洲精品aⅴ在线观看| 亚洲av综合色区一区| 亚洲成国产人片在线观看| 欧美黄色片欧美黄色片| 国产精品久久久久成人av| 黄片无遮挡物在线观看| 一级爰片在线观看| 黑人欧美特级aaaaaa片| av.在线天堂| 亚洲三区欧美一区| 亚洲一级一片aⅴ在线观看| 国产毛片在线视频| 欧美 亚洲 国产 日韩一| 国产精品一国产av| 欧美在线黄色| 高清不卡的av网站| 一级片免费观看大全| 中文字幕制服av| 亚洲综合色网址| 99国产精品免费福利视频| 不卡av一区二区三区| 免费观看a级毛片全部| 亚洲欧美色中文字幕在线| 日本av手机在线免费观看| 中文字幕人妻丝袜制服| 婷婷色av中文字幕| 在线观看三级黄色| 欧美激情极品国产一区二区三区| 中文乱码字字幕精品一区二区三区| 国产精品一国产av| 久久狼人影院| 热re99久久国产66热| 亚洲av国产av综合av卡| 在线观看人妻少妇| 精品少妇久久久久久888优播| 满18在线观看网站| 久久久久久人妻| √禁漫天堂资源中文www| 亚洲成色77777| 亚洲精品成人av观看孕妇| 亚洲免费av在线视频| 男女床上黄色一级片免费看| av不卡在线播放| 交换朋友夫妻互换小说| 亚洲三区欧美一区| 最黄视频免费看| 丰满饥渴人妻一区二区三| 男人操女人黄网站| 九九爱精品视频在线观看| 亚洲欧美成人综合另类久久久| 自线自在国产av| 99久久精品国产亚洲精品| 亚洲美女搞黄在线观看| 最近最新中文字幕免费大全7| 一区福利在线观看| 一级毛片电影观看| 亚洲成av片中文字幕在线观看| 久久 成人 亚洲| 黄频高清免费视频| 中文字幕制服av| av一本久久久久| 国产av码专区亚洲av| 黄色 视频免费看| 一边亲一边摸免费视频| 免费高清在线观看日韩| 久久久久久免费高清国产稀缺| 777米奇影视久久| 亚洲av日韩精品久久久久久密 | 亚洲欧美激情在线| 操出白浆在线播放| 美女午夜性视频免费| 亚洲欧洲精品一区二区精品久久久 | 国产亚洲精品第一综合不卡| 啦啦啦在线观看免费高清www| 国产成人一区二区在线| 嫩草影院入口| 亚洲国产欧美一区二区综合| 晚上一个人看的免费电影| 最近最新中文字幕免费大全7| 国产一卡二卡三卡精品 | 亚洲第一av免费看| 91老司机精品| 免费高清在线观看日韩| 不卡av一区二区三区| 欧美亚洲 丝袜 人妻 在线| 久久人人爽人人片av| 99热网站在线观看| 美女视频免费永久观看网站| 操出白浆在线播放| 人人妻,人人澡人人爽秒播 | 少妇被粗大猛烈的视频|