《電子技術應用》
您所在的位置:首頁 > 模擬設計 > 業界動態 > 基于 FPGA 的任意波形發生器+低通濾波器系統設計

基于 FPGA 的任意波形發生器+低通濾波器系統設計

2022-10-04
來源:FPGA技術江湖
關鍵詞: FPGA Artix-7 EGO1

  第一部分 設計概述 /Design Introduction

  1.1 設計目的

  本次設計包括基于FPGA的任意波形發生器設計實現和基于FPGA的低通濾波器設計實現。

  波形發生器是一種常見的信號源,能夠產生多種標準信號和用戶定義信號,并保證較高精度和較高穩定性,廣泛地應用于電子電路、自動控制系統和數字實驗等,諸如電話、電視、收音機、高校通信系統實驗等領域都需要用到波形發生器。本次設計基于Verilog HDL語言使用Artix-7系列芯片在EGO1開發板上實現一個簡易的任意波形發生器。該波形發生器能夠產生正弦波、方波、三角波與鋸齒波四種波形,可輸出頻率有:1MHz、100KHz、10KHz、1KHz、100Hz、10Hz、1Hz。波形類型選擇控制字與頻率設置控制字通過串口輸入。

  低通濾波可以簡單的認為:設定一個頻率點,當信號頻率高于這個頻率時不能通過,在數字信號中,這個頻率點也就是截止頻率,當頻域高于這個截止頻率時,則全部賦值為0。因為在這一處理過程中,讓低頻信號全部通過,所以稱為低通濾波。低通過濾的概念存在于各種不同的領域,諸如電子電路,數據平滑,聲學阻擋,圖像模糊等領域經常會用到。在數字圖像處理領域,從頻域看,低通濾波可以對圖像進行平滑去噪處理。

  本次設計將用兩種方式實現低通濾波器。

  方法一:利用Vivado自身具備的DDS和FIR的IP核實現;

  方法二:通過Verilog編程實現FIR的功能。

  1.2 作品展示

  

2.png

  圖1  UART串口實現

  

3.png

  圖2 波形發生器仿真

 

4.png

  圖3 低通濾波器仿真

  第二部分 系統組成及功能說明 /System Construction & Function Description

  2.1 系統概述

  任意波形發生器:計劃產生正弦波、方波、三角波與鋸齒波四種波形,可輸出頻率有:1MHz、100KHz、10KHz、1KHz、100Hz、10Hz、1Hz,利用URAT串口與串口調試助手相連進行控制字數據傳輸。已經完成全部功能,由于沒有示波器,實物調試部分不完整。

  低通濾波器:本次設計將用兩種方式實現低通濾波器。

  方法一:利用Vivado自身具備的DDS和FIR的IP核實現;

  方法二:通過Verilog編程實現FIR的功能。

  方法一使用Vivado的DDS IP核生成兩個正弦信號,頻率分別為4MHz和5MHz。同時調用乘法器IP將兩個信號進行混頻處理,乘法器將混頻信號交給FIR IP核進行過濾處理。使用Matlab進行濾波器的參數設計,這里使用窗函數法設計FIR濾波器。窗函數法設計濾波器的基本思想是:選擇一種合適的理想頻率選擇性濾波器,然后將它的沖激響應截短以獲得一個具有線性相位和因果的FIR濾波器,因此這種方法的關鍵是選區某種較好的窗函數,在濾波器性能給定的情況下,盡量選擇主瓣窄,旁瓣峰值小的窗函數。設置好以后將FIR濾波器系數導出,將生成好coe文件導入Vivado中FIR IP核中。撰寫TestBench文件進行波形仿真。

  方法二利用DDS 生成一個由三個正弦波疊加的待濾波信號,兩個正弦波的頻率分別是200Hz、800Hz。用MATLAB產生頻率為500Hz的低通濾波器抽頭系數,然后將待濾波信號送入Vivado,用Verilog編寫的FIR濾波器進行仿真,觀察濾波后的波形。

  2.2 項目系統框圖

  任意波形發生器:

 

