1 引 言

　　目前FIR濾波器的實(shí)現(xiàn)方法主要有3種：利用單片通用數(shù)字濾波器集成電路、DSP器件和可編程邏輯器件實(shí)現(xiàn)。單片通用數(shù)字濾波器使用方便，但由于字長(zhǎng)和階數(shù)的規(guī)格較少，不能完全滿足實(shí)際需要。使用DSP器件實(shí)現(xiàn)雖然簡(jiǎn)單，但由于程序順序執(zhí)行，執(zhí)行速度必然不快。

　　FPGA有著規(guī)整的內(nèi)部邏輯陣列和豐富的連線資源，特別適合于數(shù)字信號(hào)處理任務(wù)，相對(duì)于串行運(yùn)算為主導(dǎo)的通用DSP芯片來(lái)說(shuō)，其并行性和可擴(kuò)展性更好。但長(zhǎng)期以來(lái)，FPGA一直被用于系統(tǒng)邏輯或時(shí)序控制上，很少有信號(hào)處理方面的應(yīng)用，其原因主要是因?yàn)樵?a class="contentlabel" href="http://www.biyoush.com/news/listbylabel/label/FPGA">FPGA中缺乏實(shí)現(xiàn)乘法運(yùn)算的有效結(jié)構(gòu)。本文利用FPGA乘累加的快速算法，可以設(shè)計(jì)出高速的FIR數(shù)字濾波器，使FPGA在數(shù)字信號(hào)處理方面有了長(zhǎng)足的發(fā)展。

　　2 Matlab設(shè)計(jì)濾波器參數(shù)

　　以表1的濾波器參數(shù)為例，分析設(shè)計(jì)高速FIR數(shù)字濾波器的方法。

　　利用Matlab為設(shè)計(jì)FIR濾波器提供的工具箱，選擇濾波器類型為低通FIR，設(shè)計(jì)方法為窗口法，階數(shù)為16，窗口類型為Hamming，Beta為0.5，F(xiàn)s為8.6 kHz，F(xiàn)C為3.4 kHz，導(dǎo)出的濾波器系數(shù)如下：

　　3 快速FIR濾波器算法的基本原理

　　(1) 分布式算法

　　分布式算法在完成乘加功能時(shí)是通過(guò)將各輸入數(shù)據(jù)每一對(duì)應(yīng)位產(chǎn)生的部分積預(yù)先相加形成相應(yīng)的部分積，然后再對(duì)各部分積進(jìn)行累加得到最終結(jié)果。

　　對(duì)于一個(gè)N(N為偶數(shù))階線性相位FIR數(shù)字濾波器，輸出可由式(1)表示：

　　(2) 乘法器設(shè)計(jì)

　　高性能乘法器是實(shí)現(xiàn)高性能的FIR運(yùn)算的關(guān)鍵，分析乘法器的運(yùn)算過(guò)程，可以分解為部分積的產(chǎn)生和部分積的相加兩個(gè)步驟。部分積的產(chǎn)生非常簡(jiǎn)單，實(shí)現(xiàn)速度較快，而部分積相加的過(guò)程是多個(gè)二進(jìn)制數(shù)相加的加法問(wèn)題，實(shí)現(xiàn)速度通常較慢。解決乘法器速度問(wèn)題，需要分別從這兩個(gè)方面入手，減小部分積的個(gè)數(shù)，提高部分積相加運(yùn)算的速度。

　　3.1 Booth算法

　　Booth算法針對(duì)二進(jìn)制補(bǔ)碼表示的符號(hào)數(shù)之間的相乘，即可以同時(shí)處理二進(jìn)制正數(shù)/負(fù)數(shù)的乘法運(yùn)算。Booth算法乘法器可以減少乘法運(yùn)算部分積個(gè)數(shù)，提高乘法運(yùn)算的速度。

　　下面討論一個(gè)M b