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

    Measurement and analysis of Doppler shift for high-speed rail scenario①

    2015-04-17 05:33:41DingJianwen丁建文
    High Technology Letters 2015年2期
    關(guān)鍵詞:建文

    Ding Jianwen (丁建文)

    (State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University, Beijing 100044, P.R.China)

    ?

    Measurement and analysis of Doppler shift for high-speed rail scenario①

    Ding Jianwen (丁建文)②

    (State Key Laboratory of Rail Traffic Control and Safety, Beijing Jiaotong University, Beijing 100044, P.R.China)

    Based on the analysis of the impact of Doppler shift to the railway dedicated mobile communication quality under the condition of high-speed mobility, this paper presents a Doppler shift measurement approach based on the phase estimation of frequency correction channel (FCCH) which correctly reflects the time-varying characteristic of the channel with a testing system developed. The Doppler shift data under the high-speed condition is collected, processed, analyzed, and compared with the simulation results. The scientific laws of the Doppler shift distribution of radio channel under high-speed condition are obtained. These data and analysis are essential for the establishment of high-speed railway Doppler power spectrum model and the development of the key technology of anti-Doppler shift.

    Doppler shift, high-speed railway (HSR), phase estimation, software radio

    0 Introduction

    A high-speed railway (HSR) propagation channel which belongs to wireless multipath fading channel has a significant impact on the design and performance analysis of wireless train control systems[1,2]. One of the main features of the HSR channel is time-varying characteristic[3]. It is required that it should be able to track the time-varying rate and remove or diminish the fading and Doppler shift of the received signal caused by the time-varying of the channel[4].

    Compared with low-speed motion, the most obvious problem of the wireless channel under the HSR condition is the much higher Doppler shift, which will lead to signal distortion, and thus, result in the deterioration of communication quality. Doppler shift can also lead to obvious time-varying of the wireless multipath fading channel, deteriorating the performance of algorithms like channel estimation and equalization, etc.[5]. Therefore, the time-varying characteristic of the wireless multipath fading channel can be obtained by studying Doppler shift of the channel, and then its impact on the communication quality under the HSR condition can be analyzed. Correspondingly, the Doppler shift correction technology can be studied.

    In Ref.[6], the tuned free-space path loss models in high speed railway considering Doppler effect and simulation results are proposed. Previous channel Doppler shift measurement campaigns investigate mainly the car-to-car (C2C) and car-to-infrastructure (C2I) cases[7-11]. In Refs[7,8], the radio channel is investigated at the frequency band of 2.4GHz. In Ref.[9], it presents a narrowband measurement campaign at 5.2GHz. Measured and simulated Doppler spectra for a frequency selective C2I channel are reported in Ref.[10]. For the C2C and C2I cases, Ref.[11] carries out a C2X radio channel measurement campaign in the 5GHz band to investigate and analyze the joint spatio-temporal Doppler spectra.

    For the application of LTE in HSR, Ref.[12] adopts the Winner II model, the Doppler shift HSR channel model and the large-scale model based on field measurements to analyze the performance of LTE-R. In Ref.[13], the approximating function is used to estimate the Doppler frequency when the distance from the high-speed terminal to the base station is relatively close. In Ref.[14], the authors propose an algorithm of joint channel estimation and Doppler frequency offset estimation based on Ricean channel model. However, these works only conduct theoretical calculating and simulation for Doppler shift.

    To the best of the authors’ knowledge, there is no publicly reported measurement data related to Doppler shift under the HSR condition. This paper presents a Doppler shift measurement method and based on that, measurement campaigns are performed at the frequency of 930 MHz of the GSM-R system in high-speed railway. Then, simulation and quantitative analyzing of dynamic characteristics of Doppler shift are conducted. Thus, the HSR Doppler shift model could be built.

    1 Measurement campaign

    1.1 Measurement method

    Doppler shift measurement method presented by this paper realizes the measuring by using phase estimation of the FCCH.

    Since the propagation delay from the transmitter to the receiver is generally unknown, carrier synchronization must be derived at the receiver. Taking GSM-R system as an example, The base station (BS) should send signals on the broadcast control channel (BCCH) to enable the mobile station (MS) to synchronize itself to BS[15]. The precondition of the synchronization is the capture of frequency correction bursts in FCCH which prepares MS for finding and demodulating the synchronous chanel (SCH) of the cell. The frequency shift of the signal can be estimated by capturing FCCH which corresponds to a pulse sequence of the frequency correction burst (FB). In each 51 multiframe of the GSM-R control channel, the FCCH burst is mapped into No. 0, 10, 20, 30 and 40 of multiframe, which transmit frequency information to MS. All the 148 bits of FB are 0 in binary, whose structure is shown in Fig.1[16]. FCCH after GMSK modulation is a pure sine wave, whose frequency is 67.5kHz higher than center frequency of the carrier.

    Fig.1 Frequency correction burst structure

    Under the ideal condition, the phase of FCCH receiving signal is

    (1)

    where Ts=1/270.833=3.7μs,and it is a symbol period. But there is a phase error caused by background noise and frequency shift, then the phase of FCCH receiving signal is expressed as

    (2)

    where Δf is frequency shift, and εn(t) is phase error caused by background noise. Make t=iTs, then the discrete phase of FCCH is

    i=0,1,…,141 (3)

    i=1,2,…,141 (4)

    According to the above two equations, the following equation can be deduced.

    i=0,1,…,141 (5)

    The phase difference between the two adjacent symbols is determined by the frequency shift and the phase error of the adjacent symbols caused by noise. Here, suppose that the phase error of two adjacent symbols caused by background noise are equal, the following equation can be obtained:

    (6)

    Under an ideal condition, the phase difference of the FCCH’s two adjacent bursts is π/2. In fact, due to the existence of Doppler shift under the condition of high-speed movement, the actually measured phase shift is Δφ(i,i-1). Therefore, the following can be obtained:

    (7)

    where Δφi,i-1can be obtained through calculation of the actually measured values output by the wireless signal collection equipment, and the calculating procedure is as follows:

    According to the I/Q data sequence of the baseband signal output by the wireless signal collection equipment, the phase sequence of FCCH receiving symbols can be calculated:

    According to the above phase sequence, the phase difference sequence could be calculated:

    Δφ(i,i-1)={Δφ0,Δφ1,…,Δφi,}, where

    Δφi,=φi,-φi-1(9)

    Putting Δφi,i-1, into Eq.(7),Doppler shift could be calculated.

    1.2 Field measurement system

    The detailed descriptions of the measurement systems are as follows.

    Transmitter: Existing GSM-R base stations (BSs) are utilized as transmitters. The carrier frequency of broadcast control channel of these BSs is in the range of 930.2MHz to 933.8MHz with a bandwidth of 200KHz. The transmitting power is 43dBm and the dual-polarization directional antennas, with 17dBi gain, 65°horizontal, and 6.8°vertical beamwidths are mounted on towers and connected to the signal transmitter with feeder cable. The power splitter with a 3dB loss is utilized to provide multichannel signals for BS to cover the cells on both sides, as shown in Fig.2(a). The transmitting power of BS during the measurement is fixed.

    Receiving antenna: The receiving antenna used by the measurement system is an omnidirectional antenna, with 4dBi gain, 80° vertical beamwidth and located in the middle part of the train, on top of the carriage at a height of 30cm above the roof, shown in Fig.2(b). The receiving antenna above the train and the receiver inside the train are connected by a 20m feeder cable with a loss of 3dB.

    Fig.2 (a) BS Transmitting antennas (b) Receiving antenna, installed on top of the carriage

    The receiver is the Doppler shift measurement system which can measure the receiving power of the BCCH carrier and gain the FCCH information through demodulation. It is composed of air interface wireless signal collection equipment, front end processor, back end processor, global positioning system (GPS) receiver, operation terminal, GSM-R antenna, GPS antenna and data cable, as shown in Fig.3. We set the BCCH frequency of BS along the railway in the measurement system, then collect wireless signal from these frequency and calculate the frequency shift.

    (a)

    (b)

    The air interface wireless signal collection equipment used here is an 8-channels receiver, which is based on the universal software radio platform (USRP), and consists of the RF processor unit, analog to digital converter(ADC), intermediate frequency(IF) sampler, digital down converter, interface unit, and clock unit, shown in Fig.4. Four sets of the front end processor are installed in the system in order to meet the requirement of processing the signals of 8-channels simultaneously. The algorithm described in 1.1 is realized in the front end processor to analyze and process the data, as well as to obtain the coverage level distribution and the Doppler shift data. The back end processor receives these data from the front end processor, and associates them with the kilometer mark, latitude and longitude information and then stores them. The operation terminal is used for data inquiry and downloading, as well as testing configuration.

    Fig.4 Air interface wireless signal collection equipment

    1.3 Measurement environment

    Extensive measurements were carried out along the section of Zaozhuang West station to Bengbu South station of Beijing-Shanghai high-speed railway in China. Viaduct, cutting and tunnel scenario are included in this section. The train used in the measurement is a high-speed comprehensive testing train, which is 3.8m high and 203m long.

    The Doppler shift measurement system is installed in the carriage of the testing train and receives the wireless signal from BS along the track through the receiving antenna deployed on the roof of the train. Under the speed of 0~350km/h, we conducted wireless channel measurement and recorded the receiving power of the BCCH carrier and the FCCH burst in real time. Then the Doppler shift is calculated by demodulating the signal in the front end processor. Over 10 round-trip tests are carried out in this test section so as to collect Doppler shift data under different velocities and ensure the consistency and repeatability of the data.

    2 Simulation analysis

    Under the condition of HSR, supposing that the distance between BSs is 3000m, the distance between the base station and the railway track is 20m, and without consideration of fading, the Doppler shift can be given by the following equation[17]:

    fs(t)=fdcosθ(t)

    (10)

    (11)

    where Dsis half of BS to BS distance, and Ds/2 is the initial distance of the train from BS, and Dminis BS-Railway track distance, both in meters; v is the velocity of the train in m/s, t is time in seconds.

    Doppler shift and cosine angle are given by Eqs(10) and (11) respectively, where the required input parameters listed in Table 1 and the simulation result of Doppler shift in different velocities are shown in Fig.5.

    Table 1 Parameters for high speed train conditions

    Fig.5 Doppler shift trajectory under different velocities

    3 Data analysis

    Based on the collected raw frequency shift data (0~350km/h), the distribution of instantaneous Doppler shift, velocity, kilometer mark and coverage level under different velocities are obtained, shown in Fig.6.

    Fig.7 and Fig.8 show the distribution trajectory of Doppler shift and velocity in the measured results of 300km/h and 350km/h, respectively. The statistical properties of the absolute value of the instantaneous Doppler shift under speed of 300km/h and 350km/h are shown in Table 2, and the statistical distributions are shown in Fig.9(a) and (b).

    Through analyzing the collected data, the Doppler shift distribution laws under the HSR condition can be drawn:

    (1) With the increase of the train speed, there is a significant increase in the Doppler shift. When the train speed is below 80km/h, the Doppler shift is about 40Hz~80Hz. When the train speed is 300km/h, the Doppler shift is up to about 250Hz and when the train speed is 350km/h, the Doppler shift is up to 280Hz~300Hz. This is consistent with the theoretical calculation and simulation result in Chapter 2, which verifies the correctness of the Doppler shift measurement and estimation method put forward by this paper.

    (a) Doppler shift

    (b) Coverage level

    Fig.7 Field instantaneous Doppler shift and velocity distribution (350km/h)

    Fig.8 Field instantaneous Doppler shift and velocity distribution (300km/h)

    (2) Since the railway mobile communication network is of linear coverage, when the train runs toward the BS direction, the Doppler shift value is positive. When the train-BS distance is much longer than that of the BS-railway track, the Doppler shift maintains the maximum. When the train reaches the BS, the incident wave of the train and the BS is a right angle, and the Doppler shift is 0. When the train passes and runs away from BS, the Doppler shift value is negative, and when the train-BS distance is much longer than that of the BS-railway track, the Doppler shift maintains the negative maximum. In fact, since the BS is very close to the railway track (about 20~40m), and this distance is much shorter than that of the train-BS in most time, the time that Doppler shift value gradually changes from the positive maximum to 0 then to the negative maximum under speed of 350km/h lasts only 2 seconds. It can be seen that this frequency shift changes very quickly, which is consistent to the simulation result.

    (3) When inter-cell handover occurs to MS in the high speed train, the serving BS changes. Within the period of handover (around 200ms), the relationship of MS and serving BS changes from back to the incident wave to facing to the incident wave, and the Doppler shift value changes instantly from the negative maximum to positive maximum. That is, the Doppler shift doubles. It can be seen that under such condition, the Doppler shift changes fastest, which is consistent to the simulation result. The rapid change of the Doppler shift results in the increase of the bit error rate (BER) caused by the handover, which results in higher requirement for the frequency shift estimation and correction of the MS receiver.

    (4) Due to the multipath effect in the field testing, each path has experienced an obvious Doppler shift process in the relative motion of MS and BS, which leads to fluctuation of the measured instantaneous value within a certain degree.

    Table 2 Statistical properties of the absolute value of the instantaneous Doppler shift (300km/h and 350km/h)

    (a) 300km/h

    (b) 350km/h

    4 Analysis of the impact of Doppler shift on different railway communication systems

    For GSM-R system (uplink frequency 885MHz, downlink frequency 934MHz), when the train speed is 350km/h, the Doppler shift is 287Hz (uplink) and 301Hz (downlink), and the corresponding coherence time is Tc=0.423/f=1.47ms (uplink)/1.41ms (downlink)[5]which is far greater than the GSM-R signal’s symbol period (3.69s). That is, the impact of the Doppler shift on the receiver is limited under the speed of 350km/h theoretically. Practically, Doppler shift will bring certain impact on the demodulation performance of BS and MS. The simulation result in Fig.10 shows that the greater the Doppler shift is, the greater the performance loss is[18]. Moreover, under the high-speed condition, the angle of the train and BS is close to zero in the most of the time, the MS-BS distance changes rapidly, so does the frequency shift. Due to the rapid change of Doppler shift, MS may be unable to correctly receive signals from neighboring cells, which is one of the reasons for the handover failure and call dropping. To solve the above problems, it is necessary to adopt and improve the automatic frequency correction technology for BS and MS in HSR.

    Fig.10 Performance loss caused by Doppler shift

    The next generation communication system will adopt the orthogonal frequency division multiplexing (OFDM) technology, which is a high-speed data transmitting technology under the wireless condition[19]. It divides the high-speed data flow into several low-speed data flows to simultaneously transmit on the mutual orthogonal sub-carrier wave, which greatly expands the pulse width of the symbol, improves the anti-multipath fading performance and the anti-delay spread ability. However, due to widening of the symbol’s pulse width, it cannot meet the requirement that Tb is far less than coherence time Tc (Tc=0.423/f) under the high-speed condition. At this moment, the Doppler shift and Doppler spread will destroy the orthogonal characteristic between sub-carriers and bring serious inter-carrier interference (ICI) which decreases the performance of system greatly. Therefore, for the next generation railway mobile communication, further study on the broadband mobile communication system carrier frequency shift estimation and its impact on the system under HSR condition will be of great importance.

    5 Conclusion

    Doppler shift measurement and estimation are an important means to grasp the time-varying characteristic of the HSR wireless multipath fading channel. This paper presents the Doppler shift measurement based on the phase estimation of the frequency correction channel, which correctly reflects the time-varying characteristic of the channel. A group of 930 MHz Doppler shift measurements under the high-speed condition have been reported and presented for the first time. The consistency (with a difference within 15 Hz) among the actually measured, theoretically calculated and simulated Doppler shift proves the feasibility and accuracy of this method. For the measurement of the broadband wireless channel Doppler shift, this method still needs further improvement to contribute to the reasonable design for the communication system under the HSR wireless multipath fading channel.

    In the meantime, based on extensive experimental results, Doppler shift characteristics have been measured, calculated and analyzed. The scientific laws of the Doppler shift distribution of radio channel under high-speed condition are obtained. Frequency shift estimation and correction are important for ensuring the high-speed handover, reducing error rate and signal distortion, which lays a foundation for the establishment of the HSR Doppler power spectrum model and the research and development of the key technology of anti-Doppler shift.

    Reference

    [ 1] He R S, Ai B, Zhong Z D, et al. Measurement-based stochastic model for high-speed railway channels. IEEE Intelligent Transportation Systems Magazine, 2014, PP (99):1-16

    [ 2] Ai B, Cheng X, Kürner T, et al. Challenges toward wireless communications for high-speed railway. IEEE Transactions on Intelligent Transportation Systems, 2014, 15 (5): 2143-58

    [ 3] Chen D H. Carrier frequency offset estimation for OFDM in highly mobile channels. Journal of Zhengzhou University of Light Industry (Natural Science Edition), 2012, 27(5): 99-102

    [ 4] Liu L Z, Feng D S. Research on measure method of doppler shift. Radio Engineering,2009, 39(5): 51-53

    [ 5] Rappaport T S. Wireless Communications Principles and Practice. Second Edition. Beijing: Publishing House of Electronics Industry, 2013. 178-179

    [ 6] Bok J, Ryu H G. Path loss model considering Doppler shift for high speed railroad communication. In: Proceedings of the 16th International Conference on Advanced Communication Technology (ICACT), Pyeongchang, South Korea, 2014. 1-4

    [ 7] Acosta G, Tokuda K, Ingram M A. Measured joint doppler-delay power profiles for vehicle-to-vehicle communications at 2.4GHz. In: Proceedings of the Global Telecommunications Conference, Dallas, USA, 2004. 3813-3817

    [ 8] Acosta G, Ingram M A. Model development for the wideband expressway vehicle-to-vehicle 2.4GHz channel. In: Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC), Las Vegas, USA, 2006. 1283-1288

    [ 9] Maurer J, Fügen T, Wiesbeck W. Narrow-band measurement and analysis of the inter-vehicle transmission channel at 5.2GHz. In: Proceedings of the Vehicular Technology Conference (VTC), Birmingham, USA, 2002. 1274-1278

    [10] Zhao X W, Kivinen J, Vainikainen P, et al. Characterization of Doppler spectra for mobile communications at 5.3GHz. IEEE Transaction on Vehicular Technology, 2003, 52(1): 14-23

    [11] Paier A, Karedal J, Czink N, et al. First results from car-to-car and car-to-infrastructure radio channel measurements at 5.2GHz. In: Proceedings of the 18th Annual IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC’07), Athens, Greece, 2007. 1-5

    [12] Guan K, Zhong Z D, Ai B. Assessment of LTE-R using high speed railway channel model. In: Proceedings of the 3rd International Conference on Communications and Mobile Computing, Qingdao, China, 2011. 461-464

    [13] Li Z M, Xiong Z G, Wu J J, et al. Doppler frequency estimation method based on approximating function under the high-speed railway environment. Journal of Convergence Information Technology, 2012, 7(20):419-428

    [14] Yang Y Q, Fan P Y, Huang Y M. Doppler frequency offsets estimation and diversity reception scheme of high speed railway with multiple antennas on separated carriages. In: Proceedings of the 2012 International Conference on Wireless Communications & Signal Processing (WCSP 2012), Huangshan, China, 2012. 1-6

    [15] 3GPP TS 45.010 Technical Specification Group GSM/EDGE Radio Access Network, Radio subsystem synchronization (V11.0.0) http://www.3gpp.org:3GPP, 2012

    [16] 3GPP TS 45.002 Technical Specification Group GSM/EDGE Radio Access Network, Multiplexing and multiple access on the radio path (V12.1.0) http://www.3gpp.org:3GPP, 2013

    [17] 3GPP TS 36.104 Technical Specification Group Radio Access Network, Evolved Universal Terrestrial Radio Access (E-UTRA), Base Station (BS) radio transmission and reception (V12.4.0) http://www.3gpp .org:3GPP, 2014

    [18] Huawei Technologies. Method and apparatus for signal processing based on automatic frequency control. China patent: CN101207437A, 2008-06-25

    [19] Chen J W. Research on Doppler Shift Estimation Technology in OFDM System: [M.S. dissertation]. Chengdu: Southwest Jiaotong University, 2009. 1-5

    Ding Jianwen, born in 1980. He received his M.S. degree from Beijing Jiaotong University in 2005.He also received his B.S. degree from North Jiaotong University in 2002.His research focuses on dedicated mobile communication system for railway, radio propagation and channel modelling, and safety communication technology for train control system.

    10.3772/j.issn.1006-6748.2015.02.011

    ①Supported by the National Natural Science Foundation of China (No. U1334202; 61222105), the Key Grant Project of Chinese Ministry of Education (No. 313006), the Project of State Key Lab of Rail Traffic Control and Safety (No. RCS 2014ZT11).

    ②To whom correspondence should be addressed. E-mail: jwding@bjtu.edu.cn Received on Sep. 11, 2014, Jiang Wenyi, Zhong Zhangdui, Guan Ke, Ai Bo

    猜你喜歡
    建文
    Magneto-hydrodynamic simulation study of direct current multi-contact circuit breaker for equalizing breaking arc
    李建文
    冼建文
    南風(fēng)(2020年8期)2020-08-06 10:25:54
    Properties of intermediate-frequency vacuum arc in sinusoidal curved contact and butt contact
    Long-Time Dynamics of Solutions for a Class of Coupling Beam Equations with Nonlinear Boundary Conditions
    愛(ài)情豈是買(mǎi)保險(xiǎn)夢(mèng)碎20萬(wàn)“愛(ài)情借條”
    山莊
    北極光(2017年1期)2017-03-18 14:42:35
    九年級(jí)化學(xué)期末測(cè)試題
    青牛之歌
    海峽美食節(jié)
    成人av一区二区三区在线看| 法律面前人人平等表现在哪些方面| 欧美日韩成人在线一区二区| 欧美日韩乱码在线| 一夜夜www| 久久精品亚洲av国产电影网| 老司机午夜福利在线观看视频| 色尼玛亚洲综合影院| 99精品在免费线老司机午夜| 成人精品一区二区免费| 午夜福利影视在线免费观看| 妹子高潮喷水视频| 在线视频色国产色| 国产精品久久久人人做人人爽| 亚洲av成人av| 91成年电影在线观看| 中国美女看黄片| 国产欧美日韩综合在线一区二区| 性少妇av在线| 国产精品免费大片| 777米奇影视久久| www.熟女人妻精品国产| 91老司机精品| 亚洲中文日韩欧美视频| 欧美不卡视频在线免费观看 | 99久久综合精品五月天人人| 亚洲精品国产色婷婷电影| 看免费av毛片| 在线十欧美十亚洲十日本专区| 午夜两性在线视频| 午夜免费成人在线视频| 国产免费av片在线观看野外av| 午夜福利免费观看在线| 免费看十八禁软件| 久久精品人人爽人人爽视色| 亚洲第一av免费看| 国产高清激情床上av| 在线永久观看黄色视频| 男女床上黄色一级片免费看| 亚洲精品一二三| 天堂√8在线中文| 俄罗斯特黄特色一大片| 日韩人妻精品一区2区三区| 午夜精品在线福利| 亚洲欧美激情综合另类| 国产亚洲欧美精品永久| 91大片在线观看| 午夜日韩欧美国产| 九色亚洲精品在线播放| 久久这里只有精品19| 美女 人体艺术 gogo| 亚洲国产中文字幕在线视频| 我的亚洲天堂| 免费在线观看日本一区| 老鸭窝网址在线观看| 国产精品久久久人人做人人爽| 久久人人97超碰香蕉20202| 搡老岳熟女国产| 免费在线观看完整版高清| 国产成人精品在线电影| 天天操日日干夜夜撸| 国产欧美日韩一区二区精品| 两人在一起打扑克的视频| 亚洲国产中文字幕在线视频| 一级毛片精品| 欧美精品啪啪一区二区三区| 日韩 欧美 亚洲 中文字幕| 欧美av亚洲av综合av国产av| av视频免费观看在线观看| 日本精品一区二区三区蜜桃| 国产一区二区三区在线臀色熟女 | 亚洲欧美色中文字幕在线| 18禁黄网站禁片午夜丰满| 中亚洲国语对白在线视频| 黄频高清免费视频| 国产欧美亚洲国产| 欧美乱码精品一区二区三区| 真人做人爱边吃奶动态| 国产激情欧美一区二区| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三区四区五区乱码| 大片电影免费在线观看免费| 久久人妻av系列| 在线观看舔阴道视频| 亚洲成a人片在线一区二区| 国产人伦9x9x在线观看| 亚洲成人手机| 国产精品免费视频内射| 老熟妇仑乱视频hdxx| 亚洲欧美色中文字幕在线| 在线天堂中文资源库| 精品福利永久在线观看| 成人黄色视频免费在线看| 婷婷丁香在线五月| 最近最新中文字幕大全电影3 | 国产精品亚洲一级av第二区| 99久久99久久久精品蜜桃| 热re99久久国产66热| 久久中文字幕人妻熟女| 亚洲精品久久成人aⅴ小说| 精品亚洲成国产av| 丝瓜视频免费看黄片| 精品一区二区三卡| 亚洲男人天堂网一区| 狠狠婷婷综合久久久久久88av| 99riav亚洲国产免费| 色婷婷av一区二区三区视频| 最近最新中文字幕大全电影3 | 9色porny在线观看| 欧美精品啪啪一区二区三区| 亚洲全国av大片| 另类亚洲欧美激情| 黄频高清免费视频| 国产av又大| 亚洲第一欧美日韩一区二区三区| 窝窝影院91人妻| 在线观看66精品国产| av线在线观看网站| 国产精品.久久久| 夜夜夜夜夜久久久久| 久久精品亚洲熟妇少妇任你| 两个人免费观看高清视频| 日韩大码丰满熟妇| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区蜜桃| 日本a在线网址| 国产亚洲欧美精品永久| 亚洲精品中文字幕在线视频| 大陆偷拍与自拍| 婷婷丁香在线五月| 亚洲熟女毛片儿| 80岁老熟妇乱子伦牲交| 国产不卡一卡二| 丰满人妻熟妇乱又伦精品不卡| 中文字幕最新亚洲高清| 国产精品乱码一区二三区的特点 | 亚洲一区高清亚洲精品| 久久影院123| 亚洲一区中文字幕在线| 国产亚洲精品第一综合不卡| 亚洲精品粉嫩美女一区| 麻豆国产av国片精品| 日本vs欧美在线观看视频| 欧美成人免费av一区二区三区 | 女人高潮潮喷娇喘18禁视频| 久久久久精品人妻al黑| 欧美精品人与动牲交sv欧美| 国产单亲对白刺激| 王馨瑶露胸无遮挡在线观看| 悠悠久久av| 欧美黄色片欧美黄色片| av超薄肉色丝袜交足视频| 久久久国产成人精品二区 | 国产精品1区2区在线观看. | 久久精品成人免费网站| 日韩欧美在线二视频 | 国产精品九九99| 在线观看免费日韩欧美大片| 伦理电影免费视频| 美国免费a级毛片| 纯流量卡能插随身wifi吗| 亚洲国产看品久久| 黄色视频不卡| 国产精品偷伦视频观看了| √禁漫天堂资源中文www| 欧美大码av| 老司机亚洲免费影院| 波多野结衣av一区二区av| 黑丝袜美女国产一区| 久久这里只有精品19| 丁香六月欧美| 12—13女人毛片做爰片一| 精品国产美女av久久久久小说| 国产成人精品无人区| 99久久精品国产亚洲精品| 国精品久久久久久国模美| 黄色 视频免费看| 热99久久久久精品小说推荐| 女性生殖器流出的白浆| 天堂动漫精品| 久久九九热精品免费| 最新在线观看一区二区三区| 亚洲av片天天在线观看| 久久精品国产a三级三级三级| 国产精品.久久久| 波多野结衣一区麻豆| 久久久久久久久免费视频了| 国产亚洲av高清不卡| 亚洲视频免费观看视频| 最近最新免费中文字幕在线| 国产成人啪精品午夜网站| 色老头精品视频在线观看| av在线播放免费不卡| 国产乱人伦免费视频| 成人av一区二区三区在线看| 岛国毛片在线播放| 亚洲精品乱久久久久久| 不卡av一区二区三区| 久久精品熟女亚洲av麻豆精品| 电影成人av| 国产精品久久电影中文字幕 | 精品人妻1区二区| 国产精品二区激情视频| 成年人黄色毛片网站| 国产亚洲欧美98| 九色亚洲精品在线播放| 午夜免费观看网址| 亚洲熟女精品中文字幕| 大码成人一级视频| 王馨瑶露胸无遮挡在线观看| 无遮挡黄片免费观看| 黄色视频,在线免费观看| e午夜精品久久久久久久| 王馨瑶露胸无遮挡在线观看| 黑人巨大精品欧美一区二区蜜桃| 成人免费观看视频高清| 国产国语露脸激情在线看| 啦啦啦在线免费观看视频4| 欧美日韩亚洲高清精品| 国产欧美日韩精品亚洲av| 在线观看日韩欧美| 日日夜夜操网爽| 女人爽到高潮嗷嗷叫在线视频| 搡老岳熟女国产| 天堂中文最新版在线下载| 自拍欧美九色日韩亚洲蝌蚪91| 757午夜福利合集在线观看| av在线播放免费不卡| 色精品久久人妻99蜜桃| 91大片在线观看| 亚洲精品中文字幕在线视频| 欧美日本中文国产一区发布| 99热国产这里只有精品6| 国产精品电影一区二区三区 | 成人黄色视频免费在线看| 国精品久久久久久国模美| 久久人人97超碰香蕉20202| 露出奶头的视频| 在线av久久热| 国产精品久久久久成人av| 精品视频人人做人人爽| 视频在线观看一区二区三区| 久久午夜综合久久蜜桃| 久久久久久免费高清国产稀缺| 久久久水蜜桃国产精品网| 国产成人一区二区三区免费视频网站| 王馨瑶露胸无遮挡在线观看| 男女午夜视频在线观看| 黄频高清免费视频| www.精华液| 中文亚洲av片在线观看爽 | 久久热在线av| 国产高清视频在线播放一区| 女人久久www免费人成看片| 三级毛片av免费| av电影中文网址| 国产日韩欧美亚洲二区| 久久人妻熟女aⅴ| 欧美另类亚洲清纯唯美| 一边摸一边做爽爽视频免费| 高清视频免费观看一区二区| 久久久国产成人免费| 亚洲成人手机| 一级a爱片免费观看的视频| 久久亚洲真实| 亚洲情色 制服丝袜| 三上悠亚av全集在线观看| 国产一区二区激情短视频| 国产成人免费观看mmmm| 欧美人与性动交α欧美软件| 久久亚洲精品不卡| 在线观看免费日韩欧美大片| 国产1区2区3区精品| 极品少妇高潮喷水抽搐| 欧美日韩亚洲综合一区二区三区_| 大片电影免费在线观看免费| 国产亚洲精品久久久久久毛片 | 啦啦啦在线免费观看视频4| 十八禁高潮呻吟视频| 最近最新中文字幕大全免费视频| 欧美精品av麻豆av| 国产精品久久久久久精品古装| 成人亚洲精品一区在线观看| e午夜精品久久久久久久| 美女国产高潮福利片在线看| 精品乱码久久久久久99久播| bbb黄色大片| 亚洲专区字幕在线| 他把我摸到了高潮在线观看| 中文字幕人妻熟女乱码| 中文字幕精品免费在线观看视频| av视频免费观看在线观看| 久久久国产成人免费| 一本综合久久免费| 日韩欧美一区视频在线观看| 亚洲色图 男人天堂 中文字幕| 黑人操中国人逼视频| 亚洲一区中文字幕在线| 丝袜在线中文字幕| 99国产精品免费福利视频| 精品久久久久久电影网| 麻豆成人av在线观看| 人妻一区二区av| 一进一出抽搐动态| 成年版毛片免费区| 午夜福利,免费看| aaaaa片日本免费| 亚洲自偷自拍图片 自拍| 两个人免费观看高清视频| 国产熟女午夜一区二区三区| 一本综合久久免费| 久久精品aⅴ一区二区三区四区| 久久性视频一级片| 精品久久久精品久久久| 国产主播在线观看一区二区| 一区在线观看完整版| 好看av亚洲va欧美ⅴa在| 日韩精品免费视频一区二区三区| 欧美乱妇无乱码| 午夜福利影视在线免费观看| 国产91精品成人一区二区三区| 欧美日韩一级在线毛片| 老司机深夜福利视频在线观看| 国产精品免费一区二区三区在线 | 脱女人内裤的视频| 人人妻,人人澡人人爽秒播| 久久精品91无色码中文字幕| 亚洲国产精品合色在线| 人成视频在线观看免费观看| 精品福利观看| 电影成人av| 亚洲午夜精品一区,二区,三区| 日韩欧美免费精品| 美国免费a级毛片| 国产精品一区二区在线观看99| 啦啦啦视频在线资源免费观看| 国产不卡av网站在线观看| 亚洲国产中文字幕在线视频| 黄色女人牲交| 久久ye,这里只有精品| 精品一品国产午夜福利视频| 午夜成年电影在线免费观看| 久久精品aⅴ一区二区三区四区| av网站在线播放免费| 国产成人系列免费观看| 国产单亲对白刺激| 亚洲国产看品久久| 免费av中文字幕在线| 一区二区日韩欧美中文字幕| 夫妻午夜视频| 另类亚洲欧美激情| 国产在线精品亚洲第一网站| 亚洲成av片中文字幕在线观看| 欧美久久黑人一区二区| 中文字幕人妻丝袜一区二区| 精品国内亚洲2022精品成人 | 欧洲精品卡2卡3卡4卡5卡区| 久久久精品免费免费高清| 大片电影免费在线观看免费| 午夜福利一区二区在线看| 免费久久久久久久精品成人欧美视频| 午夜福利一区二区在线看| 午夜成年电影在线免费观看| 亚洲视频免费观看视频| 人人澡人人妻人| 欧美激情 高清一区二区三区| av天堂久久9| 午夜免费观看网址| 亚洲在线自拍视频| 大香蕉久久网| 亚洲国产欧美网| 成在线人永久免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲九九香蕉| 91大片在线观看| 精品亚洲成a人片在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲自偷自拍图片 自拍| 中文字幕人妻丝袜一区二区| 首页视频小说图片口味搜索| 精品少妇一区二区三区视频日本电影| 国产男女超爽视频在线观看| 制服人妻中文乱码| 精品久久久久久,| 波多野结衣一区麻豆| 欧美不卡视频在线免费观看 | 日韩 欧美 亚洲 中文字幕| 啪啪无遮挡十八禁网站| 黄色视频,在线免费观看| 国产高清国产精品国产三级| 精品人妻1区二区| 18禁黄网站禁片午夜丰满| 99国产精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 亚洲精品国产一区二区精华液| 亚洲欧美日韩另类电影网站| 亚洲av日韩在线播放| 中文字幕精品免费在线观看视频| 国产午夜精品久久久久久| 在线播放国产精品三级| 亚洲伊人色综图| 麻豆国产av国片精品| 麻豆国产av国片精品| 美女 人体艺术 gogo| 国产高清视频在线播放一区| 久久午夜亚洲精品久久| 亚洲成人免费电影在线观看| 亚洲五月天丁香| 久久精品亚洲熟妇少妇任你| 国产av精品麻豆| 一区二区三区国产精品乱码| 亚洲成av片中文字幕在线观看| 国产成人欧美在线观看 | 99精品欧美一区二区三区四区| 亚洲第一欧美日韩一区二区三区| 乱人伦中国视频| 久久午夜亚洲精品久久| 又紧又爽又黄一区二区| 亚洲国产精品一区二区三区在线| netflix在线观看网站| 国产三级黄色录像| 成人18禁高潮啪啪吃奶动态图| 久久久精品免费免费高清| e午夜精品久久久久久久| 建设人人有责人人尽责人人享有的| 色94色欧美一区二区| 亚洲男人天堂网一区| 女人精品久久久久毛片| 亚洲精品在线美女| 国产精品久久视频播放| 色播在线永久视频| 久久久水蜜桃国产精品网| 黑人欧美特级aaaaaa片| 亚洲精品一二三| 久久久国产成人精品二区 | 制服人妻中文乱码| 啦啦啦免费观看视频1| 午夜老司机福利片| 久久婷婷成人综合色麻豆| a级毛片在线看网站| 一本综合久久免费| 色94色欧美一区二区| 怎么达到女性高潮| 久久中文看片网| 国产精品久久电影中文字幕 | 国产激情久久老熟女| 欧美精品亚洲一区二区| 欧美精品一区二区免费开放| 欧美+亚洲+日韩+国产| 三级毛片av免费| 亚洲国产精品sss在线观看 | 精品国内亚洲2022精品成人 | 一级毛片精品| 欧美日韩亚洲综合一区二区三区_| 色老头精品视频在线观看| 亚洲人成电影观看| 丝瓜视频免费看黄片| 亚洲精品乱久久久久久| 1024香蕉在线观看| 夜夜躁狠狠躁天天躁| 99精品欧美一区二区三区四区| 日韩欧美在线二视频 | 极品少妇高潮喷水抽搐| 91国产中文字幕| 不卡av一区二区三区| 人人妻人人澡人人爽人人夜夜| 脱女人内裤的视频| 亚洲成av片中文字幕在线观看| 免费av中文字幕在线| 一级a爱片免费观看的视频| 久久久精品国产亚洲av高清涩受| 不卡一级毛片| 91国产中文字幕| 一区二区三区激情视频| 夜夜爽天天搞| 成年版毛片免费区| а√天堂www在线а√下载 | 亚洲av日韩精品久久久久久密| 亚洲欧美日韩高清在线视频| 国产主播在线观看一区二区| 久久草成人影院| 99国产极品粉嫩在线观看| 18禁黄网站禁片午夜丰满| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区精品| 久久婷婷成人综合色麻豆| 久久九九热精品免费| 一区二区三区激情视频| 久久九九热精品免费| 亚洲精品中文字幕在线视频| 成人免费观看视频高清| 午夜久久久在线观看| 国产真人三级小视频在线观看| 人人妻人人澡人人看| 国产精品免费视频内射| 久久国产精品大桥未久av| 欧美最黄视频在线播放免费 | 老司机亚洲免费影院| 午夜福利免费观看在线| 国产激情欧美一区二区| 国产91精品成人一区二区三区| 国产精品永久免费网站| 精品视频人人做人人爽| 欧美成人免费av一区二区三区 | 国产精品一区二区在线观看99| 久久国产亚洲av麻豆专区| 91精品三级在线观看| 久久国产精品人妻蜜桃| 另类亚洲欧美激情| 国产激情久久老熟女| 亚洲国产精品sss在线观看 | 中亚洲国语对白在线视频| 在线十欧美十亚洲十日本专区| 午夜福利免费观看在线| 久久国产亚洲av麻豆专区| 中文字幕人妻丝袜一区二区| 国产人伦9x9x在线观看| 日韩 欧美 亚洲 中文字幕| 最新在线观看一区二区三区| 身体一侧抽搐| 亚洲五月婷婷丁香| 久热爱精品视频在线9| 亚洲aⅴ乱码一区二区在线播放 | 黄片小视频在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 一进一出抽搐gif免费好疼 | 国产精华一区二区三区| 激情在线观看视频在线高清 | 真人做人爱边吃奶动态| 91精品三级在线观看| 久久影院123| 丝瓜视频免费看黄片| 久久中文字幕人妻熟女| 一级毛片女人18水好多| 美女午夜性视频免费| 亚洲人成电影观看| 老司机亚洲免费影院| 欧美日韩乱码在线| 9191精品国产免费久久| 五月开心婷婷网| 国产精品电影一区二区三区 | 黄色怎么调成土黄色| 99国产精品免费福利视频| 夜夜躁狠狠躁天天躁| 99国产精品免费福利视频| 50天的宝宝边吃奶边哭怎么回事| 一a级毛片在线观看| 可以免费在线观看a视频的电影网站| 丝瓜视频免费看黄片| 黄色片一级片一级黄色片| 久久久久久人人人人人| 久99久视频精品免费| 99久久人妻综合| 久久久国产成人免费| 婷婷成人精品国产| 一区二区三区精品91| 久久国产精品影院| 亚洲五月婷婷丁香| 丝袜在线中文字幕| 久久久精品国产亚洲av高清涩受| av中文乱码字幕在线| 最新在线观看一区二区三区| 巨乳人妻的诱惑在线观看| 久久草成人影院| 午夜福利欧美成人| 色尼玛亚洲综合影院| 在线永久观看黄色视频| 亚洲在线自拍视频| 国产一区二区激情短视频| 午夜91福利影院| 黄频高清免费视频| 黄色成人免费大全| 在线观看www视频免费| 一级作爱视频免费观看| 人人妻,人人澡人人爽秒播| 狂野欧美激情性xxxx| 久久久精品免费免费高清| 9色porny在线观看| 国产亚洲欧美在线一区二区| 9191精品国产免费久久| 久久99一区二区三区| 激情视频va一区二区三区| 在线观看免费午夜福利视频| 99国产极品粉嫩在线观看| 99久久综合精品五月天人人| 久久亚洲精品不卡| 久久影院123| 女性生殖器流出的白浆| 又黄又粗又硬又大视频| 大香蕉久久网| 亚洲人成电影免费在线| 亚洲少妇的诱惑av| 国产激情久久老熟女| 日本欧美视频一区| 夜夜躁狠狠躁天天躁| 中文字幕人妻丝袜制服| 亚洲三区欧美一区| aaaaa片日本免费| 久久精品国产亚洲av香蕉五月 | 国产在线一区二区三区精| 成年女人毛片免费观看观看9 | 69精品国产乱码久久久| 亚洲国产毛片av蜜桃av| 久久久久精品国产欧美久久久| 天堂√8在线中文| 美国免费a级毛片| 18禁裸乳无遮挡免费网站照片 | 18禁国产床啪视频网站| 国产精品一区二区在线不卡| 人人妻人人澡人人爽人人夜夜| 国产乱人伦免费视频| 日本一区二区免费在线视频| 一区福利在线观看| 免费一级毛片在线播放高清视频 | 亚洲第一av免费看|