5.png

  圖4 任意波形發生器系統框圖

  

6.png

  圖5 串口模塊設計框圖

  

7.png

  圖6 DDS原理框圖

  低通濾波器:

 

8.png

  圖7 低通濾波器設計框圖

  

9.png

  圖8 低通濾波器原理框圖

  2.3 技術說明

  (1)任意波形發生器

  本次設計使用到的是UART串口通信,主要用來接收串口數據發送的數據,字節傳輸形式:1-0-x-x-x-x-x-x-x-x-0-1,檢測到下降沿即-1-0-(起始標志)后開始進行接收字節,將接收到的一一存入寄存器中,8個bit后結束接收,就是簡單的1個字節的接收。其協議發送原理:

  

10.png

  圖9 UART數據幀格式

  根據URAT串口的數據幀格式,起始位為0,所以我們需要用一個下降沿檢測模塊,來判斷起始位,然后開始準備接收串口數據。

  然后使用一個計數器模塊,在1Bit數據傳輸的時間找到最佳時間來輸出讀取串口數據。在串口通信中波特率表示每秒能傳輸的數據位,9600就是每秒能傳輸9600位,1word=2byte、1byte=8bit,這里的9600位指的是9600bit。每個開發板上都有它的系統時鐘,本次設計用的是EGO1開發板,系統時鐘是100MHZ,我們假定波特率是9600bit/s傳輸,那么傳輸1bit就需要100M/9600個時鐘周期。那么在本次設計中傳輸1bit所需的計數值為10416,將采集數據的時刻放在波特率計數器每次循環技術的中間位置,即5208。

  利用一個控制器模塊,在下降沿檢驗通過后的信號與計數器模塊輸出的讀數據時間信號控制下,依次將8位數據放入寄存器之中,然后利用case語句,經過停止位數據后停止數據接收。

  最后在串口的總模塊中對上述三個模塊進行例化,將對應輸入輸出相連。

  下面根據下圖所示的DDS技術原理框圖,當一個參考時鐘過來 傳輸一個頻率控制字進入累加器在累加器中加入上一個參考時鐘的數據,相加后進去寄存器,然后寄存器的輸出便是正弦查表的地址,通過這個數據在波形存儲區查表,再通過D/A轉換器和低通濾波器輸出。結合項目功能要求實際分模塊實現任意波形發生器。

  波形類型控制字與頻率設置控制字是分別拿串口助手輸入的有效數據位的[7:6]與[5:2]位,我們注意到,串口傳輸過程是從低字節開始傳輸的,也就是說我們在串口助手中發送的有效數據為11010011,接收到的數據為11001011,所以需要將接收到的數據倒序然后給到兩個控制字的寄存器。

  

11.png

  圖10 DDS直接數字頻率合成技術原理

  首先我們利用MATLAB新建腳本(。m文件),鍵入相應的程序,產生了存放有位寬為8,位深為256的四種波形數據的COE文件。然后在Vivado中利用IP catalog中的Block Memory Generator核。分別新建四個Block Memory Generator核,設置Basic欄中的Memory Type為Single Port ROM,Port A Options 欄中的位寬為8,位深為256,最后在Other Options中,選擇加載本地文件選中MATLAB已經生成的正弦波、方波、鋸齒波、三角波的COE數據文件,分別將它們放入了不同的Rom中,我們就可以通過Vivado進行編程,在時鐘的驅動下,根據地址去讀取Rom中的數據然后輸出,即可產生波形數據。

  我們拿到URAT串口接收的頻率控制字后,利用case語句,根據不同的輸入,輸出不同的頻率控制字,關于頻率控制字的計算,首先系統的時鐘為100MHz,相位累加器為32位,232等于4294967296,想要輸出頻率為1kHz的信號為例,通過以下的計算公式,得到K約等于42950。

  

