基于DSP的FIR數字濾波器的實現
加入技術交流群
1．2 過渡帶及阻帶衰減
根據對過渡帶及阻帶衰減的要求，設計時可選擇窗函數的形狀，并估計窗口長度N。設待求濾波器的過渡帶用△ω表示，它近似等于窗函數的主瓣寬度。由于過渡帶△ω近似與窗口長度N成反比。即N=A／△ω，其中A決定于窗口形式，例如，矩形窗A=4π，哈明窗A=8π等。按照過渡帶及阻帶衰減情況，選擇窗函數形式。其設計原則是在保證阻帶衰減的情況下，盡量選擇主瓣比較窄的窗函數。
1．3 單位取樣響應的計算
計算濾波器的單位取樣響應h(n)時，可按正式進行：
式中，ω(n)是上面選擇好的窗函數。如果要求線性相位，則要求hd(n)和ω(n)均對(N-1)／2對稱，前面介紹的幾種窗函數已保證對(N-1)
／2偶對稱。而如要求h(n)對(N-1)／2奇對稱，只要保證hd(n)對(N-1)／2奇對稱即可。
驗算技術指標是否滿足要求時其設計出的濾波器頻率響應可采用下式進行計算：
計算上式時，可使用FFT算法。如果
不能滿足要求，那么，根據具體情況，可重復上述設計，直到滿足要求為止。2 基于DSP的FIR數字濾波器實現方案
2．1 濾波系統的差分方程
若輸入信號x(n)和濾波器的單位沖激響應h(n)在頻域分別為
和
，則其輸出信號的頻率響應為
。根據離散傅氏變換的性質，可以得到濾波系統的差分方程
。其中，h(n)的計算可依據濾波器的參數，并借助MATLAB仿真工具計算得到。2．2 實現方案
由于在卷積運算之前，h(n)的N個數值是已知的，因此，可先在程序存儲器中開辟N個單元來存放h(n)。由于輸入序列x(n)是不斷變化的，因此，在數據存儲器中可開辟N個存儲單元并對其進行移位寄存，其初始值分別為x(n)、x(n-1)……x(n-N-1)，然后采用循環尋址的方法對其進行訪問。每次輸入新的樣本時，可以以新樣本改寫滑窗中的最老的數據，而滑窗中的其他數據則不需要移動。利用片內8 kB(循環緩沖區長度)寄存器可對滑窗進行間接尋址，循環緩沖區地址首尾相鄰。8級循環緩沖區的結構如圖1所示，其中頂部為低地址，圖1中的(a)(b)(c)分別為初始狀態、輸入1個和2個樣本后的存儲情況。本文引用地址：http://www.czjhyjcfj.com/article/150927.htm
評論