• <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
    亚洲七黄色美女视频| 亚洲成国产人片在线观看| 国产有黄有色有爽视频| 国产日韩欧美在线精品| 黄色视频在线播放观看不卡| 国产一区二区 视频在线| 精品国产乱码久久久久久小说| 中文字幕制服av| 免费在线观看视频国产中文字幕亚洲| 国产区一区二久久| 男女无遮挡免费网站观看| 国产精品亚洲一级av第二区| 老熟女久久久| 少妇精品久久久久久久| 国产深夜福利视频在线观看| www日本在线高清视频| 久久精品国产99精品国产亚洲性色 | 国产91精品成人一区二区三区 | 在线天堂中文资源库| 亚洲伊人色综图| 成人永久免费在线观看视频 | 黄网站色视频无遮挡免费观看| 婷婷丁香在线五月| 交换朋友夫妻互换小说| 国产91精品成人一区二区三区 | 天天操日日干夜夜撸| 伊人久久大香线蕉亚洲五| 啦啦啦在线免费观看视频4| 一区在线观看完整版| 国产伦理片在线播放av一区| 欧美日韩精品网址| 黄色视频在线播放观看不卡| 999久久久国产精品视频| 国产xxxxx性猛交| 国产不卡一卡二| 国产无遮挡羞羞视频在线观看| 免费观看a级毛片全部| 人人妻人人澡人人看| 老司机福利观看| 丰满人妻熟妇乱又伦精品不卡| 免费久久久久久久精品成人欧美视频| 午夜成年电影在线免费观看| 99九九在线精品视频| 国产欧美日韩一区二区三区在线| 欧美精品高潮呻吟av久久| 国产欧美日韩一区二区三| av有码第一页| 最近最新免费中文字幕在线| tocl精华| 日本撒尿小便嘘嘘汇集6| 国产极品粉嫩免费观看在线| 日韩熟女老妇一区二区性免费视频| a级片在线免费高清观看视频| 涩涩av久久男人的天堂| 悠悠久久av| 老汉色∧v一级毛片| 成年动漫av网址| 午夜福利在线免费观看网站| 久久人妻熟女aⅴ| 日本黄色日本黄色录像| 狠狠婷婷综合久久久久久88av| 成在线人永久免费视频| 成年人免费黄色播放视频| 亚洲国产精品一区二区三区在线| 人人妻人人爽人人添夜夜欢视频| 大香蕉久久网| 又大又爽又粗| 麻豆国产av国片精品| 超色免费av| 大型av网站在线播放| 国产精品自产拍在线观看55亚洲 | 日韩免费高清中文字幕av| 老司机福利观看| 午夜91福利影院| 麻豆乱淫一区二区| 亚洲国产看品久久| 久久中文字幕一级| 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 国产精品免费视频内射| 亚洲第一av免费看| 天天添夜夜摸| 午夜老司机福利片| 亚洲精品久久午夜乱码| 少妇 在线观看| 国产xxxxx性猛交| 国产一区二区三区视频了| 香蕉久久夜色| 日日夜夜操网爽| 91老司机精品| 黄色a级毛片大全视频| 久久影院123| 午夜久久久在线观看| 久久久久视频综合| 日韩一卡2卡3卡4卡2021年| 国产熟女午夜一区二区三区| 久久天躁狠狠躁夜夜2o2o| 亚洲av第一区精品v没综合| 国产一区二区三区在线臀色熟女 | 一个人免费在线观看的高清视频| 精品免费久久久久久久清纯 | 亚洲专区中文字幕在线| 亚洲专区中文字幕在线| 久久精品国产亚洲av香蕉五月 | 亚洲av成人不卡在线观看播放网| 日本欧美视频一区| 亚洲欧美精品综合一区二区三区| 高清黄色对白视频在线免费看| 自线自在国产av| 女性被躁到高潮视频| 美女国产高潮福利片在线看| 欧美亚洲日本最大视频资源| 亚洲第一av免费看| 亚洲中文av在线| 久久精品亚洲av国产电影网| 操美女的视频在线观看| 久久久久国产一级毛片高清牌| 午夜激情久久久久久久| 中文字幕精品免费在线观看视频| 在线观看66精品国产| 午夜福利欧美成人| 国产一区二区激情短视频| 国产一卡二卡三卡精品| 在线观看一区二区三区激情| 国产高清视频在线播放一区| 亚洲精品一卡2卡三卡4卡5卡| 91大片在线观看| 亚洲综合色网址| 三级毛片av免费| 男女边摸边吃奶| 天天操日日干夜夜撸| 美女高潮到喷水免费观看| 亚洲成人免费av在线播放| 老汉色∧v一级毛片| 色婷婷久久久亚洲欧美| 亚洲 欧美一区二区三区| 人妻久久中文字幕网| 国产欧美日韩综合在线一区二区| 狠狠婷婷综合久久久久久88av| 宅男免费午夜| 热99久久久久精品小说推荐| 伊人久久大香线蕉亚洲五| 激情视频va一区二区三区| 女同久久另类99精品国产91| 一二三四在线观看免费中文在| 免费在线观看视频国产中文字幕亚洲| 一区二区三区国产精品乱码| 一级毛片精品| 中文字幕av电影在线播放| 国产99久久九九免费精品| 精品少妇久久久久久888优播| 亚洲精品国产色婷婷电影| 亚洲精品国产色婷婷电影| av天堂久久9| 欧美大码av| 国产日韩一区二区三区精品不卡| 午夜成年电影在线免费观看| 丰满人妻熟妇乱又伦精品不卡| 成人三级做爰电影| 亚洲精品国产一区二区精华液| 中文字幕色久视频| 大片电影免费在线观看免费| 大片电影免费在线观看免费| 法律面前人人平等表现在哪些方面| 亚洲成国产人片在线观看| 人妻久久中文字幕网| 国产片内射在线| 欧美人与性动交α欧美精品济南到| 桃红色精品国产亚洲av| 中文欧美无线码| 美女主播在线视频| 人人妻人人澡人人看| 精品一区二区三区四区五区乱码| 女警被强在线播放| 桃花免费在线播放| 自拍欧美九色日韩亚洲蝌蚪91| aaaaa片日本免费| 亚洲精品乱久久久久久| 一级片'在线观看视频| 多毛熟女@视频| 久久人妻av系列| 看免费av毛片| 99久久99久久久精品蜜桃| 午夜福利影视在线免费观看| 国产精品影院久久| 精品久久久久久久毛片微露脸| 成人手机av| xxxhd国产人妻xxx| 老熟妇仑乱视频hdxx| 国产精品.久久久| 性少妇av在线| 亚洲成人手机| 侵犯人妻中文字幕一二三四区| 99久久人妻综合| 免费黄频网站在线观看国产| 丝袜美腿诱惑在线| 国产麻豆69| 在线观看免费午夜福利视频| 可以免费在线观看a视频的电影网站| 亚洲免费av在线视频| 在线av久久热| 99热网站在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久久国产一区二区| 正在播放国产对白刺激| 精品第一国产精品| 制服人妻中文乱码| 久久99热这里只频精品6学生| 91九色精品人成在线观看| 亚洲中文日韩欧美视频| 色婷婷久久久亚洲欧美| 精品国产超薄肉色丝袜足j| 国产亚洲av高清不卡| av线在线观看网站| 免费少妇av软件| 久久精品国产亚洲av香蕉五月 | 一本色道久久久久久精品综合| 少妇的丰满在线观看| 黑人猛操日本美女一级片| 亚洲欧美日韩高清在线视频 | 亚洲国产看品久久| 精品乱码久久久久久99久播| 午夜福利免费观看在线| 国产成人精品无人区| 青草久久国产| 日本wwww免费看| 成人18禁在线播放| 高潮久久久久久久久久久不卡| 一级毛片精品| 欧美激情久久久久久爽电影 | 在线观看66精品国产| 乱人伦中国视频| 国产真人三级小视频在线观看| 国产淫语在线视频| 人人妻人人爽人人添夜夜欢视频| 亚洲男人天堂网一区| 亚洲精品中文字幕一二三四区 | 亚洲成人手机| 老司机在亚洲福利影院| 国产欧美日韩一区二区三区在线| 乱人伦中国视频| 欧美老熟妇乱子伦牲交| 老熟女久久久| 性高湖久久久久久久久免费观看| 天天添夜夜摸| 最新的欧美精品一区二区| 操美女的视频在线观看| 夫妻午夜视频| 久久精品国产99精品国产亚洲性色 | 成年版毛片免费区| 国产精品久久久久久精品古装| 国产高清视频在线播放一区| 亚洲国产看品久久| 国产在线免费精品| 男女下面插进去视频免费观看| 国产色视频综合| 国产男女内射视频| 亚洲中文日韩欧美视频| 亚洲成人免费电影在线观看| 国产免费av片在线观看野外av| 国产欧美亚洲国产| 亚洲国产欧美一区二区综合| 一本—道久久a久久精品蜜桃钙片| 熟女少妇亚洲综合色aaa.| 欧美乱码精品一区二区三区| 久久久久久久国产电影| 午夜福利一区二区在线看| 大片电影免费在线观看免费| 99精品在免费线老司机午夜| 精品国产亚洲在线| 黑人巨大精品欧美一区二区蜜桃| 少妇被粗大的猛进出69影院| 高清黄色对白视频在线免费看| 91国产中文字幕| 久久久精品国产亚洲av高清涩受| 午夜久久久在线观看| 中文字幕av电影在线播放| 一级黄色大片毛片| 成年动漫av网址| 视频区欧美日本亚洲| 中国美女看黄片| av又黄又爽大尺度在线免费看| 久久热在线av| 免费在线观看视频国产中文字幕亚洲| 久久久精品区二区三区| 啦啦啦免费观看视频1| 国产av精品麻豆| 欧美人与性动交α欧美精品济南到| 老司机午夜福利在线观看视频 | 亚洲午夜精品一区,二区,三区| 韩国精品一区二区三区| 在线天堂中文资源库| 久久久久久久久久久久大奶| 狂野欧美激情性xxxx| 国产精品电影一区二区三区 | 欧美乱码精品一区二区三区| 中文字幕制服av| 18禁国产床啪视频网站| 1024视频免费在线观看| 国产单亲对白刺激| 国产成人精品久久二区二区免费| 精品福利观看| 在线观看一区二区三区激情| 亚洲熟妇熟女久久| 国产精品秋霞免费鲁丝片| 一级毛片女人18水好多| 国产一区二区三区视频了| 亚洲成人国产一区在线观看| 久久免费观看电影| 国产亚洲午夜精品一区二区久久| 咕卡用的链子| 女人被躁到高潮嗷嗷叫费观| videosex国产| 三级毛片av免费| 男男h啪啪无遮挡| 欧美日韩亚洲高清精品| 美女高潮喷水抽搐中文字幕| 欧美亚洲 丝袜 人妻 在线| 伦理电影免费视频| 国产精品免费视频内射| 亚洲精品久久成人aⅴ小说| 又大又爽又粗| 国产亚洲午夜精品一区二区久久| 我要看黄色一级片免费的| 亚洲少妇的诱惑av| 一级片免费观看大全| 黄色毛片三级朝国网站| 日韩 欧美 亚洲 中文字幕| 久久精品成人免费网站| 国产精品一区二区在线不卡| 午夜激情av网站| 国产精品国产av在线观看| 久久精品亚洲av国产电影网| xxxhd国产人妻xxx| 不卡一级毛片| 亚洲av美国av| 久9热在线精品视频| 两个人看的免费小视频| 亚洲精品在线美女| 国产精品二区激情视频| 后天国语完整版免费观看| 一级毛片精品| 国产精品久久久av美女十八| 久久亚洲真实| 一级毛片电影观看| 激情视频va一区二区三区| 天堂动漫精品| 免费观看a级毛片全部| 亚洲色图综合在线观看| 嫩草影视91久久| 精品国产乱码久久久久久男人| 菩萨蛮人人尽说江南好唐韦庄| 久久99热这里只频精品6学生| 热re99久久精品国产66热6| 精品少妇黑人巨大在线播放| 午夜老司机福利片| 女人爽到高潮嗷嗷叫在线视频| 久久精品国产a三级三级三级| 一本一本久久a久久精品综合妖精| 一区二区三区精品91| 国产成人一区二区三区免费视频网站| 夜夜骑夜夜射夜夜干| 国产欧美亚洲国产| 国产成人啪精品午夜网站| 国产精品98久久久久久宅男小说| 一区福利在线观看| 青草久久国产| 大片免费播放器 马上看| 丁香欧美五月| 1024视频免费在线观看| 人人澡人人妻人| 国产深夜福利视频在线观看| 老司机影院毛片| e午夜精品久久久久久久| 国产精品熟女久久久久浪| 亚洲精品久久成人aⅴ小说| 一边摸一边抽搐一进一小说 | 成年女人毛片免费观看观看9 | 99久久精品国产亚洲精品| 亚洲午夜精品一区,二区,三区| 一区二区日韩欧美中文字幕| 91国产中文字幕| 一区二区三区国产精品乱码| 国产一区有黄有色的免费视频| 国产伦人伦偷精品视频| 一级毛片电影观看| 12—13女人毛片做爰片一| 精品免费久久久久久久清纯 | 国产区一区二久久| 欧美激情久久久久久爽电影 | 国产在线精品亚洲第一网站| 在线亚洲精品国产二区图片欧美| 岛国毛片在线播放| 国产成人免费观看mmmm| 黄色成人免费大全| www.999成人在线观看| 色婷婷av一区二区三区视频| 老熟妇乱子伦视频在线观看| 欧美黄色片欧美黄色片| 一进一出抽搐动态| 亚洲专区字幕在线| 另类亚洲欧美激情| 精品国产乱码久久久久久小说| 日日摸夜夜添夜夜添小说| 人人妻,人人澡人人爽秒播| 久久99一区二区三区| 国产aⅴ精品一区二区三区波| 亚洲欧美日韩另类电影网站| 女人久久www免费人成看片| 精品久久蜜臀av无| 日韩中文字幕欧美一区二区| 日韩视频一区二区在线观看| 国产男女内射视频| 少妇粗大呻吟视频| 交换朋友夫妻互换小说| 黄色怎么调成土黄色| 精品人妻1区二区| av超薄肉色丝袜交足视频| 看免费av毛片| 天天躁狠狠躁夜夜躁狠狠躁| 老鸭窝网址在线观看| 黄色毛片三级朝国网站| 午夜福利欧美成人| 久久九九热精品免费| 久久精品国产综合久久久| 精品视频人人做人人爽| 一区在线观看完整版| 国产精品熟女久久久久浪| 深夜精品福利| 欧美黄色淫秽网站| 亚洲精品一卡2卡三卡4卡5卡| a级片在线免费高清观看视频| 成人黄色视频免费在线看| 9191精品国产免费久久| 国产精品国产av在线观看| bbb黄色大片| 久久久久网色| 国产一区二区激情短视频| 欧美国产精品一级二级三级| 丝袜美足系列| 老汉色∧v一级毛片| 男女免费视频国产| 在线观看免费高清a一片| 久热这里只有精品99| 国产精品亚洲av一区麻豆| 天天躁夜夜躁狠狠躁躁| 最近最新免费中文字幕在线| 精品少妇内射三级| 在线永久观看黄色视频| 三级毛片av免费| 国产精品1区2区在线观看. | 热re99久久精品国产66热6| 久久国产精品大桥未久av| av一本久久久久| 一个人免费看片子| 一区二区三区精品91| 成人永久免费在线观看视频 | 欧美一级毛片孕妇| 午夜两性在线视频| 香蕉丝袜av| 精品一区二区三区av网在线观看 | 男女床上黄色一级片免费看| 亚洲成a人片在线一区二区| 久久国产精品影院| 精品久久久久久久毛片微露脸| 少妇的丰满在线观看| 妹子高潮喷水视频| 久久精品国产亚洲av高清一级| 国产人伦9x9x在线观看| 嫩草影视91久久| 最新的欧美精品一区二区| 搡老熟女国产l中国老女人| 国产精品二区激情视频| 日韩中文字幕视频在线看片| 不卡一级毛片| 51午夜福利影视在线观看| 亚洲精品国产一区二区精华液| 在线观看www视频免费| 一夜夜www| 国产av国产精品国产| 午夜福利欧美成人| 无人区码免费观看不卡 | 中文欧美无线码| 50天的宝宝边吃奶边哭怎么回事| 男女午夜视频在线观看| 黄色片一级片一级黄色片| 亚洲精品成人av观看孕妇| 亚洲九九香蕉| 亚洲精品久久午夜乱码| 一区二区日韩欧美中文字幕| 操出白浆在线播放| 久热爱精品视频在线9| 久久中文字幕一级| 怎么达到女性高潮| 在线观看免费视频网站a站| 成年人黄色毛片网站| 在线 av 中文字幕| 亚洲av日韩精品久久久久久密| 夫妻午夜视频| av天堂在线播放| 制服诱惑二区| 国产精品98久久久久久宅男小说| 亚洲精品国产精品久久久不卡| 国产单亲对白刺激| 亚洲全国av大片| 嫁个100分男人电影在线观看| 亚洲av片天天在线观看| 亚洲视频免费观看视频| 久久精品国产99精品国产亚洲性色 | 久久精品国产亚洲av香蕉五月 | 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美网| 99精品欧美一区二区三区四区| 18禁裸乳无遮挡动漫免费视频| 丝袜在线中文字幕| 欧美日韩成人在线一区二区| 女同久久另类99精品国产91| av超薄肉色丝袜交足视频| 岛国在线观看网站| 女人久久www免费人成看片| 国产日韩一区二区三区精品不卡| 国产野战对白在线观看| 女警被强在线播放| 美女视频免费永久观看网站| 色精品久久人妻99蜜桃| www日本在线高清视频| 久久中文字幕人妻熟女| 中文字幕av电影在线播放| 中文字幕人妻丝袜一区二区| 成人影院久久| 国产av一区二区精品久久| 午夜激情久久久久久久| 国产在线观看jvid| 国产97色在线日韩免费| 国产日韩欧美亚洲二区| 无人区码免费观看不卡 | 久久久久久久精品吃奶| 夜夜骑夜夜射夜夜干| 亚洲av电影在线进入| 99国产精品99久久久久| 亚洲精品美女久久av网站| 中文字幕最新亚洲高清| avwww免费| 午夜激情久久久久久久| 在线观看人妻少妇| 极品人妻少妇av视频| 黄色 视频免费看| 久久九九热精品免费| 成人特级黄色片久久久久久久 | 日本vs欧美在线观看视频| 男人舔女人的私密视频| 精品熟女少妇八av免费久了| 精品一区二区三区视频在线观看免费 | 久久精品人人爽人人爽视色| 国精品久久久久久国模美| 99国产精品免费福利视频| 搡老乐熟女国产| 高清视频免费观看一区二区| 国产日韩一区二区三区精品不卡| 日韩 欧美 亚洲 中文字幕| 精品午夜福利视频在线观看一区 | 天堂动漫精品| 淫妇啪啪啪对白视频| 黄片小视频在线播放| 一级毛片电影观看| 在线观看免费视频网站a站| 精品午夜福利视频在线观看一区 | 露出奶头的视频| 久久精品国产综合久久久| 俄罗斯特黄特色一大片| 我要看黄色一级片免费的| 男女之事视频高清在线观看| 在线十欧美十亚洲十日本专区| 极品教师在线免费播放| 在线观看舔阴道视频| 丝袜美腿诱惑在线| 亚洲色图综合在线观看| 国产精品自产拍在线观看55亚洲 | 999久久久国产精品视频| 两性夫妻黄色片| 成人av一区二区三区在线看| 日本wwww免费看| 国产成人系列免费观看| 亚洲综合色网址| 久久久欧美国产精品| xxxhd国产人妻xxx| 日韩视频一区二区在线观看| 女性被躁到高潮视频| 国产成人影院久久av| 国产高清视频在线播放一区| 国产成人免费观看mmmm| 午夜福利欧美成人| 国产aⅴ精品一区二区三区波| 美女高潮喷水抽搐中文字幕| 日日爽夜夜爽网站| 满18在线观看网站| 精品国产亚洲在线| 色婷婷av一区二区三区视频| 丁香欧美五月| 正在播放国产对白刺激| 亚洲精品国产精品久久久不卡| 亚洲成人手机| 免费观看av网站的网址| 亚洲欧美精品综合一区二区三区| 天天操日日干夜夜撸| 亚洲三区欧美一区| 精品熟女少妇八av免费久了| 纯流量卡能插随身wifi吗| 性少妇av在线| 精品一品国产午夜福利视频| 欧美精品av麻豆av| 欧美激情极品国产一区二区三区| 日日摸夜夜添夜夜添小说| 免费女性裸体啪啪无遮挡网站| 免费在线观看视频国产中文字幕亚洲|