12.png

  然后在相位累加器中以頻率控制字為步長進行周期的加法運算,再將相位累加器的每個值的高位作為地址去讀Rom的數據。這里輸出的為8位的DAC數據,所以用相位累加器的高八位作為Rom的地址輸入。

 ?。?)低通濾波器

  低通濾波器實現:

  方法一:利用Vivado自身具備的DDS和FIR的IP核實現

  在Vivado中新建工程,創建Block Design,用DDS IP核生成兩個正弦波信號,如圖11所示。

  

13.png

  圖11 DDS IP核

  將DDS的參數設置為圖12所示。

  

14.png

  圖12 DDS參數設置

 ?。?)調用乘法器,將兩個正弦波信號進行混頻處理,乘法器設置如下圖13

  

15.png

  圖13 乘法器IP核參數設置

  (3)調用FIR IP核,對其參數進行設置。利用Matlab對濾波器的參數進行設置。在matlab窗口輸入fdatool,在彈出的窗口中,在Response Type(濾波器響應類型)一欄中選擇Lowpass(低通),Desgned Method(實現方法)選擇Window(窗函數實現),即基于窗函數法設計FIR濾波器,我們在Options中選擇Hamming窗,Filter Order中選擇Special Order,濾波器階數選擇128(階數越高,濾波效果越好,占用的芯片資源會更多些)Frequency  Specializations中單位Units選擇MHz,采樣頻率Fs選擇100MHz(根據奈奎斯特低通信號采樣定理,采樣頻率至少高于信號的2倍,實際工程上一般得高于低通信號頻率的5-10倍,如果信號頻率過高的話,實際的AD采樣頻率核精度無法達到,所以工程上一般把低通信號經調制后搬到頻率很高的載波上,再下變頻到中頻上用帶通采樣定理進行AD采樣),之前提到采樣頻率高于信號頻率的5-10倍就可以,因為我們這里FIR濾波器時鐘是和板卡時鐘連接在一起的,沒有編寫時鐘分頻模塊,所以這里采樣頻率直接設置為100MHz(一個周期點數為100MHz/2MHz=50,板卡資源會浪費一些,但是不影響),截止頻率Fc設為2MHz(1MHz的正弦波通過,9MHz的被濾除),點擊View和Design Filter按鈕,設計的濾波器時域和頻域波形(Time Domain and Frequency domain)以及幅相響應(Magnitude and Phase Response),最后生成coe文件,加載到FIR IP核中。具體設置如下圖。

 

16.png

  圖14 濾波器參數設置

  (4)將各IP核之間進行連線,生成頂層模塊。如圖15。

  

17.png

  圖15 各IP核連線圖

  方法二:通過Verilog編程實現FIR的功能。具體設計流程如圖16。

  

18.png

  圖16 設計流程圖

  (1)利用DDS 生成一個由三個正弦波疊加的待濾波信號,兩個正弦波的頻率分別是200Hz、800Hz。

 ?。?)用MATLAB產生頻率為500Hz的低通濾波器抽頭系數,

  (3)Vivado進行FIR的verilog設計。FIR的代碼實現流程如下圖7所示。

  

19.png

  圖17 FIR verilog設計框圖

  第一級流水線的功能是將輸入信號進行延時,每到來一個時鐘信號,便將輸入信號保存到delay_pipelin1中,然后將剩下的依次移動一位。

  第二級流水線的功能是將輸入經過延時的信號和濾波器系數相乘,每到來一個時鐘便將一個新的乘積結果更新到multi_data中。

  累加積乘的主要功能是將乘積累加,累加的結果就是濾波后的信號。

 ?。?)仿真測試,創建TestBench文件進行仿真測試。

  第三部分 完成情況及性能參數 /Final Design & Performance Parameters

  任意波形發生器

  已實現功能:產生正弦波、方波、三角波與鋸齒波四種波形,可輸出頻率有:1MHz、100KHz、10KHz、1KHz、100Hz、10Hz、1Hz;利用URAT串口與串口調試助手相連進行控制字數據傳輸。由于沒有示波器,實物調試部分不完整。

  a.串口模塊仿真分析:在Testbench中依次給串口數據(data_rx)為10 11010011 011。

  

