基于matlabGUI的IIR低通數(shù)字濾波器設(shè)計

DSP 課程設(shè)計報告題目基于 matlab GUI 的 IIR 低通數(shù)字濾波器設(shè)計學(xué) 院 信息工程學(xué)院 專 業(yè) 通信工程 姓 名 學(xué) 號 指導(dǎo)教師 完成時間 2015 年 1 月 7 日 目錄1.設(shè)計的主要內(nèi)容及基本要求 .12. 總體設(shè)計思路 .13設(shè)計方案 .13.1 原始語音信號加噪處理 .13.2 濾波器設(shè)計 .23.2.1 數(shù)字濾波器的基本知識 .23.2.2 數(shù)字濾波器的選取 .33.3 濾除噪聲后的語音信號 .34.GUI 設(shè)計 .44.1 GUI 界面設(shè)計概述 .44.2 GUI 界面設(shè)計具體操作 .44.3GUI 界面演示 .45.課程設(shè)計總結(jié) .6參考文獻 .6附錄 .711.設(shè)計的主要內(nèi)容及基本要求要求錄制一段自己的語音信號后，對所采集的語音信號加入干擾噪聲，設(shè)計相應(yīng)的數(shù)字濾波器濾除噪聲并進行驗證，最后設(shè)計 GUI 界面。2. 總體設(shè)計思路1.獲取一段語音。2.對所獲取的語音信號進行加噪處理。3.設(shè)計數(shù)字濾波器。4.將加噪后的信號通過所設(shè)計的濾波器，濾除噪聲。5.驗證濾除噪聲后的語音信號。6.利用 matlab GUI 設(shè)計模板設(shè)計相應(yīng)的 GUI 界面。3設(shè)計方案3.1 原始語音信號加噪處理fs=44100; Au=1;s=audioread(C:UsersASUSDesktopyinpin.wav);sound(s,44100);t=(0:length(s)-1)/fs;f=linspace(-fs/2,fs/2,length(s);n=(Au*cos(2*pi*20000*t);x=s+n;sound(x,44100);y1=fft(s);y2=fft(x);y3=fftshift(y1);y4=fftshift(y2);subplot(221);plot(t,s);title(原語音信號);subplot(222);plot(t,x);title(加噪后的語音信號 );subplot(223);plot(f,abs(y3);title(原始語音信號頻譜);subplot(224);plot(f,abs(y4);title(加噪后的信號頻譜);20 2 4 6 8 10 12 14-1-0.500.51 低低低低低0 2 4 6 8 10 12 14-2-1012 低低低低低低低低-2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5x 104050100150200250 低低低低低低低低-2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5x 10400.511.52x 105 低低低低低低低低3.2 濾波器設(shè)計3.2.1 數(shù)字濾波器的基本知識數(shù)字濾波器分為 FIR 數(shù)字濾波器和 IIR 數(shù)字濾波器兩種，即有限沖激響應(yīng)濾波器（FIR ，F(xiàn)inite Impulse Response）濾波器和無線沖激響應(yīng)（ IIR，Infinite Impulse Response）濾波器。（一） FIR 數(shù)字濾波器 FIR 的特點： nNnzhzH10)()(不存在極點（z=0 除外） ，系統(tǒng)函數(shù)| z0|處收斂。系統(tǒng)單位沖激響應(yīng)在有限個 n 值處不為零。結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒有輸出到輸入的反饋。故只能用較高的階數(shù)達到高的選擇性。FIR數(shù)字濾波器的幅頻特性精度較之于 IIR 數(shù)字濾波器低，但是線性相位即不同頻率分量的信號經(jīng)過FIR 數(shù)字濾波器后他們的時間差不變。FIR 數(shù)字濾波器具有系統(tǒng)穩(wěn)定，易實現(xiàn)相位控制，允許設(shè)計多帶通濾波器等優(yōu)點。FIR 數(shù)字濾波器的結(jié)構(gòu)有四種：橫截性；級聯(lián)型結(jié)構(gòu)；頻率采樣型結(jié)構(gòu)；線性相位 FIR 濾波器的結(jié)構(gòu)（二）IIR 數(shù)字濾波器IIR 數(shù)字濾波器的特點：電位沖激響應(yīng) h(n)是無限長的；系統(tǒng)函數(shù) H(z)在有限 z 平面上（0 |z|）有極點存在；結(jié)構(gòu)上存在著輸出到輸入的反饋，也就是結(jié)構(gòu)上是遞歸型的。IIR 濾波器運算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運算組成。有限階 IIR 的表達式： 3)()(01 knxbknyayMNk NkkzazXYH10)(IIR 數(shù)字濾波器的結(jié)構(gòu)也有四種：直接 I 型；直接 II 型典范型；級聯(lián)型；并聯(lián)型。3.2.2 數(shù)字濾波器的選取由原始語音信號的頻譜圖可以看出，語音信號的能量集中在低頻部分，所以根據(jù)濾波器的特性和功能，所選擇的濾波器應(yīng)該具有濾除無用的高頻成分的作用，所以可以選用低通濾波器。又由于 IIR 數(shù)字濾波器的輸出不但取決于過去和現(xiàn)在的輸入，還取決于過去的輸出，所以，選擇IIR 數(shù)字濾波器比選擇 FIR 數(shù)字濾波器要好。所設(shè)計的 IIR 切比雪夫數(shù)字低通濾波器如下：T=1;Fs=1/T;Ap=1;As=10;Wp=0.55*pi/T;Ws=0.6*pi/T;wp=(2/T)*tan(Wp/2);ws=(2/T)*tan(Ws/2);N,wc=cheb1ord(wp,ws,Ap,As,s);B,A=cheby1(N,Ap,wc,s);Bz,Az=bilinear(B,A,Fs);H,W=freqz(Bz,Az,512,fs);3.3 濾除噪聲后的語音信號將加噪后的語音信號通過所設(shè)計的低通數(shù)字濾波器：r=filter(Bz,Az,x);并播放語音，聽其效果。0 1 2 3x 10400.511.5 低低低低低0 5 10 15-1-0.500.51低低低低低低低低低低低-4 -2 0 2 4x 1040100020003000 低低低低低-4 -2 0 2 4x 1040100020003000 低低低低低低低低4由濾除噪聲后的語音信號的頻譜圖可以看出，噪音已被濾除，試聽可知，原始語音信號基本沒有丟失。4.GUI 設(shè)計4.1 GUI 界面設(shè)計概述GUI 設(shè)計面板是 GUI 設(shè)計工具應(yīng)用的平面，面板上部提供了菜單和常用工具按鈕，左邊提供了多種如命令按鈕、單選按鈕、可編輯文本框、靜態(tài)文本框、彈出式菜單等。進行設(shè)計時, 首先單擊面板左邊所需的控件, 然后在右邊的圖形界面編輯區(qū)中再次單擊某一恰當位置, 這時將在該位上為圖形界面添加相應(yīng)的控件。一個圖形界面的完成, 除了設(shè)計其外觀外, 還有相當?shù)囊徊糠质峭ㄟ^屬性的設(shè)來完成的。應(yīng)用 MATLAB 制作這樣一個過程是非常方便的, 我們可以通過 GUI 操作來看到。該環(huán)境下要設(shè)計一個界面友好的仿真軟件，一般應(yīng)完成以下兩個步驟：（1）GUI 界面設(shè)計。主要是通過不同的文本框、按鈕等許多工具的使用，設(shè)計出一個圖形用戶界面。要清楚這個圖形界面的功能是什么，即在圖形界面上的操作會引發(fā)什么樣的結(jié)果。（2）回調(diào)函數(shù)的設(shè)計。用戶應(yīng)根據(jù)設(shè)計好的圖形界面的功能，針對各個不同的圖形對象來編寫出能夠?qū)崿F(xiàn)該功能的函數(shù)代碼，確保這個圖形界面能夠完成所預(yù)定的功能。4.2 GUI 界面設(shè)計具體操作在 matlab 命令行窗口中輸入 guide 或點擊 matlab 界面上面的菜單中的 圖標，打開 gui設(shè)計的圖形界面。然后拖入所要的圖形控件，按需要修改外觀和空間屬性，直至滿足要求。設(shè)置屬性點擊 guide 界面上方的 Run 按鈕，會生成一個 fig 文件，一個 m 文件，其中 fig 文件就是界面的圖形，m 文件是界面的回調(diào)函數(shù)，在 m 文件里每個控件的回調(diào)函數(shù)都已經(jīng)自動生成，控件要做的工作就是在文件框架下定義某些特殊要求的狀態(tài)并補充完整回調(diào)函數(shù)場, 使單擊控件時激活回調(diào)程序完成一定的功能。4.3GUI 界面演示51.初始界面2.輸入相應(yīng)參數(shù)并點擊相應(yīng)控件后界面65.課程設(shè)計總結(jié)這次的課程設(shè)計，讓我對數(shù)字信號處理有了更深的認識，對 Matlab 能更熟練地應(yīng)用，對matlab GUI 有了更多的了解。對于 GUI 在數(shù)字信號處理中的應(yīng)用中，數(shù)字信號處理這門學(xué)科的知識是基礎(chǔ)，要掌握數(shù)字信號處理的相關(guān)知識的原理后，并用代碼來實現(xiàn)，才能很好地結(jié)合MATLAB 進行 GUI 編程。參考文獻1陳后金數(shù)字信號處理第二版北京：高等教育出版社，2008.112劉衛(wèi)國MATLAB 程序設(shè)計與應(yīng)用第二版北京：高等教育出版社，2006.77附錄%GUI 源程序function varargout = gui_lmj(varargin)% GUI_LMJ MATLAB code for gui_lmj.fig% GUI_LMJ, by itself, creates a new GUI_LMJ or raises the existing% singleton*.% H = GUI_LMJ returns the handle to a new GUI_LMJ or the handle to% the existing singleton*.% GUI_LMJ(CALLBACK,hObject,eventData,handles,.) calls the local% function named CALLBACK in GUI_LMJ.M with the given input arguments.% GUI_LMJ(Property,Value,.) creates a new GUI_LMJ or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before gui_lmj_OpeningFcn gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to gui_lmj_OpeningFcn via varargin.% *See GUI Options on GUIDEs Tools menu. Choose “GUI allows only one% instance to run (singleton)“.% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help gui_lmj% Last Modified by GUIDE v2.5 07-Jan-2015 10:29:29% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct(gui_Name, mfilename, .gui_Singleton, gui_Singleton, .gui_OpeningFcn, gui_lmj_OpeningFcn, .gui_OutputFcn, gui_lmj_OutputFcn, .gui_LayoutFcn, , .gui_Callback, );if nargin endif nargoutvarargout1:nargout = gui_mainfcn(gui_State, varargin:);else8gui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before gui_lmj is made visible.function gui_lmj_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to gui_lmj (see VARARGIN)% Choose default command line output for gui_lmjhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes gui_lmj wait for user response (see UIRESUME)% uiwait(handles.figure1);% - Outputs from this function are returned to the command line.function varargout = gui_lmj_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure%