頻率取樣法

對理想濾波器的頻率響應取樣。由下式 $H(k)=H_{d}(e^{j\omega })|_{\omega ={\frac {2\pi }{N}}k},k\in [0,N-1]$ 即可設計不同k點(離散的頻率點) 之值。舉例，將較小的k對應的響應值設計為1，較大的k對應的響應值設計為0,可得低通濾波器響應；將較小的k對應的響應值設計為0，較大的k對應的響應值設計為1,可得高通濾波器響應。

方法介紹

給定一個理想濾波器的離散時間傅立葉轉換$H_{d}(f)$，$h[n]$則為我們要設計的有限脈衝響應濾波器的脈衝響應，在$n\in [0,N-1]$區間設計。考慮$R(f)$是$r[n]=h[n+P]$的離散時間傅立葉轉換。

此方法基本精神要求為

${\begin{aligned}R({\frac {m}{N}}f_{s})=H_{d}({\frac {m}{N}}f_{s}),&\forall m\in 0,1,2,...N-1\\\\&f_{s}:sampling\ frequency\end{aligned}}$

若以 normalized frequency $F={\frac {f}{f_{s}}}$ 對公式進行變數變換 則 $R({\frac {m}{N}})=H_{d}({\frac {m}{N}})$

步驟一

取樣 normalized過後的理想濾波器的離散時間傅立葉轉換$H_{d}({\frac {m}{N}}),\forall m\in 0,1,2,...N-1$

步驟二

求$H_{d}({\frac {m}{N}})$的逆離散傅立葉轉換(IDFT)(inverse discrete Fourier transform) $r_{1}[n]={\frac {1}{N}}\sum _{m=0}^{N-1}H_{d}({\frac {m}{N}})exp(j{\frac {2\pi m}{N}}n)$

步驟三

考慮$N$的奇偶情形

若$N$是奇數，則

$r[n]={\begin{cases}r_{1}[n]&{\text{for}}\ n=0,1,...,k\ \ k={\frac {N-1}{2}}\\r_{1}[n+N]&{\text{for}}\ n=-k,-k+1,...,-1.\end{cases}}$

若$N$是偶數，則 $r[n]={\begin{cases}r_{1}[n]&{\text{for}}\ n=1,...,k\ \ k={\frac {N}{2}}\\r_{1}[n+N]&{\text{for}}\ n=-k,-k+1,...,-1.\end{cases}}$

步驟四

我們據此法得到的有限響應濾波器$h[n]$，即

$h[n]=r[n-k]$

相關證明

若$R(F)$是$r[n]$的離散時間傅立葉轉換

$R(F)=\sum _{n=-\infty }^{\infty }{r[n]e^{-j2\pi Fn}}=\sum _{n=0}^{N-1}{r_{1}[n]e^{-j2\pi Fn}}$

令$F={\frac {m}{N}}$,

$R({\frac {m}{N}})=\sum _{n=0}^{N-1}{r_{1}[n]e^{-j{\frac {2\pi m}{N}}n}}$

又$r_{1}[n]$是$H_{d}({\frac {m}{N}})$的IDFT

∴$R({\frac {m}{N}})=H_{d}({\frac {m}{N}})$

證明了照上述方法設計可得到在特定頻率上都會與理想濾波器響應吻合的有限響應近似濾波器

優點

• 簡單且直觀

缺點

• 所得到的有限響應濾波器並非最優化的，有限響應造成的漣波 大小介於用MSE和MINIMAX方法設計的濾波器之間

• 難以跟權重函數結合

應用實例

考慮理想的離散希爾伯特轉換濾波器，我們打算據此設計15點有限響應濾波器

希爾伯特轉換的頻率響應虛部

步驟一

先對此理想濾波器進行15個點的採樣

對此理想濾波器進行取樣

步驟二

$r_{1}[n]=ifft([0,-0.9,-1,-1,-1,-1,-1,-0.7,0.7,1,1,1,1,1,0.9])$

ifft為(逆離散傅立葉轉換)

步驟三

${\begin{aligned}r[n]=&[-0.0315,0.0182,-0.0693,0.0132,-0.1690,0.0078,-0.6206,0.0000,0.6206,\\&-0.0078,0.1690,-0.0132,0.0693,-0.0182,0.0315]\\\end{aligned}}$

步驟四

${\begin{aligned}h[n]=r[n-8]&=[-0.0315,0.0182,-0.0693,0.0132,-0.1690,0.0078,-0.6206,0.0000,\\&\ \ 0.6206,-0.0078,0.1690,-0.0132,0.0693,-0.0182,0.0315]\end{aligned}}$

步驟2~3的脈衝響應

最後設計出的濾波器頻率響應

紅線為頻率響應

參考文獻

• Jian-Jiun Ding (2024), Advanced Digital Signal Processing