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

    Design of a 20-Gsps 12-bit time-interleaved analog-to-digital conversion system

    2021-04-01 01:31:32RuoShiDongLeiZhaoJiaJunQinWenTaoZhongYiChunFanShuBinLiuQiAn
    Nuclear Science and Techniques 2021年3期

    Ruo-Shi Dong· Lei Zhao· Jia-Jun Qin · Wen-Tao Zhong ·Yi-Chun Fan· Shu-Bin Liu · Qi An

    Abstract The time-interleaved analog-to-digital conversion (TIADC) technique is an effective method for increasing the sampling rate in a waveform digitization system. In this study, a 20-Gsps TIADC system was designed.A wide-bandwidth performance was achieved by optimizing the analog circuits, and a sufficient effective number of bits (ENOB) performance guaranteed using the perfect reconstruction algorithm for mismatch error correction.The proposed system was verified by tests,and the results indicated that a - 3 dB bandwidth of 6 GHz and the ENOB performance of 8.7 bits at 1 GHz and 7.6 bits at 6 GHz were successfully achieved.

    Keywordst Time-interleaved technique · High-speed A/D conversion · High bandwidth · Mismatch error correction

    1 Introduction

    Waveform digitization is a significant technique, especially when high-speed signals are involved in physics experiments,such as SiPM detectors[1],and many circuits for waveform digitization are developed [2–4]. In the traditional signal measurement method, the readout electronics are usually based on amplification and shaping circuits[5,6],and some other techniques are also involved,such as the charge-to-time conversion[7].Comparing with the traditional method, the input signal is directly sampled in a waveform digitization system. The original signal waveform is reconstructed by a series of sampling points,by which further details regarding charge and time can be obtained by analyzing it. Moreover, there is no baseline pile-up in a waveform digitization system, and the dead time is significantly shorter than that of the traditional method.Owing to these advantages,waveform digitization has recently become a popular research topic. In a highspeed waveform digitization system, the sampling rate is one of the most important characteristics, which directly determines how detailed the input waveform can be depicted on the time axis. The TIADC technique is an effective method for increasing the sampling rate. In a TIADC system, multiple analog-to-digital converters(ADCs) are employed. Assuming the number of ADCs is denoted as M, all having a sampling rate of Fswhile the sampling clock phase of each is shifted by a step size of 2πFs/M, an equivalent sampling rate of M × Fscan be achieved by interleaving the entire sampling data.

    Because there is more than one ADC in a TIADC system,the mismatch error introduced by different ADCs will inevitably deteriorate system performance. Several correction methods have been studied and developed to neutralize the influence of the mismatch error, and both foreground methods [8–10] and background methods[11–13]are involved.The perfect reconstruction correction algorithm is an effective foreground correction method,where the finite impulse response (FIR) filters are used to correct the ADC sampling data. In addition, an improved perfect reconstruction correction algorithm was developed,where a broadband correction is achieved using only one group of filter coefficients [10]. This method is especially suitable for error correction of the pulse waveform measurement owing to the broad range of the frequency spectrum.

    In recent years,several high-speed TIADC systems have been developed[10,14].In addition to the requirement of a high sampling rate, many other features are required for a high performance of a TIADC system; a few of the key parameters are susceptible to the analog circuits used in the TIADC design.For instance,the bandwidth performance is an important parameter that not only refers to the ADC but is also susceptible to the analog buffers, as well as the signal transmission line. Thus, the analog transmission circuit is critical for the system. Moreover, the sampling clock circuit should also be precisely designed owing to the aperture jitter directly influencing the signal-to-noise ratio(SNR) performance, especially in the high-frequency range.Therefore,the analog circuit design is significant for developing a TIADC system with high-speed performance.

    In this study, a 20-Gsps 12-bit TIADC system is designed, which mainly focuses on the analog signal transmission and clock generation circuits. Simulations were conducted on these circuits to guarantee their high performance.The sampled data were received by an FPGA and transferred to a PC for further analysis. The system performance was tested, and the results indicated that the- 3 dB bandwidth was up to 6 GHz, and an ENOB performance of 7.6 bits at 6 GHz was achieved after performing the perfect reconstruction correction.

    The remainder of this paper is organized as follows.The proposed design regarding the hardware circuit of a 20-Gsps TIADC system is introduced in Sect. 2, which is followed by the presentation of the error correction algorithm in Sect. 3.The test results are provided and discussed in Sect. 4. Finally, the conclusions are summarized in Sect. 5.

    2 20-Gsps TIADC system design

    The structure of the proposed 20-Gsps TIADC system is shown in Fig. 1a, which employs two 10-Gsps 12-bit ADCs for interleaved sampling. There are four sub-ADC banks in each ADC, and thus, the entire system can be regarded as an eight-channel 2.5-Gsps TIADC system.Because ADCs function in a dual-edge sampling mode,the sampling clock of each ADC is required to be 5 GHz for a 10-Gsps sampling,which indicates that the sampling clock period is 200 ps; both the rising edge and falling edge are used for sampling.For an equivalent 20-Gsps sampling,the time skew between the two 5-GHz sampling clocks should be 50 ps. In this case, there are four clock edges in total,which are spread in an equal time interleave of 50 ps.Hence, the phase skew between the two sampling clocks originating from the clock generating circuit should be 90°.

    The proposed system can be divided into the following three parts: (1) an analog signal circuit, responsible for transmitting a high-speed input signal to the ADC; (2) a clock generation circuit, responsible for generating a highperformance sampling clock for ADCs and other clocks for data interface and digital circuits; and (3) a digital circuit,responsible for high-speed data receiving, buffering, and uploading to a PC, and is mostly implemented through a field-programmable gate array (FPGA).

    2.1 High-bandwidth analog signal transmission circuit

    The structure of the analog signal transmission circuit is shown in Fig. 1b. The input signal is split into two channels by a power splitter,and a transformer is employed for each channel to convert a single-ended signal to a differential signal to meet the input requirement of the ADC.The analog signal transmission circuit plays a significant role in the determination of system bandwidth and signal imbalance, which is a critical source of mismatch error. Therefore,the power splitter and transformer should be carefully chosen.

    The most concerned parameters of a power splitter include the bandwidth,insertion loss,amplitude imbalance,and return loss. Generally, power splitters can be divided into the following two types: Wilkinson and resistive. The Wilkinson power splitter achieves the lowest insertion loss by splitting the input signal power into multi-channel outputs, thus providing an ideal - 3 dB output for a 1:2 Wilkinson power splitter. However, a high-performance Wilkinson power splitter, having both a wide frequency range down to DC frequencies and a low-amplitude imbalance,has a complex circuit structure and is used as a block module rather than a surface mount device, which is not suitable for integrated systems. In contrast, a resistive power splitter splits the input signal based on the signal amplitude. It usually has a simple circuit structure and can be easily integrated into a surface mount device. It simultaneously provides a broadband response down to the DC range and a low-amplitude imbalance.The main drawback of the resistive power splitter is that the ideal insertion loss is- 6 dB,which is larger than that of the Wilkinson power splitter.

    Fig. 1 a Structure of the proposed 20-Gsps 12-bit TIADC system. b Structure of an analog signal transmission circuit

    Based on these characteristics,the type of power splitter chosen should be thoroughly considered. Several parameters of a few typical surface mount power splitters are listed in Table 1, where column ‘‘Type’’ R and W indicate a resistive and Wilkinson power splitter, respectively; the corresponding S-parameters of the insertion and return losses are shown in Fig. 2a and b, respectively.

    Owing to the requirement for a wide frequency range for a pulse waveform measurement,the bandwidth range of an analog device should not only encompass the DC frequencies but also a certain part of the high-frequency band.A lower insertion loss can ensure a larger dynamic range of the measurement,which should not be too large.Moreover,the return loss can cause a harmonic distortion of the waveform; therefore, it should be minimized. The amplitude imbalance of a power splitter is the main source of the gain error between the two ADC channels,and a large gain error can decrease the efficiency of the mismatch error correction. Considering the aforementioned, the PD-0030SM was selected for use as the front-end power splitter in this study.

    A transformer is employed to transform a single-ended signal into a differential signal. Owing to the existence of an internal differential 100-Ω termination resistor between the ADC input pins, and considering the fact that the characteristic impedances of the power splitter and transmission line are equal to 50 Ω,a transformer with a resistor ratio of 1:2 is required for impedance matching.The BAL-0009SMG was chosen to be used as the coupling transformer owing to its wide bandwidth ranging from 0.5 kHz to 9 GHz. The key parameters are listed in Table 2.

    Considering the transmission requirement of high-frequency signals up to the frequency of the GHz-order, the dielectric of a printed circuit board(PCB)and the trace of a transmission line can significantly influence the signal quality. In this system, Megtron 6 is used as the PCB dielectric, which can support the transmission of signals with a frequency of up to 25 GHz. The transmission lines are designed to have the same lengths, and similar corner shapes for the two channels to meet the requirement for signal consistency.

    Table 1 Parameters of several typical power splitters

    The S-parameters of the transmission line sections are extracted using the ANSYS software.A simulation setup of the proposed analog signal transmission circuit was constructed, as illustrated in Fig. 3a. In the proposed system,there is a short transmission line from the input SMA connector to the power splitter, and the transmission lines from the power splitter to the two transformers are the same. Finally,there is a differential transmission line from each transformer to each ADC.The simulation result of the insertion loss is shown in Fig. 3b, where the dashed lines refer to the insertion loss without considering the transmission lines, and the solid lines refer to the transmission lines. At low frequencies, the insertion loss is not significantly influenced by the transmission lines. However, at high frequencies, the transmission lines introduce a nonnegligible insertion loss. According to the simulation results, the - 3 dB bandwidth of the analog signal transmission circuit is approximately 7.5 GHz.

    Fig. 2 (Color online) a Insertion losses of a few typical power splitters. b Return losses of a few typical power splitters

    Table 2 Transformer parameters

    Fig. 3 (Color online) a Simulation circuit of the analog signal transmission circuit.b Insertion loss of the analog signal transmission circuit

    2.2 High-performance clock generation circuit

    In an ADC sampling system, the sampling clock jitter can deteriorate the SNR performance. In addition, the aperture jitter caused by the time uncertainty of the ADC conversion can contribute to the deterioration of the SNR performance.The total jitter results from the square root of the quadratic sum of the sampling clock jitter and the aperture jitter. The jitter leads to an uncertainty in the sampling position and is reflected in the form of an amplitude noise of sampling points. This can cause a deterioration of the SNR, which is given by the following:

    where findenotes the input signal frequency, and tjitterdenotes the total jitter. To achieve an improved system bandwidth performance, a high-performance clock generation circuit is highly recommended.For instance,if a 7.6-bit ENOB is expected at a 6-GHz input signal, the maximum tolerable jitter is calculated to be approximately 112 fs.

    The structure of the clock generation circuit is shown in Fig. 4a. A high-performance phase-locked loop (PLL)(LMX2852) is employed as a 5-GHz sampling clock generator, and an ultra-low phase noise oscillator is used to generate the reference clock for the PLL.

    The sources of the phase noise mainly include the reference oscillator, PLL phase detector, voltage-controlled oscillator (VCO), and PLL output divider. The frequency range of the VCO in LMX2582 is 3550–7100 MHz,which indicates that a 5-GHz sampling clock can be generated from the VCO directly while the PLL output divider is bypassed.Hence,the white noise of the PLL output divider can be avoided,which mainly aggravates the phase noise at frequency points distant from the center clock frequency.The simulation results of the phase noise of different parts of the PLL are shown in Fig. 4b. The phase noise of the OSC contributes the most at the near-end offset frequencies,while PLL and VCO contribute the same at the middle offset frequency range.At far-end offset frequencies,VCO is the only source of the phase noise.

    Fig.4 (Color online)a Structure of the analog signal transmission circuit.b Simulation results of phase noise.c Phase noise simulation result of the ADC sampling clock

    The clock jitter can be calculated as follows:

    where A denotes the area of the integrated phase noise power and fodenotes the output clock frequency of the PLL.The phase noise power integration region is set to be 0.1 kHz–5 GHz, and the output clock jitter has an estimated value of 81.47 fs.

    Additional clock transmission circuits have been employed to ensure an improved performance for a larger margin of clock jitter. Band-pass filters (BPFs) are commonly used for noise filtering. The parameters of the two BPF types are listed in Table 3.The insertion losses in the stop-band of both BPFs are sufficiently high to achieve effective filtering. Furthermore, a narrow pass-band of 4900–5200 MHz can be achieved by cascading two BPFs after the PLL output, which eliminates the phase noise.

    After the sampling clock is filtered, a power splitter splits the clock for two ADC channels, and a phase shifter of each clock path shifts the clock phase to adjust the phase difference between the two clocks to 90°.The phase shifter requires a small phase step size for the phase adjustment to ensure high precision because it directly affects the skew error,which should be minimal.Finally,a transformer with a resistor ratio of 1:2 couples the sampling clock from a single-ended to a differential pair. All the devices need to have a low insertion loss at 5 GHz. Meanwhile, an RF amplifier is inserted prior to the BPFs to compensate for the insertion loss of the BPFs and other devices to guarantee a sufficiently large clock amplitude at the ADC input end.

    The phase noise simulation was conducted on the entire clock transmission circuit;the results are shown in Fig. 4c,which indicates that the phase noise of the sampling clock is effectively eliminated by the clock transmission circuit.In particular, when the offset frequency is greater than 300 MHz, the phase noise decreases sharply owing to the existence of the BPFs.Using Eq. (2),the expected jitter of the sampling clock at the ADC input end is calculated to be 60.71 fs. Note, the estimated jitter value does not consider aperture jitter,which is one of the internal characteristics of the ADC and cannot be directly obtained. The best SNR and ENOB performance achieved in the simulations without considering any other noise factors, such as aperture jitter, quantization error, and harmonic distortion, is shown in Table 4.

    Table 3 Parameters of two types of BPF

    2.3 Data interface and digital circuit

    The block diagram of the digital circuit is shown in Fig. 5a. The digital circuit is responsible for data transformation, buffering, and uploading. The line rate of the original data of a single 10-Gsps, 12-bit ADC is up to 128 Gbps after an internal 60b/64b coding in the ADC.To achieve this high-speed data transformation,the JESD204B interface,which is a high-speed serial interface protocol,is employed for the data transformation from ADC to FPGA.Two JESD204B links are set for each ADC, and there are eight lanes per link. The line rate of each JESD204B lane was calculated to be 10 Gbps after the 8b/10b coding. An FPGA with greater than 32 Gigabit transceivers is responsible for receiving and deserializing the data. A system reference clock of 250 MHz is provided by the clock generation circuit for the JESD frame alignment.The FPGA is also responsible for data caching and transferring to the PC through the SFP using the TCP protocol.

    To estimate the data transfer performance of the JESD204B lanes,the S-parameters of the PCB traces of the JESD lanes were extracted using the ANSYS software,and eye diagram simulations were conducted. The simulation result of one of the JESD lanes is shown in Fig. 5b, where the eye diagram is widely open, which indicates that the data interface can ensure effective data transfer at a line rate of 10 Gbps.

    Table 4 The best SNR and ENOB performances achieved in the clock simulation

    Fig. 5 (Color online) a Block diagram of a digital circuit. b Eye diagram simulation result of the 10-Gbps JESD204B lane

    3 Perfect reconstruction correction of 20-Gsps TIADC system

    3.1 Perfect reconstruction correction algorithm

    As previously indicated, TIADC systems suffer from mismatch errors, including the gain error, skew error, and offset error. In the proposed 20-Gsps TIADC system, a perfect reconstruction algorithm[10]was implemented in a wide range of the frequency band to mitigate the mismatch error and enhance system performance. The correction procedure of this algorithm is based on an FIR filter and can be achieved by either software or hardware modification. The block diagram of the system mismatch error and the corresponding correction process are presented in Fig. 6a. Each of the channels suffers from a certain group of mismatch error parameters; after sampling, the data of each channel are reorganized by an up-sampling of Mtimes,and an FIR filter is applied for each channel.Finally,the data array of the M channels are summed as a timeinterleaved data array after correction.

    Fig. 6 (Color online) a Mismatch error of sampling and its correction. b Correction performance dependence on the FFT point number.c Correction performance dependence on the filter order

    Considering the correction process,if the input signal in the time domain is denoted as x(t), the mismatch error can be expressed as follows:

    1. The signal sampled at the ADC input end is denoted as gmx(t), where gmdenotes the gain error parameter of the mth channel.The expression of gmin the frequency domain is denoted as~gm(ω) and is obtained by the Fourier transform, which is correlated with the frequency. The gain error is mainly caused by the signal splitting circuit inconsistency and ADC inconsistency. For the zeroth channel, the gain error can be defined as g0=1.

    3. The offset stage added to the sampling signal for the mth ADC channel is defined as the offset error, and it changes the signal expression to x(t) + Δom, where Δomusually denotes a constant independent of the signal frequency.

    In an ideal TIADC system, where gm=1, Δt0=0, and Δom=0,Fm( jω)equals e-jwmTs,which corresponds to the frequency response of an mTstime delay for the mth channel,an ideal expression( jω) can be derived as the correction target.In a non-ideal TIADC system,the offset error is usually corrected by measuring the offset of the noise waveform or by random chopping [15, 16]; therefore, this part of the error can be omitted during the derivation.Considering the other two types of mismatch errors,Fm( jω)(m=0,1...M-1) acts as an FIR filter for each ADC channel to obtain an ideal frequency domain expression( jω)from Xs( jω).Further details regarding this algorithm can be found in a previous study[10].

    During the correction process, the gain error and skew error of each channel at the frequency calibrated by the sine fitting are calculated first. Then, a series of mismatch error parameters at all frequencies are obtained by interpolation, and the numerical expression of Fmjω( ) in the Nyquist zoom is calculated. Finally, an N-order FIR parameter in the time domain for each channel is calculated by the inverse fast Fourier transform (IFFT) and windowing, and thus, in total, M × N FIR filter parameters are obtained. The obtained parameters are used to generate an FIR filter, which is exerted on the sampling array to conduct the mismatch error correction.

    3.2 Simulations

    To evaluate the effectiveness of the mismatch error correction algorithm and determine suitable values of the key algorithm parameters, simulations of the 20-Gsps TIADC were conducted. There were four sub-ADC banks in a single ADC, and the mismatch errors between them could not be ignored. Thus, it was more effective to apply the correction process to eight interleaved channels of the proposed system, rather than only two separated ADCs.

    According to the calibration results of the gain error and skew error, the simulations were conducted to study the key parameters of the algorithm. The first parameter studied was the FFT point number, which determined the intensity of the numerical expression of Fmjω( ). The correction performance for a different FFT point number is presented in Fig. 6b,where the FFT point number changed from 512 to 16,384. This shows that the correction effect was similar when the FFT point number was set within the indicated range.

    Next, the FIR filter order N was studied. Because the IFFT result is an infinite array and there must be a truncation when an FIR filter is used, the filter order is a significant parameter that determines the precision of the frequency response of the FIR filter relative to the frequency response of an ideal filter. Hence, it is directly correlated with the correction effect.The simulation results of the ENOB performance for different FIR filter orders are presented in Fig. 6c, where N is swept from 10 to 150. As presented in Fig. 6c,the correction effect was significantly sensitive to the filter order N. Particularly, the ENOB performance improved remarkably with the value of N when N was less than 100.Consequently,the simulation results indicate that the proposed correction algorithm is effective and provides an apparent improvement in the ENOB performance.

    4 Test result

    Tests were conducted to verify the system design performance. The results are summarized below.

    Fig. 7 (Color online) a The S21 parameter of the analog signal transmission circuit. b Bandwidth performance of the 20-Gsps TIADC system. c Frequency spectrum of the single ADC before correction. d Frequency spectrum of ADC1 and ADC2 before correction. e Frequency spectrum of the single ADC after correction.f Frequency spectrum of the ADC1 and ADC2 after correction.g The ENOB performance of the 20-Gsps TIADC system

    4.1 Analog signal transmission circuit

    First, the test was conducted on an analog signal transmission circuit to verify its performance. A test board,including a power splitter and transformer, was constructed. There were four single-ended outputs of the transformers. The S21parameter from the input signal to each transformer output was tested using a network analyzer.

    Theoretically, the differential output power can be calculated as the summation of both single-ended powers. As a result, the S21parameter from the input signal to each transformer’s differential output was calculated. The test results are shown in Fig. 7a, where the - 3 dB bandwidth of the analog signal transmission circuit was wider than 7.5 GHz.

    4.2 System performance

    In the next phase of the test,the dynamic performance of the 20-Gsps 12-bit TIADC system was tested according to the IEEE standard. The input signal was generated by an RF signal generator with an output frequency ranging from DC frequencies to 6 GHz.

    First,the bandwidth performance was tested.The output power of the signal generator was a constant value,and the output frequency was changed to measure the output amplitude at different frequencies. The test results are shown in Fig. 7b, where the - 3 dB bandwidth of the entire system was up to 6 GHz. The bandwidth performance curve had a similar shape to that of the simulation result on the insertion loss of the analog transmission circuit when the transmission lines were included,as shown in Fig. 3b. Hence, the insertion loss of the transmission lines can be assumed to be the primary cause of the change in gain between 5 and 6 GHz.

    Next, the ENOB performance was tested. The output frequency of the signal generator was adjusted to ensure that the signal amplitude sampled by the ADC was - 1 dBFS approximately. By exerting the FFT on the test result, the frequency spectrum was obtained, and the signal-to-noise and distortion ratio (SINAD), as well as the ENOB performance, were calculated as follows:

    Next, the correction effect was tested and verified. A four-channel correction was conducted for a single ADC,and an eight-channel correction was performed for the interleaved ADC1 and ADC2. The frequency spectrum at the 347-MHz input after the error correction is shown in Fig. 7e (single ADC) and Fig. 7f (interleaved ADC1 and ADC2),where the mismatch error peaks are observed to be significantly eliminated by the correction algorithm. The ENOB was approximately 8.7 bits at 347 MHz after the error correction.

    Finally, to verify the system performance in a wide frequency band, we also performed a frequency sweep from 50 MHz to 6 GHz.By adjusting the amplitude of the signal generator, it was ensured that the signal amplitude sampled by the ADC was approximately- 1 dBFS at each frequency in the frequency band during the test.The ENOB performance before and after correction is shown in Fig. 7g, where the dashed line indicates the ENOB performance prior to correction. During the correction process,the FIR filter order was changed between 100 and 150 for a correction performance comparison. Before the error correction, the ENOB performance was approximately 8 bits at frequencies lower than 1 GHz. However, as the signal frequency increased, the performance deteriorated immediately and was up to 5 bits at frequencies greater than 3.2 GHz. Following the error correction, the ENOB performance was 8.7 bits at 347 MHz,8 bits at 4 GHz,and 7.6 bits at 6 GHz, which sufficiently coincided with the results of a single ADC.

    Consequently, the proposed system has a wide-bandwidth performance and sufficient ENOB performance.The perfect reconstruction mismatch error correction algorithm can be considered to be effective in the proposed system.

    5 Conclusion

    In this study, a 20-Gsps 12-bit TIADC system was designed and tested. A high-quality analog signal transmission circuit and clock generating circuit were designed to enhance the system performance. A mismatch error correction method based on the perfect reconstruction correction algorithm was developed,and its efficiency was verified by tests. The test result showed that the sampling system bandwidth was up to 6 GHz, while the ENOB performance was 8.7 bits at 1 GHz, 8 bits at 4 GHz, and 7.6 bits at 6 GHz, and was significantly improved by the mismatch error correction process.

    亚洲国产看品久久| 天堂中文最新版在线下载| 欧美日韩亚洲综合一区二区三区_| 日韩大码丰满熟妇| 国产亚洲欧美精品永久| 高清欧美精品videossex| 三上悠亚av全集在线观看| 午夜福利,免费看| cao死你这个sao货| 日韩av免费高清视频| 欧美乱码精品一区二区三区| 丝瓜视频免费看黄片| 国产91精品成人一区二区三区 | 中文字幕人妻丝袜一区二区| 后天国语完整版免费观看| 久久久久精品国产欧美久久久 | 在线观看www视频免费| 国产xxxxx性猛交| 国产精品av久久久久免费| 18禁黄网站禁片午夜丰满| 久久人人爽av亚洲精品天堂| 日韩中文字幕欧美一区二区 | 精品一区二区三区av网在线观看 | 成人三级做爰电影| av不卡在线播放| av在线app专区| 亚洲成色77777| 日韩伦理黄色片| 91老司机精品| 国产高清视频在线播放一区 | 亚洲国产最新在线播放| 午夜福利视频在线观看免费| 久久精品国产亚洲av高清一级| 亚洲伊人久久精品综合| 91九色精品人成在线观看| 女人被躁到高潮嗷嗷叫费观| 午夜福利免费观看在线| 亚洲av男天堂| 99久久99久久久精品蜜桃| 18禁观看日本| 丝袜美腿诱惑在线| 久久久久国产一级毛片高清牌| 国产成人精品久久二区二区免费| 精品人妻1区二区| 男女床上黄色一级片免费看| 美女脱内裤让男人舔精品视频| 国产av国产精品国产| 制服人妻中文乱码| 美女午夜性视频免费| 午夜91福利影院| 啦啦啦中文免费视频观看日本| 精品一区二区三区av网在线观看 | 国产爽快片一区二区三区| 女性被躁到高潮视频| 亚洲精品日韩在线中文字幕| 少妇人妻 视频| 中文字幕色久视频| 一级毛片女人18水好多 | 午夜福利视频在线观看免费| 一区福利在线观看| 亚洲一码二码三码区别大吗| 久久久久视频综合| 精品福利观看| 国产不卡av网站在线观看| 亚洲av国产av综合av卡| 搡老乐熟女国产| 美国免费a级毛片| 99久久人妻综合| 精品人妻熟女毛片av久久网站| 亚洲精品久久久久久婷婷小说| 男女国产视频网站| 侵犯人妻中文字幕一二三四区| 免费少妇av软件| 亚洲欧美一区二区三区黑人| 国产亚洲欧美精品永久| 亚洲精品一卡2卡三卡4卡5卡 | 宅男免费午夜| 两性夫妻黄色片| xxxhd国产人妻xxx| 男女午夜视频在线观看| 交换朋友夫妻互换小说| 黄片播放在线免费| 老司机靠b影院| 19禁男女啪啪无遮挡网站| 欧美人与性动交α欧美精品济南到| 看免费av毛片| videos熟女内射| 黄频高清免费视频| 亚洲av电影在线观看一区二区三区| 另类精品久久| 大香蕉久久网| 婷婷丁香在线五月| 手机成人av网站| 日韩 欧美 亚洲 中文字幕| 美女大奶头黄色视频| 久久ye,这里只有精品| 国产精品一二三区在线看| 亚洲成av片中文字幕在线观看| 欧美日韩视频精品一区| 女人爽到高潮嗷嗷叫在线视频| 在线观看人妻少妇| 免费观看人在逋| 成人三级做爰电影| 国产精品久久久久久精品古装| 亚洲精品美女久久av网站| 国产国语露脸激情在线看| 搡老岳熟女国产| 亚洲一码二码三码区别大吗| 校园人妻丝袜中文字幕| videos熟女内射| 高清黄色对白视频在线免费看| 久久久精品94久久精品| 亚洲精品在线美女| 国产精品av久久久久免费| 精品久久蜜臀av无| 韩国高清视频一区二区三区| 黑丝袜美女国产一区| 日韩中文字幕视频在线看片| 一区二区三区激情视频| 最近中文字幕2019免费版| 日韩中文字幕视频在线看片| 另类精品久久| 亚洲av日韩精品久久久久久密 | 国产成人av教育| 国产成人啪精品午夜网站| 少妇 在线观看| 一本一本久久a久久精品综合妖精| 婷婷色综合大香蕉| 色婷婷久久久亚洲欧美| 免费在线观看完整版高清| 日韩伦理黄色片| 最近手机中文字幕大全| 黄片小视频在线播放| 精品熟女少妇八av免费久了| 欧美精品啪啪一区二区三区 | 精品国产一区二区三区四区第35| 婷婷色麻豆天堂久久| 午夜免费鲁丝| 免费黄频网站在线观看国产| 久久精品亚洲av国产电影网| 日韩免费高清中文字幕av| 国产精品麻豆人妻色哟哟久久| 日韩av不卡免费在线播放| 亚洲av男天堂| 无限看片的www在线观看| 国产又爽黄色视频| 免费女性裸体啪啪无遮挡网站| 亚洲精品乱久久久久久| 热re99久久精品国产66热6| 女人久久www免费人成看片| 首页视频小说图片口味搜索 | cao死你这个sao货| 久久热在线av| 亚洲欧美色中文字幕在线| 婷婷色麻豆天堂久久| 日日夜夜操网爽| 国产日韩欧美亚洲二区| 欧美亚洲 丝袜 人妻 在线| 亚洲,一卡二卡三卡| 欧美av亚洲av综合av国产av| 男人操女人黄网站| 色婷婷久久久亚洲欧美| 人妻一区二区av| 国产午夜精品一二区理论片| 亚洲精品久久久久久婷婷小说| 国产av国产精品国产| 精品亚洲成国产av| 国产色视频综合| 欧美日韩国产mv在线观看视频| 午夜福利乱码中文字幕| 在线观看免费午夜福利视频| 97人妻天天添夜夜摸| 高清视频免费观看一区二区| 亚洲国产av影院在线观看| 下体分泌物呈黄色| 90打野战视频偷拍视频| 亚洲国产精品一区三区| 日本欧美视频一区| 日韩伦理黄色片| 国产片特级美女逼逼视频| 一本一本久久a久久精品综合妖精| 日韩制服骚丝袜av| 中文乱码字字幕精品一区二区三区| 宅男免费午夜| 成年女人毛片免费观看观看9 | 中文乱码字字幕精品一区二区三区| 日韩 亚洲 欧美在线| 国产一区有黄有色的免费视频| 美女脱内裤让男人舔精品视频| 少妇的丰满在线观看| 黄色毛片三级朝国网站| 欧美日韩精品网址| 大香蕉久久成人网| 成人国产一区最新在线观看 | 我的亚洲天堂| 久久久久国产一级毛片高清牌| 好男人视频免费观看在线| 999精品在线视频| 国产精品二区激情视频| 欧美人与性动交α欧美精品济南到| 啦啦啦中文免费视频观看日本| 日韩制服丝袜自拍偷拍| av片东京热男人的天堂| 一本久久精品| 国产精品国产av在线观看| 久久精品成人免费网站| 青青草视频在线视频观看| 少妇粗大呻吟视频| 性色av乱码一区二区三区2| 久久ye,这里只有精品| 精品少妇久久久久久888优播| 九草在线视频观看| 国产片特级美女逼逼视频| 侵犯人妻中文字幕一二三四区| 别揉我奶头~嗯~啊~动态视频 | 亚洲熟女精品中文字幕| 黄色片一级片一级黄色片| 桃花免费在线播放| 两个人看的免费小视频| 午夜福利影视在线免费观看| 丝袜在线中文字幕| 国产欧美日韩精品亚洲av| 男女边摸边吃奶| 国产精品欧美亚洲77777| 亚洲自偷自拍图片 自拍| 欧美成人午夜精品| 久久久精品免费免费高清| 精品福利永久在线观看| 超色免费av| 99久久精品国产亚洲精品| 夫妻午夜视频| 美女福利国产在线| 2021少妇久久久久久久久久久| 久久久精品区二区三区| 自线自在国产av| 免费看不卡的av| 成年美女黄网站色视频大全免费| 国产精品一国产av| www.精华液| av国产久精品久网站免费入址| 在线观看www视频免费| avwww免费| 成在线人永久免费视频| 91国产中文字幕| 国产精品免费大片| 手机成人av网站| 19禁男女啪啪无遮挡网站| 啦啦啦在线观看免费高清www| 亚洲精品中文字幕在线视频| 久久久久久久久久久久大奶| 最近手机中文字幕大全| www.av在线官网国产| 中文精品一卡2卡3卡4更新| 少妇粗大呻吟视频| 最新的欧美精品一区二区| 精品少妇内射三级| 国产精品二区激情视频| 免费av中文字幕在线| www.自偷自拍.com| 香蕉国产在线看| 免费在线观看视频国产中文字幕亚洲 | 免费日韩欧美在线观看| videosex国产| 亚洲欧美日韩高清在线视频 | 91字幕亚洲| 日韩一本色道免费dvd| 人妻人人澡人人爽人人| 国产一区二区激情短视频 | 一级黄片播放器| 最近中文字幕2019免费版| 久久热在线av| 久久精品aⅴ一区二区三区四区| 99久久精品国产亚洲精品| 亚洲av综合色区一区| 国产伦人伦偷精品视频| 亚洲中文日韩欧美视频| 国产在线视频一区二区| 中国国产av一级| 欧美精品高潮呻吟av久久| 69精品国产乱码久久久| 午夜两性在线视频| 天天躁夜夜躁狠狠久久av| 王馨瑶露胸无遮挡在线观看| 亚洲国产日韩一区二区| 久久国产精品人妻蜜桃| av视频免费观看在线观看| 亚洲av男天堂| 日韩制服丝袜自拍偷拍| 亚洲熟女精品中文字幕| 妹子高潮喷水视频| 亚洲精品久久久久久婷婷小说| 女人久久www免费人成看片| 欧美另类一区| 男女无遮挡免费网站观看| 久久这里只有精品19| 久热爱精品视频在线9| 久久久欧美国产精品| 大型av网站在线播放| 亚洲色图 男人天堂 中文字幕| 脱女人内裤的视频| 免费黄频网站在线观看国产| 精品人妻熟女毛片av久久网站| 日韩一本色道免费dvd| 精品久久久久久电影网| 高清不卡的av网站| 国产爽快片一区二区三区| 岛国毛片在线播放| 午夜老司机福利片| 日韩免费高清中文字幕av| 久久久国产欧美日韩av| 精品免费久久久久久久清纯 | 黄色毛片三级朝国网站| 丝瓜视频免费看黄片| 亚洲人成77777在线视频| 人体艺术视频欧美日本| 永久免费av网站大全| 久久午夜综合久久蜜桃| 建设人人有责人人尽责人人享有的| www.999成人在线观看| 视频区欧美日本亚洲| 国产一区二区三区av在线| 黑人欧美特级aaaaaa片| 亚洲欧美激情在线| www.999成人在线观看| 亚洲av成人精品一二三区| 国产日韩欧美视频二区| 久久免费观看电影| 在现免费观看毛片| 满18在线观看网站| 欧美性长视频在线观看| 香蕉丝袜av| 91精品伊人久久大香线蕉| 亚洲国产日韩一区二区| 日日摸夜夜添夜夜爱| 人妻一区二区av| 欧美国产精品一级二级三级| 赤兔流量卡办理| 午夜激情久久久久久久| 巨乳人妻的诱惑在线观看| 日本vs欧美在线观看视频| 成人国产一区最新在线观看 | 999久久久国产精品视频| 国产一区二区三区av在线| 黄网站色视频无遮挡免费观看| 国产精品久久久久久人妻精品电影 | 国产熟女午夜一区二区三区| 国产精品久久久久久精品古装| 国产黄色免费在线视频| 亚洲精品成人av观看孕妇| 精品福利永久在线观看| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 中文乱码字字幕精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 男男h啪啪无遮挡| 在线天堂中文资源库| 久久青草综合色| 亚洲av美国av| 人人妻人人添人人爽欧美一区卜| 欧美日本中文国产一区发布| 国精品久久久久久国模美| 91精品三级在线观看| 我的亚洲天堂| 十八禁网站网址无遮挡| 国产1区2区3区精品| 十八禁网站网址无遮挡| 久久国产亚洲av麻豆专区| 大香蕉久久网| 精品高清国产在线一区| 夫妻性生交免费视频一级片| 亚洲欧美一区二区三区久久| 久久中文字幕一级| 一区二区三区乱码不卡18| 亚洲av电影在线进入| 在线观看免费高清a一片| www.av在线官网国产| 中文字幕av电影在线播放| 捣出白浆h1v1| 亚洲中文日韩欧美视频| 一级毛片黄色毛片免费观看视频| 夜夜骑夜夜射夜夜干| 最近中文字幕2019免费版| 欧美日韩国产mv在线观看视频| 看十八女毛片水多多多| 精品国产国语对白av| 午夜福利在线免费观看网站| 精品久久久精品久久久| 最近最新中文字幕大全免费视频 | 大香蕉久久成人网| 蜜桃国产av成人99| 欧美日韩福利视频一区二区| 伊人久久大香线蕉亚洲五| 亚洲熟女精品中文字幕| 日韩制服丝袜自拍偷拍| 首页视频小说图片口味搜索 | 精品一区在线观看国产| a级片在线免费高清观看视频| 国产精品一区二区在线不卡| 深夜精品福利| 国精品久久久久久国模美| 女人久久www免费人成看片| 精品高清国产在线一区| 免费观看av网站的网址| 欧美精品亚洲一区二区| 国产成人影院久久av| 国产精品国产三级国产专区5o| 国产精品国产三级专区第一集| 建设人人有责人人尽责人人享有的| 香蕉国产在线看| 中文字幕亚洲精品专区| 亚洲国产av新网站| 一个人免费看片子| 多毛熟女@视频| 1024香蕉在线观看| 亚洲,一卡二卡三卡| av有码第一页| 又黄又粗又硬又大视频| 欧美在线一区亚洲| 国产视频一区二区在线看| 只有这里有精品99| 亚洲国产欧美日韩在线播放| 亚洲伊人久久精品综合| 欧美精品高潮呻吟av久久| 免费久久久久久久精品成人欧美视频| 久久毛片免费看一区二区三区| 久久国产精品人妻蜜桃| 五月开心婷婷网| 欧美乱码精品一区二区三区| 在线看a的网站| 9热在线视频观看99| 一级,二级,三级黄色视频| 激情视频va一区二区三区| 久久九九热精品免费| 午夜精品国产一区二区电影| 日本vs欧美在线观看视频| 在线观看一区二区三区激情| √禁漫天堂资源中文www| 午夜免费男女啪啪视频观看| 亚洲 欧美一区二区三区| 少妇裸体淫交视频免费看高清 | 国产精品国产av在线观看| 69精品国产乱码久久久| 精品亚洲乱码少妇综合久久| 2018国产大陆天天弄谢| 亚洲精品一二三| 青青草视频在线视频观看| 国产精品.久久久| 无限看片的www在线观看| 性色av乱码一区二区三区2| 亚洲一区中文字幕在线| 国产精品国产三级专区第一集| 国产无遮挡羞羞视频在线观看| 99香蕉大伊视频| 国产深夜福利视频在线观看| 欧美黑人精品巨大| 国产成人欧美在线观看 | 制服人妻中文乱码| 日韩,欧美,国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 交换朋友夫妻互换小说| 亚洲精品中文字幕在线视频| 七月丁香在线播放| 精品亚洲成国产av| 亚洲欧洲日产国产| 国产精品二区激情视频| 欧美乱码精品一区二区三区| 日韩一本色道免费dvd| 国产成人欧美| 极品人妻少妇av视频| 国产在线观看jvid| 亚洲精品国产色婷婷电影| 精品少妇黑人巨大在线播放| 9191精品国产免费久久| 国产成人精品久久二区二区91| 亚洲精品成人av观看孕妇| av福利片在线| 天天躁狠狠躁夜夜躁狠狠躁| 精品一区在线观看国产| 视频区欧美日本亚洲| 色婷婷久久久亚洲欧美| 免费人妻精品一区二区三区视频| 久久人人97超碰香蕉20202| 亚洲精品乱久久久久久| 亚洲人成网站在线观看播放| 久久国产亚洲av麻豆专区| 亚洲欧美一区二区三区黑人| 午夜两性在线视频| 极品人妻少妇av视频| 久久亚洲精品不卡| 亚洲精品久久午夜乱码| 色94色欧美一区二区| 美女国产高潮福利片在线看| 日本色播在线视频| 在线观看人妻少妇| 精品一区在线观看国产| 色94色欧美一区二区| 手机成人av网站| 波多野结衣一区麻豆| 亚洲精品一区蜜桃| 午夜久久久在线观看| 日本黄色日本黄色录像| 婷婷色综合www| 日韩av在线免费看完整版不卡| 国产一级毛片在线| 一本—道久久a久久精品蜜桃钙片| 久久久精品免费免费高清| 亚洲精品av麻豆狂野| 黄频高清免费视频| 日韩中文字幕欧美一区二区 | 80岁老熟妇乱子伦牲交| 在线观看免费高清a一片| 欧美在线一区亚洲| 欧美日韩亚洲高清精品| 啦啦啦在线观看免费高清www| 欧美黑人欧美精品刺激| 中文字幕精品免费在线观看视频| 丝袜人妻中文字幕| 国产深夜福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久鲁丝午夜福利片| 女性被躁到高潮视频| 成人手机av| 免费黄频网站在线观看国产| 亚洲情色 制服丝袜| 久久人妻福利社区极品人妻图片 | 18禁黄网站禁片午夜丰满| 亚洲第一av免费看| 9色porny在线观看| 免费观看av网站的网址| 在线观看www视频免费| 欧美黄色淫秽网站| 亚洲精品一二三| 91九色精品人成在线观看| a级毛片在线看网站| 欧美+亚洲+日韩+国产| 免费久久久久久久精品成人欧美视频| 日本黄色日本黄色录像| 一边摸一边抽搐一进一出视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美色中文字幕在线| 丰满迷人的少妇在线观看| 久久中文字幕一级| 亚洲精品在线美女| 18禁裸乳无遮挡动漫免费视频| kizo精华| netflix在线观看网站| 色精品久久人妻99蜜桃| 成人国产一区最新在线观看 | av有码第一页| 狂野欧美激情性bbbbbb| 国产成人一区二区三区免费视频网站 | 亚洲九九香蕉| 精品福利观看| 日本一区二区免费在线视频| 亚洲精品久久成人aⅴ小说| 亚洲,一卡二卡三卡| 母亲3免费完整高清在线观看| 久久人人97超碰香蕉20202| 一级毛片电影观看| 成人亚洲欧美一区二区av| 侵犯人妻中文字幕一二三四区| 日本vs欧美在线观看视频| 尾随美女入室| av不卡在线播放| 美女主播在线视频| 人妻 亚洲 视频| 精品国产国语对白av| 亚洲五月婷婷丁香| 天天影视国产精品| 国产成人免费无遮挡视频| 美女大奶头黄色视频| 好男人电影高清在线观看| 国产高清videossex| 国产熟女欧美一区二区| 在线观看www视频免费| 啦啦啦在线观看免费高清www| 大香蕉久久网| 亚洲精品成人av观看孕妇| 十分钟在线观看高清视频www| 欧美黑人欧美精品刺激| 美女国产高潮福利片在线看| 久久久精品国产亚洲av高清涩受| 少妇被粗大的猛进出69影院| 国语对白做爰xxxⅹ性视频网站| 免费人妻精品一区二区三区视频| 国产又爽黄色视频| 91精品伊人久久大香线蕉| 亚洲精品久久久久久婷婷小说| 亚洲国产毛片av蜜桃av| 日日爽夜夜爽网站| 亚洲精品日韩在线中文字幕| 精品亚洲成a人片在线观看| 亚洲情色 制服丝袜| 国产成人精品无人区| 天天躁夜夜躁狠狠躁躁| 欧美另类一区| 少妇裸体淫交视频免费看高清 | 欧美精品高潮呻吟av久久| 啦啦啦在线免费观看视频4| 亚洲国产欧美一区二区综合| 精品亚洲成a人片在线观看| xxxhd国产人妻xxx| 天天操日日干夜夜撸| 日本五十路高清| 久久人人爽av亚洲精品天堂| 欧美在线一区亚洲| 午夜免费成人在线视频| 久久综合国产亚洲精品| 久久人妻熟女aⅴ| 男人舔女人的私密视频| 国产一区二区 视频在线| 国产成人精品久久久久久| 纯流量卡能插随身wifi吗| 制服人妻中文乱码| 国产女主播在线喷水免费视频网站|