0 引言
在進(jìn)行高沖擊信號(hào)采集試驗(yàn)時(shí)，雖然在采集板中設(shè)計(jì)了硬件濾波單元，但依然存在著采集信號(hào)毛刺較多的現(xiàn)象，這為后續(xù)信號(hào)的識(shí)別和分析帶來(lái)了不便，因此需要采用軟件濾波的方法對(duì)信號(hào)進(jìn)行再次處理。
常用的數(shù)字濾波器有IIR濾波器和FIR濾波器。從性能上看，IIR濾波器所用的系數(shù)少，可以提高運(yùn)算速度，但其相位響應(yīng)，特別是帶沿為非線形，穩(wěn)定性不能一直得到保證。FIR數(shù)字濾波器可以得到嚴(yán)格的線性相位，但對(duì)同一幅頻特性要求，實(shí)現(xiàn)階數(shù)高，運(yùn)算量大。因此一般的IIR和FIR數(shù)字濾波器難以做到嚴(yán)格的線性相位與小運(yùn)算量兼顧，再加上采用浮點(diǎn)運(yùn)算，影響了對(duì)信號(hào)的實(shí)時(shí)處理速度。
在對(duì)采集的沖擊信號(hào)進(jìn)行軟件濾波處理時(shí)，為了加快計(jì)算速度，對(duì)濾波器的性能要求并不是很高，濾波器設(shè)計(jì)應(yīng)簡(jiǎn)單易行，因此整系數(shù)濾波器是很好的選擇。這里詳細(xì)說(shuō)明了建立在極零點(diǎn)抵消基礎(chǔ)上的簡(jiǎn)單整系數(shù)濾波器的遞歸實(shí)現(xiàn)方法，并在硬件自帶開(kāi)發(fā)環(huán)境IAREmbedded Workbench中編寫(xiě)濾波算法，將其應(yīng)用在所設(shè)計(jì)的信號(hào)采集系統(tǒng)中，取得了較好的試驗(yàn)效果。

1 建立在極零點(diǎn)抵消基礎(chǔ)上的簡(jiǎn)單整系數(shù)濾波器
建立在極零點(diǎn)抵消基礎(chǔ)上的整系數(shù)濾波器以“梳狀濾波器”為基礎(chǔ)，通過(guò)設(shè)置不同的極點(diǎn)位置及抵消相應(yīng)位置處的零點(diǎn)，分別設(shè)計(jì)出具有低通、高通、帶通和帶阻特性的整系數(shù)濾波器。
若一系統(tǒng)的傳遞函數(shù)為：

可以看出，H(z)有M個(gè)零點(diǎn)均勻分布在單位圓上，它的頻率響應(yīng)為：

這構(gòu)成了一個(gè)最簡(jiǎn)單的“梳妝濾波器”，它由其在2π內(nèi)有M個(gè)等分零點(diǎn)，并且幅度為梳妝而得名。
若z=1處設(shè)置一個(gè)極點(diǎn)，用以抵消H(z)處的零點(diǎn)，則系統(tǒng)的傳遞函數(shù)為：

當(dāng)M=10時(shí)，其極一零分布以及幅頻特性分別如圖1和圖2所示。可以看出，該濾波器具有低通特性，系數(shù)為整數(shù)。

由圖2可以看出，濾波器的阻帶性能不好，這是因sin函數(shù)的較大邊瓣引起的。為了壓縮這些邊瓣，取得更好的性能，需將濾波器的傳遞函數(shù)HLP(z)取k次方，得到：

為了進(jìn)行歸一化處理，需將式(5)除以N，濾波器的系數(shù)就是非整數(shù)，無(wú)法體現(xiàn)進(jìn)行濾波器整系數(shù)的優(yōu)勢(shì)。因此可以先將信號(hào)數(shù)據(jù)進(jìn)行濾波，濾波后的數(shù)據(jù)再除以N，這樣既利用了整系數(shù)濾波器的良好性能，又達(dá)到了歸一化的目的。通過(guò)這種方法設(shè)計(jì)出的濾波器，性能基本可滿(mǎn)足要求，在設(shè)計(jì)帶通和帶阻濾波器時(shí)，頻率的選擇會(huì)受到一定限制。

2 低通數(shù)字濾波器的軟件實(shí)現(xiàn)
2．1 整系數(shù)低通濾波器的實(shí)現(xiàn)
分析所要應(yīng)用的硬件電路和采集信號(hào)的參數(shù)時(shí)，結(jié)合上一節(jié)所介紹的濾波器設(shè)計(jì)方法，進(jìn)行簡(jiǎn)單整系數(shù)低通濾波器的設(shè)計(jì)。要求截止頻率fp=1 kHz，采樣頻率fs=60 kHz，通帶衰減ap=3 dB，阻帶衰減as=40 dB。由式(4)及式(5)，有：

由已知條件可知，相應(yīng)數(shù)字濾波器3 dB的通帶邊界頻率為：