av一区二区在线观看_亚洲男人的天堂网站_日韩亚洲视频_在线成人免费_欧美日韩精品免费观看视频_久草视

您的位置:首頁技術文章
文章詳情頁

Python利用FFT進行簡單濾波的實現

瀏覽:85日期:2022-08-06 11:15:38

1、流程

大體流程如下,無論圖像、聲音、ADC數據都是如下流程:

(1)將原信號進行FFT;

(2)將進行FFT得到的數據去掉需要濾波的頻率;

(3)進行FFT逆變換得到信號數據;

2、算法仿真

2.1 生成數據:

#采樣點選擇1400個,因為設置的信號頻率分量最高為600Hz,根據采樣定理知采樣頻率要大于信號頻率2倍,所以這里設置采樣頻率為1400Hz(即一秒內有1400個采樣點)x=np.linspace(0,1,1400)#設置需要采樣的信號,頻率分量有180,390和600y=2*np.sin(2*np.pi*180*x) + 3*np.sin(2*np.pi*390*x)+4*np.sin(2*np.pi*600*x)

2.2 對生成的數據進行FFT變換

yy=fft(y) #快速傅里葉變換yf=abs(fft(y))# 取模yf1=abs(fft(y))/((len(x)/2)) #歸一化處理yf2 = yf1[range(int(len(x)/2))] #由于對稱性,只取一半區間

2.3顯示轉換結果:

顯示原始FFT模值:

#混合波的FFT(雙邊頻率范圍)plt.figure(2)plt.plot(xf,yf,’r’) #顯示原始信號的FFT模值plt.title(’FFT of Mixed wave(two sides frequency range)’,fontsize=7,color=’#7A378B’) #注意這里的顏色可以查詢顏色代碼表

Python利用FFT進行簡單濾波的實現

顯示原始FFT歸一化后的模值:

#混合波的FFT(歸一化)plt.figure(3)plt.plot(xf1,yf1,’g’)plt.title(’FFT of Mixed wave(normalization)’,fontsize=9,color=’r’)

Python利用FFT進行簡單濾波的實現

由于對稱,只取一半區間進行顯示

plt.figure(4)plt.plot(xf2,yf2,’b’)plt.title(’FFT of Mixed wave)’,fontsize=10,color=’#F08080’)

Python利用FFT進行簡單濾波的實現

3、利用FFT進行濾波

例如將頻率為600HZ的噪聲濾掉,這里直接將該頻段的數據置零:

yy=fft(y) #快速傅里葉變換yreal = yy.real# 獲取實數部分yimag = yy.imag# 獲取虛數部分test_y =yyfor i in range(len(yy)): if i <=900 and i>=500: test_y[i]=0

對置零后的數據進行逆變換:

test = np.fft.ifft(test_y) #對變換后的結果應用ifft函數,應該可以近似地還原初始信號。

對還原的數據進行FFT變換的結果:

Python利用FFT進行簡單濾波的實現

濾波后的數據和原數據相對比:

藍色的為原數據,橙色的為濾波后的數據

Python利用FFT進行簡單濾波的實現

假設將400Hz和600Hz的信號都濾掉得到的信號圖像如下:

Python利用FFT進行簡單濾波的實現

4、對隨機噪聲進行濾波

源碼:

noise_size = 1400noise_array = np.random.normal(0, 2, noise_size) adc_value=[] for i in range(noise_size): adc_value.append(0) y= np.array(adc_value) + noise_arrayyy=fft(y) #快速傅里葉變換yf=abs(fft(y))# 取模yf1=abs(fft(y))/((len(y)/2)) #歸一化處理yf2 = yf1[range(int(len(y)/2))] #由于對稱性,只取一半區間#混合波的FFT(雙邊頻率范圍)xf = np.arange(len(y)) plt.figure(1)plt.plot(xf,yf,’r’) #顯示原始信號的FFT模值plt.title(’FFT of Mixed wave(two sides frequency range)’,fontsize=7,color=’#7A378B’) #注意這里的顏色可以查詢顏色代碼表yy=fft(y) #快速傅里葉變換yreal = yy.real# 獲取實數部分yimag = yy.imag# 獲取虛數部分test_y =yyfor i in range(len(yy)): if i <=1200 and i>=200: test_y[i]=0test = np.fft.ifft(test_y) #對變換后的結果應用ifft函數,應該可以近似地還原初始信號。y=testyy=fft(y) #快速傅里葉變換yf=abs(fft(y))# 取模yf1=abs(fft(y))/((len(y)/2)) #歸一化處理yf2 = yf1[range(int(len(y)/2))] #由于對稱性,只取一半區間#混合波的FFT(雙邊頻率范圍)xf = np.arange(len(y)) plt.figure(2)plt.plot(xf,yf,’r’) #顯示原始信號的FFT模值plt.title(’FFT of Mixed wave(two sides frequency range)’,fontsize=7,color=’#7A378B’) #注意這里的顏色可以查詢顏色代碼表

運行結果:

原數據頻譜圖:

Python利用FFT進行簡單濾波的實現

濾波后的頻譜圖:

Python利用FFT進行簡單濾波的實現

濾波后(藍色線)與原數據(紅色線)對比:

Python利用FFT進行簡單濾波的實現

以上這篇Python利用FFT進行簡單濾波的實現就是小編分享給大家的全部內容了,希望能給大家一個參考,也希望大家多多支持好吧啦網。

標簽: Python 編程
相關文章:
主站蜘蛛池模板: 一级毛片在线播放 | 日韩视频一区二区 | 日韩免费1区二区电影 | 国产成人99久久亚洲综合精品 | 欧美电影一区 | 天天操操操操操 | 国产高清精品在线 | 国产精品一区二区久久精品爱微奶 | 青青99| 久久躁日日躁aaaaxxxx | 91在线精品视频 | 黄色毛片在线观看 | 久草网址 | 一级毛片视频免费观看 | 成人影院在线视频 | 日韩一二三区视频 | 久久久一| 韩日一区 | 国产1区 | 国产区精品 | 99精品国产一区二区三区 | 欧美日韩淫片 | 天堂一区在线观看 | av影音资源 | 毛片在线视频 | 久久三级av| 久久久精品视频免费看 | 午夜精品久久久 | 久久精品亚洲欧美日韩久久 | 成人午夜性成交 | 91久久久精品国产一区二区蜜臀 | 雨宫琴音一区二区在线 | 中文在线一区二区 | 久草视频观看 | 免费中文字幕日韩欧美 | 伊人青青久久 | 免费一级片 | 欧美网址在线观看 | 一本在线 | 久久看看 | 久久久精品久久 |