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

    基于改進的DA算法實現(xiàn)高階FIR濾波器

    2014-03-05 17:56王哲周前能
    現(xiàn)代電子技術(shù) 2014年4期
    關(guān)鍵詞:乘法器桌面消耗

    王哲+周前能

    摘 要: 介紹了幾種應(yīng)用于高階FIR濾波器的算法,著重介紹DA算法,該算法主要基于LUT形式,從根本上減少了乘法器資源的消耗,提高了系統(tǒng)的實時性。并以此為基礎(chǔ)介紹了幾種改進的算法,分別是基于LUT分解、基于OBC編碼以及基于雙向選擇器的DA改進算法,這些算法都從不同程度上降低了資源的使用,提高了信號處理速度,降低了系統(tǒng)復(fù)雜性。最后,對FIR濾波器的改進算法進行了全面的分析比較和仿真,仿真結(jié)果均正確并且符合設(shè)計要求。

    關(guān)鍵字: DA算法; FIR濾波器; LUT分解; OBC編碼; 雙向選擇器

    中圖分類號: TN911?34 文獻標識碼: A 文章編號: 1004?373X(2014)04?0008?05

    Realization of high?order FIR filter based on improved DA algorithm

    WANG Zhe, ZHOU Qian?neng

    (College of Photoelectric engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China)

    Abstract: Several kinds of algorithms applied to high?order FIR filter are introduced, among which the DA algorithm is emphatically introduced. The algorithm mainly based on the form of LUT reduced the resource consumption of multiplier radically and improved the real?time performance of the system. On this basis, several kinds of improved algorithms are presented in this paper. They are the improved DA algorithms, respectively based on LUT decomposition, OBC coding and two?way selector. They improved the speed of signal processing and reduced the system complexity. The comprehensive analysis, comparison and simulation for the improved DA algorithms applied to FIR filter are carried out. The results of simulation are correct and meet the design requirements.

    Keywords: DA algorithm; FIR filter; LUT decomposition; OBC coding; two?way selector

    數(shù)字濾波器是很多數(shù)字信號處理系統(tǒng)的基本單元,它們廣泛的應(yīng)用于通信、信號處理、圖像處理、模式識別等很多領(lǐng)域。目前有多種實現(xiàn)FIR濾波器的方法,諸如使用DSP芯片、專用集成電路ASIC以及FPGA等,利用DSP來實現(xiàn)將會受限于運行速度以及較漫長的開發(fā)時間,ASIC將會受限于開發(fā)周期、成本以及通用性等問題,而FPGA的實現(xiàn)將會平衡這些因素,由于其開發(fā)周期短、實現(xiàn)方式容易、通用性強等特點,在很多領(lǐng)域都有著廣泛的應(yīng)用,并且發(fā)展?jié)摿薮蟆?/p>

    FIR濾波器的主要的設(shè)計方法就是利用乘累加來實現(xiàn)。這個方法的實現(xiàn)主要有兩種方式,一種是直接乘加的結(jié)構(gòu),N階的FIR要完成N次乘法和N-1次加法操作,這種結(jié)構(gòu)將會消耗大量的硬件資源和面積;另一種則是利用DA算法來實現(xiàn),它將乘累加的計算形式轉(zhuǎn)換為對查找表的尋址以及移位求和來實現(xiàn),這種解決辦法基于FPGA的查找表結(jié)構(gòu)并且可以顯著提高系統(tǒng)的運行效率,但是在減少乘法器使用的同時會增加大量LUT的使用,同樣是一種面積換速度思想的體現(xiàn)[1]。

    為了達到更高的運行速度,優(yōu)化系統(tǒng)資源消耗,目前對FIR濾波器算法以及DA算法的改進有很多方法,包括部分串行結(jié)構(gòu)、分割查找表、OBC編碼、多路選擇分解、CSD編碼等。此外對于FIR濾波的運行效率和資源優(yōu)化還可以通過使用流水線、狀態(tài)控制機等方式來滿足。本文分別介紹了DA算法、基于查找表分割的DA算法、基于OBC編碼的DA算法、基于多路選擇的分解算法,并在ISE中進行Verilog硬件代碼的編寫,并編寫Testbench進行波形仿真,最后與傳統(tǒng)的FIR濾波器算法進行了性能、資源以及實時性處理的比較和分析。

    1 DA算法簡介

    分布式算法是一種應(yīng)用廣泛的FIR濾波器設(shè)計方法,主要基于查找表尋址來實現(xiàn)乘累加。分布式算法在完成乘加功能時是通過將各輸入數(shù)據(jù)每一比特位對LUT進行尋址,LUT本質(zhì)上就是一個RAM。目前FPGA中多使用4輸入的LUT,即為四位地址線,所以每一個LUT可以看成一個有4位地址線的16×1的RAM。如表1所示,當用戶通過原理圖或HDL語言描述了一個邏輯電路以后,并計算好邏輯電路所有可能的結(jié)果,并把結(jié)果事先寫入RAM,這樣,每輸入一個信號進行邏輯運算就等于輸入一個地址進行查表,找出地址對應(yīng)的內(nèi)容,然后輸出即可。

    圖1 實際邏輯與LUT實現(xiàn)方式對比

    因此,LUT中需要預(yù)先存儲FIR所有系數(shù)的所有可能組合之和,單個LUT的結(jié)構(gòu)則需要構(gòu)建[2N,]N為系數(shù)個數(shù)。每一次輸入數(shù)據(jù)得到部分積后,再對各部分積進行累加得到最終結(jié)果。設(shè)[wi]是濾波器設(shè)計的系數(shù),[xi(n)]是變量,可以看作是n時刻的第i個采樣輸入數(shù)據(jù),y(n)代表n時刻的系統(tǒng)響應(yīng)。那么它們的內(nèi)積為:

    [y(n)=i=0N-1wixi(n)] (1)

    用二進制補碼來表示輸入數(shù)據(jù)[xi(n)]:

    [xi=-xi,m-1+j=0m-1xi,m-1-j2-j] (2)

    式中[xi,m-1]是[xi]的最高位,也就是符號位。將式(2)代入式(1)可得:

    [y(n)=i=0N-1wi[-xi,m-1+j=0m-1xi,m-1-j2-j] =i=0N-1j=0m-1wixi,m-1-j2-j-i=0N-1wixi,m-1] (3)

    因為[xi,m-1-j]的值只有兩種可能性:1和0,所以[i=0N-1wixi,m-1-j]的計算結(jié)果有[2N]個不同的結(jié)果,單個LUT存儲下所有[i=0N-1wixi,m-1-j]的可能結(jié)果,并且利用來自位移存儲器[xi,m-1-j]作為LUT的地址信號從LUT中取出每一個比特位與系數(shù)的乘積。因此y(n)的值可以通過不斷的移位和累加操作得到[2]。四個系數(shù)的基于LUT的DA算法結(jié)構(gòu)如圖2所示。

    圖2 DA算法基本結(jié)構(gòu)

    以上簡要介紹了DA算法的基本原理,就是實現(xiàn)方式而言,分布式算法又可以分為并行方式、串行方式和串并結(jié)合的方式。并行方式首先要建立一個n×m的輸入緩沖陣列,將數(shù)據(jù)讀入,并對其進行排列,用所有n個輸入數(shù)據(jù)的相同位對分布式查找表進行尋址,全并行方式實現(xiàn)需要消耗大量的分布式查找表,將會消耗很多的系統(tǒng)資源,但是可以很容易采用流水線方式來實現(xiàn),所以可以得到很高的運算速度,它是一種資源換取效率的方法[3]。如圖3所示。

    圖3 并行方式實現(xiàn)DA算法

    串行方式只使用一個分布式查找表,先從最低位開始處理,用所有的n個輸入變量的最低位對分布式查找表進行尋址,得到部分和后右移一位,放入寄存器中,與此同時,n個輸入變量的次低位已經(jīng)開始對分布式查找表進行二次尋址,等得到另一個部分和,并將其同樣右移一位,與前一個部分和相加,重復(fù)以上步驟,直至得到最終的和,并用前面累加的結(jié)果減去符號位的部分和。與全并行方式相比,速度降低了很多,但是消耗的系統(tǒng)資源較少,設(shè)計思想簡單,是一種效率換資源的方法。圖2所示的即為串行形式。串并結(jié)合的方式則是對前面兩種方式的平衡。

    2 DA算法的改進

    當濾波器階數(shù)較小的時候,DA算法的查找表規(guī)模對資源的消耗可以接受,但是隨著濾波器階數(shù)的增加,查找表的規(guī)模會隨著濾波器階數(shù)的增大而呈指數(shù)增長,如果還是使用傳統(tǒng)的DA算法來實現(xiàn)FIR,其資源將會消耗巨大,所以在高階濾波器的設(shè)計中必須對DA算法進行改進。

    為了縮小LUT,對LUT的地址線的數(shù)量進行分割,可以將規(guī)模較大的LUT分割為一些規(guī)模較小的LUT,由于FIR濾波器是線性濾波器,因此低階濾波器的輸出可以相加,并且由此實現(xiàn)一個高階濾波器的輸出。如果再額外的加上流水線寄存器,這一結(jié)構(gòu)的改變不僅不會降低運行速度,而且可以極大的減小設(shè)計規(guī)模。

    如圖4所示,實現(xiàn)一個4N的DA設(shè)計需要3個輔助的加法器,LUT表格的規(guī)模從1個[B×24N]的查找表減低到4個[B×2N]的查找表,查找表規(guī)模[4]減小到了原來的[14]。圖4即為對查找表的分割,每一級加法器前還可以加流水線,縮短延時路徑,提高時鐘頻率和數(shù)據(jù)吞吐率。

    圖4 基于查找表分割的DA算法

    OBC編碼即偏移二進制編碼,基本原理即在分布算法中,將輸入數(shù)據(jù)通過OBC編碼由{0,1}映射到{1,-1},使得ROM表的上下兩部分具有鏡像對稱性關(guān)系,利用這種對稱性可以將ROM表的大小從原來的壓[2N]縮到[2N-1],即壓縮一半[5]。下面來推導(dǎo)基于OBC編碼的DA算法:輸入數(shù)據(jù)可以表示為 :

    [xi=12[xi-(-xi)]] (4)

    式中:[-xi]可以表示為式(5):

    [-x1,m-1+j=0m-1x1,m-1-j2-j+2-(m-1)] (5)

    將式(5)和式(2)代入式(4)可得:[xi=12[xi-(-xi)] =12[-xi,m-1+j=0m-1xi,m-1-j2-j+xi,m-1-j=0m-1xi,m-1-j2-j-2-(m-1)] =12[-(xi,m-1-xi,m-1)+j=0m-1(xi,m-1-j-xi,m-1-j)2-j-2-(m-1)]] (6)

    為了讓式子更簡潔,將兩部分整合為一個式子,定義[di,j],其取值范圍為[-1,1]。

    [di,j=xi,j-xi,j, j≠m-1-(xi,m-1-xi,m-1), j=m-1] (7)

    所以得到:

    [xi=12j=0m-1di,m-1-j2-j+2-(m-1)] (8)

    所以:

    [y=j=0m-1i=0N-112widi,m-1-j2-j-12i=0N-1wi2-(m-1)] (9)

    為了讓式子更簡潔,再次定義[Cj=i=0N-112widi,j],[Cex=i=0N-112wi],所以輸出數(shù)據(jù)可以寫成:

    [y=j=0m-1Cj2-j-Cex2-(m-1)] (10)

    基于OBC查找表的LUT表如圖5所示,優(yōu)化后的OBC查找表如圖6所示。

    圖5 基于OBC查找表的LUT表

    圖6 優(yōu)化后的OBC查找表

    可以觀察到LUT的系數(shù)組合成鏡像,可以通過對[x0j]的控制來減小LUT的規(guī)模,當[x0j]=0時,其他階數(shù)相同比特位上的數(shù)據(jù)與其異或后的值不變;當[x0j=1]時,其他階數(shù)相同比特位上的數(shù)據(jù)異或后的值將取反[6]。S1為最高位標志位,當S1=1時,LUT的輸出取反,當S1=0時,LUT表的值不變。所以LUT表的最終輸出由[x0j]和S1共同決定;S2雙向選擇器控制著輸入移位寄存器的方向和[Cex]。上述原理結(jié)構(gòu)如圖7所示。

    最后要介紹的方法的核心思想就是使用多個雙向選擇器,利用系數(shù)的對稱性逐步代替掉LUT。該結(jié)構(gòu)的推導(dǎo)如式(11)所示:

    [ LUT(bk-1,bk-2,…,b1,b0)=LUT(bk-2,…,b1,b0)+wk-1,0= LUT(bk-3,...,b1,b0)+wk-1,0+wk-2,0=…= wk-1,0+wk-2,0+…+w0,0]

    從式(11)可以看出,該算法就是把寄存器的輸出的最高位不斷反復(fù)的從LUT分離出來,由一個雙向選擇器和一個全加器來代替這部分LUT,從而達到減小LUT的目的,如圖8所示;然后再次通過使用選擇開關(guān)逐步替換掉次高位,再次減小LUT的容量,在對LUT進行反復(fù)的簡化之后,原來LUT結(jié)構(gòu)被雙向選擇器和加法器完全代替掉,最后得到不包含LUT的DA算法結(jié)構(gòu)[7],如圖9所示。

    圖7 基于OBC編碼的DA算法的實現(xiàn)

    圖8 最高位被替換掉的LUT

    3 各算法實現(xiàn)比較和性能分析

    根據(jù)項目需求,該濾波器主要應(yīng)用于基站中頻處理模塊中,F(xiàn)IR濾波器放置于低采樣率的一端,并完成多個窄帶子載波的整形濾波,用于補償前端濾波器濾波特性的不足,同時完成一定倍數(shù)的插值(抽?。V波。由于該濾波器設(shè)計抽取插值倍數(shù)很高,通帶頻率也只有十幾kHz,此外,該模塊用于多個子載波的處理,因此對資源占用的要求也非??量?。綜上所述,為了滿足通帶的性能指標以及盡量減少資源消耗的要求,本文需要設(shè)計一個基于DA算法的120階FIR濾波器才能滿足參數(shù)要求。首先使用FDATool導(dǎo)出120階,定點數(shù)據(jù)類型為(16,-2,18)的FIR濾波器系數(shù),并全部轉(zhuǎn)換為16進制數(shù),使用Verilog進行代碼編寫,在ISE環(huán)境下分別綜合仿真以上幾種算法。

    圖9 改進的最終DA結(jié)構(gòu)

    得到資源消耗如表1所示。

    表1 各算法資源消耗對比

    從表1的資源分析可以看出,并行FIR濾波器主要通過大量乘法器來并行處理數(shù)據(jù),消耗大量邏輯資源,但是處理速度最快;串行FIR濾波器復(fù)用單個乘法器來完成120次乘法運算,輸入一個數(shù)據(jù)需要8個時鐘周期,輸出一個數(shù)據(jù)需要120個時鐘周期,快慢時鐘相差120倍,但是滿足了15倍下采樣的時鐘需求,處理速度最慢但消耗LUT資源最少;部分串行結(jié)構(gòu)平衡了二者的利弊,通過對時序的控制復(fù)用了8個乘法器,但是依然需要消耗乘法器資源;DA算法完全不在使用乘法器資源,但是會消耗大量的LUT資源,處理速度處于中等水平;LUT分解將大規(guī)模的查找表分割為了多個小規(guī)模的查找表,增加了流水線,處理速度要快于傳統(tǒng)的DA算法,而且還減少了LUT的資源消耗;基于雙向選擇器的改進算法進一步減小了查找表消耗,但是卻增加了寄存器的使用,處理速度趨于傳統(tǒng)的DA算法和LUT分解改進算法之間。如圖10~圖13所示。

    在有相同控制信號的Testbench的仿真下,對比各波形傳輸延遲可以看出,各算法的處理速度大致相互關(guān)系為:

    全并行>DA改進(LUT分解)≥DA分解(雙向選擇器)>傳統(tǒng)DA算法>全串行。

    圖10 并行FIR濾波器

    圖11 串行FIR濾波器

    圖12 DA算法的FIR濾波器

    圖13 改進的DA算法的FIR濾波器

    4 結(jié) 論

    通過對FIR濾波器傳統(tǒng)算法、DA算法及其改進的算法進行比較,當FPGA乘法器資源充足的情況下及濾波器階數(shù)較低的情況下,使用傳統(tǒng)算法比較理想,實現(xiàn)起來也較為容易;而本文提到的三種改進的DA算法很適用于在高階FIR濾波器的設(shè)計情況,基于OBC編碼的DA算法可以減少一半的LUT使用;在濾波器階數(shù)較大時(16階以上),LUT規(guī)模將會呈指數(shù)增加,基于查找表分解的改進算法可以從很大程度上減小LUT的規(guī)模,對而基于選擇器的DA算法可以完全不再占用更多的LUT和乘法器,可以節(jié)省LUT的資源到達50%以上,顯著降低了系統(tǒng)的復(fù)雜程度并提高了運行的速度。三種算法在資源緊缺和濾波器階數(shù)很高的情況下還可以結(jié)合使用,再加上流水線寄存器,可以在保證資源消耗不多的情況下,整合資源利用效率,提高整個信號處理的速度,減少傳輸延遲,減小系統(tǒng)復(fù)雜度,因此改進后的算法將會優(yōu)于傳統(tǒng)算法。目前,改進DA算法的FIR濾波器已經(jīng)應(yīng)用于系統(tǒng)設(shè)計中。

    參考文獻

    [1] 田耕.無線通信FPGA設(shè)計[M].北京:電子工業(yè)出版社,2007.

    [2] 程佩青.數(shù)字信號處理教程[M].北京:清華大學(xué)出版社,2010.

    [3] 劉朋全.基于FPGA的FIR數(shù)字濾波器的設(shè)計與實現(xiàn)[D].西安:西北工業(yè)大學(xué),2006.

    [4] 孫友軍.基于FPGA的FIR數(shù)字濾波器算法研究與設(shè)計實現(xiàn)[D].北京:冶金自動化研究設(shè)計院,2010.

    [5] 顧文奕.FIR數(shù)字濾波器的優(yōu)化與驗證[J].電子測量技術(shù),2008(5):183?189.

    [6] GUO Rui, DEBRUNNER L S. A novel adaptive filter implementation scheme using distributed arithmetic [C]// Proceedings of the Forty Fifth Asilomar Conference on Signals, Systems and Computers. Pacific Grove, CA: [s.n.], 2011: 160?164.

    [7] YOO Heejong, ANDERSIN D V. Hardware?efficient distributed arithmetic architecture for high?order digital filters [C]// Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing. [S.l.]: IEEE. 2005:125?128.

    可以觀察到LUT的系數(shù)組合成鏡像,可以通過對[x0j]的控制來減小LUT的規(guī)模,當[x0j]=0時,其他階數(shù)相同比特位上的數(shù)據(jù)與其異或后的值不變;當[x0j=1]時,其他階數(shù)相同比特位上的數(shù)據(jù)異或后的值將取反[6]。S1為最高位標志位,當S1=1時,LUT的輸出取反,當S1=0時,LUT表的值不變。所以LUT表的最終輸出由[x0j]和S1共同決定;S2雙向選擇器控制著輸入移位寄存器的方向和[Cex]。上述原理結(jié)構(gòu)如圖7所示。

    最后要介紹的方法的核心思想就是使用多個雙向選擇器,利用系數(shù)的對稱性逐步代替掉LUT。該結(jié)構(gòu)的推導(dǎo)如式(11)所示:

    [ LUT(bk-1,bk-2,…,b1,b0)=LUT(bk-2,…,b1,b0)+wk-1,0= LUT(bk-3,...,b1,b0)+wk-1,0+wk-2,0=…= wk-1,0+wk-2,0+…+w0,0]

    從式(11)可以看出,該算法就是把寄存器的輸出的最高位不斷反復(fù)的從LUT分離出來,由一個雙向選擇器和一個全加器來代替這部分LUT,從而達到減小LUT的目的,如圖8所示;然后再次通過使用選擇開關(guān)逐步替換掉次高位,再次減小LUT的容量,在對LUT進行反復(fù)的簡化之后,原來LUT結(jié)構(gòu)被雙向選擇器和加法器完全代替掉,最后得到不包含LUT的DA算法結(jié)構(gòu)[7],如圖9所示。

    圖7 基于OBC編碼的DA算法的實現(xiàn)

    圖8 最高位被替換掉的LUT

    3 各算法實現(xiàn)比較和性能分析

    根據(jù)項目需求,該濾波器主要應(yīng)用于基站中頻處理模塊中,F(xiàn)IR濾波器放置于低采樣率的一端,并完成多個窄帶子載波的整形濾波,用于補償前端濾波器濾波特性的不足,同時完成一定倍數(shù)的插值(抽取)濾波。由于該濾波器設(shè)計抽取插值倍數(shù)很高,通帶頻率也只有十幾kHz,此外,該模塊用于多個子載波的處理,因此對資源占用的要求也非??量獭>C上所述,為了滿足通帶的性能指標以及盡量減少資源消耗的要求,本文需要設(shè)計一個基于DA算法的120階FIR濾波器才能滿足參數(shù)要求。首先使用FDATool導(dǎo)出120階,定點數(shù)據(jù)類型為(16,-2,18)的FIR濾波器系數(shù),并全部轉(zhuǎn)換為16進制數(shù),使用Verilog進行代碼編寫,在ISE環(huán)境下分別綜合仿真以上幾種算法。

    圖9 改進的最終DA結(jié)構(gòu)

    得到資源消耗如表1所示。

    表1 各算法資源消耗對比

    從表1的資源分析可以看出,并行FIR濾波器主要通過大量乘法器來并行處理數(shù)據(jù),消耗大量邏輯資源,但是處理速度最快;串行FIR濾波器復(fù)用單個乘法器來完成120次乘法運算,輸入一個數(shù)據(jù)需要8個時鐘周期,輸出一個數(shù)據(jù)需要120個時鐘周期,快慢時鐘相差120倍,但是滿足了15倍下采樣的時鐘需求,處理速度最慢但消耗LUT資源最少;部分串行結(jié)構(gòu)平衡了二者的利弊,通過對時序的控制復(fù)用了8個乘法器,但是依然需要消耗乘法器資源;DA算法完全不在使用乘法器資源,但是會消耗大量的LUT資源,處理速度處于中等水平;LUT分解將大規(guī)模的查找表分割為了多個小規(guī)模的查找表,增加了流水線,處理速度要快于傳統(tǒng)的DA算法,而且還減少了LUT的資源消耗;基于雙向選擇器的改進算法進一步減小了查找表消耗,但是卻增加了寄存器的使用,處理速度趨于傳統(tǒng)的DA算法和LUT分解改進算法之間。如圖10~圖13所示。

    在有相同控制信號的Testbench的仿真下,對比各波形傳輸延遲可以看出,各算法的處理速度大致相互關(guān)系為:

    全并行>DA改進(LUT分解)≥DA分解(雙向選擇器)>傳統(tǒng)DA算法>全串行。

    圖10 并行FIR濾波器

    圖11 串行FIR濾波器

    圖12 DA算法的FIR濾波器

    圖13 改進的DA算法的FIR濾波器

    4 結(jié) 論

    通過對FIR濾波器傳統(tǒng)算法、DA算法及其改進的算法進行比較,當FPGA乘法器資源充足的情況下及濾波器階數(shù)較低的情況下,使用傳統(tǒng)算法比較理想,實現(xiàn)起來也較為容易;而本文提到的三種改進的DA算法很適用于在高階FIR濾波器的設(shè)計情況,基于OBC編碼的DA算法可以減少一半的LUT使用;在濾波器階數(shù)較大時(16階以上),LUT規(guī)模將會呈指數(shù)增加,基于查找表分解的改進算法可以從很大程度上減小LUT的規(guī)模,對而基于選擇器的DA算法可以完全不再占用更多的LUT和乘法器,可以節(jié)省LUT的資源到達50%以上,顯著降低了系統(tǒng)的復(fù)雜程度并提高了運行的速度。三種算法在資源緊缺和濾波器階數(shù)很高的情況下還可以結(jié)合使用,再加上流水線寄存器,可以在保證資源消耗不多的情況下,整合資源利用效率,提高整個信號處理的速度,減少傳輸延遲,減小系統(tǒng)復(fù)雜度,因此改進后的算法將會優(yōu)于傳統(tǒng)算法。目前,改進DA算法的FIR濾波器已經(jīng)應(yīng)用于系統(tǒng)設(shè)計中。

    參考文獻

    [1] 田耕.無線通信FPGA設(shè)計[M].北京:電子工業(yè)出版社,2007.

    [2] 程佩青.數(shù)字信號處理教程[M].北京:清華大學(xué)出版社,2010.

    [3] 劉朋全.基于FPGA的FIR數(shù)字濾波器的設(shè)計與實現(xiàn)[D].西安:西北工業(yè)大學(xué),2006.

    [4] 孫友軍.基于FPGA的FIR數(shù)字濾波器算法研究與設(shè)計實現(xiàn)[D].北京:冶金自動化研究設(shè)計院,2010.

    [5] 顧文奕.FIR數(shù)字濾波器的優(yōu)化與驗證[J].電子測量技術(shù),2008(5):183?189.

    [6] GUO Rui, DEBRUNNER L S. A novel adaptive filter implementation scheme using distributed arithmetic [C]// Proceedings of the Forty Fifth Asilomar Conference on Signals, Systems and Computers. Pacific Grove, CA: [s.n.], 2011: 160?164.

    [7] YOO Heejong, ANDERSIN D V. Hardware?efficient distributed arithmetic architecture for high?order digital filters [C]// Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing. [S.l.]: IEEE. 2005:125?128.

    可以觀察到LUT的系數(shù)組合成鏡像,可以通過對[x0j]的控制來減小LUT的規(guī)模,當[x0j]=0時,其他階數(shù)相同比特位上的數(shù)據(jù)與其異或后的值不變;當[x0j=1]時,其他階數(shù)相同比特位上的數(shù)據(jù)異或后的值將取反[6]。S1為最高位標志位,當S1=1時,LUT的輸出取反,當S1=0時,LUT表的值不變。所以LUT表的最終輸出由[x0j]和S1共同決定;S2雙向選擇器控制著輸入移位寄存器的方向和[Cex]。上述原理結(jié)構(gòu)如圖7所示。

    最后要介紹的方法的核心思想就是使用多個雙向選擇器,利用系數(shù)的對稱性逐步代替掉LUT。該結(jié)構(gòu)的推導(dǎo)如式(11)所示:

    [ LUT(bk-1,bk-2,…,b1,b0)=LUT(bk-2,…,b1,b0)+wk-1,0= LUT(bk-3,...,b1,b0)+wk-1,0+wk-2,0=…= wk-1,0+wk-2,0+…+w0,0]

    從式(11)可以看出,該算法就是把寄存器的輸出的最高位不斷反復(fù)的從LUT分離出來,由一個雙向選擇器和一個全加器來代替這部分LUT,從而達到減小LUT的目的,如圖8所示;然后再次通過使用選擇開關(guān)逐步替換掉次高位,再次減小LUT的容量,在對LUT進行反復(fù)的簡化之后,原來LUT結(jié)構(gòu)被雙向選擇器和加法器完全代替掉,最后得到不包含LUT的DA算法結(jié)構(gòu)[7],如圖9所示。

    圖7 基于OBC編碼的DA算法的實現(xiàn)

    圖8 最高位被替換掉的LUT

    3 各算法實現(xiàn)比較和性能分析

    根據(jù)項目需求,該濾波器主要應(yīng)用于基站中頻處理模塊中,F(xiàn)IR濾波器放置于低采樣率的一端,并完成多個窄帶子載波的整形濾波,用于補償前端濾波器濾波特性的不足,同時完成一定倍數(shù)的插值(抽?。V波。由于該濾波器設(shè)計抽取插值倍數(shù)很高,通帶頻率也只有十幾kHz,此外,該模塊用于多個子載波的處理,因此對資源占用的要求也非??量?。綜上所述,為了滿足通帶的性能指標以及盡量減少資源消耗的要求,本文需要設(shè)計一個基于DA算法的120階FIR濾波器才能滿足參數(shù)要求。首先使用FDATool導(dǎo)出120階,定點數(shù)據(jù)類型為(16,-2,18)的FIR濾波器系數(shù),并全部轉(zhuǎn)換為16進制數(shù),使用Verilog進行代碼編寫,在ISE環(huán)境下分別綜合仿真以上幾種算法。

    圖9 改進的最終DA結(jié)構(gòu)

    得到資源消耗如表1所示。

    表1 各算法資源消耗對比

    從表1的資源分析可以看出,并行FIR濾波器主要通過大量乘法器來并行處理數(shù)據(jù),消耗大量邏輯資源,但是處理速度最快;串行FIR濾波器復(fù)用單個乘法器來完成120次乘法運算,輸入一個數(shù)據(jù)需要8個時鐘周期,輸出一個數(shù)據(jù)需要120個時鐘周期,快慢時鐘相差120倍,但是滿足了15倍下采樣的時鐘需求,處理速度最慢但消耗LUT資源最少;部分串行結(jié)構(gòu)平衡了二者的利弊,通過對時序的控制復(fù)用了8個乘法器,但是依然需要消耗乘法器資源;DA算法完全不在使用乘法器資源,但是會消耗大量的LUT資源,處理速度處于中等水平;LUT分解將大規(guī)模的查找表分割為了多個小規(guī)模的查找表,增加了流水線,處理速度要快于傳統(tǒng)的DA算法,而且還減少了LUT的資源消耗;基于雙向選擇器的改進算法進一步減小了查找表消耗,但是卻增加了寄存器的使用,處理速度趨于傳統(tǒng)的DA算法和LUT分解改進算法之間。如圖10~圖13所示。

    在有相同控制信號的Testbench的仿真下,對比各波形傳輸延遲可以看出,各算法的處理速度大致相互關(guān)系為:

    全并行>DA改進(LUT分解)≥DA分解(雙向選擇器)>傳統(tǒng)DA算法>全串行。

    圖10 并行FIR濾波器

    圖11 串行FIR濾波器

    圖12 DA算法的FIR濾波器

    圖13 改進的DA算法的FIR濾波器

    4 結(jié) 論

    通過對FIR濾波器傳統(tǒng)算法、DA算法及其改進的算法進行比較,當FPGA乘法器資源充足的情況下及濾波器階數(shù)較低的情況下,使用傳統(tǒng)算法比較理想,實現(xiàn)起來也較為容易;而本文提到的三種改進的DA算法很適用于在高階FIR濾波器的設(shè)計情況,基于OBC編碼的DA算法可以減少一半的LUT使用;在濾波器階數(shù)較大時(16階以上),LUT規(guī)模將會呈指數(shù)增加,基于查找表分解的改進算法可以從很大程度上減小LUT的規(guī)模,對而基于選擇器的DA算法可以完全不再占用更多的LUT和乘法器,可以節(jié)省LUT的資源到達50%以上,顯著降低了系統(tǒng)的復(fù)雜程度并提高了運行的速度。三種算法在資源緊缺和濾波器階數(shù)很高的情況下還可以結(jié)合使用,再加上流水線寄存器,可以在保證資源消耗不多的情況下,整合資源利用效率,提高整個信號處理的速度,減少傳輸延遲,減小系統(tǒng)復(fù)雜度,因此改進后的算法將會優(yōu)于傳統(tǒng)算法。目前,改進DA算法的FIR濾波器已經(jīng)應(yīng)用于系統(tǒng)設(shè)計中。

    參考文獻

    [1] 田耕.無線通信FPGA設(shè)計[M].北京:電子工業(yè)出版社,2007.

    [2] 程佩青.數(shù)字信號處理教程[M].北京:清華大學(xué)出版社,2010.

    [3] 劉朋全.基于FPGA的FIR數(shù)字濾波器的設(shè)計與實現(xiàn)[D].西安:西北工業(yè)大學(xué),2006.

    [4] 孫友軍.基于FPGA的FIR數(shù)字濾波器算法研究與設(shè)計實現(xiàn)[D].北京:冶金自動化研究設(shè)計院,2010.

    [5] 顧文奕.FIR數(shù)字濾波器的優(yōu)化與驗證[J].電子測量技術(shù),2008(5):183?189.

    [6] GUO Rui, DEBRUNNER L S. A novel adaptive filter implementation scheme using distributed arithmetic [C]// Proceedings of the Forty Fifth Asilomar Conference on Signals, Systems and Computers. Pacific Grove, CA: [s.n.], 2011: 160?164.

    [7] YOO Heejong, ANDERSIN D V. Hardware?efficient distributed arithmetic architecture for high?order digital filters [C]// Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing. [S.l.]: IEEE. 2005:125?128.

    猜你喜歡
    乘法器桌面消耗
    如此消耗卡路里
    玉鋼燒結(jié)降低固體燃料消耗實踐
    桌面舞龍
    一種低開銷的近似乘法器設(shè)計
    降低鋼鐵料消耗的生產(chǎn)實踐
    基于APP在線控制雙擠出頭FDM桌面3D打印機的研制
    桌面云技術(shù)在鐵路行業(yè)中的應(yīng)用
    我們消耗很多能源
    基于FPGA的流水線單精度浮點數(shù)乘法器設(shè)計*
    當灰塵厚厚地落滿了桌面
    freevideosex欧美| 精品一区二区三区视频在线| av在线app专区| 一个人看视频在线观看www免费| 制服丝袜香蕉在线| 国产精品麻豆人妻色哟哟久久| kizo精华| 久久久久久久久久久久大奶| 欧美国产精品一级二级三级| 狂野欧美激情性bbbbbb| 国产精品人妻久久久久久| 国产免费一级a男人的天堂| 国产 一区精品| 狂野欧美激情性xxxx在线观看| 欧美日韩av久久| 少妇精品久久久久久久| 永久网站在线| 国产男人的电影天堂91| 亚洲国产精品成人久久小说| 午夜福利视频在线观看免费| 亚洲成色77777| 哪个播放器可以免费观看大片| 成人18禁高潮啪啪吃奶动态图 | 国产深夜福利视频在线观看| .国产精品久久| 亚洲欧美成人综合另类久久久| 国产精品女同一区二区软件| 999精品在线视频| 国产精品不卡视频一区二区| 欧美日韩在线观看h| 青春草亚洲视频在线观看| 亚洲欧美一区二区三区国产| 日日摸夜夜添夜夜添av毛片| 亚洲熟女精品中文字幕| 国产视频首页在线观看| 热re99久久精品国产66热6| 午夜免费男女啪啪视频观看| 欧美精品人与动牲交sv欧美| 在线观看国产h片| 欧美xxxx性猛交bbbb| 在线观看免费视频网站a站| 我的老师免费观看完整版| 爱豆传媒免费全集在线观看| a 毛片基地| h视频一区二区三区| 我的老师免费观看完整版| 国产男人的电影天堂91| 性色avwww在线观看| 国产黄色视频一区二区在线观看| 欧美精品亚洲一区二区| 亚洲五月色婷婷综合| 99热6这里只有精品| 久久久久久久大尺度免费视频| 午夜福利视频在线观看免费| 亚洲精品久久成人aⅴ小说 | 亚洲精品乱久久久久久| 国产精品一区www在线观看| 青春草视频在线免费观看| 亚洲欧美一区二区三区黑人 | 又大又黄又爽视频免费| 国产在视频线精品| 国产免费又黄又爽又色| 一级,二级,三级黄色视频| 高清av免费在线| 成年人午夜在线观看视频| 中文字幕久久专区| 18在线观看网站| 久久午夜综合久久蜜桃| 麻豆乱淫一区二区| 欧美激情 高清一区二区三区| 纵有疾风起免费观看全集完整版| 成人漫画全彩无遮挡| 国产精品一二三区在线看| 色5月婷婷丁香| 中文欧美无线码| 黄片播放在线免费| 黑人巨大精品欧美一区二区蜜桃 | 国产伦精品一区二区三区视频9| 丝袜美足系列| 精品久久久噜噜| 亚洲欧美精品自产自拍| 中文字幕最新亚洲高清| 嘟嘟电影网在线观看| 久久久久国产网址| 久久久久视频综合| 一级黄片播放器| 日韩在线高清观看一区二区三区| 只有这里有精品99| 久久99热这里只频精品6学生| 日本黄色片子视频| 韩国av在线不卡| 亚洲伊人久久精品综合| 最黄视频免费看| 天天躁夜夜躁狠狠久久av| 成人亚洲精品一区在线观看| 免费观看在线日韩| 国产精品无大码| 国产男人的电影天堂91| 国产 精品1| 久久毛片免费看一区二区三区| 天天躁夜夜躁狠狠久久av| 制服人妻中文乱码| 亚洲不卡免费看| 大陆偷拍与自拍| xxx大片免费视频| 亚洲欧美中文字幕日韩二区| 欧美人与善性xxx| 久久久久精品久久久久真实原创| 女性被躁到高潮视频| 日韩一区二区视频免费看| 成人毛片a级毛片在线播放| 97在线视频观看| 麻豆乱淫一区二区| 色网站视频免费| 欧美精品国产亚洲| 国产精品久久久久久久久免| 插逼视频在线观看| 国产片内射在线| 免费看av在线观看网站| 免费观看在线日韩| 综合色丁香网| 国产一区有黄有色的免费视频| 男的添女的下面高潮视频| 国产精品一区二区在线不卡| 成人二区视频| 永久网站在线| 18禁在线播放成人免费| av在线app专区| 永久免费av网站大全| 一级毛片电影观看| 黑人巨大精品欧美一区二区蜜桃 | 一区二区三区精品91| 国产黄片视频在线免费观看| 国产视频首页在线观看| 婷婷成人精品国产| 色94色欧美一区二区| 免费不卡的大黄色大毛片视频在线观看| 丝袜脚勾引网站| 国产精品不卡视频一区二区| 七月丁香在线播放| 日本黄色片子视频| 成年美女黄网站色视频大全免费 | 免费看光身美女| 极品人妻少妇av视频| 一区二区三区乱码不卡18| 欧美97在线视频| 国产精品国产三级专区第一集| 亚洲色图综合在线观看| 一级毛片aaaaaa免费看小| 国产黄频视频在线观看| 精品少妇内射三级| 蜜桃在线观看..| av福利片在线| 欧美日韩av久久| 亚洲久久久国产精品| 国产一区二区在线观看av| 在线 av 中文字幕| 欧美激情 高清一区二区三区| 亚洲欧洲国产日韩| 一本一本综合久久| 人妻制服诱惑在线中文字幕| 免费久久久久久久精品成人欧美视频 | 高清不卡的av网站| 日韩欧美一区视频在线观看| 国产精品无大码| 美女xxoo啪啪120秒动态图| 春色校园在线视频观看| 色婷婷av一区二区三区视频| 午夜福利视频精品| 我要看黄色一级片免费的| 在线观看www视频免费| 中国美白少妇内射xxxbb| 日韩av不卡免费在线播放| av一本久久久久| av不卡在线播放| 在线观看一区二区三区激情| 在线看a的网站| 色吧在线观看| av国产精品久久久久影院| 国产精品一区www在线观看| 久热这里只有精品99| 久久国产精品大桥未久av| 菩萨蛮人人尽说江南好唐韦庄| 日韩强制内射视频| 春色校园在线视频观看| 简卡轻食公司| 色婷婷久久久亚洲欧美| 曰老女人黄片| 久久久精品免费免费高清| 成人毛片60女人毛片免费| 高清av免费在线| 亚洲欧美日韩另类电影网站| 精品久久久噜噜| 久久久久久久久久人人人人人人| 日本黄大片高清| 日韩不卡一区二区三区视频在线| 黄色怎么调成土黄色| 久久精品夜色国产| 免费看不卡的av| 母亲3免费完整高清在线观看 | 麻豆精品久久久久久蜜桃| 中国美白少妇内射xxxbb| 欧美日韩精品成人综合77777| 观看美女的网站| 欧美成人精品欧美一级黄| 老司机影院毛片| 亚洲一级一片aⅴ在线观看| 午夜老司机福利剧场| 一级爰片在线观看| 永久网站在线| 九色成人免费人妻av| 国产精品嫩草影院av在线观看| av在线老鸭窝| 国产有黄有色有爽视频| 亚洲国产av新网站| 久久 成人 亚洲| 999精品在线视频| 国产伦理片在线播放av一区| av网站免费在线观看视频| 激情五月婷婷亚洲| 久久精品久久久久久久性| 久久久欧美国产精品| 婷婷色麻豆天堂久久| 97在线视频观看| 狂野欧美白嫩少妇大欣赏| 两个人免费观看高清视频| 午夜福利网站1000一区二区三区| 国产男人的电影天堂91| a级毛片黄视频| 久久婷婷青草| 久久久午夜欧美精品| 中国国产av一级| 欧美xxxx性猛交bbbb| 中国美白少妇内射xxxbb| 狠狠婷婷综合久久久久久88av| 久久久久视频综合| 天堂8中文在线网| 成年人午夜在线观看视频| 丝瓜视频免费看黄片| 亚洲在久久综合| 最新的欧美精品一区二区| 久久人人爽人人爽人人片va| 熟女人妻精品中文字幕| 午夜福利,免费看| 免费播放大片免费观看视频在线观看| 国产高清有码在线观看视频| 日本91视频免费播放| 日韩人妻高清精品专区| 视频中文字幕在线观看| 九九爱精品视频在线观看| 国产有黄有色有爽视频| 亚洲精品乱久久久久久| 在线观看免费高清a一片| 妹子高潮喷水视频| av视频免费观看在线观看| 免费观看在线日韩| 国产国拍精品亚洲av在线观看| 99久久精品国产国产毛片| 九九爱精品视频在线观看| 亚洲,欧美,日韩| 国产极品天堂在线| 久久人人爽av亚洲精品天堂| 婷婷色av中文字幕| 国产片内射在线| 在线 av 中文字幕| 久久精品熟女亚洲av麻豆精品| 婷婷色麻豆天堂久久| 一级二级三级毛片免费看| 欧美成人精品欧美一级黄| 成人免费观看视频高清| 最黄视频免费看| 狂野欧美激情性xxxx在线观看| 国产精品国产三级国产专区5o| 免费看不卡的av| 国产成人精品福利久久| 一区在线观看完整版| 大片免费播放器 马上看| 伊人久久国产一区二区| 亚洲av二区三区四区| 3wmmmm亚洲av在线观看| xxx大片免费视频| 日韩欧美精品免费久久| 国产深夜福利视频在线观看| 国产不卡av网站在线观看| 满18在线观看网站| 男女高潮啪啪啪动态图| 嘟嘟电影网在线观看| .国产精品久久| 国产一区亚洲一区在线观看| 久久久久精品性色| 亚洲美女搞黄在线观看| 精品熟女少妇av免费看| 亚洲av电影在线观看一区二区三区| 欧美日韩亚洲高清精品| 黄色一级大片看看| 国产伦精品一区二区三区视频9| 高清欧美精品videossex| 国产高清三级在线| 免费高清在线观看日韩| 国产 一区精品| 人人妻人人添人人爽欧美一区卜| 国产精品99久久99久久久不卡 | 精品一品国产午夜福利视频| 久久99一区二区三区| 搡女人真爽免费视频火全软件| 国产精品女同一区二区软件| 麻豆精品久久久久久蜜桃| 男人添女人高潮全过程视频| 中文字幕人妻丝袜制服| 啦啦啦在线观看免费高清www| 久久国内精品自在自线图片| 国产日韩欧美亚洲二区| 久久女婷五月综合色啪小说| 少妇丰满av| 免费观看无遮挡的男女| 99热这里只有是精品在线观看| 另类亚洲欧美激情| 一级,二级,三级黄色视频| 一级a做视频免费观看| 国产精品久久久久成人av| 午夜影院在线不卡| 交换朋友夫妻互换小说| 精品久久国产蜜桃| 日韩亚洲欧美综合| 夜夜爽夜夜爽视频| 国产有黄有色有爽视频| 日日撸夜夜添| 色婷婷久久久亚洲欧美| 精品一区二区三卡| 久久99一区二区三区| 国产男女内射视频| 国产精品一区www在线观看| 亚洲丝袜综合中文字幕| av网站免费在线观看视频| 欧美另类一区| 免费黄色在线免费观看| 热99久久久久精品小说推荐| 一级,二级,三级黄色视频| 97在线视频观看| 欧美xxxx性猛交bbbb| 国产不卡av网站在线观看| 最新的欧美精品一区二区| 国产亚洲午夜精品一区二区久久| 精品久久久精品久久久| 又黄又爽又刺激的免费视频.| 国产高清不卡午夜福利| 久久久久久伊人网av| 国内精品宾馆在线| 又黄又爽又刺激的免费视频.| 亚洲第一区二区三区不卡| 精品午夜福利在线看| 久久久久网色| 夜夜爽夜夜爽视频| 亚洲精品亚洲一区二区| 久久久国产一区二区| 午夜激情av网站| 18+在线观看网站| 国产精品偷伦视频观看了| 天堂8中文在线网| 女的被弄到高潮叫床怎么办| www.av在线官网国产| 亚洲美女视频黄频| 亚洲国产精品成人久久小说| 精品99又大又爽又粗少妇毛片| 亚洲精品aⅴ在线观看| 91精品伊人久久大香线蕉| 丁香六月天网| 夜夜看夜夜爽夜夜摸| 精品久久久久久久久av| 三级国产精品片| 国产在线视频一区二区| 一区在线观看完整版| 国产高清有码在线观看视频| 久久久久久人妻| 国产av码专区亚洲av| xxx大片免费视频| 简卡轻食公司| kizo精华| 国产极品粉嫩免费观看在线 | 蜜桃国产av成人99| 免费av不卡在线播放| 2018国产大陆天天弄谢| 少妇人妻精品综合一区二区| 人妻制服诱惑在线中文字幕| 色吧在线观看| a 毛片基地| 国产在线视频一区二区| 成人国产av品久久久| 久久亚洲国产成人精品v| 亚洲激情五月婷婷啪啪| 中文乱码字字幕精品一区二区三区| 久久久a久久爽久久v久久| 边亲边吃奶的免费视频| 99九九线精品视频在线观看视频| 人妻系列 视频| 看免费成人av毛片| 精品一区二区免费观看| 国产在视频线精品| 夜夜爽夜夜爽视频| 18在线观看网站| av国产精品久久久久影院| 涩涩av久久男人的天堂| 欧美精品人与动牲交sv欧美| 9色porny在线观看| 国产爽快片一区二区三区| 一级毛片黄色毛片免费观看视频| 少妇的逼好多水| 在线观看www视频免费| 国产免费现黄频在线看| 精品午夜福利在线看| 亚洲人与动物交配视频| 亚洲欧美成人精品一区二区| 久久热精品热| 免费大片18禁| 丝瓜视频免费看黄片| 免费观看av网站的网址| 成人国语在线视频| 亚洲五月色婷婷综合| 九草在线视频观看| 日日撸夜夜添| 免费黄色在线免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 我的女老师完整版在线观看| 精品一区二区三区视频在线| 欧美激情 高清一区二区三区| 亚洲欧美成人综合另类久久久| 一本一本综合久久| 亚洲欧美成人精品一区二区| 日韩不卡一区二区三区视频在线| 男人爽女人下面视频在线观看| 伦精品一区二区三区| 我的女老师完整版在线观看| 欧美日韩av久久| 国产成人精品福利久久| 成人国产av品久久久| 女性生殖器流出的白浆| 成人国产av品久久久| 久久久久久久久久人人人人人人| 亚洲国产精品专区欧美| av不卡在线播放| a级片在线免费高清观看视频| 少妇丰满av| 免费日韩欧美在线观看| 欧美老熟妇乱子伦牲交| 男人操女人黄网站| 亚洲四区av| 久久av网站| 国产一区亚洲一区在线观看| 建设人人有责人人尽责人人享有的| 免费大片黄手机在线观看| xxx大片免费视频| 亚洲精品456在线播放app| 狠狠婷婷综合久久久久久88av| av天堂久久9| 国产一级毛片在线| 制服丝袜香蕉在线| 精品少妇久久久久久888优播| 亚洲第一av免费看| 亚洲熟女精品中文字幕| 亚洲精品国产av成人精品| 18+在线观看网站| av免费在线看不卡| 国产视频内射| 在线天堂最新版资源| 一级二级三级毛片免费看| 极品少妇高潮喷水抽搐| 亚洲av中文av极速乱| 少妇猛男粗大的猛烈进出视频| 日韩电影二区| 69精品国产乱码久久久| 人妻少妇偷人精品九色| 看十八女毛片水多多多| 国产精品欧美亚洲77777| 99热这里只有精品一区| 国产永久视频网站| 国产成人免费无遮挡视频| 最新中文字幕久久久久| 天堂中文最新版在线下载| 综合色丁香网| 亚洲欧美色中文字幕在线| 亚洲熟女精品中文字幕| 男人爽女人下面视频在线观看| 久久鲁丝午夜福利片| 九色成人免费人妻av| 蜜桃久久精品国产亚洲av| 天天躁夜夜躁狠狠久久av| videossex国产| 91精品一卡2卡3卡4卡| av免费在线看不卡| 成人二区视频| 中文天堂在线官网| 大片免费播放器 马上看| 熟女av电影| 又粗又硬又长又爽又黄的视频| 在线亚洲精品国产二区图片欧美 | 日日摸夜夜添夜夜爱| 在线观看免费视频网站a站| 亚洲美女视频黄频| 国产亚洲一区二区精品| 高清黄色对白视频在线免费看| 国产成人a∨麻豆精品| 国产精品不卡视频一区二区| 亚洲国产成人一精品久久久| 中文天堂在线官网| 亚洲成人手机| 91精品一卡2卡3卡4卡| 亚州av有码| 国产精品蜜桃在线观看| 亚洲一区二区三区欧美精品| 国产亚洲精品久久久com| 亚洲中文av在线| 国产高清国产精品国产三级| 日本欧美视频一区| 日韩精品有码人妻一区| 麻豆乱淫一区二区| 高清黄色对白视频在线免费看| 男男h啪啪无遮挡| 日韩精品免费视频一区二区三区 | 久久久国产精品麻豆| 日韩强制内射视频| 黄色欧美视频在线观看| 草草在线视频免费看| 亚洲成人一二三区av| 男的添女的下面高潮视频| 性色avwww在线观看| 国产在线免费精品| 久久影院123| 久久狼人影院| 成人国产av品久久久| 国产精品.久久久| 人人澡人人妻人| 欧美变态另类bdsm刘玥| 在线观看三级黄色| 国产在线一区二区三区精| 国产 一区精品| 日韩一区二区三区影片| 欧美精品一区二区大全| 高清欧美精品videossex| 欧美日韩综合久久久久久| xxxhd国产人妻xxx| 熟女av电影| 在线观看国产h片| 午夜免费观看性视频| 最近最新中文字幕免费大全7| 国产爽快片一区二区三区| 国产熟女午夜一区二区三区 | 22中文网久久字幕| 五月天丁香电影| 精品熟女少妇av免费看| 夜夜骑夜夜射夜夜干| 最近中文字幕2019免费版| 毛片一级片免费看久久久久| 一个人看视频在线观看www免费| 国产日韩一区二区三区精品不卡 | 亚洲性久久影院| 久久ye,这里只有精品| 一区二区日韩欧美中文字幕 | 免费观看在线日韩| av电影中文网址| 国产成人精品在线电影| 亚洲精品av麻豆狂野| 国产一区二区三区综合在线观看 | 日韩av免费高清视频| 久久久国产欧美日韩av| www.色视频.com| 丝瓜视频免费看黄片| 看非洲黑人一级黄片| 日本黄色片子视频| 只有这里有精品99| 麻豆乱淫一区二区| 蜜臀久久99精品久久宅男| 欧美xxxx性猛交bbbb| 欧美精品一区二区免费开放| 成年av动漫网址| 女人精品久久久久毛片| 国产亚洲av片在线观看秒播厂| 一级片'在线观看视频| 久久毛片免费看一区二区三区| 丰满迷人的少妇在线观看| 高清视频免费观看一区二区| 亚洲色图综合在线观看| 波野结衣二区三区在线| 国产亚洲精品久久久com| 纵有疾风起免费观看全集完整版| a级毛片黄视频| 国产极品天堂在线| 久久精品夜色国产| 男人操女人黄网站| 老熟女久久久| av又黄又爽大尺度在线免费看| h视频一区二区三区| 成人漫画全彩无遮挡| av有码第一页| 最近手机中文字幕大全| 永久网站在线| 国产精品99久久99久久久不卡 | 高清黄色对白视频在线免费看| 欧美日韩av久久| 欧美亚洲 丝袜 人妻 在线| 在线观看www视频免费| av免费观看日本| 久久精品国产亚洲av天美| 亚洲综合色惰| 久久国产精品男人的天堂亚洲 | 久久久久视频综合| 在线观看免费视频网站a站| 欧美性感艳星| 看十八女毛片水多多多| 最后的刺客免费高清国语| 22中文网久久字幕| 国产无遮挡羞羞视频在线观看| 精品国产一区二区久久| 少妇精品久久久久久久| 国国产精品蜜臀av免费| 亚洲国产精品999|