20.png

  圖18  URAT串口仿真波形

  如上圖所示,FPGA接收到的數據為11001011,根據仿真結果,URAT串口測試成功。

  b.任意波形發生模塊仿真分析:

  

21.png

  圖19 波形發生模塊仿真1

 

22.png

  圖20 波形發生模塊仿真2

  如上圖所示,在仿真中分別給波形類型控制字為0、1、2、3,對應頻率設置控制字為5、6、3、2。分別產生了頻率為10kHz的正弦波、100kHz的三角波、100Hz的三角波、10Hz的鋸齒波。其中100kHz的三角波為方波,詢問老師后確定是系統顯示問題,其波形數據圖    所示。根據分析知波形發生器仿真測試成功。

  c.系統仿真測試與分析

  最后通過頂層模塊,例化DDS模塊與UART模塊,將相應輸入輸出連接起來,編寫仿真激勵文件,給串口助手發送值,驗證波形輸出。

  

23.png

  圖21 系統仿真測試1

  

24.png

  圖22 系統仿真測試2

  

25.png

  圖23 系統仿真測試3

  

26.png

  圖24  系統仿真測試4

  如上所示,根據仿真波形可知,當仿真輸入1011010011011,FPGA收到11001011,類型控制字11(3),頻率設置控制字0100 (4),輸出周期為 1ms(即頻率為1KHz)的 鋸齒波;

  當仿真輸入1000011111011,FPGA收到11111000,類型控制字00(0), 頻率設置控制字0111(7),輸出周期為1000ns頻率 (即頻率為1MHz )的正弦波;

  當仿真輸入1001010111011,FPGA收到11101010,類型控制字01(1),頻率設置控制字0101(5),輸出周期為100us( 即頻率為10KHz)的方波  ;

  當仿真輸入1010011011011,FPGA收到11011001,類型控制字10(2),頻率設置控制字0110(6),輸出周期為10us(即頻率為100KHz)的三角波。

  根據分析可知,系統仿真成功,接下來Run Synthesis→Run Implementation,然后在Windows的下拉選項中點擊I/O Port,根據EGO1硬件手冊上的引腳說明綁定引腳,選擇邏輯電平為LVCMOS33,然后保存即可生成引腳約束文件。最后Generate Bitstream→Open Hardware Manager→Open Target→Auto Connect→選中生成的bit文件→Program即可進行實物測試。

  

27.png

  圖25 系統實物調試

  

28.png

  圖26 系統實物調試

  將BIT文件下入EGO1開發板后,打開串口助手,發送“de”,即二進制的“11011110”,可以觀察到代表FPGA串口接收數據的LED1中的八個燈的數據(亮為“1”,滅為“0”)為11011110;發送“95”,即二進制的“10010101”時,LED1燈群展示數據為10010101。由于沒有示波器,不能看到波形,僅能通過LED2中的LED7—LED2顯示可以看出DDS波形發生器的波形類型控制字與頻率設置控制字數據都是對的。至此,正確實現任意波形發生器功能。

  低通濾波器

  低通濾波器通過IP核設計得以實現,仿真波形如下圖8所示。

  

29.png


  圖27 FIR—IP核實現仿真波形

  通過上面仿真波形可以看出4MHz和5MHz的正弦混頻信號經過FIR低通濾波器后將9MHz的高頻正弦信號濾除掉了,輸出為我們需要的1MHz低頻信號。

  

30.png

  圖27 FIR-Verilog代碼實現仿真波形

  如上圖所示,DDS產生頻率為200Hz與800Hz輸入,DDS_out為疊加信號,低通濾波器輸出波形不符合功能要求,未能定位到具體原因,仿真錯誤、系統功能未正確實現。




更多信息可以來這里獲取==>>電子技術應用-AET<<

圖片.jpg


