數(shù)字信號處理實訓(xùn)論文-模擬帶通濾波器和數(shù)字濾波器的仿真

摘要數(shù)字濾波在數(shù)字信號處理中，占有重要的地位。數(shù)字濾波包括FIR和IIR兩種濾波方式，其中FIR濾波具有很多優(yōu)點，可以在幅度特性隨意設(shè)計的同時，保證精確、嚴(yán)格的線性相位，濾波穩(wěn)定，不會出現(xiàn)遞歸型結(jié)構(gòu)中的極限振蕩等不穩(wěn)定現(xiàn)象，且誤差較小，可采用FFT算法實現(xiàn)，因此運算效率高。設(shè)計FIR濾波器常用的方法有窗函數(shù)法與頻率抽樣法，但是這兩種方法均不易精確控制通帶與阻帶的邊界頻率，所以在實際應(yīng)用中有一定的局限性。本次實訓(xùn)包含硬件跟軟件兩個局部。硬件局部先通過mutime仿真軟件仿真，得到差不多的中心頻率，帶寬，再畫原理圖，做電路板，做出的電路板通過調(diào)節(jié)相關(guān)的電位器，那么比擬容易調(diào)到所需的中心頻率跟帶寬。軟件局部做的是FIR濾波器，通過matlab軟件設(shè)計FIR濾波器，包含了高通，低通，帶通，帶阻，通過matlab軟件輸入不同的頻率及選擇相關(guān)函數(shù)設(shè)計FIR濾波器的窗函數(shù)。關(guān)鍵字：數(shù)字濾波器；Matlab；FIR；Abstractcharacteristicsrandomdesignatthesametime,ensureaccurate,strictlinearphase,filterstability,won'tappearrecursivetypestructureoftheextremeoscillationandunstablephenomenon,andlesserror,canuseFFTalgorithm,sotheoperationefficiency.DesignFIRfiltersthecommonlyusedmethodwithawindowfunctionmethodandfrequencysamplingmethod,butthetwomethodsarenoteasyprecisecontrolandstopbandbandpassboundaryfrequency,soinpracticalapplicationaresubjecttosomelimitations.Thehardwareandsoftwaretrainingcontainstwoparts.Firstmutimehardwarepartthroughsimulationsoftwaresimulation,getaboutofcenterfreq,bandwidth,drawprinciplediagrams,makecircuitboards,makebyadjustingthecircuitboardsofrelatedpotentiometer,itiseasiertoadjusttotherequiredcenterfrequencywithbandwidth.PartofthesoftwaretodoisFIRfilters,throughthematlabsoftwaredesignFIRfilters,containsahigh,lowpass,band-pass,withresistance,throughthematlabsoftwareinputdifferentfrequencyandchooserelevantfunctiondesignofFIRfilterswindowfunction.Keyword:digitalfilter;Matlab;FIR目錄引言 6１硬件局部 6帶通濾波器的設(shè)計及參數(shù)計算 6電路板PCB圖及mutime仿真軟件的電路仿真圖見附錄 7硬件設(shè)計所需相關(guān)器件 7芯片5532介紹 7硬件的調(diào)試與檢測 8檢測硬件電路 8模塊調(diào)試檢測 8整體調(diào)試檢測 82軟件局部 8濾波器概述 8窗函數(shù)的介紹 9FIR數(shù)字濾波器 9數(shù)字濾波器設(shè)計原理 9濾波器設(shè)計方法 103MATLABGUI的簡介 104基于MatlabGUI的數(shù)字濾波器設(shè)計 11GUI界面設(shè)計概述 11GUI頁面設(shè)計 12“濾波器設(shè)計軟件〞回調(diào)函數(shù)編寫 14軟件運行和結(jié)果 155總結(jié) 16謝辭 18參考文獻(xiàn) 19附錄 20引言數(shù)字濾波器是一種用來過濾時間離散信號的數(shù)字系統(tǒng)，通過對抽樣數(shù)據(jù)進(jìn)行數(shù)學(xué)處理來到達(dá)頻域濾波的目的。根據(jù)其單位沖激響應(yīng)函數(shù)的時域特性可分為兩類：無限沖激響應(yīng)〔IIR〕濾波器和有限沖激響應(yīng)〔FIR〕濾波器。與IIR濾波器相比，F(xiàn)IR的實現(xiàn)是非遞歸的，總是穩(wěn)定的；更重要的是，F(xiàn)IR濾波器在滿足幅頻響應(yīng)要求的同時，可以獲得嚴(yán)格的線性相位特性。因此，它在高保真的信號處理，如數(shù)字音頻、圖像處理、數(shù)據(jù)傳輸、生物醫(yī)學(xué)等領(lǐng)域得到廣泛應(yīng)用。數(shù)字濾波在數(shù)字信號處理中，占有重要的地位。數(shù)字濾波包括FIR和IIR兩種濾波方式，其中FIR濾波具有很多優(yōu)點，可以在幅度特性隨意設(shè)計的同時，保證精確、嚴(yán)格的線性相位，濾波穩(wěn)定，不會出現(xiàn)遞歸型結(jié)構(gòu)中的極限振蕩等不穩(wěn)定現(xiàn)象，且誤差較小，可采用FFT算法實現(xiàn)，因此運算效率高。設(shè)計FIR濾波器常用的方法有窗函數(shù)法與頻率抽樣法，但是這兩種方法均不易精確控制通帶與阻帶的邊界頻率，所以在實際應(yīng)用中有一定的局限性。１硬件局部1.1帶通濾波器的設(shè)計及參數(shù)計算圖1.1所示是一個多路負(fù)反應(yīng)二階有源帶通濾波器的原理圖，它使用單個通用運算放大〔通用運放〕接成單電源供電模式，易于實現(xiàn)。它的上限截止頻率和下限截止頻率可以非常近，具有非常很強(qiáng)的頻率選擇性。令C1=C2=C，Req是R1和R2并聯(lián)的值。品質(zhì)因數(shù)Q等于中心頻率除以帶寬，Q

