仍然使用FDATOOL設(shè)計(jì)濾波器，當(dāng)前設(shè)計(jì)一個(gè)數(shù)字帶通濾波器。至于用的是沖擊響應(yīng)不變法，還是其它的方法。暫時(shí)不考慮。FIR 需要的階數(shù)太多，也不考慮。使用IIR濾波，線性相位就不要想了。可以選巴特沃茲（最大平整度），或切比雪夫（最大陡降特性。）發(fā)現(xiàn)在相同的性能下切比雪夫需要的階數(shù)少。

生成的傳遞函數(shù)是按照多個(gè)二階單元級(jí)聯(lián)。系統(tǒng)提供 SOS（Second Order Section）也可以稱為“救命”矩陣。其思路是將高階傳遞函數(shù)分解為多個(gè)穩(wěn)定的二階函數(shù)級(jí)聯(lián)，保持系統(tǒng)穩(wěn)定。因?yàn)槭褂玫氖荌IR （無(wú)限沖擊響應(yīng)）構(gòu)成，注定其是非線性相位。會(huì)導(dǎo)致濾波后的信號(hào)波形畸變。導(dǎo)出濾波器參數(shù)文件.

每個(gè)SOS參數(shù)可以使用直接II型，實(shí)現(xiàn)如下：

編寫一段代碼測(cè)試一下。

待濾波信號(hào)的頻譜特性

可以看到org 信號(hào)上有噪聲，有直流偏置和低頻交流干擾。經(jīng)過(guò)帶通濾波后，去除了直流，去除了噪聲，基本還原原始信號(hào)signal。剛開始時(shí)并不穩(wěn)定需要等一段時(shí)間才能達(dá)到穩(wěn)定。

可見，此段代碼可以將一個(gè)輸入序列x 轉(zhuǎn)換為一個(gè)輸出序列 y，完成對(duì)x序列的濾波。

按照生成的結(jié)構(gòu)圖，編寫m文件的解釋：

每采集到一個(gè)數(shù)據(jù)時(shí)，以輸入數(shù)據(jù)作為參數(shù)，調(diào)用BandPassFilter（），得到一個(gè)濾波輸出。相當(dāng)于實(shí)時(shí)處理。

經(jīng)過(guò)濾波后的信號(hào)

經(jīng)過(guò)濾波+陷波后的信號(hào)