本站內容除特別聲明的原創文章之外,轉載內容只為傳遞更多信息,并不代表本網站贊同其觀點。轉載的所有的文章、圖片、音/視頻文件等資料的版權歸版權所有權人所有。本站采用的非本站原創文章及圖片等內容無法一一聯系確認版權者。如涉及作品內容、版權和其它問題,請及時通過電子郵件或電話通知我們,以便迅速采取適當措施,避免給雙方造成不必要的經濟損失。聯系電話:010-82306118;郵箱:aet@chinaaet.com。
热re99久久精品国产66热_欧美小视频在线观看_日韩成人激情影院_庆余年2免费日韩剧观看大牛_91久久久久久国产精品_国产原创欧美精品_美女999久久久精品视频_欧美大成色www永久网站婷_国产色婷婷国产综合在线理论片a_国产精品电影在线观看_日韩精品视频在线观看网址_97在线观看免费_性欧美亚洲xxxx乳在线观看_久久精品美女视频网站_777国产偷窥盗摄精品视频_在线日韩第一页
  • <strike id="ygamy"></strike>
  • 
    
      • <del id="ygamy"></del>
        <tfoot id="ygamy"></tfoot>
          <strike id="ygamy"></strike>
          国产日韩欧美夫妻视频在线观看| 亚洲老司机av| 亚洲精品一品区二品区三品区| 一片黄亚洲嫩模| 精品成人国产在线观看男人呻吟| 亚洲精品黄网在线观看| 欧美区视频在线观看| 欧美国产亚洲视频| 久久精品女人| 国产日韩欧美电影在线观看| 免费在线看一区| 久久综合国产精品| 在线观看亚洲专区| 新片速递亚洲合集欧美合集| 亚洲伦伦在线| 91久久夜色精品国产网站| 国产日产欧美a一级在线| 欧美精品在欧美一区二区少妇| 国产精品乱人伦中文| 亚洲激情六月丁香| 一本色道久久综合亚洲91| 久久久噜噜噜久久中文字幕色伊伊| 亚洲黄色三级| 欧美日韩成人免费| 亚洲第一网站免费视频| 久久亚洲精品视频| 欧美亚洲色图校园春色| 国产精品夜夜夜一区二区三区尤| 国产精品久久久亚洲一区| 亚洲国产黄色| 国产一区二区三区在线播放免费观看| 久久激情五月婷婷| 国内精品国产成人| 曰韩精品一区二区| 欧美日韩精品国产| 亚洲精品综合| 99av国产精品欲麻豆| 亚洲美女色禁图| 国产精品自在线| 欧美在线观看一区二区三区| 欧美一区二视频| 亚洲精品孕妇| 老司机午夜精品| 欧美在线亚洲| 亚洲欧美国产精品va在线观看| 欧美午夜视频一区二区| 一区二区三区国产精华| 欧美视频一区在线| 亚洲性视频h| 午夜亚洲精品| 亚洲一区二区久久| 国产亚洲一区二区精品| 国产一区二区三区四区在线观看| 亚洲人体1000| 欧美精品99| 久久av一区二区三区亚洲| 麻豆成人综合网| 国语自产精品视频在线看抢先版结局| 欧美国产日韩免费| 国产老女人精品毛片久久| 欧美精品一区二区久久婷婷| 久久久久久9| 亚洲高清资源综合久久精品| 国产视频亚洲| 日韩视频免费看| 欧美激情按摩在线| 一区二区三区日韩精品视频| 亚洲精品自在在线观看| 狠狠入ady亚洲精品| 亚洲精品美女免费| 亚洲综合第一| 欧美日韩在线免费视频| 亚洲在线一区二区三区| 亚洲一区二区三区乱码aⅴ蜜桃女| 亚洲永久视频| 韩国av一区二区三区| 国产乱子伦一区二区三区国色天香| 午夜精品久久久久久久久久久| 亚洲免费观看在线视频| 久久精品国产96久久久香蕉| 久久久精品tv| 亚洲一区二区三区高清| 亚洲激情午夜| 亚洲人成精品久久久久| 欧美aa在线视频| 免费在线看成人av| 香蕉久久a毛片| 欧美伊人影院| 久久综合久久美利坚合众国| 亚洲三级色网| 国产精品一二三四| 最新69国产成人精品视频免费| 久久国产精品亚洲va麻豆| 亚洲人午夜精品| 欧美制服丝袜| 亚洲摸下面视频| 国产精品蜜臀在线观看| 欧美区日韩区| 性视频1819p久久| 亚洲欧美久久久| 欧美第一黄网免费网站| 在线观看不卡| 美女在线一区二区| 国产欧美三级| 亚洲二区精品| 国产精品羞羞答答| 亚洲一区一卡| 另类人畜视频在线| 激情综合视频| 亚洲人成在线观看网站高清| 欧美日韩国产在线一区| 欧美日韩国产亚洲一区| 亚洲经典视频在线观看| 亚洲国产精品成人va在线观看| 亚洲视频www| 亚洲激情一区二区| 欧美性猛交xxxx乱大交蜜桃| 欧美日韩一区二区三| 欧美激情影院| 久久综合免费视频影院| 一本色道久久综合亚洲精品婷婷| 久久久久久久波多野高潮日日| 亚洲风情亚aⅴ在线发布| 久久se精品一区精品二区| 老司机免费视频一区二区三区| 99re热这里只有精品免费视频| 欧美日本高清| 欧美日韩国产bt| 国产婷婷色一区二区三区在线| 影音先锋另类| 欧美区亚洲区| 欧美黄色视屏| 一区二区三区免费观看| 欧美日韩福利视频| 午夜精品久久久久久久| 亚洲欧洲另类| 欧美日韩久久不卡| 黄色成人av网站| 伊人男人综合视频网| 欧美视频一区二区三区四区| 亚洲第一黄色| 久久嫩草精品久久久精品| 欧美在线一二三| 国产精品成人一区二区网站软件| 欧美久久久久久蜜桃| 亚洲大胆人体视频| 欧美三级特黄| 亚洲欧洲日夜超级视频| 欧美专区在线播放| 国产精品久久久久久久久久三级| 国产精品中文字幕欧美| 在线观看视频一区二区| 在线观看视频一区二区欧美日韩| 亚洲国产欧美一区| 免费一级欧美在线大片| 亚洲在线观看视频网站| 国产一区导航| 美女黄网久久| 欧美人与禽性xxxxx杂性| 欧美电影免费观看大全| 国产一区亚洲| 久久精品国产综合| 先锋影院在线亚洲| 久久国产日本精品| 一本久久知道综合久久| 国产亚洲欧美一区| 美女视频一区免费观看| 国产精品一区二区在线观看| 国产亚洲一本大道中文在线| 亚洲欧洲美洲综合色网| 国产精品jvid在线观看蜜臀| 欧美日韩亚洲高清| 亚洲国产精品99久久久久久久久| 国语自产精品视频在线看8查询8| 先锋资源久久| 欧美激情视频网站| 99精品99久久久久久宅男| 国产精品国产三级国产aⅴ浪潮| 亚洲乱码日产精品bd| 美女任你摸久久| 国产美女扒开尿口久久久| 午夜欧美大尺度福利影院在线看| 艳女tv在线观看国产一区| 亚洲国产精品va在线观看黑人| 欧美天堂在线观看| 欧美尤物一区| 欧美性猛交99久久久久99按摩| 亚洲精品在线三区| 亚洲第一中文字幕在线观看| 国产精品成人一区二区网站软件| 久久亚洲色图| 欧美视频一区在线| 亚洲激情在线观看| 欧美专区在线观看一区| 亚洲精品在线免费观看视频| 免费欧美电影| 亚洲男女自偷自拍图片另类| 国产午夜精品理论片a级大结局| 亚洲天堂第二页| 亚洲精品欧美日韩| 亚洲伊人一本大道中文字幕| 欧美日韩美女一区二区| 亚洲午夜羞羞片| 久久av最新网址| 国产欧美综合一区二区三区| 国产精品一区二区黑丝| 欧美伊人久久| 亚洲欧美激情视频| 国产精品女人毛片| 亚洲欧美激情视频在线观看一区二区三区| 国产精品久久久久9999高清| 欧美电影在线观看完整版| 一区二区电影免费观看| 欧美在线综合视频| 国产精品美女久久久| 正在播放日韩| 欧美日本簧片| 亚洲国产精品一区二区www| 日韩视频中午一区| 国产日韩在线看片| 久久综合图片| 亚洲欧美在线网| 一区二区三区日韩在线观看| 亚洲一区二区三区乱码aⅴ| 日韩一二在线观看| 99re亚洲国产精品| 你懂的一区二区| 欧美天天视频| 国产精品一二三四区| 91久久在线观看| 欧美aaaaaaaa牛牛影院| 亚洲一区二区久久| 在线成人免费观看| 亚洲欧洲日本国产| 亚洲国产成人porn| 欧美日韩亚洲一区二| 亚洲精品一区二区网址| 久久不射网站| 99视频精品| 欧美日韩和欧美的一区二区| 国产精品成人观看视频国产奇米| 日韩视频第一页| 久久久精品午夜少妇| 国产一区二区剧情av在线| 国产一区二区三区四区老人| 国产精品女主播| 国产农村妇女精品一二区| 国产精品成人av性教育| 国产精品美女| 国产亚洲一区在线| 亚洲视频一区二区| 国产精品福利在线| 欧美亚洲在线观看| 中文一区二区在线观看| 欧美日韩在线播放三区四区| 国产精品裸体一区二区三区| 国产综合久久久久久| 国产精品自拍三区| 亚洲韩国青草视频| 一本大道久久精品懂色aⅴ| 久久久久久有精品国产| 国产一区二区三区在线观看免费视频| 欧美日韩精品二区| 午夜亚洲视频| 日韩亚洲欧美成人一区| 国产日韩欧美一区| 国产欧美日韩精品在线| 欧美精品www| 欧美新色视频| 夜夜嗨av一区二区三区网页| 久久精品一区二区三区四区| 国产日韩专区在线| 亚洲已满18点击进入久久| 99国产精品99久久久久久| 国产乱码精品一区二区三区av| 亚洲精品美女在线观看| 国产精品一二三四| 91久久久精品| 日韩亚洲在线观看| 亚洲一区二区毛片| 亚洲在线免费视频| 91久久久久久国产精品| 一区二区在线不卡| 亚洲黄色免费网站| 国产乱理伦片在线观看夜一区| 国产精品国产精品国产专区不蜜| 亚洲国产精品一区二区www在线| 午夜国产精品视频免费体验区| 美女国内精品自产拍在线播放| 麻豆精品一区二区av白丝在线| 国产在线不卡精品| 亚洲高清视频在线观看| 亚洲国产精品久久久久秋霞蜜臀| 夜夜嗨av一区二区三区四区| 国产欧美在线观看| 亚洲小说区图片区| 国产精品家庭影院| 欧美亚洲成人网| …久久精品99久久香蕉国产| 国产免费亚洲高清| 日韩亚洲欧美精品| 蜜臀久久99精品久久久画质超高清| 亚洲高清在线观看| 午夜国产不卡在线观看视频| 亚洲电影中文字幕| 精品福利免费观看| 亚洲一区二区三区四区中文| 男人插女人欧美| 久久精品91久久香蕉加勒比| 亚洲一级高清| 国产精品乱码一区二区三区| 午夜精品久久久久久久白皮肤| 亚洲一区二区视频在线观看| 亚洲激情校园春色| 欧美国产三级| 亚洲一区制服诱惑| 亚洲美女黄色| 一本色道久久综合亚洲精品小说| 国产欧亚日韩视频| 久久久久久久一区二区三区| 国产欧美日韩另类视频免费观看| 国产香蕉久久精品综合网| 性娇小13――14欧美| 久久9热精品视频| 亚洲一区二区四区| 国产伦精品一区二区三区照片91| 久久这里只有精品视频首页| 国产美女诱惑一区二区|