=

fC/BW。由式可以看出可以通過讓R3的值遠(yuǎn)大于Req來獲得大的Q值。Q值越大，頻率選擇性越好，帶寬越小。反之那么反。令中心頻率為fc，那么計算公式如下：其中圖1.1電路原理圖1.2電路板PCB及mutime仿真軟件的電路仿真圖見附錄1.3硬件設(shè)計所需相關(guān)器件芯片5532介紹NE5532/SE5532/SA5532/NE5532A/SE5532A/SA5532A是一種雙運放高性能低噪聲運算放大器。相比擬大多數(shù)標(biāo)準(zhǔn)運算放大器，如1458，它顯示出更好的噪聲性能，提高輸出驅(qū)動能力和相當(dāng)高的小信號和電源帶寬。這使該器件特別適合應(yīng)用在高品質(zhì)和專業(yè)音響設(shè)備，儀器和控制電路和通道放大器。如果噪音非常最重要的，因此建議使用5532A版，因為它能保證噪聲電壓指標(biāo)。NE5532是高性能低噪聲運放，與很多標(biāo)準(zhǔn)運放〔如1458〕相似，它具有較好的噪聲性能，優(yōu)良的輸出驅(qū)動能力及相當(dāng)高的小信號與電源帶寬。特點如下：〔1〕小信號帶寬：10MHz；〔2〕輸出驅(qū)動能力：600Ω，10V；〔3〕輸入噪聲電壓：5nV/√HZ〔典型值〕；〔4〕DC電壓增益：50000；〔5〕AC電壓增益：10KHz時2200；〔6〕電源帶寬：140KHz；〔7〕轉(zhuǎn)換速率：9V/μS；〔8〕大電源電壓范圍：±3～±20V。管腳圖如下列圖1.4所示圖1.4芯片5532管腳圖硬件的調(diào)試與檢測.1檢測硬件電路電路板制作完后，先檢查了電路線路，由于線寬設(shè)為40mil，所以沒發(fā)現(xiàn)斷線現(xiàn)象。同時為了防止氧化，涂上了松香助焊劑。將元器件焊接完畢后，檢測各焊點，沒發(fā)現(xiàn)問題。.2模塊調(diào)試檢測經(jīng)過調(diào)試濾波器，能正常顯示波形，但在截止頻率時，幅值不能準(zhǔn)確對應(yīng)衰減。調(diào)整可調(diào)電阻，將其值對應(yīng)。.3整體調(diào)試檢測接線完畢，調(diào)整示波器，截止頻率對應(yīng)衰減幅值，正常。2軟件局部2.1濾波器概述隨著信息時代和數(shù)字世界的到來,數(shù)字信號處理已成為當(dāng)今一門極其重要的學(xué)科和技術(shù)領(lǐng)域。數(shù)字信號處理在通信、語音、圖像,自動控制、雷達(dá)、軍事、航空航天、醫(yī)療和家用電器等眾多領(lǐng)域得到了廣泛的應(yīng)用。在數(shù)字信號處理中,數(shù)字濾波器占有極其重要的地位?，F(xiàn)代數(shù)字濾波器可以用軟件或設(shè)計專用的數(shù)字處理硬件兩種方式來實現(xiàn),用軟件來實現(xiàn)數(shù)字濾波器優(yōu)點是隨著濾波器參數(shù)的改變,很容易改變?yōu)V波器的性能。根據(jù)數(shù)字濾波器單脈沖響應(yīng)的時域特性可將數(shù)字濾波器分為兩種,即IIR〔InfiniteImpulseResponse〕無限長脈沖響應(yīng)數(shù)字濾波器和FIR(FiniteImpulseResponse)有限長脈沖響應(yīng)數(shù)字濾波器。從功能上分類,可分為低通、高通、帶通、帶阻濾波器。2.2窗函數(shù)的介紹設(shè)計濾波器盡量要求窗函數(shù)滿足以下兩項要求：(1)窗譜主瓣盡可能地窄，以獲取較陡的過渡帶。(2)盡量減少窗譜的最大旁瓣的相對幅度。也就是能量盡量集中于主瓣，這樣使尖峰和波紋減小，就可增大阻帶的衰減。但是這兩項要求是不能同時滿足的。中選用主瓣寬度較窄時，雖然得到陡峭的過渡帶，但通帶和阻帶的波動明顯增加；中選用最小的旁瓣幅度時，雖能得到平坦的幅度響應(yīng)和較小的阻帶波紋，但過渡帶加寬，即主瓣會加寬。因此，實際所選用的窗函數(shù)往往是它們的折中。設(shè)計FIR濾波器常用的窗函數(shù)有：矩形窗、三角窗、漢寧窗、海明窗、布拉克曼窗、凱塞窗等。FIR數(shù)字濾波器.1數(shù)字濾波器設(shè)計原理根據(jù)數(shù)字濾波器沖激響應(yīng)的時域特征,可將數(shù)字濾波器分為2種,即無限長沖激響應(yīng)濾波器(IIRDF)和有限長沖激響應(yīng)濾波器(FIRDF)。FIRDF具有突出的優(yōu)點:系統(tǒng)總是穩(wěn)定的、易于實現(xiàn)線性相位、允許設(shè)計多通帶(或多阻帶)濾波器。因此FIRDF在數(shù)字信號處理中得到廣泛的應(yīng)用。但與IIRDF相比,在滿足同樣的阻帶衰減的情況下需要較高的階數(shù)。濾波器階數(shù)越高將占用更多的DSP運算時間。因此,對FIRDF的設(shè)計目標(biāo)是在滿足指標(biāo)要求的情況下盡量減少濾波器的階數(shù)。數(shù)字濾波器可以理解為是一個計算程序或算法,將代表輸入信號的數(shù)字時間序列轉(zhuǎn)化為代表輸出信號的數(shù)字時間序列,并在轉(zhuǎn)化過程中,使信號按預(yù)定的形式變化。FIRDF的沖激響應(yīng)h(k)是有限長的M階FIRDF系統(tǒng)函數(shù)可表示為濾波器的輸出：它的設(shè)計問題實質(zhì)上是確定能滿足所要求的轉(zhuǎn)移序列或脈沖響應(yīng)的常數(shù)問題,設(shè)計方法主要有窗函數(shù)法、頻率采樣法和等波紋最正確逼近法等。假設(shè)要逼近的理想濾波器的頻率響應(yīng)為常用的有巴特沃思濾波器、切比雪夫型濾波器、橢圓濾波器和巴塞爾濾波器。.2FIR濾波器設(shè)計方法在這里我們利用窗函數(shù)法設(shè)計FIR濾波器。窗函數(shù)法的根本想法是選取某一種適宜的理想頻率選擇性濾波器(這種濾波器總是有一個非因果，無限長的脈沖響應(yīng)),然后將它的脈沖響應(yīng)截斷(或加窗)以得到一個線形相位和因果的FIR濾波器。因此,這種方法的重點在于選擇某種恰當(dāng)?shù)拇昂瘮?shù)和一種適宜的理想濾波器。由于數(shù)字濾波器傳輸函數(shù)只與頻域的相對值有關(guān)，故在設(shè)計時可先將濾波器設(shè)計指標(biāo)進(jìn)行歸一化處理，設(shè)采樣頻率為Fs，歸一化頻率的計算公式是：利用完全設(shè)計法設(shè)計數(shù)字濾波器的步驟：〔1〕將設(shè)計指標(biāo)歸一化處理?！?〕根據(jù)歸一化頻率，選擇函數(shù)確定最小階數(shù)N和頻率參數(shù)Wn?！?〕確定窗口值。Windows指定窗函數(shù)類型，默認(rèn)為boxcar窗；可選Hanning、Hamming、Blackman、kaiser、bartlett等窗，每種窗都可以由Matlab的相應(yīng)函數(shù)生成。〔4〕確定傳遞函數(shù)的分母系數(shù)。函數(shù)fir1的調(diào)用格式為B=fir1(n,Wn,'ftype',Windows)，同時選擇在此函數(shù)中選擇要設(shè)計的濾波器的類型。其中，n為濾波器階數(shù)，Wc為截止頻率ftype決定濾波器類型，ftype=high，設(shè)計高通FIR濾波器，ftype=stop，設(shè)計帶阻FIR濾波器?！?〕再用freqz函數(shù)驗證設(shè)計結(jié)果。3MATLABGUI的簡介MATLAB是美國MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計算的高級技術(shù)計算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大局部。MATLAB是矩陣實驗室〔MatrixLaboratory〕的簡稱，和Mathematica、Maple并稱為三大數(shù)學(xué)軟件。它在數(shù)學(xué)類科技應(yīng)用軟件中在數(shù)值計算方面首屈一指。MATLAB可以進(jìn)行矩陣運算、繪制函數(shù)和數(shù)據(jù)、實現(xiàn)算法、創(chuàng)立用戶界面、連接其他編程語言的程序等，主要應(yīng)用于工程計算、控制設(shè)計、信號處理與通訊、圖像處理、信號檢測、金融建模設(shè)計與分析等領(lǐng)域。MATLAB的根本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分相似，故用MATLAB來解算問題要比用C，F(xiàn)ORTRAN等語言完相同的事情簡捷得多，并且mathwork也吸收了像Maple等軟件的優(yōu)點,使MATLAB成為一個強(qiáng)大的數(shù)學(xué)軟件。在新的版本中也參加了對C，F(xiàn)ORTRAN，C++，JAVA的支持?？梢灾苯诱{(diào)用,用戶也可以將自己編寫的實用程序?qū)氲組ATLAB函數(shù)庫中方便自己以后調(diào)用。MATALB可以創(chuàng)立圖形用戶界面GUI(GraphicalUserInterface),它是用戶和計算機(jī)之間交流的工具。MATLAB將所有GUl支持的用戶控件都集成在這個環(huán)境中并提供界面外觀、屬性和行為響應(yīng)方式的設(shè)置方法,隨著版本的提高,這種能力還會不斷加強(qiáng)。而且具有強(qiáng)大的繪圖功能,可以輕松的獲得更高質(zhì)量的曲線圖。4基于MatlabGUI的數(shù)字濾波器設(shè)計4.1GUI界面設(shè)計概述GUI設(shè)計面板是GUI設(shè)計工具應(yīng)用的平面，面板上部提供了菜單和常用工具按鈕，左邊提供了多種如命令按鈕、單項選擇按鈕、可編輯文本框、靜態(tài)文本框、彈出式菜單等。進(jìn)行設(shè)計時,首先單擊面板左邊所需的控件,然后在右邊的圖形界面編輯區(qū)中再次單擊某一恰當(dāng)位置,這時將在該位上為圖形界面添加相應(yīng)的控件。一個圖形界面的完成,除了設(shè)計其外觀外,還有相當(dāng)?shù)囊痪植渴峭ㄟ^屬性的設(shè)來完成的。因此在設(shè)置這些屬性時,要注意下面幾個常用又很重要的屬性設(shè)置：1、控件風(fēng)格和外觀〔1〕BackgroundColor：設(shè)置控件背景顏色，使用[RGB]或顏色定義?！?〕CData：在控件上顯示的真彩色圖像，使用矩陣表示?！?〕ForegroundColor：文本顏色?！?〕String屬性：控件上的文本，以及列表框和彈出菜單的選項?！?〕Visible：控件是否可見。2、對象的常規(guī)信息〔1〕Enable屬性：表示此控件的使能狀態(tài)，設(shè)置為on〞，表示可選，為“off〞時那么表示不可選?！?〕Style：控件對象類型。〔3〕Tag：控件表示〔用戶定義〕?！?〕TooltipString屬性：提示信息顯示。當(dāng)鼠標(biāo)指針位于此控件上時，顯示提示信息?！?〕UserData：用戶指定數(shù)據(jù)?！?〕Position：控件對象的尺寸和位置?！?〕Units：設(shè)置控件的位置及大小的單位〔8〕有關(guān)字體的屬性，如FontAngle，F(xiàn)ontName等。3、控件回調(diào)函數(shù)的執(zhí)行〔1〕BusyAction：處理回調(diào)函數(shù)的中斷。有兩種選項：即Cancel：取消中斷事件，queue：排隊〔默認(rèn)設(shè)置〕?！?〕ButtonDownFcn屬性：按鈕按下時的處理函數(shù)?！?〕CallBack屬性：是連接程序界面整個程序系統(tǒng)的實質(zhì)性功能的紐帶。該屬性值應(yīng)該為一個可以直接求值的字符串，在該對象被選中和改變時，系統(tǒng)將自動地對字符串進(jìn)行求值。〔4〕CreateFcn：在對象產(chǎn)生過程中執(zhí)行的回調(diào)函數(shù)?！?〕DeleteFcn：刪除對象過程中執(zhí)行的回調(diào)函數(shù)。〔6〕Interruptible屬性：指定當(dāng)前的回調(diào)函數(shù)在執(zhí)行時是否允許中斷，去執(zhí)行其他的函數(shù)。4、控件當(dāng)前狀態(tài)信息〔1〕ListboxTop：在列表框中顯示的最頂層的字符串的索引?！?〕Max：最大值?！?〕Min：最小值?！?〕Value：控件的當(dāng)前值。應(yīng)用MATLAB制作這樣一個過程是非常方便的,我們可以通過GUI操作來看到。該環(huán)境下要設(shè)計一個界面友好的仿真軟件，一般應(yīng)完成以下兩個步驟：〔1〕GUI界面設(shè)計。主要是通過不同的文本框、按鈕等許多工具的使用，設(shè)計出一個圖形用戶界面。要清楚這個圖形界面的功能是什么，即在圖形界面上的操作會引發(fā)什么樣的結(jié)果?！?〕回調(diào)函數(shù)的設(shè)計。用戶應(yīng)根據(jù)設(shè)計好的圖形界面的功能，針對各個不同的圖形對象來編寫出能夠?qū)崿F(xiàn)該功能的函數(shù)代碼，確保這個圖形界面能夠完成所預(yù)定的功能。4.2GUI頁面設(shè)計在matlab命令行窗口中輸入guide或點擊matlab界面上面的菜單中的圖標(biāo)，翻開gui設(shè)計的圖形界面如圖所要的圖形控件，按需要修改外觀和空間屬性，直至滿足要求，如圖。guide界面上方的Run按鈕，會生成一個fig文件，一個m文件，其中fig文件就是界面的圖形，m文件是界面的回調(diào)函數(shù)，在m文件里每個控件的回調(diào)函數(shù)都已經(jīng)自動生成，控件要做的工作就是在文件框架下定義某些特殊要求的狀態(tài)并補(bǔ)充完整回調(diào)函數(shù)場,使單擊控件時激活回調(diào)程序完成一定的功能。圖4.1GUI設(shè)計的圖形界面圖“濾波器設(shè)計軟件〞GUI界面設(shè)計4.3“濾波器設(shè)計軟件〞回調(diào)函數(shù)編寫1、DigitalFilter下拉菜單在下拉DigitalFilter菜單里選擇FIR濾波器設(shè)計中選擇FIR選項時，使IIR的Windows的選擇窗口不能使用，程序和效果圖〔圖4.3〕functionDigitalFilter_Callback(hObject,eventdata,handles)%讀取此時選擇的濾波器選擇，“IIR、FIR〞DigitalFilter_value=get(handles.DigitalFilter,'Value');%中選擇了FIR時，set(handles.Windows,'enable','on')end圖4.3濾波器選擇2、FilterType下拉菜單在下拉菜單FilterType中選擇濾波器類型,選擇高通濾波器時,〔圖〕functionFilterType_Callback(hObject,eventdata,handles)%讀取此時設(shè)計的濾波器的類型，“Highpass〞FilterType_value=get(handles.FilterType,'Value')圖4.4頻率參數(shù)輸入圖4.6“濾波器設(shè)計軟件〞主面板4.4軟件運行和結(jié)果1、運行和結(jié)果顯示在圖形界面下,按“運行〞,出現(xiàn)如下圖界面，選擇要設(shè)計的濾波器選項，編輯框中輸人要求設(shè)計的數(shù)字濾波器的性能指標(biāo),,選擇“Run〞命令按鈕,出現(xiàn)如下圖的模擬低通原型濾波器幅頻響應(yīng)和相頻響應(yīng)曲線,最小階次。2、以FIR高通濾波器為例,Fs=2200HZ,p1=500HZ,Fst1=300HZ,Rp=3dB,Rs=20dB,n=4,利用六種不同濾波器進(jìn)行設(shè)計，顯示Boxar窗高通濾波器的運行結(jié)果。圖4.8Boxar窗高通濾波器5總結(jié)參考文獻(xiàn)[1]胡光書.數(shù)字信號處理—理論、算法與實現(xiàn).北京清華大學(xué)出版社1997[2]吳湘淇.信號、系統(tǒng)與信號處理.北京電子工業(yè)出版社1996[3]余卞章．?dāng)?shù)字信號處理教程〔第二版〕．西安西北工業(yè)大學(xué)出版社2002[4]趙紅怡張常年．?dāng)?shù)字信號處理及其MATLAB實現(xiàn)．北京化學(xué)工業(yè)出版社2002[5]黃文梅熊桂林楊勇．信號分析與處理—MATLAB語言及應(yīng)用．長沙國防科技大學(xué)出版社2000[6]陳懷琛吳大正高西全．MATLAB在電子信息課程中的應(yīng)用．北京電子工業(yè)出版2002[7]網(wǎng)福昌.通信原理[M].清華大學(xué)出版社2005[8].通信系統(tǒng)仿真原理[M].機(jī)械工業(yè)出版社，2005[9]鐘麟王峰.MATLAB仿真技術(shù)與應(yīng)用教程[M].國防工業(yè)出版社2004[10]劉樹棠.現(xiàn)代通信系統(tǒng)(Matlab版)〔第二版〕電子工業(yè)出版社2007[11]徐平平等.數(shù)字通信－根底與應(yīng)用(第二版)電子工業(yè)出版社2021附錄1、電路板PCB及mutime仿真軟件的電路仿真圖圖1.2電路板PCB圖1.3電路仿真圖2、局部程序設(shè)計%選擇FIR濾波器case1%FIR中的Windows選擇switchWindows_value%選擇設(shè)計boxcar濾波器case1%選擇濾波器類型switchFilterType_value%高通濾波器case2w=boxcar(n+1);b=fir1(n,Wn,'high',w);[h,f]=freqz(b,1,Nn,Fs_value);axes(handles.Magnitude);if(DisplayType_value==1)plot(f,abs(h))elseplot(f,20*log10(abs(h)))endgridon;axes(handles.Phase);plot(f,angle(h));gridon;%選擇設(shè)計bartlett濾波器case2%選擇濾波器類型switchFilterType_value%高通濾波器case2w=bartlett(n+1);b=fir1(n,Wn,'high',w);[h,f]=freqz(b,1,Nn,Fs_value);axes(handles.Magnitude);if(DisplayType_value==1)plot(f,abs(h))elseplot(f,20*log10(abs(h)))endgridon;axes(handles.Phase);plot(f,angle(h));gridon;%選擇設(shè)計Blackman濾波器case3%選擇濾波器類型switchFilterType_value%高通濾波器case2w=blackman(n+1);b=fir1(n,Wn,'high',w);[h,f]=freqz(b,1,Nn,Fs_value);axes(handles.Magnitude);if(DisplayType_value==1)plot(f,abs(h))elseplot(f,20*log10(abs(h)))endgridon;axes(handles.Phase);plot(f,angle(h));gridon;%選擇設(shè)計hanning濾波器case4%選擇濾波器類型switchFilterType_value%高通濾波器case2w=hanning(n+1);b=fir1(n,Wn,'high',w);[h,f]=freqz(b,1,Nn,Fs_value);axes(handles.Ma