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

    基于多相濾波的多通道數字終端設計*

    2015-03-22 11:31:31金乘進岳友嶺吳俊霖
    天文研究與技術 2015年4期
    關鍵詞:傅里葉頻域時域

    趙 欣,金乘進,朱 巖,岳友嶺,吳俊霖

    (1. 中國科學院國家天文臺,北京 100012;2. 中國科學院大學,北京 100049)

    基于多相濾波的多通道數字終端設計*

    趙 欣1,2,金乘進1,朱 巖1,岳友嶺1,吳俊霖1,2

    (1. 中國科學院國家天文臺,北京 100012;2. 中國科學院大學,北京 100049)

    數字終端技術在射電天文數據處理方面已成為主流,在國內外都有廣泛的應用?;仡櫫藶V波的基本原理,多相濾波技術在多通道濾波中的應用,探討了對通道帶通性能的優(yōu)化方法。使用美國UC Berkeley CASPER團組研發(fā)的終端開發(fā)平臺實現了一個觀測帶寬為400 MHz,通道數16 k,頻譜分辨率約為25 kHz的數字終端,并進行了掃頻測試。在設計中,模擬信號首先通過模擬數字轉換器轉換成數字信號,利用多相濾波器組算法實現多通道輸出,并對每個通道內的數據進行取平方運算并累加得到功率譜。累加后的數據可通過10 GbE傳輸到其他設備存儲或繼續(xù)處理。該終端可在FAST河外中性氫頻譜巡天觀測中使用。

    多相濾波;多通道濾波;10 GbE;數字終端

    早期的射電天文數據處理采用模擬數據終端,結構簡單且造價相對便宜,如1951年Ewen和Purcell首次發(fā)現來自銀河系中心的中性氫譜線所用的模擬終端[1];烏魯木齊天文站1998年研制的用于脈沖星觀測的模擬多通道消色散終端等。與模擬終端相比,數字終端具有處理方式靈活、性能穩(wěn)定、可移植性好、易于升級的特點,且更容易存儲大量數據。隨著電子以及計算機技術的發(fā)展,數字終端已經成為射電望遠鏡終端的主流技術。如2007年安裝在青海德令哈觀測基地的由紫金山天文臺研制的 “數字FFT頻譜及配套中頻系統(tǒng)”;云南天文臺與北京大學合作研制的0.7~1.4 GHz分米波自相關型寬帶頻譜儀[2];上海天文臺研制的數字濾波和數字基帶轉換器(Digital Base-band Converter, DBBC)[3];美國國立射電天文臺GreenBank射電望遠鏡的VEGAS多用途數字終端等。

    本文對濾波的基本原理和濾波器性能優(yōu)化的方法做了回顧,并對多相濾波技術在多通道頻譜分析中的應用做了說明。同時對計算量進行了分析,由于多相濾波進行了降采樣操作,故節(jié)省了計算量,在射電天文頻譜觀測中得到了廣泛應用。本文采用多相濾波技術在CASPER研制的終端開發(fā)平臺上設計和實現了一個觀測帶寬為400 MHz,通道數16 k的數字頻譜終端,該終端可應用于FAST河外中性氫巡天觀測*Nan Rendong, et al. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) project[J]. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) project, 2011.。

    1 濾波和多通道頻譜分析基本原理

    1.1 濾波基本原理回顧

    濾波是在頻域上對數據進行選擇,帶通濾波器形狀為矩形的濾波函數使得通帶內的信號不失真地通過濾波器,通帶外的信號則被濾掉。當濾波器作用于信號時,在頻域上實現的操作為信號的頻域響應和帶通濾波器的 “矩形” 頻域響應相乘,對應在時域上的操作為信號和帶通濾波器時域函數的卷積。在實際操作中,往往對 “矩形” 濾波器的時域sinc函數做一定長度的截取(見圖1(1)至圖1(3)),并用這個截取后的sinc函數和輸入信號卷積。此時由于對頻域 “矩形” 濾波器的時域函數進行了截取,頻域的濾波器性能發(fā)生了變化。根據卷積定理,時域上無限長的sinc函數和截斷函數(這里為矩形窗)的相乘對應到頻域的操作為二者對應的頻域響應的卷積,這里即為 “矩形” 濾波器和一個sinc函數(截斷函數矩形窗對應的頻率響應)的卷積。與 “矩形” 濾波器相比,此時得到的濾波函數的頻域響應不再為矩形,出現過渡帶和震蕩,存在頻譜泄露現象*The Polyphase Filter Bank Technique[OL]. CASPER Memo 42, 2011.(見圖1(4))。

    為了改善該濾波函數的頻域帶通性能,可以加長其時域函數截斷長度(即矩形窗長度)。此時截斷函數(矩形窗)的頻率響應sinc函數的主瓣變窄,使得濾波函數通帶形狀更接近矩形,但吉布斯現象增強;可以選擇其他合適的窗函數(如漢明窗)改善,得到較好的濾波性能(見圖1(5)至1(8))。漢明窗和矩形窗的頻率響應對比如圖2,頻域采用歸一化頻率,可以看出漢明窗的頻域響應旁瓣抑制(約-50 dB)較矩形窗的旁瓣抑制(約-15 dB)要好,但其主瓣寬度較寬。

    圖1 (1)為一個矩形窗函數,起到對序列進行截斷的作用;(2)為一個無限長的sinc函數,其頻域響應為一個頻域矩形函數;(3)對圖1(2)的無限長的sinc函數采用圖1(1)的矩形窗相乘進行截斷操作后得到的序列;(4)為圖1(3)截斷序列所對應的頻域響應,有較大震蕩和較寬過渡帶,存在頻譜泄露現象。(5)為加長的窗函數(此處為漢明窗);(6)為無限長的sinc函數;(7)為加長的窗函數和無限長sinc函數相乘得到的序列;(8)為圖1(7)對應序列的頻域響應,經過加長序列和采用窗函數后,震蕩和過渡帶減小,頻域響應性能得到改善

    圖2 矩形窗和漢明窗的頻率響應對比

    1.2 多通道輸出和多相濾波基本原理

    射電天文中頻譜觀測往往需要在一定的頻率范圍內做多個通道,并獲取每個通道的功率。根據帕斯瓦爾定律,要獲取每個通道的功率,不需要得到每個通道內的時域函數,而是通過得到通道內的頻域函數的均值后,直接進行平方運算得到功率[4]。取通道帶寬對應的 “矩形” 濾波器,進行傅里葉逆變換后得到時域的sinc函數,輸入信號與該函數相乘后進行傅里葉變換得到頻域 “矩形” 濾波器和輸入信號頻譜的卷積,得到的卷積值即為 “矩形” 濾波器通帶內的輸入信號頻譜的平均值。將得到的頻域上的卷積值以通帶帶寬為間隔取樣,取平方后,得到多通道輸出的功率譜。

    對于滿足奈奎斯特采樣定理的采樣頻率為fs的離散序列,取N個時域數據進行快速傅里葉轉換后可以得到N個頻域數據。N個頻域數據即為在原頻譜與一個主瓣寬度為2fs/N的sinc函數卷積上的間距為fs/N的采樣,通過取平方運算后可以得到包括各個頻點的帶通為主瓣寬度2fs/Nsinc函數內的功率。為了改善通道帶通性能形狀,可以在時域乘以頻率通道帶寬為fs/N的 “矩形” 濾波器函數對應到時域上的sinc函數。但因時域上矩形窗長度的限制,通道形狀難以 “矩形”。為了獲得上述 “矩形” 通道的性能,需要增加時域矩形窗的長度,對時域上的sinc函數取更長的序列(如8倍的N/fs,見圖1(7)),再進行快速傅里葉變換運算后得到更多的數據點(8倍, 見圖1(8))。此時的通道函數頻域響應為加長的矩形窗函數和sinc函數各自的頻域響應的卷積結果。通過改變時域矩形窗的長度,選擇不同的窗函數可以改善通道形狀,從而影響每個通道的平坦度和矩形系數。由于只加長了時域sinc函數截斷長度,其主瓣寬度并未改變,所以通道帶寬仍為fs/N。但此時的頻域數據量為原來的8倍,相鄰數據點之間間隔為fs/8N,因此只需要進行8倍降采樣即可得到通道帶寬為fs/N頻域數據。多相濾波[5-6]即為實現該計算的一種手段。與對加長窗長度和改變窗函數的序列直接進行8N個點的快速傅里葉變換相比,多相濾波計算量大大減少;而與采用矩形窗只進行N個點的快速傅里葉變換相比,計算量增加也不多(見表1),其硬件實現結構如圖3。

    表1 多相結構與直接離散傅里葉變換的乘法次數對比

    1.3 濾波器性能優(yōu)化分析

    如上兩節(jié)所述,通過改變時域sinc函數截取長度和選擇不同的窗函數可以改善通道帶通性能。圖4(1)展示了采用相同的時域漢明窗和相同的時域sinc函數的情況下,只改變sinc函數的截取長度,在符合奈奎斯特采樣定理要求下用相同采樣率分別取一個主瓣長度的時域sinc函數(虛線)和16倍主瓣長度的時域sinc函數(實線)所對應的多相濾波的通道形狀。圖4(2)是對兩個相鄰通道交疊處的放大。當sinc函數截取長度越長時,對應的通道函數越趨近于矩形。時域sinc函數相同,故時域主瓣寬度保持不變,因此對應的頻域通道的寬度也相同。當時域的矩形窗足夠長時,相鄰通道的頻譜響應的交疊點大致相同(在-3 dB左右)。

    圖3 多相濾波硬件實現結構圖

    圖4 不同截取長度和不同主瓣寬度的sinc函數對頻域通帶形狀的影響

    2 多通道數字終端的實現與測試

    和直接快速傅里葉變換相比,由于多相濾波的降采樣處理,計算量大大減少,在射電天文數據處理方面得到了廣泛的應用。國內云南天文臺、上海天文臺、國家天文臺太陽射電團組等都對該技術進行了研究和應用。美國UC Berkeley的CASPER團組開發(fā)了一系列的射電天文數字終端軟硬件平臺,對于不需了解底層硬件的天文學家是一個不錯的選擇,FAST終端開發(fā)將采用該平臺。本文在CASPER終端平臺上實現了一個觀測帶寬為400 MHz、通道數16 k,可通過10 GbE網卡傳輸數據的數字終端。該終端可在FAST多波束河外中性氫巡天中得到應用。河外中性氫巡天觀測的需求如下:

    帶寬:400 MHz,觀測中性氫紅移從0到0.26,頻率范圍為1 050 MHz~1 450 MHz。

    采樣精度:8 bits。

    通道數:12 k。

    河外中性氫觀測速度分辨率應達到10 km/s,結合帶寬,通道數為

    (1)

    在設計中,多相結構(PolyphaseFilterBank,PFB)的通道形狀如圖4中實線,直接快速傅里葉變換的通道形狀為圖5中虛線。設計中多相結構的窗長度是直接進行快速傅里葉變換的矩形窗長度的8倍。多相結構第一旁瓣抑制約-60dB,快速傅里葉變換第一旁瓣抑制約為-15dB。取矩形系數為3dB帶寬和60dB帶寬的比值(越接近1越好),則多相結構的矩形系數為0.63,快速傅里葉變換的矩形系數為0.42,多相結構通帶形狀更優(yōu)。圖6展示了設計中的幾個通道形狀。

    圖5 多相結構和快速傅里葉變換的單個通道頻域響應

    圖6 設計中多通道形狀展示

    該多通道數字終端在MATLAB/Simulink環(huán)境下進行仿真設計,利用MATLAB和Xilinx開發(fā)的SystemGenerator[7]將mdl*https://casper.berkeley.edu/wiki/Wideband_Spectrometer文件 “翻譯” 成硬件語言運行,通過python語言操作硬件模塊實現顯示和存儲。mdl文件截圖見圖7。

    圖7 硬件設計圖

    信號處理流程為:400MHz的模擬信號首先經過模擬數字轉換器采樣輸出量化比特數為8bit的數字信號;之后數字信號經過所設計的多相濾波和傅里葉變換相結合的多通道結構實現多通道輸出,最后將輸出頻譜數據經過累加后顯示并通過10GbE網卡傳送到其他設備上存儲。

    現場可編程門陣列時鐘頻率為200MHz,而10GbE網卡工作頻率為156.25MHz,最大緩存空間為8 704Bytes。接口工作在異步模式下,若不加控制直接將現場可編程門陣列處理數據不斷地傳送給10GbE模塊,緩存最終將溢出?,F場可編程門陣列的頻譜數據是以突發(fā)數據包(burst)的形式向10GbE模塊傳送并且每條頻譜的大小(213×8B=216B)超過最大緩存空間。利用兩個突發(fā)數據包之間的間隔的累加時間,可以采用先進先出模塊(First In First Out, FIFO)[8]解決這個問題。每個突發(fā)數據包先存放在先進先出模塊中,再由先進先出模塊向10 GbE模塊發(fā)送傳送請求,傳送數據大小由先進先出模塊控制,可以將一條頻譜數據分成8次(23次)傳輸給10 GbE,邏輯控制模塊如圖8。先進先出模塊每次向10 GbE模塊傳送一個突發(fā)數據包中的210個數據,等待(212-210)個時鐘后(給10 GbE模塊足夠的傳送時間),傳送第2次210個數據,再等待(212-210)個時鐘,以此類推,直到傳送第1個突發(fā)數據包中的最后一組210個數據。結合輸入先進先出模塊的有效信號valid以及FIFO本身的空信號進行其他的邏輯控制,最終形成了圖8的傳送邏輯控制模塊。

    圖8 FIFO傳送邏輯控制模塊

    10 GbE模塊[9]采用UDP協(xié)議[10]通過以太網傳送數據包給指定服務器。輸入數據比特數為64,設定好目的ip和目的端口后,通過tx_valid置1輸入數據,在輸入數據的最后一個時鐘同時將tx_end_of_frame置1完成一次數據輸入,并進行UDP封裝后傳送向目的ip和端口。在目的ip處,可通過Wireshark Network Analyzer監(jiān)測和查看收到的UDP包,并通過程序將UDP封裝的頭部文件去除后以二進制文件格式存儲在指定目錄下。

    2.1 多通道頻譜數字終端的測試

    測試實驗中使用安捷倫N5181A信號發(fā)生器進行掃頻來測試該終端的性能。信號發(fā)生器設置為步進頻率自動單次掃描模式,輸出信號功率為0 dBm,每個單頻信號持續(xù)時間500 ms。信號源輸出的信號傳給模擬數字轉換器,經模擬數字轉換器采樣后傳送給數字信號處理板卡ROACH進行多相濾波處理。圖9從下至上依次為安捷倫N5181A信號發(fā)生器、控制服務器和終端開發(fā)平臺。圖10掃頻頻段為110 MHz~115 MHz,積分時間0.025 3 s。圖11為當輸入單頻信號從49.95 MHz遞增至50.1 MHz(步進為0.001 MHz),第2049和2050兩個通道(由左及右)的輸出值變化情況,由此可以掃出通道的形狀。該設計采用8段數據,每段數據的長度為214個點,由于縱坐標采用功率值,故交疊點在-6 dB位置。

    圖9 掃頻測試系統(tǒng)

    圖10 掃頻測試

    圖11 通道形狀圖

    3 結束語

    本文回顧了濾波的基本原理,分析了濾波器性能優(yōu)化方法,并對多相濾波技術在多通道頻譜分析中的應用做了說明。通過增加時域矩形窗的長度和額外加窗函數的方法,可對濾波器帶通性能進行優(yōu)化。與直接快速傅里葉變換相比,多相濾波器由于在時域上采用了更長的矩形窗,其濾波通道性能得到很大改善,但計算量增加不多。在UC Berkeley的CASPER團組研制的終端開發(fā)平臺上利用多相濾波技術實現了一個觀測帶寬為400 MHz,通道數16 k的數字終端,并完成了頻譜的生成顯示和10 GbE的傳輸。該終端可應用于FAST多波束河外中性氫巡天觀測。

    [1] Ewen H I, Purcell E M. Observation of a line in the galactic radio spectrum: radiation from galactic hydrogen at 1 420 Mc./sec [J]. Nature, 1951, 168: 356.

    [2] 陳敬英, 王書浩, 夏志國, 等. 云南天文臺太陽射電頻譜儀的網絡系統(tǒng)[J]. 云南天文臺臺刊, 2000(1): 41-48. Cheng Jingying, Wang Shuhao, Xia Zhiguo, et al. The network system of solar radio spectrographs at Yunnan Observatory[J]. Publications of Yunnan Observatory, 2000(1): 41-48.

    [3] 項英, 張秀忠. 數字濾波技術在射電天文測量中的應用[J]. 天文學進展, 2004, 22(2): 95-103. Xiang Ying, Zhang Xiuzhong. Application of digital filtering in radio astronomy[J]. Progress in Astronomy, 2004, 22(2): 95-103.

    [4] 董亮, 汪敏, 苗愛敏, 等. 基于FPGA嵌入式的多通道功率監(jiān)測處理器設計[J]. 電子測量技術, 2010, 33(3): 58-63+105. Dong Liang, Wang Min, Miao Aimin, et al. Multi-channel power detecting processor’s designing based on FPGA embedded system[J]. Electronic Measurement Technology, 2010, 33(3): 58-63+105.

    [5] 陳嵐, 張秀忠. 用于VLBI數字基帶轉換的多相濾波技術研究[J]. 天文學進展, 2008, 26(1): 87-94. Chen Lan, Zhang Xiuzhong. The study of DBBC based on poly-phase filter banks and FFT in VLBI[J].Progress in Astronomy, 2008, 26(1): 87-94.

    [6] 朱凱, 甘恒謙, 朱巖, 等. 多相濾波器組譜分析方法的性能討論與天文觀測應用[J]. 天文研究與技術——國家天文臺臺刊, 2011, 8(1): 81-90. Zhukai, Gan Hengqian, Zhu Yan, et al. On the performance of spectral analysis with a polyphase filter bank and its application in radio astronomy[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2011, 8(1): 81-90.

    [7] 紀志成. FPGA數字信號處理設計教程[M]. 西安: 西安電子科技大學出版社, 2008.

    [8] 張西洋, 何樂生, 董亮, 等. 800~975 MHz太陽射電數字觀測終端的設計與實現[J]. 天文研究與技術——國家天文臺臺刊, 2014, 11(2): 118-126. Zhang Xiyang, He lesheng, Dong Liang, et al. Design and implementation of a digital observation terminal for solar radio observation within the 800MHz-975MHz band[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2014, 11(2): 118-126.

    [9] 劉東亮, John Ford, Glen Langston, 等. 基于IBOB的射電望遠鏡寬帶數字頻譜儀系統(tǒng)的設計與測試[J]. 天文研究與技術——國家天文臺臺刊, 2010, 7(1): 27-32. Liu Dongliang, John Ford, Glen Langston, et al. Design and test of an IBOB based wideband digital spectrometer system for a radio telescope[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2010, 7(1): 27-32.

    [10]崔鶴, 劉云清, 盛家進. 基于FPGA的UDP/IP協(xié)議棧的研究與實現[J]. 長春理工大學學報: 自然科學版, 2014, 37(2): 133-137. Cui He, Liu Yunqing, Sheng Jiajin. Research and implementation on UDP/IP protocol stack base on FPGA[J]. Journal of Changchun University of Science and Technology: Natural Science Edition, 2014, 37(2): 133-137.

    A Design of a Multi-Channel Digital Backend Based on a Polyphase Filter Bank

    Zhao Xin1,2, Jin Chengjin1, Zhu Yan1, Yue Youling1, Wu Junlin1,2

    (1. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China, Email: cjjin@bao.ac.cn; 2. University of Chinese Academy of Sciences, Beijing 100049, China)

    The digital-backend technique has become the main-stream signal-processing technology of radio astronomy and has been applied in China and abroad. In this paper we review the principles of signal filtering and application of polyphase filter banks (PFBs) in multi-channel signal filtering. We particularly discuss approaches to optimize the bandpass performance using a PFB. Using the backend-developing platform provided by the CASPER of the UC Berkeley, we have designed and constructed a digital backend with a bandwidth of 400MHz divided into 16k channels (corresponding to a frequency resolution of 25kHz). We have also carried out a frequency-sweep test of the backend. The processing of the digital backend has been designed to be as follows. Input analog signals are converted into digital signals by an ADC. The digital signals are filtered by a PFB in multiple frequency channels, and the filtered signals in each channel are transformed by the FFT block to yield a power spectrum. The power spectra in multiple channels are summed and sent to other devices through a 10GbE interface for further processing or storage. The backend as designed can be applied in extragalactic HI surveys using the FAST radio telescope.

    Polyphase filter bank; Multi-channel fitering; 10GbE; Digital backend

    國家重點基礎研究發(fā)展計劃 (973計劃) (2012CB821800);國家自然科學基金 (11373045);國家重點基礎研究發(fā)展計劃 (973計劃) (2013CB837900);國家自然科學基金國際合作與交流項目 (11261140641) 和中國科學院重點部署項目 (KJZD-EW-T01) 資助.

    2014-12-28;修定日期:2015-02-11

    趙 欣,女,碩士. 研究方向:天文數據后端. Email: zhaoxin@bao ac cn

    金乘進,男,研究員. 研究方向:接收機及后端. Email: cjjin@bao ac cn

    P111.44

    A

    1672-7673(2015)04-0495-08

    CN 53-1189/P ISSN 1672-7673

    猜你喜歡
    傅里葉頻域時域
    雙線性傅里葉乘子算子的量化加權估計
    基于小波降噪的稀疏傅里葉變換時延估計
    測控技術(2018年7期)2018-12-09 08:58:26
    基于時域信號的三電平逆變器復合故障診斷
    測控技術(2018年11期)2018-12-07 05:49:02
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設計
    雷達學報(2018年3期)2018-07-18 02:41:34
    基于極大似然準則與滾動時域估計的自適應UKF算法
    基于傅里葉變換的快速TAMVDR算法
    基于改進Radon-Wigner變換的目標和拖曳式誘餌頻域分離
    基于時域逆濾波的寬帶脈沖聲生成技術
    一種基于頻域的QPSK窄帶干擾抑制算法
    基于頻域伸縮的改進DFT算法
    電測與儀表(2015年3期)2015-04-09 11:37:24
    亚洲欧美日韩东京热| 国产精品福利在线免费观看| 亚洲国产精品专区欧美| 欧美zozozo另类| 三级毛片av免费| av黄色大香蕉| 成人高潮视频无遮挡免费网站| 久久久久国产网址| 中文天堂在线官网| 日韩av不卡免费在线播放| 99久久中文字幕三级久久日本| 一级毛片久久久久久久久女| 麻豆av噜噜一区二区三区| 精品一区二区三卡| 欧美三级亚洲精品| 久久这里只有精品中国| 一个人观看的视频www高清免费观看| 久久久久精品性色| 精品一区二区三卡| 国产精品一区二区性色av| 美女主播在线视频| 偷拍熟女少妇极品色| 简卡轻食公司| 日韩一区二区视频免费看| 高清午夜精品一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 美女黄网站色视频| 欧美zozozo另类| 女人被狂操c到高潮| a级一级毛片免费在线观看| 亚洲成色77777| 日韩欧美精品v在线| 国国产精品蜜臀av免费| 在线免费十八禁| 又黄又爽又刺激的免费视频.| 亚洲经典国产精华液单| 久久精品久久久久久噜噜老黄| 青春草亚洲视频在线观看| 久99久视频精品免费| 国产免费又黄又爽又色| 久久久久久九九精品二区国产| 精品久久久久久电影网| 国产高清有码在线观看视频| 一个人观看的视频www高清免费观看| 久久久国产一区二区| 2018国产大陆天天弄谢| 亚洲电影在线观看av| 国产黄片美女视频| 丝袜喷水一区| 99热这里只有精品一区| 免费看av在线观看网站| 男人狂女人下面高潮的视频| 欧美极品一区二区三区四区| 九九爱精品视频在线观看| 成人亚洲精品av一区二区| 99热这里只有精品一区| 高清在线视频一区二区三区| 亚洲国产成人一精品久久久| 亚洲精品乱码久久久v下载方式| 99热这里只有精品一区| 国产久久久一区二区三区| 老师上课跳d突然被开到最大视频| 天堂av国产一区二区熟女人妻| 91精品一卡2卡3卡4卡| 免费观看精品视频网站| 寂寞人妻少妇视频99o| 日韩制服骚丝袜av| 午夜福利在线观看吧| 日本欧美国产在线视频| 成人美女网站在线观看视频| 一个人看的www免费观看视频| 国产精品综合久久久久久久免费| 激情 狠狠 欧美| 99热6这里只有精品| 秋霞伦理黄片| 91狼人影院| 男人舔女人下体高潮全视频| 如何舔出高潮| 国产精品美女特级片免费视频播放器| 美女高潮的动态| 一级av片app| 如何舔出高潮| 天天躁日日操中文字幕| 亚洲av.av天堂| 国产黄频视频在线观看| 国产一级毛片在线| 啦啦啦啦在线视频资源| 天堂俺去俺来也www色官网 | 能在线免费看毛片的网站| 嫩草影院精品99| 91久久精品国产一区二区三区| 亚洲精品影视一区二区三区av| 成人鲁丝片一二三区免费| 国产精品美女特级片免费视频播放器| 成年女人在线观看亚洲视频 | 久久精品人妻少妇| 亚洲欧美成人综合另类久久久| 久久韩国三级中文字幕| 亚洲精品中文字幕在线视频 | 亚洲国产日韩欧美精品在线观看| 国产白丝娇喘喷水9色精品| 最近最新中文字幕免费大全7| 欧美丝袜亚洲另类| 免费观看精品视频网站| 精品一区二区三区视频在线| 国产精品久久久久久av不卡| 午夜精品在线福利| 精品久久久精品久久久| 国产乱来视频区| 男插女下体视频免费在线播放| 亚洲欧美清纯卡通| 18禁裸乳无遮挡免费网站照片| av在线天堂中文字幕| 久久久色成人| 黄色欧美视频在线观看| 精品国产露脸久久av麻豆 | 啦啦啦啦在线视频资源| 国产黄色小视频在线观看| 国产黄a三级三级三级人| 一级a做视频免费观看| 夫妻午夜视频| or卡值多少钱| 国产精品精品国产色婷婷| 嫩草影院精品99| or卡值多少钱| 最近手机中文字幕大全| 亚洲精品影视一区二区三区av| 建设人人有责人人尽责人人享有的 | 免费看光身美女| 久久精品国产亚洲av涩爱| 波多野结衣巨乳人妻| 天堂av国产一区二区熟女人妻| 国产免费又黄又爽又色| 伊人久久精品亚洲午夜| 黄色一级大片看看| 久久亚洲国产成人精品v| 免费av观看视频| 日韩欧美一区视频在线观看 | 亚洲综合精品二区| 狂野欧美激情性xxxx在线观看| 最近手机中文字幕大全| 午夜视频国产福利| 日日啪夜夜撸| 免费观看性生交大片5| 亚洲精品国产av成人精品| 国产成人精品一,二区| 成人亚洲欧美一区二区av| 日韩av免费高清视频| 嫩草影院入口| 国产成人a区在线观看| 色综合色国产| 国产精品爽爽va在线观看网站| 高清欧美精品videossex| 天堂网av新在线| 亚洲激情五月婷婷啪啪| 国国产精品蜜臀av免费| 久久精品国产亚洲av天美| 免费看不卡的av| 亚洲成人精品中文字幕电影| av在线天堂中文字幕| 成人二区视频| 日本三级黄在线观看| 十八禁国产超污无遮挡网站| 亚洲性久久影院| 国产毛片a区久久久久| 国产色婷婷99| 人人妻人人看人人澡| 三级经典国产精品| 人人妻人人看人人澡| 亚洲国产高清在线一区二区三| 国产 一区精品| 91久久精品国产一区二区三区| 成年人午夜在线观看视频 | 极品少妇高潮喷水抽搐| 欧美+日韩+精品| 精品不卡国产一区二区三区| 哪个播放器可以免费观看大片| 日韩中字成人| 国产av在哪里看| 搡女人真爽免费视频火全软件| 99久久精品热视频| 久久99热这里只频精品6学生| 能在线免费看毛片的网站| 亚洲人成网站在线播| 免费观看精品视频网站| 大陆偷拍与自拍| 国产精品国产三级专区第一集| 亚洲av成人av| 免费观看性生交大片5| 干丝袜人妻中文字幕| 熟妇人妻不卡中文字幕| 欧美+日韩+精品| av在线蜜桃| eeuss影院久久| 婷婷六月久久综合丁香| 午夜福利在线观看免费完整高清在| 美女国产视频在线观看| 99久久中文字幕三级久久日本| 国产成人精品婷婷| kizo精华| 欧美日韩亚洲高清精品| 国产一区有黄有色的免费视频 | 最近的中文字幕免费完整| 最后的刺客免费高清国语| 啦啦啦中文免费视频观看日本| 舔av片在线| 亚洲欧洲日产国产| 中文字幕av在线有码专区| 国产一区二区三区综合在线观看 | 日韩视频在线欧美| 亚洲在久久综合| 国产白丝娇喘喷水9色精品| 久久久久精品久久久久真实原创| 插阴视频在线观看视频| 国产在线男女| 久久99精品国语久久久| 国产精品久久久久久精品电影小说 | 国产精品无大码| 久久久久九九精品影院| 一级毛片aaaaaa免费看小| 最近中文字幕2019免费版| 亚洲人成网站在线观看播放| 久久这里只有精品中国| 国产黄色免费在线视频| av福利片在线观看| 男女啪啪激烈高潮av片| 久久精品国产自在天天线| 黑人高潮一二区| 免费观看无遮挡的男女| 卡戴珊不雅视频在线播放| 精品久久久久久久末码| 国产精品熟女久久久久浪| 久久草成人影院| 日日啪夜夜爽| 国产免费又黄又爽又色| 伦精品一区二区三区| 黄色一级大片看看| 午夜福利高清视频| 国产亚洲午夜精品一区二区久久 | 性插视频无遮挡在线免费观看| 国产激情偷乱视频一区二区| 中文字幕av成人在线电影| 久久午夜福利片| 欧美97在线视频| 亚洲欧洲国产日韩| 最后的刺客免费高清国语| 欧美激情国产日韩精品一区| 国产成人免费观看mmmm| 欧美成人午夜免费资源| 午夜福利在线观看免费完整高清在| 69人妻影院| 99热这里只有精品一区| 亚洲精品日韩av片在线观看| 嫩草影院入口| 久久精品国产亚洲av涩爱| www.色视频.com| 成年免费大片在线观看| 日本免费在线观看一区| 免费观看无遮挡的男女| 午夜免费男女啪啪视频观看| 亚洲av免费在线观看| 在线免费观看不下载黄p国产| 欧美最新免费一区二区三区| 2018国产大陆天天弄谢| 亚洲av电影不卡..在线观看| 男女边摸边吃奶| 黄片无遮挡物在线观看| 国产精品一区www在线观看| 色网站视频免费| 国产亚洲91精品色在线| 中文在线观看免费www的网站| 大陆偷拍与自拍| 国产黄色小视频在线观看| 国产乱人偷精品视频| 人妻夜夜爽99麻豆av| 日韩三级伦理在线观看| 蜜桃久久精品国产亚洲av| 日本黄大片高清| 听说在线观看完整版免费高清| 青春草视频在线免费观看| 国内揄拍国产精品人妻在线| 老司机影院成人| 看十八女毛片水多多多| 一级毛片久久久久久久久女| 久久久久免费精品人妻一区二区| 麻豆国产97在线/欧美| 久久精品夜色国产| 国产精品一区二区在线观看99 | 亚洲怡红院男人天堂| 乱码一卡2卡4卡精品| 男女边吃奶边做爰视频| 日韩大片免费观看网站| 亚洲欧洲日产国产| 成人亚洲精品一区在线观看 | 欧美成人精品欧美一级黄| 日本三级黄在线观看| 不卡视频在线观看欧美| 黄色配什么色好看| 国产爱豆传媒在线观看| 黑人高潮一二区| 一级毛片我不卡| 欧美极品一区二区三区四区| 午夜视频国产福利| 亚洲激情五月婷婷啪啪| 久久久久久久久大av| 能在线免费看毛片的网站| 亚洲aⅴ乱码一区二区在线播放| 精品一区二区三区视频在线| 久久久久久久久久久丰满| 亚洲av中文字字幕乱码综合| 91狼人影院| 国产爱豆传媒在线观看| 精品一区二区三区视频在线| 亚洲欧美日韩东京热| 丰满乱子伦码专区| 欧美xxxx性猛交bbbb| 老司机影院毛片| 美女主播在线视频| 秋霞在线观看毛片| 国产一区有黄有色的免费视频 | 亚洲国产欧美在线一区| 亚洲人成网站在线播| 中文欧美无线码| 一个人观看的视频www高清免费观看| 欧美成人一区二区免费高清观看| 欧美+日韩+精品| 亚洲av免费在线观看| 精品久久久久久成人av| 国产伦一二天堂av在线观看| 久久久久久久午夜电影| 国产亚洲av片在线观看秒播厂 | 日韩成人av中文字幕在线观看| 久久久久久久国产电影| 国产一区二区三区av在线| 亚洲va在线va天堂va国产| 国产精品.久久久| 一级毛片 在线播放| 免费看日本二区| 久久热精品热| 2021少妇久久久久久久久久久| 免费av不卡在线播放| 国产精品一区二区三区四区久久| 老女人水多毛片| 久久鲁丝午夜福利片| 久久精品综合一区二区三区| 亚洲精品第二区| 国产黄片美女视频| 免费av毛片视频| 中文欧美无线码| 午夜视频国产福利| 黄色配什么色好看| 久久这里只有精品中国| 免费大片18禁| 国产人妻一区二区三区在| 国产 一区 欧美 日韩| 欧美一级a爱片免费观看看| 夫妻性生交免费视频一级片| 天天躁夜夜躁狠狠久久av| 偷拍熟女少妇极品色| 欧美日韩国产mv在线观看视频 | 欧美丝袜亚洲另类| 91精品一卡2卡3卡4卡| 天天躁夜夜躁狠狠久久av| 欧美性感艳星| 别揉我奶头 嗯啊视频| 免费观看a级毛片全部| 人妻少妇偷人精品九色| 亚洲av中文av极速乱| 国产精品久久久久久av不卡| 九九在线视频观看精品| 少妇的逼好多水| 中文字幕免费在线视频6| 国内少妇人妻偷人精品xxx网站| 春色校园在线视频观看| 一个人看视频在线观看www免费| 精品午夜福利在线看| 男人和女人高潮做爰伦理| 熟女电影av网| 最近最新中文字幕免费大全7| 天堂俺去俺来也www色官网 | 精品一区在线观看国产| 亚洲av电影不卡..在线观看| 久久99精品国语久久久| 禁无遮挡网站| 伦理电影大哥的女人| 床上黄色一级片| 观看免费一级毛片| 久久鲁丝午夜福利片| 一级毛片电影观看| 高清视频免费观看一区二区 | 欧美极品一区二区三区四区| 亚洲av成人精品一二三区| 在线 av 中文字幕| 看免费成人av毛片| 哪个播放器可以免费观看大片| 亚洲欧美成人精品一区二区| 人人妻人人澡人人爽人人夜夜 | 内射极品少妇av片p| 国产av在哪里看| 一级a做视频免费观看| 秋霞在线观看毛片| 免费人成在线观看视频色| 欧美+日韩+精品| 又大又黄又爽视频免费| 精品久久久久久久久av| 日韩 亚洲 欧美在线| 大陆偷拍与自拍| 国产视频首页在线观看| 亚洲国产成人一精品久久久| 22中文网久久字幕| 一区二区三区免费毛片| 最近中文字幕2019免费版| av在线天堂中文字幕| 国产探花极品一区二区| av在线播放精品| 欧美日本视频| 久久精品久久精品一区二区三区| 熟妇人妻久久中文字幕3abv| 亚洲婷婷狠狠爱综合网| 成年女人看的毛片在线观看| 18禁裸乳无遮挡免费网站照片| 国产成人精品久久久久久| 国产片特级美女逼逼视频| 全区人妻精品视频| 六月丁香七月| 午夜视频国产福利| 欧美日韩亚洲高清精品| 日韩欧美国产在线观看| av专区在线播放| 国产白丝娇喘喷水9色精品| 日韩精品有码人妻一区| 欧美成人a在线观看| 亚洲av中文字字幕乱码综合| 日韩精品青青久久久久久| 免费观看性生交大片5| 大又大粗又爽又黄少妇毛片口| av网站免费在线观看视频 | 一区二区三区高清视频在线| 高清av免费在线| a级毛色黄片| 欧美xxⅹ黑人| 噜噜噜噜噜久久久久久91| 日韩 亚洲 欧美在线| 亚洲精品自拍成人| 国产精品熟女久久久久浪| 日韩av在线免费看完整版不卡| 国产片特级美女逼逼视频| 亚洲精品久久久久久婷婷小说| 色播亚洲综合网| 一个人看的www免费观看视频| 国产男人的电影天堂91| 一级毛片aaaaaa免费看小| av免费观看日本| 久久精品久久久久久久性| 国产一区二区在线观看日韩| 18禁动态无遮挡网站| 我的女老师完整版在线观看| 乱码一卡2卡4卡精品| 精品久久久久久久久av| 建设人人有责人人尽责人人享有的 | 日韩欧美 国产精品| 成人午夜高清在线视频| 久久久久久久久久成人| www.色视频.com| 一级爰片在线观看| 免费观看a级毛片全部| 亚洲不卡免费看| 日韩精品青青久久久久久| 美女cb高潮喷水在线观看| 中文天堂在线官网| 国产精品国产三级国产专区5o| 欧美bdsm另类| 99热这里只有精品一区| 午夜激情久久久久久久| 久久精品国产鲁丝片午夜精品| 国产成人精品婷婷| 内射极品少妇av片p| 在线观看美女被高潮喷水网站| 波野结衣二区三区在线| 超碰av人人做人人爽久久| 女人久久www免费人成看片| 日本欧美国产在线视频| 大话2 男鬼变身卡| 日韩欧美精品免费久久| 一级毛片aaaaaa免费看小| av免费观看日本| freevideosex欧美| av在线天堂中文字幕| 国产色婷婷99| 五月伊人婷婷丁香| 人妻系列 视频| 成年女人看的毛片在线观看| 亚洲高清免费不卡视频| 免费人成在线观看视频色| 亚洲精品色激情综合| 亚洲一区高清亚洲精品| 一个人看视频在线观看www免费| 亚洲综合色惰| 久久久久久久国产电影| 成年女人看的毛片在线观看| 一级二级三级毛片免费看| 久久久久久久久中文| 欧美日韩精品成人综合77777| 亚洲精品中文字幕在线视频 | 亚洲欧美日韩卡通动漫| 亚洲欧美成人精品一区二区| 国模一区二区三区四区视频| 91久久精品电影网| 成年人午夜在线观看视频 | 欧美bdsm另类| 国产色婷婷99| 一区二区三区四区激情视频| 午夜精品国产一区二区电影 | 五月伊人婷婷丁香| 一级爰片在线观看| xxx大片免费视频| 伊人久久国产一区二区| 七月丁香在线播放| 男人狂女人下面高潮的视频| 国产成年人精品一区二区| 国产伦在线观看视频一区| 搡老乐熟女国产| 亚洲色图av天堂| 久久精品综合一区二区三区| 97超碰精品成人国产| 国产精品久久久久久av不卡| 亚洲精华国产精华液的使用体验| 少妇裸体淫交视频免费看高清| 你懂的网址亚洲精品在线观看| 最近手机中文字幕大全| 亚洲自拍偷在线| 亚洲国产精品国产精品| 免费播放大片免费观看视频在线观看| av专区在线播放| 亚洲av成人av| 国产精品一区二区三区四区免费观看| 精品欧美国产一区二区三| 国产精品一区二区性色av| 欧美成人精品欧美一级黄| 日韩三级伦理在线观看| 九九爱精品视频在线观看| 男人舔女人下体高潮全视频| 又爽又黄无遮挡网站| 欧美成人精品欧美一级黄| 亚洲成色77777| 国产精品女同一区二区软件| 伦精品一区二区三区| 色综合站精品国产| 波野结衣二区三区在线| 性色avwww在线观看| 街头女战士在线观看网站| 亚洲欧美一区二区三区国产| 最新中文字幕久久久久| 欧美潮喷喷水| 成人漫画全彩无遮挡| 国产白丝娇喘喷水9色精品| 午夜福利在线观看吧| 亚洲真实伦在线观看| 日本一本二区三区精品| 高清av免费在线| 国产白丝娇喘喷水9色精品| 亚洲av电影在线观看一区二区三区 | 卡戴珊不雅视频在线播放| 国产高清三级在线| 青春草亚洲视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产一区二区三区av在线| 亚洲国产精品成人综合色| 中文天堂在线官网| 色尼玛亚洲综合影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 大又大粗又爽又黄少妇毛片口| 婷婷六月久久综合丁香| 少妇的逼水好多| 好男人视频免费观看在线| 欧美日韩综合久久久久久| 欧美丝袜亚洲另类| 日本与韩国留学比较| 亚洲不卡免费看| 乱人视频在线观看| 欧美一级a爱片免费观看看| 丝袜美腿在线中文| 日韩一本色道免费dvd| 在线免费观看不下载黄p国产| 我的女老师完整版在线观看| 精品少妇黑人巨大在线播放| 久久午夜福利片| 国产精品无大码| 成人欧美大片| 中国国产av一级| 亚洲经典国产精华液单| a级毛色黄片| 国产亚洲91精品色在线| 美女被艹到高潮喷水动态| 性插视频无遮挡在线免费观看| 人人妻人人澡欧美一区二区| 欧美 日韩 精品 国产| 亚洲综合色惰| 国产成人freesex在线| 在现免费观看毛片| 亚洲精品自拍成人| 亚洲人成网站高清观看| 亚洲国产精品成人综合色| 看十八女毛片水多多多| 亚洲精品成人av观看孕妇| 亚洲自偷自拍三级| 日韩,欧美,国产一区二区三区| 五月天丁香电影| 九九久久精品国产亚洲av麻豆| 伦精品一区二区三区| 99九九线精品视频在线观看视频| 美女脱内裤让男人舔精品视频| 国产精品久久久久久久久免| 国产极品天堂在线| 嘟嘟电影网在线观看|