• <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é)
    免费久久久久久久精品成人欧美视频 | 97热精品久久久久久| 人妻少妇偷人精品九色| 极品教师在线视频| av国产精品久久久久影院| 亚洲av成人精品一区久久| 欧美日韩亚洲高清精品| 又大又黄又爽视频免费| 高清午夜精品一区二区三区| 久久精品国产a三级三级三级| 亚洲色图综合在线观看| 亚洲激情五月婷婷啪啪| 人妻 亚洲 视频| 高清欧美精品videossex| 久久精品久久久久久噜噜老黄| h视频一区二区三区| 在线 av 中文字幕| 99久久精品国产国产毛片| av在线播放精品| 成人影院久久| 久久综合国产亚洲精品| 国产精品爽爽va在线观看网站| 国产精品精品国产色婷婷| 99九九线精品视频在线观看视频| 我要看日韩黄色一级片| 日日摸夜夜添夜夜爱| 毛片一级片免费看久久久久| 国模一区二区三区四区视频| 日韩免费高清中文字幕av| 免费黄色在线免费观看| 午夜免费男女啪啪视频观看| 国产永久视频网站| 26uuu在线亚洲综合色| 多毛熟女@视频| 老司机影院成人| 免费黄频网站在线观看国产| 黄色配什么色好看| 欧美bdsm另类| 少妇的逼好多水| 成年人午夜在线观看视频| 久久久久久久亚洲中文字幕| 亚洲第一区二区三区不卡| 少妇人妻久久综合中文| 丰满少妇做爰视频| 新久久久久国产一级毛片| 免费av不卡在线播放| 精品国产三级普通话版| 观看美女的网站| 亚洲av免费高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 久久韩国三级中文字幕| 在线观看免费日韩欧美大片 | 国产亚洲精品久久久com| 亚洲欧美清纯卡通| 一区二区三区乱码不卡18| 国产老妇伦熟女老妇高清| 少妇人妻精品综合一区二区| 蜜桃久久精品国产亚洲av| 免费少妇av软件| videos熟女内射| 国产 一区 欧美 日韩| 99久久精品一区二区三区| 日韩,欧美,国产一区二区三区| 黑丝袜美女国产一区| 欧美xxⅹ黑人| 国产黄色免费在线视频| 久久久久精品久久久久真实原创| 久久精品久久久久久噜噜老黄| 最新中文字幕久久久久| 久久人妻熟女aⅴ| 亚洲第一区二区三区不卡| 欧美激情极品国产一区二区三区 | 伦精品一区二区三区| 日韩成人伦理影院| 国产成人精品福利久久| 国产在线视频一区二区| 丝瓜视频免费看黄片| 国产精品熟女久久久久浪| 日日啪夜夜撸| 欧美一级a爱片免费观看看| 国产日韩欧美在线精品| 欧美xxxx黑人xx丫x性爽| 在线 av 中文字幕| 免费黄网站久久成人精品| kizo精华| 国产成人精品久久久久久| 少妇人妻一区二区三区视频| www.色视频.com| 男女国产视频网站| 在线观看免费高清a一片| 久久99热6这里只有精品| 国产高清三级在线| 国产又色又爽无遮挡免| 国产免费福利视频在线观看| 精品国产三级普通话版| 嫩草影院新地址| 国产精品99久久99久久久不卡 | 亚洲真实伦在线观看| 日本wwww免费看| 国产精品欧美亚洲77777| 精品人妻偷拍中文字幕| 国产精品.久久久| 新久久久久国产一级毛片| 乱系列少妇在线播放| 久久综合国产亚洲精品| 天天躁日日操中文字幕| 能在线免费看毛片的网站| 天堂8中文在线网| 日韩av免费高清视频| 成人特级av手机在线观看| 成人漫画全彩无遮挡| 一级毛片久久久久久久久女| 女人十人毛片免费观看3o分钟| 午夜福利视频精品| 九色成人免费人妻av| 国产av精品麻豆| 欧美 日韩 精品 国产| 国产免费一区二区三区四区乱码| 七月丁香在线播放| 男女边摸边吃奶| 成人亚洲精品一区在线观看 | 高清毛片免费看| 欧美zozozo另类| 久久99热这里只有精品18| 青青草视频在线视频观看| 国产国拍精品亚洲av在线观看| 一本—道久久a久久精品蜜桃钙片| 国产男人的电影天堂91| 噜噜噜噜噜久久久久久91| 亚洲va在线va天堂va国产| 美女视频免费永久观看网站| 亚洲精品一二三| 三级经典国产精品| 欧美区成人在线视频| 高清在线视频一区二区三区| 六月丁香七月| 亚洲av.av天堂| 99精国产麻豆久久婷婷| 国产精品人妻久久久久久| 亚洲av中文字字幕乱码综合| 国产黄片美女视频| 内地一区二区视频在线| 美女内射精品一级片tv| a级一级毛片免费在线观看| 麻豆乱淫一区二区| 国产欧美日韩一区二区三区在线 | 国产伦理片在线播放av一区| 丝袜喷水一区| av专区在线播放| 日韩国内少妇激情av| 亚洲精品aⅴ在线观看| 九九久久精品国产亚洲av麻豆| 亚洲精品色激情综合| 久久97久久精品| 久久精品国产亚洲网站| 亚洲av成人精品一区久久| 九九在线视频观看精品| 26uuu在线亚洲综合色| 久久人人爽av亚洲精品天堂 | 久久女婷五月综合色啪小说| 少妇人妻一区二区三区视频| 搡老乐熟女国产| 中文乱码字字幕精品一区二区三区| 极品少妇高潮喷水抽搐| 少妇高潮的动态图| freevideosex欧美| 男女国产视频网站| 久久精品夜色国产| 亚洲丝袜综合中文字幕| 3wmmmm亚洲av在线观看| 一本一本综合久久| 99re6热这里在线精品视频| 久久99热这里只有精品18| 国产精品久久久久久精品电影小说 | 黄色日韩在线| 在线观看一区二区三区| 国产淫片久久久久久久久| 91精品一卡2卡3卡4卡| 少妇被粗大猛烈的视频| 国产在视频线精品| 18+在线观看网站| 在线观看一区二区三区| 久久久久网色| 亚洲最大成人中文| 一级片'在线观看视频| 国产中年淑女户外野战色| 97在线人人人人妻| 黑人猛操日本美女一级片| 在线观看美女被高潮喷水网站| 日本wwww免费看| 少妇被粗大猛烈的视频| 国产精品国产三级专区第一集| 在线天堂最新版资源| 色综合色国产| 秋霞伦理黄片| 国产又色又爽无遮挡免| 日日啪夜夜爽| 亚洲美女搞黄在线观看| 免费看日本二区| 麻豆精品久久久久久蜜桃| 精品亚洲成国产av| 小蜜桃在线观看免费完整版高清| 夜夜看夜夜爽夜夜摸| 免费播放大片免费观看视频在线观看| 深爱激情五月婷婷| 免费黄色在线免费观看| 人人妻人人看人人澡| 国产亚洲5aaaaa淫片| 欧美一区二区亚洲| 久久99热6这里只有精品| 国产伦理片在线播放av一区| av播播在线观看一区| 在线观看美女被高潮喷水网站| 超碰97精品在线观看| 亚洲性久久影院| 深夜a级毛片| 精品少妇久久久久久888优播| 亚洲不卡免费看| 国产高清不卡午夜福利| 99久久中文字幕三级久久日本| 18禁在线无遮挡免费观看视频| 亚洲国产精品国产精品| 亚洲无线观看免费| 中文字幕亚洲精品专区| 老师上课跳d突然被开到最大视频| 三级经典国产精品| 精华霜和精华液先用哪个| 97超视频在线观看视频| 亚洲va在线va天堂va国产| 国产久久久一区二区三区| 色综合色国产| 视频区图区小说| 一级二级三级毛片免费看| 国产片特级美女逼逼视频| 丝袜喷水一区| 久久毛片免费看一区二区三区| 国产成人一区二区在线| 汤姆久久久久久久影院中文字幕| 午夜激情久久久久久久| 王馨瑶露胸无遮挡在线观看| 极品教师在线视频| 欧美激情极品国产一区二区三区 | 中文资源天堂在线| 伊人久久国产一区二区| 国产成人一区二区在线| av国产免费在线观看| 超碰av人人做人人爽久久| 久久久久久久大尺度免费视频| 亚洲经典国产精华液单| 亚洲av男天堂| 97热精品久久久久久| 亚洲第一av免费看| 国产爱豆传媒在线观看| 97在线视频观看| 国产视频内射| 人人妻人人爽人人添夜夜欢视频 | 国产成人免费观看mmmm| 亚洲av福利一区| 高清不卡的av网站| 干丝袜人妻中文字幕| 内射极品少妇av片p| 草草在线视频免费看| 婷婷色综合www| 小蜜桃在线观看免费完整版高清| 蜜臀久久99精品久久宅男| 偷拍熟女少妇极品色| 香蕉精品网在线| 大陆偷拍与自拍| 国产精品人妻久久久影院| 乱系列少妇在线播放| 80岁老熟妇乱子伦牲交| av在线观看视频网站免费| 亚洲精品国产色婷婷电影| 一级黄片播放器| av福利片在线观看| 一区二区三区四区激情视频| 色5月婷婷丁香| 午夜激情久久久久久久| 老女人水多毛片| 少妇丰满av| 久久精品国产亚洲网站| 高清黄色对白视频在线免费看 | 婷婷色av中文字幕| 男女免费视频国产| 午夜老司机福利剧场| 内射极品少妇av片p| 五月玫瑰六月丁香| 麻豆乱淫一区二区| 亚洲不卡免费看| 肉色欧美久久久久久久蜜桃| 高清日韩中文字幕在线| 狂野欧美激情性bbbbbb| 精品国产露脸久久av麻豆| 少妇丰满av| 亚洲四区av| 免费看不卡的av| 日本av免费视频播放| 人人妻人人添人人爽欧美一区卜 | 亚洲精品第二区| av视频免费观看在线观看| 久久久亚洲精品成人影院| 天堂8中文在线网| 精品人妻一区二区三区麻豆| 男女免费视频国产| 狂野欧美激情性xxxx在线观看| 嘟嘟电影网在线观看| 夫妻性生交免费视频一级片| 国产国拍精品亚洲av在线观看| 少妇熟女欧美另类| 五月开心婷婷网| 欧美精品国产亚洲| 99热这里只有是精品在线观看| 看十八女毛片水多多多| 爱豆传媒免费全集在线观看| 国产无遮挡羞羞视频在线观看| 精品国产一区二区三区久久久樱花 | 自拍欧美九色日韩亚洲蝌蚪91 | .国产精品久久| 在线观看美女被高潮喷水网站| 国产午夜精品一二区理论片| 夫妻午夜视频| 久久婷婷青草| 一本—道久久a久久精品蜜桃钙片| 黄色视频在线播放观看不卡| 久久av网站| 精品国产乱码久久久久久小说| 综合色丁香网| 国产日韩欧美亚洲二区| 美女脱内裤让男人舔精品视频| 美女国产视频在线观看| 婷婷色综合www| 在线观看免费视频网站a站| 日日啪夜夜撸| 一级毛片我不卡| 欧美老熟妇乱子伦牲交| 日本wwww免费看| 狂野欧美激情性bbbbbb| 久久女婷五月综合色啪小说| 97在线视频观看| 熟妇人妻不卡中文字幕| 国产免费一级a男人的天堂| 观看美女的网站| 99热这里只有是精品50| 在线观看免费日韩欧美大片 | 一区二区三区精品91| 国产人妻一区二区三区在| 日本-黄色视频高清免费观看| av国产久精品久网站免费入址| 久久久久国产网址| 高清不卡的av网站| 性高湖久久久久久久久免费观看| 欧美成人a在线观看| 亚洲中文av在线| 啦啦啦中文免费视频观看日本| 内射极品少妇av片p| 在线观看美女被高潮喷水网站| 国产精品人妻久久久久久| av一本久久久久| 赤兔流量卡办理| 国产精品免费大片| 80岁老熟妇乱子伦牲交| 一级毛片电影观看| 国产精品不卡视频一区二区| 超碰av人人做人人爽久久| 尾随美女入室| 亚洲精品久久午夜乱码| 少妇高潮的动态图| 精品久久久久久久末码| 高清午夜精品一区二区三区| 在现免费观看毛片| 天天躁夜夜躁狠狠久久av| 日韩 亚洲 欧美在线| 中文字幕av成人在线电影| 国产美女午夜福利| 亚洲精品亚洲一区二区| 色视频www国产| 91久久精品电影网| 国产中年淑女户外野战色| 亚洲一区二区三区欧美精品| 韩国高清视频一区二区三区| 亚洲成人中文字幕在线播放| 免费黄网站久久成人精品| 搡女人真爽免费视频火全软件| 久久精品国产a三级三级三级| 亚洲伊人久久精品综合| 观看免费一级毛片| 亚洲内射少妇av| 99久久精品一区二区三区| 欧美另类一区| 欧美+日韩+精品| 日本爱情动作片www.在线观看| 91狼人影院| 视频中文字幕在线观看| 国产精品.久久久| 在线观看美女被高潮喷水网站| 最近最新中文字幕大全电影3| 国产精品一区二区性色av| 欧美性感艳星| 国产色爽女视频免费观看| 久久国产乱子免费精品| 国产v大片淫在线免费观看| 一区在线观看完整版| 在线亚洲精品国产二区图片欧美 | 国产精品偷伦视频观看了| 干丝袜人妻中文字幕| 美女主播在线视频| 国产深夜福利视频在线观看| 国产精品人妻久久久影院| 大又大粗又爽又黄少妇毛片口| 男女边吃奶边做爰视频| 日日摸夜夜添夜夜添av毛片| 免费观看性生交大片5| 久久久久国产网址| 噜噜噜噜噜久久久久久91| 又黄又爽又刺激的免费视频.| 国内少妇人妻偷人精品xxx网站| 国产精品女同一区二区软件| 免费黄频网站在线观看国产| 2022亚洲国产成人精品| 色视频在线一区二区三区| 免费看不卡的av| 精品一区在线观看国产| 国产毛片在线视频| 国精品久久久久久国模美| 亚洲最大成人中文| 18禁裸乳无遮挡免费网站照片| av.在线天堂| 国产精品嫩草影院av在线观看| 边亲边吃奶的免费视频| 久久久久国产精品人妻一区二区| 亚洲欧美精品专区久久| 久久精品国产亚洲av天美| 搡女人真爽免费视频火全软件| 乱码一卡2卡4卡精品| 久久久国产一区二区| 日韩在线高清观看一区二区三区| 国产av国产精品国产| 久久女婷五月综合色啪小说| 中文乱码字字幕精品一区二区三区| 亚洲精品视频女| 性色av一级| 国内揄拍国产精品人妻在线| 国产日韩欧美在线精品| 久久久久人妻精品一区果冻| 女性生殖器流出的白浆| 国产精品嫩草影院av在线观看| 久久精品人妻少妇| 欧美日韩亚洲高清精品| 一本—道久久a久久精品蜜桃钙片| 一级毛片aaaaaa免费看小| 青春草亚洲视频在线观看| 久久久久国产网址| 99久久精品国产国产毛片| 国产视频内射| 免费观看a级毛片全部| 日韩中文字幕视频在线看片 | 美女xxoo啪啪120秒动态图| 欧美激情极品国产一区二区三区 | 日日啪夜夜撸| 熟妇人妻不卡中文字幕| 日韩三级伦理在线观看| 久久毛片免费看一区二区三区| 久久99热这里只频精品6学生| 亚洲av二区三区四区| 老司机影院成人| 亚洲精品,欧美精品| 九九爱精品视频在线观看| 久久久久视频综合| 蜜桃在线观看..| 激情 狠狠 欧美| 国产在线免费精品| 99久国产av精品国产电影| 日韩人妻高清精品专区| 午夜精品国产一区二区电影| 五月玫瑰六月丁香| 91狼人影院| 王馨瑶露胸无遮挡在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲精品中文字幕在线视频 | 日韩强制内射视频| 国产白丝娇喘喷水9色精品| 日韩亚洲欧美综合| 国产精品成人在线| 又爽又黄a免费视频| 国产亚洲午夜精品一区二区久久| 99国产精品免费福利视频| 亚洲国产最新在线播放| 亚洲av欧美aⅴ国产| 美女视频免费永久观看网站| 成人国产麻豆网| 国产在线一区二区三区精| 亚洲av男天堂| 亚洲人成网站高清观看| 国产一区亚洲一区在线观看| 国产色爽女视频免费观看| 亚洲伊人久久精品综合| 久久精品国产亚洲av天美| 日韩伦理黄色片| 欧美3d第一页| 黄片无遮挡物在线观看| 亚洲精品成人av观看孕妇| 人妻夜夜爽99麻豆av| 亚洲欧美清纯卡通| 黑丝袜美女国产一区| 毛片一级片免费看久久久久| 大香蕉久久网| 欧美zozozo另类| 麻豆成人av视频| 欧美性感艳星| 中国美白少妇内射xxxbb| 97在线视频观看| 亚洲丝袜综合中文字幕| 大片电影免费在线观看免费| 国产成人a区在线观看| 一本—道久久a久久精品蜜桃钙片| 高清在线视频一区二区三区| 五月伊人婷婷丁香| 国产精品蜜桃在线观看| 久久精品国产自在天天线| 国产一级毛片在线| 欧美日韩在线观看h| 色视频www国产| 一区在线观看完整版| 又黄又爽又刺激的免费视频.| 美女xxoo啪啪120秒动态图| 美女脱内裤让男人舔精品视频| av不卡在线播放| 亚洲av欧美aⅴ国产| 少妇人妻精品综合一区二区| 国产极品天堂在线| 国产一区二区在线观看日韩| 妹子高潮喷水视频| 亚洲美女黄色视频免费看| 嘟嘟电影网在线观看| 欧美bdsm另类| 老熟女久久久| 日韩三级伦理在线观看| 美女中出高潮动态图| 久久久精品94久久精品| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| 欧美人与善性xxx| 综合色丁香网| 黄色一级大片看看| 美女福利国产在线 | 亚洲激情五月婷婷啪啪| 免费黄色在线免费观看| 亚洲精品,欧美精品| 国产69精品久久久久777片| 日本av免费视频播放| 蜜桃亚洲精品一区二区三区| 国产真实伦视频高清在线观看| 啦啦啦啦在线视频资源| 只有这里有精品99| 在线亚洲精品国产二区图片欧美 | 亚洲欧美日韩东京热| 午夜福利视频精品| 日本免费在线观看一区| 2021少妇久久久久久久久久久| 亚洲av日韩在线播放| 少妇 在线观看| 亚洲美女黄色视频免费看| 国产v大片淫在线免费观看| 少妇人妻 视频| 中文字幕人妻熟人妻熟丝袜美| 十分钟在线观看高清视频www | 久久人人爽人人片av| 在线观看免费日韩欧美大片 | 国国产精品蜜臀av免费| 中文字幕久久专区| 亚洲精品国产色婷婷电影| 国产伦理片在线播放av一区| 啦啦啦中文免费视频观看日本| 黄色一级大片看看| 成人无遮挡网站| 免费观看在线日韩| 狂野欧美白嫩少妇大欣赏| freevideosex欧美| 成人二区视频| 五月伊人婷婷丁香| 99久久人妻综合| xxx大片免费视频| 久久这里有精品视频免费| 成人影院久久| 人人妻人人爽人人添夜夜欢视频 | 婷婷色综合大香蕉| 日日啪夜夜爽| 一级片'在线观看视频| 97超碰精品成人国产| 国产综合精华液| 成人亚洲精品一区在线观看 | 一级爰片在线观看| 精品久久久久久久久av| 日韩成人伦理影院| 欧美成人a在线观看| 日本-黄色视频高清免费观看| 欧美精品亚洲一区二区| 亚洲精品亚洲一区二区| 久久精品夜色国产| av免费观看日本| 国产 一区 欧美 日韩| 欧美国产精品一级二级三级 | 夫妻性生交免费视频一级片| 美女cb高潮喷水在线观看| 欧美精品人与动牲交sv欧美| 久久精品久久久久久久性| 亚洲av成人精品一二三区| 国产高清三级在线| 搡老乐熟女国产| 国产在线视频一区二区| 欧美成人一区二区免费高清观看| 欧美日韩一区二区视频在线观看视频在线| 精品人妻熟女av久视频| 国产精品成人在线| 国产爽快片一区二区三区|