【文章內(nèi)容簡(jiǎn)介】 ，在輸入抽樣頻率Fs的前提下，判斷濾波器通帶臨界頻率（FpFp2）、濾波器阻帶臨界頻率（FstFst2）的歸一化頻率wp wpwsws1是否在[0,1]之間，如不正確顯示錯(cuò)誤對(duì)話框，其程序和效果圖（圖14）如下：function Fp1_Callback(hObject, eventdata, handles)%檢查輸入的通帶臨界頻率Fp1的歸一化頻率是否在[0,1]之間Fs_value=str2double(get(,39。String39。))。Fp1_value=str2double(get(,39。String39。))。wp1=2*Fp1_value/Fs_value。%如果不在[0,1]之間，顯示輸入錯(cuò)誤對(duì)話框if(wp1=1) errordlg(‘wp1=2*Fp1/Fs,歸一化頻率不在【0，1】之間，請(qǐng)輸入正確的參數(shù)’,‘錯(cuò)誤信息’) end圖 14 輸入?yún)?shù)錯(cuò)誤信息框顯示FpFstFst2的回調(diào)函數(shù)同理編程Run按鈕根據(jù)輸入的參數(shù)顯示設(shè)計(jì)的濾波器的幅頻特性和相頻特性圖，其程序如下：function Run_Callback(hObject, eventdata, handles)%點(diǎn)擊Run立即運(yùn)行AutoChoose,m文件，實(shí)現(xiàn)濾波器設(shè)計(jì)程序AutoChoose(handles)Auto Run按鈕當(dāng)選擇了Auto Run按鈕時(shí)，能夠根據(jù)所選的IIR的濾波器（Butterworth、 ChebyshevⅠ等）、FIR的窗口選項(xiàng)、圖形顯示類型（Linear、Logarithmic）立即顯示濾波器的幅頻特性和相頻特性圖，其中IIR的濾波器（Butterworth、 ChebyshevⅠ等）的回調(diào)函數(shù)程序如下：function FilterDesign_Callback(hObject, eventdata, handles)AutoRun_value=get(,39。Value39。)。%當(dāng)選擇了Filter（Butterworth、 ChebyshevⅠ等）中一種Filter運(yùn)行AutoChoose,m文件 if(AutoRun_value==1) AutoChoose(handles) endQuit按鈕退出濾波器設(shè)計(jì)窗口，其程序如下：function Quit_Callback(hObject, eventdata, handles)%點(diǎn)擊Quti按鈕退出Closefunction AutoChoose(handles)Nn=128。DigitalFilter_value=get(,39。Value39。)。FilterDesign_value=get(,39。Value39。)。Windows_value=get(,39。Value39。)。FilterType_value=get(,39。Value39。)。DisplayType_value=get(,39。Value39。)。Order_value=get(,39。Value39。)。Rp_value=str2double(get(,39。String39。))。Rs_value=str2double(get(,39。String39。))。Fs_value=str2double(get(,39。String39。))。Fp1_value=str2double(get(,39。String39。))。Fp2_value=str2double(get(,39。String39。))。Fs1_value=str2double(get(,39。String39。))。Fs2_value=str2double(get(,39。String39。))。wp1=2*Fp1_value/Fs_value。wp2=2*Fp2_value/Fs_value。ws1=2*Fs1_value/Fs_value。ws2=2*Fs2_value/Fs_value。wp=[wp1,wp2]。ws=[ws1,ws2]。if(DigitalFilter_value==1) if(FilterDesign_value==1) 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 else if(FilterDesign_value==2) if((FilterType_value==1)||(FilterType_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 %選擇設(shè)計(jì)Chebyshev1濾波器 ……case 3 %選擇設(shè)計(jì)Chebyshev2濾波器……case 4 %選擇設(shè)計(jì)Ellipse濾波器……endcase 2 %選擇FIR濾波器 switch Windows_value %FIR中的Windows選擇case 1 %選擇設(shè)計(jì)boxar濾波器……case 2 %選擇設(shè)計(jì)Bartlett濾波器……case 3 %選擇設(shè)計(jì)Blackman濾波器……case 4 %選擇設(shè)計(jì)hanning濾波器……case 5 %選擇設(shè)計(jì)hamming濾波器……case 6 %選擇設(shè)計(jì)kaiser濾波器……end end 運(yùn)行和結(jié)果顯示1 運(yùn)行和結(jié)果顯示在圖形界面下, 按“ 運(yùn)行” ,出現(xiàn)如圖15 所示界面，選擇要設(shè)計(jì)的濾波器選項(xiàng)，編輯框中輸人要求設(shè)計(jì)的數(shù)字濾波器的性能指標(biāo), 選擇“Run”命令按鈕, 出現(xiàn)如圖所示的模擬低通原型濾波器幅頻響應(yīng)和相頻響應(yīng)曲線, 最小階次。當(dāng)選擇了“Auto Run”按鈕，在