2016-08-18 82 views
0

我試圖估計多次應用簡單1-2-1濾波器的效果,並確定了殘差比例。具體而言,我想重現這一情節:重複1-2-1濾波器的頻率響應

from Small et al., 2013

我用scipy.signal.freqz如下

filt = np.array([0.25,0.5,0.25]) 
w, h=signal.freqz(filt) 

而且我認爲,對於一個重複的過濾器,我只需要自己多次乘以h(因爲它在頻域中,而濾波就是卷積)。

但是,我無法得到和他們在論文中一樣的情節。我有三個主要問題,

  1. 我以爲1-2-1濾波器只是三角濾波器,有沒有其他方法來檢查其頻域響應?

  2. 如何在python中檢查重複1-2-1過濾器的頻率響應?本身不是多次?

  3. 我很難理解freqz輸出中的w(標準化頻率)單位。有人可以向我解釋如何在劇情中轉換爲波數?

謝謝。

+0

是包含情節可在網上的文件? –

+0

鏈接在這裏,但我不認爲它可以在線。 http://link.springer.com/article/10.1007/s00382-013-1980-9 你覺得在這裏發佈這樣的情節並不好嗎?這是論文中的一個小問題,他們試圖說他們的平滑方法是有效的。 –

+0

你問過關於理解單位。本文爲這些單位提供了背景。 –

回答

0

原來我沒有錯。通過繪製傳遞函數的絕對值並將歸一化頻率除以2π得到完全相同的曲線圖,並且應用多次濾波器的mutilpe時間恰好將頻率響應與其自身相比較。

filt = np.array([0.25,0.5,0.25]) 
w, h=signal.freqz(filt) 
plt.plot(w/(2*pi), abs(h**400), label='400 pass') 

Comparison between frequency response of repeating 1-2-1 filter