【正文】 *sin(wd*n).*cos(w0*n)./(pi*n)。H=20*log10(abs(fft(h,1024)))。subplot(221),stem(nn,hd,39。)。n39。title(39。)。subplot(222),stem(nn,w,39。)。title(39。)。n39。subplot(223),stem(nn,h,39。)。xlabel(39。)。實(shí)際沖激響應(yīng)39。w=(512:511)/511。k39。axis([ 140 20])。\omega/\pi39。title(39。)。color39。w39。運(yùn)行結(jié)果如下：圖37 用布萊克曼窗設(shè)計(jì)FIR帶通濾波器的運(yùn)行結(jié)果該低通濾波器的設(shè)計(jì)參數(shù)中過渡帶寬度為，由運(yùn)行結(jié)果（圖37）可以看出，設(shè)計(jì)出的濾波器實(shí)際特性過渡帶為，所以該設(shè)計(jì)滿足要求。close all。wpl=*pi。wsl=*pi。wd1=(wpl+wsl)/2。Bt=min(abs(wslwpl),abs(wshwph))。wn=[wd1/pi,wd2/pi]。stop39。wd11=(wpl+wsl)/2。M=(N1)/2。n=nn+eps。 % 理想沖激響應(yīng)w=hanning(N)39。 % 實(shí)際濾波器的幅頻特性HH=[H(513:1024) H(1:512)]。k39。xlabel(39。)。理想沖激響應(yīng)39。axis([17 17 ])。k39。axis([17 17 ])。漢寧窗39。xlabel(39。)。k39。axis([17 17 ])。n39。title(39。)。subplot(224),plot(w,HH,39。)。xlabel(39。)。濾波器幅頻特性39。set(gcf,39。,39。)。 頻率采樣法工程上，常給出頻域上的技術(shù)指標(biāo)，所以有時(shí)會(huì)采用更直接的方法頻率采樣法設(shè)計(jì)FIR數(shù)字濾波器。clear all。N=33。N1=fix(wc/(2*pi/N))。A=[ones(1,N1+1),zeros(1,N2),ones(1,N1)]。 H=A.*exp(j*theta)。v=1:N。stem (v ,h,39。)。n39。ylabel(39。),title(39。)。M=500。w=linspace(0,pi,M)。*w)。plot(w./pi,abs(X),39。)。幅頻特性39。xlabel(39。)。Hd(\omega)39。axis([0,1,])。plot(w./pi,20*log10(abs(X)),39。)。幅頻特性39。xlabel(39。)。dB39。axis([0,1,50,10])。2. 用頻率采樣法設(shè)計(jì)高通濾波器，其理想特性為： 0 N=33 1 程序如下：close all。clc。wc=pi/3。N2=N2*N11。 A(N1+2)=。theta=pi*[0:N1]*(N1)/N。 h=real(ifft(H))。subplot(2,2,1),stem(v ,h,39。)。沖激響應(yīng)39。ylabel(39。)。M=500。w=linspace(0,pi,M)。*w)。k39。xlabel(39。)。Hd(w)39。title(39。)。subplot(2,2,3),plot(w./pi,20*log10(abs(X)),39。)。幅頻特性39。\omega/\pi39。ylabel(39。)。運(yùn)行結(jié)果如下：311 用頻率采樣法設(shè)計(jì)FIR高通濾波器的運(yùn)行結(jié)果該低通濾波器的理想特性為阻帶頻率，通帶頻率，由運(yùn)行結(jié)果可以看出，設(shè)計(jì)出的濾波器實(shí)際特性為截止頻率，阻帶起始頻率，過渡帶寬度。clear all。N=33。N1=fix(wc/(2*pi/N))。theta=pi*[0:N1]*(N1)/N。h=real(ifft(H))。subplot(2,2,1),stem (v ,h,39。)。沖激響應(yīng)39。ylabel(39。)。M=500。w=linspace(0,pi,M)。*w)。k39。xlabel(39。)。Hd(\omega)39。axis([0,1,])。幅頻特性39。subplot(2,2,3),plot(w./pi,20*log10(abs(X)),39。)。幅頻特性39。xlabel(39。)。dB39。axis([0,1,80,10])。 等波紋最佳逼近法等波紋最佳逼近法是一種優(yōu)化設(shè)計(jì)法，它克服了窗函數(shù)法和頻率采樣法的缺點(diǎn)，使最大誤差（即波紋的峰值）最小化，并在整個(gè)逼近頻段上均勻分布。1. 下面用remez函數(shù)設(shè)計(jì)一個(gè)27階的FIR低通濾波器，程序如下：%n=27，通帶截止頻率為，阻帶截止頻率為，clear。clc。 f=[0 1]。 w=[ ]。 freqz(b)。close all。n=30。 a=[0 0 1 1]。 b=remez(n,f,a,w)。運(yùn)行結(jié)果如下：圖314 濾波器的幅頻特性曲線和相頻特性曲線 本章小結(jié)通過采用三種不同的方法設(shè)計(jì)濾波器，得出結(jié)論：窗口設(shè)計(jì)法概念清楚，但臨界頻率難以控制；頻率采樣法在取樣點(diǎn)處精確保證頻響要求，克服了窗函數(shù)設(shè)計(jì)法臨界頻率不易控制的缺點(diǎn)，但不能確保截止頻率的自由取值；等波紋最佳逼近法在相同階數(shù)下，可以獲得更好的頻率特性和衰減特性，但通帶不平滑。第4章 FIR濾波器的設(shè)計(jì)和仿真 設(shè)計(jì)濾波系統(tǒng)Simulink是一種用來實(shí)現(xiàn)計(jì)算機(jī)仿真的軟件工具。Simulink中的信號(hào)處理模塊（Signal Processing Blockset）提供了200多個(gè)高級(jí)數(shù)字信號(hào)處理（DSP）和數(shù)學(xué)函數(shù)，包括變換、矩陣運(yùn)算、FIR、IIR、自適應(yīng)濾波器等的實(shí)現(xiàn)。通過調(diào)用Simulink中的功能模塊構(gòu)建數(shù)字濾波器的仿真模型，在仿真過程中，通過雙擊各功能模塊改變參數(shù)，以獲得不同狀態(tài)下的仿真結(jié)果。圖41 濾波器仿真模型 仿真結(jié)果分析使用圖41所示的系統(tǒng)模型進(jìn)行濾波器仿真，信號(hào)源為100Hz的正弦波，疊加一個(gè)高頻噪聲后送入低通濾波器中。該低通濾波器的單位沖激響應(yīng)如圖43所示。仿真模塊運(yùn)行后，示波器顯示的波形如圖4圖45和圖46所示。圖45 Scope2（疊加噪聲后的信號(hào)）圖45所示的受噪聲干擾的信號(hào)是由信號(hào)源和一個(gè)經(jīng)過高通濾波器的隨機(jī)信號(hào)疊加而成。 圖46 Scope3（經(jīng)過濾波器濾除噪聲后的信號(hào)）疊加了高頻噪聲的信號(hào)經(jīng)過利用FDA Tool工具箱及窗函數(shù)法設(shè)計(jì)的FIR低通濾波器，濾除高頻噪聲后的信號(hào)如圖46所示，可以看出，經(jīng)過低通濾波器處理后的信號(hào)與信號(hào)源波形一致，（Scope1與Scope3波形對(duì)比），說明此低通濾波器的設(shè)計(jì)還是比較成功的?；贛atlab/Simulink的仿真技術(shù)為能在一個(gè)建模環(huán)境中進(jìn)行可視化操作和與動(dòng)態(tài)系統(tǒng)（用Matlab Simulink描述）進(jìn)行交互提供了一種有效的解決方案，拓展了Matlab和Simulink在數(shù)字信號(hào)處理系統(tǒng)中的應(yīng)用能力。在同樣的階數(shù)下，等波紋最佳逼近法可以獲得最佳的頻率特性和衰耗特性，具有通帶和阻帶平坦，過渡帶窄等優(yōu)點(diǎn)；頻率取樣法設(shè)計(jì)過程簡單，但阻帶衰減明顯，若適當(dāng)選取過渡帶樣點(diǎn)值，會(huì)取得較窗函數(shù)法略好的衰耗特性；窗函數(shù)設(shè)計(jì)法在階數(shù)較低時(shí)，阻帶特性不滿足設(shè)計(jì)要求，只有當(dāng)濾波器階數(shù)較高時(shí)，使用海明窗和凱塞窗可以達(dá)到阻帶衰耗要求。FDA Tool是一個(gè)很實(shí)用的濾波器設(shè)計(jì)工具，避免了繁瑣的編程過程，可以很直觀的觀察所設(shè)計(jì)濾波器的各項(xiàng)性能指標(biāo)。致 謝誠摯地感謝XXX老師的熱情幫助和精心指導(dǎo)！如果沒有老師的指導(dǎo)，我想我是很難在這么短的時(shí)間內(nèi)順利完成論文的設(shè)計(jì)。李老師的教學(xué)態(tài)度、敬業(yè)精神對(duì)我產(chǎn)生了很大影響。最后，再次對(duì)在大學(xué)四年里關(guān)心和幫助我的老師們表示衷心地感謝！參考文獻(xiàn)[1]高西全、丁玉美，數(shù)字信號(hào)處理。北京：航空航天大學(xué)出版社 [3]賀超英，MATLAB應(yīng)用與實(shí)驗(yàn)教程。北京: 科學(xué)出版社[5]張葛祥，李娜，MATLAB仿真技術(shù)與應(yīng)用。西安：西安電子科技大學(xué)出版社[7]董長虹等，MATLAB信號(hào)處理與應(yīng)用。北京：科學(xué)出版社畢業(yè)設(shè)計(jì)小結(jié) 經(jīng)過了三個(gè)月的學(xué)習(xí)和實(shí)踐，畢業(yè)設(shè)計(jì)已經(jīng)結(jié)束。我此次畢業(yè)設(shè)計(jì)所選的課題是基于MATLAB的FIR數(shù)字濾波器的設(shè)計(jì)與仿真。通過這次的畢業(yè)設(shè)計(jì)，使我對(duì)使用MATLAB軟件進(jìn)行數(shù)字信號(hào)處理有了更深的理解和掌握，在理論和實(shí)踐上都得到了很大的提高，所以這次任務(wù)的完成使我學(xué)到了很多東西。此次的畢業(yè)設(shè)計(jì)使我學(xué)習(xí)并熟悉了MATLAB的編程方法和數(shù)字濾波器的設(shè)計(jì)思路，了解了數(shù)字濾波器的工作原理，鞏固了MATLAB語言的使用能力，提高了自己動(dòng)手的能力，學(xué)到了很多知識(shí)，并且提高了自己分析問題的能力。撰寫畢業(yè)論文的過程是訓(xùn)練我們獨(dú)立地進(jìn)行科學(xué)研究的過程。撰寫畢業(yè)論文的過程，同時(shí)也是對(duì)專業(yè)知識(shí)學(xué)習(xí)的過程，而且是更生動(dòng)、更切實(shí)、更深入的專業(yè)知識(shí)的