<![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]><![CDATA[<dfn id="yhprb"><delect id="yhprb"></delect></dfn>]]><![CDATA[<dfn id="yhprb"></dfn>]]><![CDATA[<dfn id="yhprb"><delect id="yhprb"></delect></dfn>]]><![CDATA[<dfn id="yhprb"></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"><strike id="yhprb"></strike></s></dfn>]]><![CDATA[<small id="yhprb"></small>]]><![CDATA[<small id="yhprb"><delect id="yhprb"></delect></small>]]><![CDATA[<small id="yhprb"></small>]]><![CDATA[<small id="yhprb"></small>]]><![CDATA[<dfn id="yhprb"><delect id="yhprb"></delect></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]>
IIR濾波器零相位數字濾波及其應用
作者：
時(shí)間：2011-01-13
來(lái)源：網(wǎng)絡(luò )
收藏
3.1零相位數字濾波的算法
為了克服差分數字濾波中存在的上述兩個(gè)問(wèn)題，可以采用一種零相位濾波的方法，該方法的基本思路是：先確定出濾波器的初始條件，然后將原序列的首尾進(jìn)行擴展，把擴展后序列通過(guò)濾波器，將所得結果反轉后再次通過(guò)濾波器，最后將所得結果再反轉，并去掉首尾的擴展部分，即可得到零相位濾波后的輸出序列[3]。本文采用四次差分濾波方式給出一種便于實(shí)現的詳細算法，假設輸入信號為x(n)，n=0,1,…,P。
(1) 編寫(xiě)差分濾波函數DiffFilter(b,a,x)，調用參數如前所述。
(2) 對于IIR濾波器通常有M=N，求濾波器的初始條件e(q)，q=0,1,…,N-1,這里記E=［e(0),e(1),…,e(N-1)］T，則有
E=1+a(1)1-1…0
a(2)01…0
MMMOM
a(N)00…-11× b(1)-b(0)a(0)
b(2)-b(0)a(1)
M
b(Nn)-b(0)a(N)(3)
(3) 對原輸入序列x(n)的首尾進(jìn)行擴展，在x(n)的前面添加3N個(gè)數，在后面添加3N個(gè)數，擴展后的序列記為x′(n)，n=0，1，…，P+6N。
x′(n)=2x(0)-x(3N-n)
x(n-3N)
2x(P)-x(P-1-(n-P-3N-n))
n＜3N
3Nn＜P+3N
P+3NnP+6N(4)
(4) 進(jìn)行正向濾波：先用b，a和x′(n)調用DiffFilter(b,a,x)函數，計算x′(n)通過(guò)
H(z)=∑Mm=0b(m)z-m/(1+∑Mm=1a(m)z-m)(5)
的輸出y′(n)，然后生成一長(cháng)度為P+6N+1的單位脈沖信號h(n)，用z，a和h(n)調用DiffFilter(b,a,x)函數，計算h(n)通過(guò)
H(z)=∑M-1m=0y′(0)e(m)z-m/(1+∑Mm=1a(m)z-m)(6)
的輸出y′(n)，并記y(n)=y′(n)+y″(n)。
(5) 進(jìn)行反向濾波：將正向濾波結果反轉，記x′(n)=y(P+6N-n)，重復步驟4°，得到反向濾波后的結果。
(6) 將反向濾波的結果反轉，并刪除首尾的擴展部分，即可得到最后的濾波結果。
3.2零相位數字濾波的實(shí)現
Borland公司推出的Delphi編程語(yǔ)言，具有合理的單元化結構、優(yōu)化的編譯環(huán)境，開(kāi)發(fā)速度快、編程效率高。在實(shí)現同樣功能的情況下與其它語(yǔ)言相比，不僅編寫(xiě)的代碼量少、程序可移植性強，而且還有許多優(yōu)秀的組件包可供使用。最為方便的是可以使用動(dòng)態(tài)數組，隨時(shí)能夠改變數組的長(cháng)度，這一點(diǎn)非常適合數字信號處理。
作者利用Delphi7編寫(xiě)了利用雙線(xiàn)性Z變換法設計Butterworth型IIR濾波器和利用窗函數法設計FIR濾波器，以及差分數字濾波算法和零相位數字濾波算法應用程序，并作為集成測試軟件平臺的一個(gè)虛擬儀器。其中零相位數字濾波算法程序設計流程如圖6所示。
圖6零相位數字濾波算法的流程圖
圖7所示為圖2所示信號零相位數字濾波后的波形，通過(guò)與原信號(圖2)和差分濾波后的信號(圖4)對比不難看出：零相位數字濾波后的輸出與原信號中的相位基本一致，并且起始部分沒(méi)有畸變。但不可否認，零相位數字濾波算法相對于普通差分濾波算法計算量要大許多，但以目前計算機的運算能力，計算量稍大并不是什么主要問(wèn)題。
圖7零相位濾波后的波形
4結束語(yǔ)
本文介紹了一種利用四次差分濾波算法，實(shí)現零相位數字濾波的方法，并利用Delphi7編寫(xiě)了應用軟件。通過(guò)與普通差分濾波器的實(shí)例對比分析，說(shuō)明零相位數字濾波不僅能夠避免相移，而且還能改善差分濾波起始部分的波形畸變。這一點(diǎn)在數字信號處理中具有重要的應用價(jià)值。
本文引用地址：http://dyxdggzs.com/article/151093.htm
濾波器相關(guān)文章:濾波器原理
濾波器相關(guān)文章:濾波器原理
低通濾波器相關(guān)文章:低通濾波器原理
電源濾波器相關(guān)文章:電源濾波器原理
數字濾波器相關(guān)文章:數字濾波器原理
技術(shù)專(zhuān)區
国产精品自在自线亚洲|国产精品无圣光一区二区|国产日产欧洲无码视频|久久久一本精品99久久K精品66|欧美人与动牲交片免费播放
<![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]><![CDATA[<dfn id="yhprb"><delect id="yhprb"></delect></dfn>]]><![CDATA[<dfn id="yhprb"></dfn>]]><![CDATA[<dfn id="yhprb"><delect id="yhprb"></delect></dfn>]]><![CDATA[<dfn id="yhprb"></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"><strike id="yhprb"></strike></s></dfn>]]><![CDATA[<small id="yhprb"></small>]]><![CDATA[<small id="yhprb"><delect id="yhprb"></delect></small>]]><![CDATA[<small id="yhprb"></small>]]><![CDATA[<small id="yhprb"></small>]]><![CDATA[<dfn id="yhprb"><delect id="yhprb"></delect></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]><![CDATA[<dfn id="yhprb"><s id="yhprb"></s></dfn>]]>
加入技術(shù)交流群
評論