基于matlab信號(hào)處理工具箱的數(shù)字濾波器設(shè)計(jì)與仿真

1、論文名稱：論文名稱：基于MATLAB信號(hào)處理工具箱的數(shù)字濾波器設(shè)計(jì)與仿真作者：丁磊潘貞存叢偉摘要：摘要：傳統(tǒng)的數(shù)字濾波器的設(shè)計(jì)過程復(fù)雜，計(jì)算工作量大，濾波特性調(diào)整困難，影響了它的應(yīng)用。本文介紹了一種利用MATLAB信號(hào)處理工具箱（SignalProcessingToolbox）快速有效的設(shè)計(jì)由軟件組成的常規(guī)數(shù)字濾波器的設(shè)計(jì)方法。給出了使用MATLAB語言進(jìn)行程序設(shè)計(jì)和利用信號(hào)處理工具箱的FDATool工具進(jìn)行界面設(shè)計(jì)的詳細(xì)步驟。利用MA

2、TLAB設(shè)計(jì)濾波器，可以隨時(shí)對(duì)比設(shè)計(jì)要求和濾波器特性調(diào)整參數(shù)，直觀簡便，極大的減輕了工作量，有利于濾波器設(shè)計(jì)的最優(yōu)化。本文還介紹了如何利用MATLAB環(huán)境下的仿真軟件Simulink對(duì)所設(shè)計(jì)的濾波器進(jìn)行模擬仿真。關(guān)鍵字：關(guān)鍵字：數(shù)字濾波器MATLABFIRIIR引言：引言：在電力系統(tǒng)微機(jī)保護(hù)和二次控制中，很多信號(hào)的處理與分析都是基于對(duì)正弦基波和某些整次諧波的分析，而系統(tǒng)電壓電流信號(hào)（尤其是故障瞬變過程）中混有各種復(fù)雜成分，所以濾波器一直

3、是電力系統(tǒng)二次裝置的關(guān)鍵部件【1】。目前微機(jī)保護(hù)和二次信號(hào)處理軟件主要采用數(shù)字濾波器。傳統(tǒng)的數(shù)字濾波器設(shè)計(jì)使用繁瑣的公式計(jì)算，改變參數(shù)后需要重新計(jì)算，在設(shè)計(jì)濾波器尤其是高階濾波器時(shí)工作量很大。利用MATLAB信號(hào)處理工具箱（SignalProcessingToolbox）可以快速有效的實(shí)現(xiàn)數(shù)字濾波器的設(shè)計(jì)與仿真。1數(shù)字濾波器及傳統(tǒng)設(shè)計(jì)方法數(shù)字濾波器及傳統(tǒng)設(shè)計(jì)方法數(shù)字濾波器可以理解為是一個(gè)計(jì)算程序或算法，將代表輸入信號(hào)的數(shù)字時(shí)間序列轉(zhuǎn)化為

4、代表輸出信號(hào)的數(shù)字時(shí)間序列，并在轉(zhuǎn)化過程中，使信號(hào)按預(yù)定的形式變化。數(shù)字濾波器有多種分類，根據(jù)數(shù)字濾波器沖激響應(yīng)的時(shí)域特征，可將數(shù)字濾波器分為兩種，即無限長沖激響應(yīng)（IIR）濾波器和有限長沖激響應(yīng)（FIR）濾波器。IIR數(shù)字濾波器具有無限寬的沖激響應(yīng)，與模擬濾波器相匹配。所以IIR濾波器的設(shè)計(jì)可以采取在模擬濾波器設(shè)計(jì)的基礎(chǔ)上進(jìn)一步變換的方法。FIR數(shù)字濾波器的單位脈沖響應(yīng)是有限長序列。它的設(shè)計(jì)問題實(shí)質(zhì)上是確定能滿足所要求的轉(zhuǎn)移序列或脈沖

5、響應(yīng)的常數(shù)問題，設(shè)計(jì)方法主要有窗函數(shù)法、頻率采樣法和等波紋最佳逼近法等。在對(duì)濾波器實(shí)際設(shè)計(jì)時(shí)，整個(gè)過程的運(yùn)算量是很大的。例如利用窗函數(shù)法【2】設(shè)計(jì)M階FIR低通濾波器時(shí)，首先要根據(jù)（1）式計(jì)算出理想低通濾波器的單位沖激響應(yīng)序列，然后根據(jù)（2）式計(jì)算出M個(gè)濾波器系數(shù)。當(dāng)濾波器階數(shù)比較高時(shí)，計(jì)算量比較大，設(shè)計(jì)過程中改變參數(shù)或?yàn)V波器類型時(shí)都要重新計(jì)算。（1）（2）本文將以一個(gè)FIR濾波器的設(shè)計(jì)為例來說明如何使用MATLAB設(shè)計(jì)數(shù)字濾波器：在小

6、電流接地系統(tǒng)中注入83.3Hz的正弦信號(hào)，對(duì)其進(jìn)行跟蹤分析，要求設(shè)計(jì)一帶通數(shù)字濾波器，濾除工頻及整次諧波，以便在非常復(fù)雜的信號(hào)中分離出該注入信號(hào)。參數(shù)要求：96階FIR數(shù)字濾波器，采樣頻率1000Hz，采用Hamming窗函數(shù)設(shè)計(jì)。本例中，首先在FilterType中選擇Bpass（帶通濾波器）；在DesignMethod選項(xiàng)中選擇FIRWindow（FIR濾波器窗函數(shù)法），接著在WindowSpecifications選項(xiàng)中選取Ham

7、ming；指定Filterder項(xiàng)中的Specifyder＝95；由于采用窗函數(shù)法設(shè)計(jì)，只要給出通帶下限截止頻率Fc1和通帶上限截止頻率Fc2，選取Fc1＝70Hz，F(xiàn)c2＝84Hz。設(shè)置完以后點(diǎn)擊DesignFilter即可得到所設(shè)計(jì)的FIR濾波器。通過菜單選項(xiàng)Analysis可以在特性區(qū)看到所設(shè)計(jì)濾波器的幅頻響應(yīng)、相頻響應(yīng)、零極點(diǎn)配置和濾波器系數(shù)等各種特性。設(shè)計(jì)完成后將結(jié)果保存為1.fda文件。在設(shè)計(jì)過程中，可以對(duì)比濾波器幅頻相頻特