• <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*
    亚洲自拍偷在线| 亚洲久久久久久中文字幕| a级毛片免费高清观看在线播放| 欧美一区二区精品小视频在线| 国产真实乱freesex| 久久99热这里只有精品18| 嫩草影视91久久| 国产精品亚洲一级av第二区| 夜夜看夜夜爽夜夜摸| 日本黄大片高清| 久久欧美精品欧美久久欧美| 色哟哟哟哟哟哟| 伊人久久精品亚洲午夜| 久久天躁狠狠躁夜夜2o2o| 亚洲18禁久久av| 久久精品夜夜夜夜夜久久蜜豆| 日本黄大片高清| 免费av不卡在线播放| av专区在线播放| 亚洲av不卡在线观看| 亚洲精品456在线播放app | 亚洲欧美日韩东京热| 久久久久久九九精品二区国产| 久久国产乱子伦精品免费另类| 国产免费av片在线观看野外av| av在线观看视频网站免费| 欧美不卡视频在线免费观看| 成人午夜高清在线视频| 亚洲av美国av| 午夜a级毛片| 国产极品精品免费视频能看的| 国产免费一级a男人的天堂| 最好的美女福利视频网| avwww免费| 亚洲美女视频黄频| 久久国产乱子免费精品| 天天一区二区日本电影三级| 日韩欧美国产在线观看| 亚洲精品456在线播放app | 女生性感内裤真人,穿戴方法视频| 深夜a级毛片| 黄色丝袜av网址大全| 日韩欧美在线二视频| 日韩欧美在线二视频| 欧美乱色亚洲激情| 午夜日韩欧美国产| 伦理电影大哥的女人| 此物有八面人人有两片| 美女免费视频网站| 美女高潮喷水抽搐中文字幕| 女人十人毛片免费观看3o分钟| 国产成人av教育| 蜜桃亚洲精品一区二区三区| 中文字幕av在线有码专区| 91狼人影院| 欧美日韩综合久久久久久 | 禁无遮挡网站| 日韩 亚洲 欧美在线| 免费搜索国产男女视频| 欧美性猛交╳xxx乱大交人| 琪琪午夜伦伦电影理论片6080| 久久精品国产亚洲av香蕉五月| 久久国产乱子伦精品免费另类| 好男人电影高清在线观看| 日本黄色视频三级网站网址| 搞女人的毛片| 在线播放无遮挡| 免费一级毛片在线播放高清视频| 老女人水多毛片| 中文字幕免费在线视频6| 欧美丝袜亚洲另类 | 日韩有码中文字幕| 毛片女人毛片| 最后的刺客免费高清国语| 欧美成人性av电影在线观看| 国产精品爽爽va在线观看网站| 国产成人福利小说| 嫩草影院入口| 男插女下体视频免费在线播放| 亚洲欧美日韩高清专用| 免费电影在线观看免费观看| 脱女人内裤的视频| 欧美国产日韩亚洲一区| 一级黄片播放器| 看片在线看免费视频| 国产免费男女视频| 一区二区三区激情视频| 国产精品不卡视频一区二区 | 欧美日本视频| 国产一区二区三区视频了| 欧美+日韩+精品| 色综合欧美亚洲国产小说| 久久九九热精品免费| 亚洲欧美日韩卡通动漫| 国产精品爽爽va在线观看网站| 国产高清激情床上av| 亚洲第一区二区三区不卡| 成人美女网站在线观看视频| 欧美日韩乱码在线| 国产伦人伦偷精品视频| 中国美女看黄片| 麻豆国产97在线/欧美| 日韩欧美一区二区三区在线观看| 黄色日韩在线| 国产熟女xx| 美女高潮的动态| 亚洲在线自拍视频| 人妻久久中文字幕网| 国产真实伦视频高清在线观看 | www.999成人在线观看| 欧美bdsm另类| 国产色爽女视频免费观看| 亚洲最大成人av| 天堂动漫精品| 国产乱人视频| 欧美在线黄色| 午夜福利18| 欧美一区二区精品小视频在线| 久久香蕉精品热| 91字幕亚洲| 久久精品久久久久久噜噜老黄 | av欧美777| 亚洲电影在线观看av| 亚洲中文日韩欧美视频| 可以在线观看毛片的网站| 日本三级黄在线观看| 亚洲在线观看片| 亚洲天堂国产精品一区在线| 看黄色毛片网站| 99国产精品一区二区三区| 久久精品国产亚洲av涩爱 | 国产在线精品亚洲第一网站| 国产免费男女视频| 国产一级毛片七仙女欲春2| 91午夜精品亚洲一区二区三区 | 又紧又爽又黄一区二区| 好看av亚洲va欧美ⅴa在| av天堂中文字幕网| 中文亚洲av片在线观看爽| 欧美最新免费一区二区三区 | 动漫黄色视频在线观看| 欧美最新免费一区二区三区 | 男女做爰动态图高潮gif福利片| 国产精品久久久久久精品电影| 宅男免费午夜| 欧美一区二区国产精品久久精品| 国产精品爽爽va在线观看网站| 岛国在线免费视频观看| 久99久视频精品免费| 免费在线观看亚洲国产| 亚洲第一欧美日韩一区二区三区| 一边摸一边抽搐一进一小说| 非洲黑人性xxxx精品又粗又长| 色视频www国产| 国产色爽女视频免费观看| 99精品在免费线老司机午夜| 又黄又爽又刺激的免费视频.| 久久伊人香网站| 日韩国内少妇激情av| netflix在线观看网站| 啦啦啦韩国在线观看视频| 亚洲无线观看免费| 国产一级毛片七仙女欲春2| 欧美三级亚洲精品| 国产一区二区三区视频了| 此物有八面人人有两片| 国产精品久久久久久精品电影| 国产成年人精品一区二区| 男女床上黄色一级片免费看| 久久婷婷人人爽人人干人人爱| 亚洲av.av天堂| 欧美日本亚洲视频在线播放| 午夜两性在线视频| 人妻制服诱惑在线中文字幕| 国内精品美女久久久久久| 日韩精品青青久久久久久| 精品免费久久久久久久清纯| 国产av不卡久久| 国产亚洲精品久久久com| 久9热在线精品视频| 国产亚洲欧美98| 久9热在线精品视频| 免费在线观看成人毛片| 99久久精品一区二区三区| 欧美乱色亚洲激情| 一进一出抽搐gif免费好疼| 亚洲国产欧洲综合997久久,| 天美传媒精品一区二区| 观看美女的网站| 搡女人真爽免费视频火全软件 | 男女下面进入的视频免费午夜| 国产成年人精品一区二区| 亚洲精品影视一区二区三区av| 欧美激情在线99| 久9热在线精品视频| 午夜精品久久久久久毛片777| 又爽又黄a免费视频| 9191精品国产免费久久| 国产精品1区2区在线观看.| 日韩人妻高清精品专区| 亚洲av免费在线观看| 午夜福利免费观看在线| 老鸭窝网址在线观看| 一级毛片久久久久久久久女| 美女xxoo啪啪120秒动态图 | 国产成人啪精品午夜网站| 看黄色毛片网站| 午夜日韩欧美国产| 少妇的逼水好多| av专区在线播放| h日本视频在线播放| 欧美+亚洲+日韩+国产| 亚洲真实伦在线观看| 亚洲av免费在线观看| 听说在线观看完整版免费高清| 在现免费观看毛片| 日本一本二区三区精品| 国产精品不卡视频一区二区 | 丰满人妻一区二区三区视频av| 久久久久久大精品| 美女cb高潮喷水在线观看| 一本精品99久久精品77| 国产伦精品一区二区三区四那| 一本综合久久免费| 一个人看的www免费观看视频| 搡女人真爽免费视频火全软件 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 色视频www国产| 国产精品久久电影中文字幕| 亚洲午夜理论影院| 亚洲成av人片在线播放无| 国产视频一区二区在线看| 国产一区二区在线观看日韩| 免费高清视频大片| 国产人妻一区二区三区在| 国产在线男女| 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 国产亚洲精品综合一区在线观看| 男人舔女人下体高潮全视频| 欧美日韩乱码在线| 亚洲午夜理论影院| 亚洲片人在线观看| 能在线免费观看的黄片| 午夜两性在线视频| 国产三级黄色录像| 成熟少妇高潮喷水视频| 婷婷精品国产亚洲av在线| 99国产综合亚洲精品| 精品无人区乱码1区二区| 日韩中文字幕欧美一区二区| 国产精品三级大全| 精品久久久久久久久久免费视频| 亚洲成人免费电影在线观看| 成人午夜高清在线视频| 久久久精品大字幕| 日韩欧美一区二区三区在线观看| 精品人妻熟女av久视频| 国产精品1区2区在线观看.| 免费av毛片视频| 精品欧美国产一区二区三| 精品福利观看| 天堂√8在线中文| 国产美女午夜福利| 一级黄色大片毛片| eeuss影院久久| 亚洲成人久久爱视频| 欧美日韩国产亚洲二区| 淫秽高清视频在线观看| 久久精品国产99精品国产亚洲性色| 亚洲无线观看免费| 色综合亚洲欧美另类图片| 成人av在线播放网站| 一个人看的www免费观看视频| 亚洲,欧美精品.| 国产真实乱freesex| 女同久久另类99精品国产91| 成年女人永久免费观看视频| 欧美一级a爱片免费观看看| 午夜福利高清视频| 欧美xxxx黑人xx丫x性爽| 精品国内亚洲2022精品成人| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 日韩欧美在线乱码| 变态另类丝袜制服| 亚洲不卡免费看| 成人三级黄色视频| 日韩av在线大香蕉| 在线播放无遮挡| 久久亚洲精品不卡| 成人鲁丝片一二三区免费| 欧美区成人在线视频| 88av欧美| av黄色大香蕉| 波多野结衣高清无吗| 毛片女人毛片| 国产亚洲精品久久久com| 天堂av国产一区二区熟女人妻| 成年女人看的毛片在线观看| 久久久久久久久大av| 18禁在线播放成人免费| 免费在线观看成人毛片| 欧美高清成人免费视频www| 国产真实伦视频高清在线观看 | 在线十欧美十亚洲十日本专区| 久久久久久久久久成人| 男人舔奶头视频| 久久精品影院6| 日本撒尿小便嘘嘘汇集6| 亚洲天堂国产精品一区在线| 特级一级黄色大片| 国内少妇人妻偷人精品xxx网站| 我要看日韩黄色一级片| 色精品久久人妻99蜜桃| 成人鲁丝片一二三区免费| 亚洲av不卡在线观看| 国产精品电影一区二区三区| 大型黄色视频在线免费观看| 久久久久国内视频| 久久久久国产精品人妻aⅴ院| 亚洲精品色激情综合| 级片在线观看| 精品久久久久久,| 51国产日韩欧美| 老司机深夜福利视频在线观看| 人人妻人人澡欧美一区二区| 51午夜福利影视在线观看| 国产成人av教育| 黄色日韩在线| 久久午夜福利片| 久久久久国产精品人妻aⅴ院| 国产精品一区二区三区四区久久| 婷婷丁香在线五月| 一区福利在线观看| 国产麻豆成人av免费视频| 亚洲中文字幕一区二区三区有码在线看| or卡值多少钱| 麻豆一二三区av精品| 在线观看一区二区三区| 亚洲美女搞黄在线观看 | 男女那种视频在线观看| 国产免费一级a男人的天堂| 18美女黄网站色大片免费观看| 人妻丰满熟妇av一区二区三区| 亚洲av.av天堂| 国产视频内射| 久久人人精品亚洲av| 嫁个100分男人电影在线观看| 少妇的逼水好多| 欧美日韩中文字幕国产精品一区二区三区| 国产在视频线在精品| 精品国产三级普通话版| 99久久精品国产亚洲精品| 又粗又爽又猛毛片免费看| 欧美黑人巨大hd| 黄色一级大片看看| 中文字幕熟女人妻在线| netflix在线观看网站| 久久精品国产清高在天天线| 午夜激情欧美在线| 欧美一区二区精品小视频在线| a级毛片a级免费在线| 人人妻人人看人人澡| 男女做爰动态图高潮gif福利片| 亚洲美女搞黄在线观看 | 九色国产91popny在线| 久久这里只有精品中国| 国产69精品久久久久777片| 日本五十路高清| 高清毛片免费观看视频网站| 欧美乱色亚洲激情| 欧美在线一区亚洲| 中文资源天堂在线| 91在线精品国自产拍蜜月| 日韩欧美国产在线观看| 最新中文字幕久久久久| 特大巨黑吊av在线直播| 日本与韩国留学比较| 两人在一起打扑克的视频| 亚洲精品在线美女| 一个人免费在线观看电影| 亚洲无线观看免费| 欧美精品国产亚洲| 看黄色毛片网站| 欧美bdsm另类| 国产成人欧美在线观看| 亚洲 欧美 日韩 在线 免费| 国产真实乱freesex| 亚洲最大成人中文| 18美女黄网站色大片免费观看| 91九色精品人成在线观看| 黄色视频,在线免费观看| 亚洲av成人精品一区久久| 十八禁网站免费在线| 欧美丝袜亚洲另类 | 少妇人妻一区二区三区视频| 首页视频小说图片口味搜索| 99久久久亚洲精品蜜臀av| 亚洲国产色片| 我的老师免费观看完整版| 精品人妻熟女av久视频| 亚洲第一区二区三区不卡| 国产欧美日韩一区二区精品| 九九在线视频观看精品| 简卡轻食公司| 一级黄色大片毛片| 久久久国产成人免费| 99国产精品一区二区三区| 午夜福利免费观看在线| 级片在线观看| 亚洲内射少妇av| 国产人妻一区二区三区在| 波多野结衣高清无吗| 日韩欧美精品免费久久 | 少妇丰满av| 亚洲av电影不卡..在线观看| 男人舔奶头视频| 制服丝袜大香蕉在线| 窝窝影院91人妻| 美女高潮喷水抽搐中文字幕| 蜜桃久久精品国产亚洲av| 精品久久久久久,| 欧美日韩亚洲国产一区二区在线观看| 亚洲乱码一区二区免费版| 免费看日本二区| 婷婷亚洲欧美| 亚洲成av人片免费观看| 国产av麻豆久久久久久久| 国产精品精品国产色婷婷| 丰满人妻熟妇乱又伦精品不卡| 亚洲av美国av| 黄色丝袜av网址大全| 亚洲18禁久久av| 色尼玛亚洲综合影院| 久久中文看片网| 亚洲av不卡在线观看| 国产欧美日韩一区二区三| 国产精品伦人一区二区| 一区福利在线观看| 特大巨黑吊av在线直播| 国产欧美日韩一区二区三| 成人亚洲精品av一区二区| 黄色女人牲交| 校园春色视频在线观看| 此物有八面人人有两片| 别揉我奶头~嗯~啊~动态视频| av在线老鸭窝| 高清日韩中文字幕在线| 真人做人爱边吃奶动态| 哪里可以看免费的av片| 好看av亚洲va欧美ⅴa在| 99在线人妻在线中文字幕| 亚洲成av人片在线播放无| 亚洲国产欧洲综合997久久,| 国产熟女xx| 精品久久久久久久久av| 久久久久久久久久黄片| 成人亚洲精品av一区二区| 成人特级黄色片久久久久久久| 一区二区三区免费毛片| 国产乱人伦免费视频| 久久精品影院6| 精品午夜福利视频在线观看一区| 久久人人精品亚洲av| www.www免费av| 天堂av国产一区二区熟女人妻| 久久人人爽人人爽人人片va | 免费在线观看成人毛片| 老司机深夜福利视频在线观看| 色播亚洲综合网| 欧美激情国产日韩精品一区| 免费搜索国产男女视频| 国产白丝娇喘喷水9色精品| 性色av乱码一区二区三区2| 级片在线观看| 精品无人区乱码1区二区| 美女被艹到高潮喷水动态| 在线免费观看的www视频| 男插女下体视频免费在线播放| 波多野结衣巨乳人妻| 亚洲自拍偷在线| 欧美xxxx性猛交bbbb| 伊人久久精品亚洲午夜| 男女床上黄色一级片免费看| 男女那种视频在线观看| 亚洲,欧美精品.| 国产亚洲精品综合一区在线观看| 淫秽高清视频在线观看| 看十八女毛片水多多多| а√天堂www在线а√下载| 一区二区三区激情视频| 老司机午夜十八禁免费视频| 人人妻人人澡欧美一区二区| 在线播放国产精品三级| 亚洲美女视频黄频| a级一级毛片免费在线观看| 久久人人爽人人爽人人片va | 一区二区三区激情视频| 热99re8久久精品国产| 日韩中文字幕欧美一区二区| 两人在一起打扑克的视频| 亚洲专区中文字幕在线| 成熟少妇高潮喷水视频| 亚洲精品日韩av片在线观看| 淫妇啪啪啪对白视频| 国产一区二区在线av高清观看| 脱女人内裤的视频| 男女之事视频高清在线观看| 又黄又爽又免费观看的视频| 青草久久国产| 三级国产精品欧美在线观看| 女人被狂操c到高潮| 91字幕亚洲| 中文字幕免费在线视频6| 亚洲精品在线观看二区| 国产久久久一区二区三区| 自拍偷自拍亚洲精品老妇| 一个人观看的视频www高清免费观看| 亚洲一区二区三区不卡视频| 亚洲七黄色美女视频| 999久久久精品免费观看国产| 51午夜福利影视在线观看| 999久久久精品免费观看国产| 极品教师在线免费播放| 国产精品久久久久久久电影| 免费在线观看成人毛片| 村上凉子中文字幕在线| 精品久久久久久久久久免费视频| 欧美午夜高清在线| 欧美日本视频| 午夜福利免费观看在线| 又紧又爽又黄一区二区| 俺也久久电影网| 91字幕亚洲| 免费黄网站久久成人精品 | 男女那种视频在线观看| 黄片小视频在线播放| 一本综合久久免费| 一本久久中文字幕| 成人av一区二区三区在线看| 99热这里只有精品一区| 3wmmmm亚洲av在线观看| 午夜福利在线观看吧| 国产免费一级a男人的天堂| 精品日产1卡2卡| 亚洲精品色激情综合| 欧美绝顶高潮抽搐喷水| 国产欧美日韩一区二区精品| 久久精品国产自在天天线| 窝窝影院91人妻| 夜夜爽天天搞| 人人妻,人人澡人人爽秒播| 直男gayav资源| 最近视频中文字幕2019在线8| 中文亚洲av片在线观看爽| xxxwww97欧美| 亚洲,欧美精品.| 日韩国内少妇激情av| 亚洲国产色片| 国产主播在线观看一区二区| 亚洲av中文字字幕乱码综合| 国模一区二区三区四区视频| 欧美日本亚洲视频在线播放| 夜夜躁狠狠躁天天躁| 男女之事视频高清在线观看| 日韩欧美 国产精品| 国产久久久一区二区三区| 中文字幕熟女人妻在线| 动漫黄色视频在线观看| 老司机福利观看| 亚洲av日韩精品久久久久久密| 麻豆成人av在线观看| 国产精品亚洲av一区麻豆| 久久人人爽人人爽人人片va | 亚洲国产欧美人成| 又紧又爽又黄一区二区| 久久久国产成人精品二区| 天堂√8在线中文| 国产精品国产高清国产av| 美女 人体艺术 gogo| 亚洲,欧美精品.| 国产av麻豆久久久久久久| 一区二区三区四区激情视频 | 亚洲欧美激情综合另类| 88av欧美| 成年版毛片免费区| 非洲黑人性xxxx精品又粗又长| 亚洲,欧美精品.| 欧美日韩瑟瑟在线播放| 国产美女午夜福利| 我的女老师完整版在线观看| 欧美+日韩+精品| 成人av在线播放网站| 亚洲综合色惰| 免费看光身美女| 亚洲人成电影免费在线| 亚洲激情在线av| 国产精品乱码一区二三区的特点| 99精品久久久久人妻精品| 免费av毛片视频| 99国产精品一区二区三区| 亚洲人与动物交配视频| 国产欧美日韩一区二区三| 亚洲,欧美,日韩| 国产成人aa在线观看| 欧美午夜高清在线| 亚洲最大成人中文| 一级黄片播放器| 亚洲av成人av| 国语自产精品视频在线第100页| 美女高潮的动态| 女生性感内裤真人,穿戴方法视频| 99视频精品全部免费 在线| 18禁黄网站禁片午夜丰满|