【正文】 erType_value==2)) [n,Wn]=cheb1ord(wp1,ws1,Rp_value,Rs_value)； set(,39。string39。,num2str(n)) else if((FilterType_value==3)||(FilterType_value==4)) [n,Wn]=cheb1ord(wp,ws,Rp_value,Rs_value)； set(,39。string39。,num2str(n)) end end else if(FilterDesign_value==3) if((FilterType_value==1)||(FilterType_value==2)) [n,Wn]=cheb2ord(wp1,ws1,Rp_value,Rs_value)； set(,39。string39。,num2str(n)) else if((FilterType_value==3)||(FilterType_value==4)) [n,Wn]=cheb2ord(wp,ws,Rp_value,Rs_value)； set(,39。string39。,num2str(n)) end end else if(FilterDesign_value==4) if((FilterType_value==1)||(FilterType_value==2)) [n,Wn]=ellipord(wp1,ws1,Rp_value,Rs_value)； set(,39。string39。,num2str(n)) else if((FilterType_value==3)||(FilterType_value==4)) [n,Wn]=ellipord(wp,ws,Rp_value,Rs_value)； set(,39。string39。,num2str(n)) end end end end end endelse if(DigitalFilter_value==2) if((FilterType_value==1)||(FilterType_value==2)) [n,Wn]=buttord(wp1,ws1,Rp_value,Rs_value)。 set(,39。string39。,num2str(n)) else if((FilterType_value==3)||(FilterType_value==4)) [n,Wn]=buttord(wp,ws,Rp_value,Rs_value)。 set(,39。string39。,num2str(n)) end end endend MinOrder_value=get(,39。Value39。)。 if(MinOrder_value==0) n=str2double(get(,39。String39。)) endswitch DigitalFilter_value %數(shù)字濾波器IIR或者FIR選擇case 1 %選擇IIR濾波器switch FilterDesign_value%IIR中Butterworth、ChebyshevChebyshevEllipise濾波器 case 1 %選擇Butterworth濾波器 switch FilterType_value%選擇濾波器類型 case 1 %低通濾波器 [b,a]=butter(n,Wn)。[h,f]=freqz(b,a,Nn,Fs_value)。 axes()。 %把下面程序得到的圖畫在Magazineif(DisplayType_value==1) %選擇Linear畫幅頻圖plot(f,abs(h)) %畫幅頻圖 else plot(f,20*log10(abs(h))) %選擇Logarithmic畫幅頻圖 end grid on。 axes()。 %把下面程序得到的圖畫在Phaseplot(f,angle(h))。 %畫相頻圖grid on。 case 2 %高通濾波器 [b,a]=butter(n,Wn,39。high39。)。 [h,f]=freqz(b,a,Nn,Fs_value)。axes()。 if(DisplayType_value==1) plot(f,abs(h)) else plot(f,20*log10(abs(h))) end grid on。 axes()。plot(f,angle(h))。grid on。 case 3 %帶通濾波器[b,a]=butter(n,Wn)。[h,f]=freqz(b,a,Nn,Fs_value)。axes()。 if(DisplayType_value==1) plot(f,abs(h)) else plot(f,20*log10(abs(h))) end grid on。 axes()。plot(f,angle(h))。grid on。 case 4 %帶阻濾波器[b,a]=butter(n,Wn,39。stop39。)。 [h,f]=freqz(b,a,Nn,Fs_value)。axes()。 if(DisplayType_value==1) plot(f,abs(h)) else plot(f,20*log10(abs(h))) end grid on。 axes()。plot(f,angle(h))。grid on。 endcase 2 %選擇設計Chebyshev1濾波器 ……case 3 %選擇設計Chebyshev2濾波器……case 4 %選擇設計Ellipse濾波器……endcase 2 %選擇FIR濾波器 switch Windows_value %FIR中的Windows選擇case 1 %選擇設計boxar濾波器……case 2 %選擇設計Bartlett濾波器……case 3 %選擇設計Blackman濾波器……case 4 %選擇設計hanning濾波器……case 5 %選擇設計hamming濾波器……case 6 %選擇設計kaiser濾波器……end end 運行和結果顯示1 運行和結果顯示在圖形界面下, 按“ 運行” ,出現(xiàn)如圖15 所示界面，選擇要設計的濾波器選項，編輯框中輸人要求設計的數(shù)字濾波器的性能指標, 選擇“Run”命令按鈕, 出現(xiàn)如圖所示的模擬低通原型濾波器幅頻響應和相頻響應曲線, 最小階次。當選擇了“Auto Run”按鈕，在“濾波器類型選擇”旁的下拉菜單框中選擇其他類型的濾波器會立刻出現(xiàn)此類型的濾波器的圖形。例如選擇設計IIR的Butterworth低通濾波器，其Fs=1000Hz，F(xiàn)p1=100Hz，F(xiàn)s1=300Hz,Rp=3dB,Rs=20dB,n=4,顯示如圖16 圖15 運行顯示界面 圖 16 Butterworth低通濾波器2 顯示不同濾波器的顯示結果選擇設計IIR低通濾波器，其Fs=1000Hz，F(xiàn)p1=100Hz，F(xiàn)s1=300Hz,Rp=3dB,Rs=20dB,n=4,利用四種不同濾波器進行設計，如圖11120顯示： 圖17 Butterworth低通濾波器 圖18 ChebyshevⅠ低通濾波器 圖19 ChebyshevⅡ低通濾波器 圖20 Ellipse低通濾波器選擇設計FIR帶通濾波器，其Fs=4000Hz，F(xiàn)p1=900Hz，F(xiàn)p2=1300Hz，F(xiàn)s1=600Hz, Fs2=1500Hz, Rp=1dB,Rs=40dB,n=20,Rp=3dB,Rs=20dB,n=4,利用四種不同濾波器進行設計，如圖2222226顯示： 圖21 Boxar窗帶通濾波器 圖22 Bartlett窗帶通濾波器 圖22 Blackman窗帶通濾波器 圖23 Hanning窗帶通濾波器 圖24 Hamming窗帶通濾波器 圖25 Kaiser窗帶通濾波器5 設計總結和心得 設計總結 GUI是實現(xiàn)人機交互的中介，具有強大的功能，可以完成許多復雜的程序模塊。使用它，需要具有一定的知識儲備和必要的經(jīng)驗技巧。并且要充分利用好MATLAB的幫助文檔，仔細研讀HELP是最好的辦法。需要了解函數(shù)句柄等必要基礎知識，熟悉各控件對象的基本屬性和方法操作，知曉不同控件的合適使用條件及其特有的功能，并會采用不同的使用手段來實現(xiàn)相同功能的設計。同時還需要詳細掌握菜單和控件。菜單很簡單，就是弄清除菜單之間的關系和如何調(diào)用就可以?？丶氖褂弥饕怯煤肅reateF和Callback屬性。CreateF中的語句就是在程序運行時，就立即執(zhí)行腳本。如果希望界面可控，那么最好用Callback屬性。在相應控件下，添加相應的腳本就可以實現(xiàn)比較復雜計算繪圖等功能。 在設計GUI的時候，要注意一定的原則和步驟，分析界面所要求實現(xiàn)的主要功能，明確設計任務，構思草圖，設計界面和屬性，編寫對象的相應代碼，實現(xiàn)控件的交互調(diào)用。 對于GUI在數(shù)字信號處理中的應用中，數(shù)字信號處理這門學科的知識是基礎，要掌握數(shù)字信號處理的相關知識的原理后，并用代碼來實現(xiàn)，才能很好地結合MATLAB進行GUI編程。參考文獻1．程佩青編著．數(shù)字信號處理教程．第三版．北京：清華大學出版社，2．徐成波　陶紅艷　楊菁　楊如明編著．數(shù)字信號處理及MATLAB實現(xiàn)．第二版．北京：清華大學出版社，3．王宏．．北京：4