1、設計IIR濾波器                                       使用Matlab的FDAtool工具可以很方便地IIR設計濾波器。使用File菜單中的export選項可以把濾波器的參數(shù)輸出到

2、Matlab的工作空間中。若濾波器為IIR型，則輸出的變量名為G和SOS。它把高階的IIR濾波器轉(zhuǎn)換為一系列二階IIR濾波器的級聯(lián)。SOS為二階IIR濾波器的參數(shù)，G為縮放系數(shù)。計算IIR濾波器的輸出的時候，先將輸入數(shù)據(jù)乘以G, 然后一一通過SOS中的每個IIR濾波器，就得到最終的結(jié)果了。但是如果二階IIR濾波器的程序只支持定點小數(shù)運算的話，那么就需要調(diào)整一下系數(shù)了。因為設計出來的G有時候很小（不知道會不會很大），這樣一旦把輸入數(shù)據(jù)乘以G后，就會使得二階IIR濾波器的輸入過小，造成很大的誤差。我們需要把G分配到每個二階IIR濾波器的系數(shù)中，使得每次得到的中間結(jié)果都剛好不會溢出。若以IIRi表示

3、第I個二階IIR濾波器的話，則我們需要保證IIR1, IIR1 * IIR2, IIR1 * IIR2 * IIR3, 的頻率響應都不超過0dB，*號表示級聯(lián)。為了使得輸出的精度最好，我們還需要讓上面這些濾波器的頻率響應中的最大值最接近0dB。我用Python寫了個小程序來實現(xiàn)這樣的計算，其中最重要的部分就是計算二階IIR濾波器的頻率響應，還好在網(wǎng)上找到了直接計算的公式，省去了不少麻煩。【作者: 紅巖若愚】【訪問統(tǒng)計: 1562】【2006年01月6日 星期五 12:45】【 加入博采】【打印】 Trackback博客手拉手回復- 評論人：anonymous   2006-

4、09-27 17:02:07  請問一下應該怎么乘呀？- 評論人：Jimmeny   2006-04-30 12:35:37  問題解決了,在調(diào)它的庫的時候,IIR的系數(shù)一定要做到雙精度字對齊- 評論人：Jimmeny   2006-04-26 15:12:24  我用的是 -O2的優(yōu)化級別,146條指令是經(jīng)過C編譯器優(yōu)化之后的結(jié)果的了, 之前我也用過DSPF_sp_biquad,2階的時候是可以的,但記得級聯(lián)之后的結(jié)果就完全不一樣了,不知道是不是我的操作有誤,我的系數(shù)也是用tf2sos轉(zhuǎn)出來的,不過我會再試一下的,謝謝你啦

5、:)到時我把結(jié)果告訴你 另外我試了移頻的算法,對抑制嘯叫的效果確實不明顯,我試過6HZ,12hz和24hz,基本上都不能避免的一個問題就是有聲音的時候難免有一種波浪聲,好象是被調(diào)制過一樣,.- 評論人：紅巖若愚   2006-04-26 14:19:32  我曾做過6713的C語言程序的優(yōu)化，按照TMS320C6000 Optimizing Compiler User&s Guide中介紹的辦法優(yōu)化一下速度能快不少。不知道你設置的編譯優(yōu)化級別多少。另外我看了一下6713的DSPlib, 其中單精度二次IIR, DSPF_sp_biquad， 計算一個取樣數(shù)

6、據(jù)只需要4個時鐘，因此減少到80條指令應該沒有問題。 另外，利用matlab的tf2sos函數(shù)可以把6階的IIR轉(zhuǎn)換為3個2階的IIR級聯(lián)，這樣就可以直接到用DSPlib的函數(shù)了。- 評論人：Jimmeny   2006-04-26 13:36:43  剛才寫錯了,是15萬條,而不是15條啊. 我想了一下;我是48K采樣的,兩通道,數(shù)據(jù)量1S就是96K,問題是我現(xiàn)在處理1個6階的IIR需要110條指令,另外每一短還要乘因子和搬數(shù)據(jù),就是附加大概36條指令,96000*(110+36)*10=140M,而對220M的DSP就已經(jīng)超過一半的了.我現(xiàn)在是想辦法把146條

7、指令那里減少到80條左右,這個可能要用到優(yōu)化匯編和多級流水線,這個我就感覺比較棘手的了.- 評論人：紅巖若愚   2006-04-26 12:26:08  可能是我理解有誤，不過按照你的數(shù)據(jù)，10個6階IIR處理一個取樣數(shù)據(jù)需要150個指令周期，不知道取樣頻率是多少，如果是44100的話，那么一秒鐘就是：44100*150 = 6.615M個指令，這好像不存在資源不夠的問題啊。- 評論人：Jimmeny   2006-04-26 10:39:15  我現(xiàn)在做的是10個6階的IIR濾波器，左右聲道同時處理，用的是1300MIPS的671

8、3，我是這樣計算的，一個DMA中斷可以處理1024點的32位數(shù)據(jù)，共可以執(zhí)行230萬條單周期指令，做一段的6階EQ一共耗費了我15條指令的資源（左右聲道同時處理時），10段下來就接近160萬了，加上UI的資源開支，我的程序就有點提襟見肘的感覺了，我的IIR是用C來寫的，優(yōu)化得不好，我一直在尋找一個做多階IIR的浮點處理的匯編程序，也看了一下54X和55X里面的LIB SRC文件，不適合用在6000上，現(xiàn)在還在頭痛著，我想減少階數(shù)，用2階來做，這樣就可以調(diào)6000LIB庫的IIR了，但2階的邊沿不陡，做好的EQ很難,所以我現(xiàn)在還在猶豫中。- 評論人：紅巖若愚   2006-04-24 14:54:35  你的意思是說10個二次IIR濾波器級聯(lián)所需要的MIPS么，不知道你的程序具體用了多少MIPS，我現(xiàn)在做的一個項目芯片是5509A，32位的IIR濾波器，一個（左有聲道一共）大概是3.3MIPS左右。程序是在DSPLIB的源程序基礎上修改的。不太清楚 “多級流水線“ 是什么意思。- 評論人：Jimmeny   2006-04-24 10:37:04  我在SPtool里面做的IIR濾波器,再導出到工作區(qū),然后用tf2sos轉(zhuǎn)換得出系數(shù),這樣得出來的系數(shù)在級聯(lián)直接2型的每級之間不需要做sca