[開發技巧]·Python實現信號濾波（基於scipy）

[開發技巧]·Python實現信號濾波（基於scipy）

 

我的網站--> http://www.yansongsong.cn

GitHub主頁--> https://github.com/xiaosongshine

​

1.背景介紹

 

在深度學習中，有時會使用Matlab進行濾波處理，再將處理過的數據送入神經網絡中。這樣是通常的處理方法，可是處理起來卻有些繁瑣，而且有時系統難以運行Matlab。Python做爲一種十分強大的語言，是支持信號濾波濾波處理的。

本文將以實戰的形式基於scipy模塊使用Python實現簡單濾波處理，包括內容有1.低通濾波，2.高通濾波，3.帶通濾波，4.帶阻濾波器。具體的含義你們能夠查閱大學課程，信號與系統。簡單的理解就是低通濾波指的是去除高於某一閾值頻率的信號；高通濾波去除低於某一頻率的信號；帶通濾波指的是相似低通高通的結合保留中間頻率信號；帶阻濾波也是低通高通的結合只是過濾掉的是中間部分。上面所說的內容會在實戰部分加以介紹，能夠對比理解一下。

如何實現的呢？個人理解，是經過時域轉換爲頻域，在頻域信號中去除相應頻域信號，最後在逆轉換還原爲時域型號。具體的內容仍是要查閱大學課程，信號與系統。本身學的很通常就不班門弄斧了。

有什麼做用呢？My Opinions，能夠消除一些干擾信號，以低通濾波爲例，例如咱們若是隻是統計脈搏信號波形，應該在1Hz左右，卻發現波形信號上有不少噪音，這些噪音都是成百上千Hz的，這些對於脈搏信號波形就屬於無用的噪音，咱們就能夠經過低通濾波器將超出某一閾值的信號過濾掉，此時獲得的波形就會比較平滑了。

 

2.實戰演練

首先咱們使用到了scipy模塊，能夠經過下述命令進行安裝：（我使用的Python==3.6）

pip install scipy

1).低通濾波

這裏假設採樣頻率爲1000hz,信號自己最大的頻率爲500hz，要濾除400hz以上頻率成分，即截至頻率爲400hz,則wn=2*400/1000=0.8。Wn=0.8

from scipy import signal

b, a = signal.butter(8, 0.8, 'lowpass')   #配置濾波器 8 表示濾波器的階數
filtedData = signal.filtfilt(b, a, data)  #data爲要過濾的信號

2).高通濾波

這裏假設採樣頻率爲1000hz,信號自己最大的頻率爲500hz，要濾除100hz如下頻率成分，即截至頻率爲100hz,則wn=2*100/1000=0.2。Wn=0.2

from scipy import signal

b, a = signal.butter(8, 0.2, 'highpass')   #配置濾波器 8 表示濾波器的階數
filtedData = signal.filtfilt(b, a, data)  #data爲要過濾的信號

3).帶通濾波

這裏假設採樣頻率爲1000hz,信號自己最大的頻率爲500hz，要濾除100hz如下，400hz以上頻率成分，即截至頻率爲100，400hz,則wn1=2*100/1000=0.2，Wn1=0.2； wn2=2*400/1000=0.8，Wn2=0.8。Wn=[0.02,0.8]

from scipy import signal

b, a = signal.butter(8, [0.2,0.8], 'bandpass')   #配置濾波器 8 表示濾波器的階數
filtedData = signal.filtfilt(b, a, data)  #data爲要過濾的信號

4).帶阻濾波

這裏假設採樣頻率爲1000hz,信號自己最大的頻率爲500hz，要濾除100hz以上，400hz如下頻率成分，即截至頻率爲100，400hz,則wn1=2*100/1000=0.2，Wn1=0.2； wn2=2*400/1000=0.8，Wn2=0.8。Wn=[0.02,0.8]，和帶通類似，可是帶通是保留中間，而帶阻是去除。

from scipy import signal

b, a = signal.butter(8, [0.2,0.8], 'bandstop')   #配置濾波器 8 表示濾波器的階數
filtedData = signal.filtfilt(b, a, data)  #data爲要過濾的信號

3.函數介紹 

1.函數的介紹

(1).濾波函數

scipy.signal.filtfilt(b, a, x, axis=-1, padtype='odd', padlen=None, method='pad', irlen=None)

輸入參數：

b: 濾波器的分子係數向量

a: 濾波器的分母系數向量

x: 要過濾的數據數組。（array型）

axis: 指定要過濾的數據數組x的軸

padtype: 必須是「奇數」、「偶數」、「常數」或「無」。這決定了用於過濾器應用的填充信號的擴展類型。{‘odd’, ‘even’, ‘constant’, None}

padlen：在應用濾波器以前在軸兩端延伸X的元素數目。此值必須小於要濾波元素個數- 1。（int型或None）

method：肯定處理信號邊緣的方法。當method爲「pad」時，填充信號；填充類型padtype和padlen決定，irlen被忽略。當method爲「gust」時，使用古斯塔夫森方法，而忽略padtype和padlen。{「pad」 ，「gust」}

irlen：當method爲「gust」時，irlen指定濾波器的脈衝響應的長度。若是irlen是None，則脈衝響應的任何部分都被忽略。對於長信號，指定irlen能夠顯著改善濾波器的性能。（int型或None）

輸出參數：

y:濾波後的數據數組

（2）.濾波器構造函數(僅介紹Butterworth濾波器)

scipy.signal.butter(N, Wn, btype='low', analog=False, output='ba')

輸入參數：

N:濾波器的階數

Wn：歸一化截止頻率。計算公式Wn=2*截止頻率/採樣頻率。（注意：根據採樣定理，採樣頻率要大於兩倍的信號自己最大的頻率，才能還原信號。截止頻率必定小於信號自己最大的頻率，因此Wn必定在0和1之間）。當構造帶通濾波器或者帶阻濾波器時，Wn爲長度爲2的列表。

btype : 濾波器類型{‘lowpass’, ‘highpass’, ‘bandpass’, ‘bandstop’},

output : 輸出類型{‘ba’, ‘zpk’, ‘sos’},

輸出參數：

b，a: IIR濾波器的分子（b）和分母（a）多項式係數向量。output='ba'

z,p,k: IIR濾波器傳遞函數的零點、極點和系統增益. output= 'zpk'

sos: IIR濾波器的二階截面表示。output= 'sos'

4.參考

https://blog.csdn.net/weixin_37996604/article/details/82864680 

https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.signal.filtfilt.html

https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.signal.butter.html#scipy.signal.butter