FIR低通濾波器的設(shè)計(jì)

1、技術(shù)要求基本原理2.1 FIR低通濾波器2.2圖形用戶界面GUI建立模型描述3.1單位沖激響應(yīng)及其幅頻響應(yīng)模塊目錄3.2不同窗函數(shù)低通幅頻響應(yīng)對(duì)比模塊3.3截止頻率調(diào)整模塊3.4界面的美化3.5功能模塊圖模塊功能分析或源程序代碼4.1單位沖激響應(yīng)及其幅頻響應(yīng)及其代碼4.2不同窗函數(shù)低通幅頻響應(yīng)對(duì)比模塊4.3截止頻率調(diào)整模塊調(diào)試過程及結(jié)論5.1調(diào)試過程5.2設(shè)計(jì)結(jié)果思考題心得體會(huì)參考文獻(xiàn)附錄2.2.4.4,5.8.8.8.9.,1.4.1.516.18.武漢理工大學(xué)專業(yè)課程設(shè)計(jì) 4 (數(shù)字信號(hào)處理)課程設(shè)計(jì)說明書FIR低通濾波器的設(shè)計(jì)1技術(shù)要求用窗函數(shù)法設(shè)計(jì)線性相位FIR低通濾波器。要求通帶截止

2、頻率3= n /4單位脈沖響應(yīng) h(n)的長度N=23。繪出h(n)及其幅頻響應(yīng)特性曲線。2基本原理2.1 FIR低通濾波器數(shù)字濾波器是數(shù)字信號(hào)處理學(xué)科的重要組成部分，應(yīng)用非常廣泛。數(shù)字濾波器，通常 是一種算法，或是一種數(shù)字處理設(shè)備。它的功能是將一組輸入的數(shù)字序列經(jīng)過一定的運(yùn)算 后變成為另一組輸出的數(shù)字序列。它的主要功能是對(duì)數(shù)字信號(hào)進(jìn)行處理，保留數(shù)字信號(hào)中 的有用成分，去除信號(hào)中的無用成分。數(shù)字濾波器是在模擬濾波器的基礎(chǔ)上發(fā)展起來的，但它們之間存在著一些重要差別。與模擬濾波器相比，數(shù)字濾波器具有精度高、穩(wěn)定性好、設(shè)計(jì)靈活等優(yōu)點(diǎn)。一般情況下數(shù) 字濾波器是一個(gè)線性非移變系統(tǒng)。與模擬濾波器相同的是，

3、數(shù)字濾波器也有低通、高通、 帶通、帶阻之分。本次實(shí)驗(yàn)中，就是要設(shè)計(jì)數(shù)字低通濾波器。低通濾波器是容許低于截止頻率的信號(hào)通過，但高于截止頻率的信號(hào)不能通過的濾 波器。從結(jié)構(gòu)上，數(shù)字濾波器可以分為遞歸型(IIR )數(shù)字濾波器和非遞歸型(FIR)數(shù)字濾 波器。本次實(shí)驗(yàn)要求利用FIR設(shè)計(jì)線性相位的低通濾波器。FIR最大的特點(diǎn)之一就是能夠 做成嚴(yán)格的線性相位關(guān)系。所為線性相位，就是指濾波器對(duì)不同頻率的正弦波產(chǎn)生的相位 延遲與正弦波的頻率呈線性關(guān)系。因而，在通過該濾波器后在濾波器通帶內(nèi)的所有信號(hào)頻 率成分，除了由相頻特性決定的延遲外，可以全部保留。窗函數(shù)法設(shè)計(jì)低通濾波器，窗函數(shù)法也稱為傅里葉級(jí)數(shù)法。從單位取

4、樣響應(yīng)的觀點(diǎn)來看，就是使設(shè)計(jì)的濾波器的h( n)逼近理想的濾波器單位取樣相應(yīng)hd( n)。13設(shè)計(jì)思想:先給定理想filter的頻響Hd(ejw)，設(shè)計(jì)一個(gè)FIR的filter的頻響為H(ejw )，使H(ejw ) 逼近 Hd(ejw)。設(shè)計(jì)過程:先用傅氏反變換求出理想filter的單位抽樣響應(yīng)hd(n),然后加時(shí)間窗w(n)對(duì)hd(n)截?cái)?以求得FIR filter的單位抽樣響應(yīng)h(n)。由此來看，窗函數(shù)的形狀及長度的選擇就尤為關(guān)鍵。由于長度已給定，利用矩形窗、三角窗、漢明窗、漢寧窗、布萊克曼窗和凱澤窗幾種窗函數(shù)對(duì)窗函數(shù)形狀不同進(jìn)行對(duì)比。2.2圖形用戶界面GUI圖形用戶界面(GUI)是指

5、采用圖形方式顯示的計(jì)算機(jī)操作用戶界面。與早期計(jì)算機(jī)使用的命令界面相比，圖形界面對(duì)于用戶來說在視覺上更易于接受。本次課程設(shè)計(jì)就是應(yīng)用matlab的GUI制作簡單的用戶界面，并且基于該界面仿真常見的幾種信號(hào)的模型。界面的引入，使用戶的操作更加便捷，也使仿真的結(jié)果呈現(xiàn)起來較為簡便。GUI本次實(shí)驗(yàn)中，由于對(duì)比窗函數(shù)形狀對(duì)低通濾波器幅頻響應(yīng)的影響，總共會(huì)有8個(gè)圖形的顯示。運(yùn)行后，由于圖形窗口較多，觀察比較不方便，所以決定用圖形用戶界面GUI來實(shí)現(xiàn)圖形的顯示。3建立模型描述3.1單位沖激響應(yīng)及其幅頻響應(yīng)模塊根據(jù)實(shí)驗(yàn)指導(dǎo)書要求，不僅要求畫出最終低通濾波器的幅頻響應(yīng)，還要畫出單位沖激響應(yīng)h(n)的圖形，還有單

6、位沖激響應(yīng)在頻域的圖形。對(duì)單位沖激序列做離散傅里葉變換,得到其幅頻響應(yīng)。其中，離散傅里葉變換要編寫相關(guān)函數(shù)m文件。理想單位沖激響應(yīng)的幅頻響應(yīng)的曲線形狀為抽樣信號(hào)的形狀。3.2不同窗函數(shù)低通幅頻響應(yīng)對(duì)比模塊在加窗處理的環(huán)節(jié)，窗函數(shù)長度已經(jīng)確定，選擇何種類型的窗函數(shù)就是影響最終低通 濾波器幅頻響應(yīng)的性能的關(guān)鍵。本實(shí)驗(yàn)中，未規(guī)定窗函數(shù)的類型，所以在程序中，我決定 將不同用窗函數(shù)截?cái)嗪蟮牡屯V波器的形狀比較作為一個(gè)拓展功能。Matlab中提供的窗函數(shù)有:(1)矩形窗boxcar( n)(2)三角窗tnang(n)(3)漢明窗ham ming(n)(4)漢寧窗hanning(n)(5)布萊克曼窗blac

7、kma n(n)(6)凱澤窗kaiser( n)各窗函數(shù)參數(shù),旁瓣峰值幅度，過渡帶寬，阻帶最小衰減如圖一。窗函數(shù)旁瓣峰值過渡阻帶雖小幅度/dB帶寬哀減/dB矩形窗-13E V-12漢寧窗-31Sjt /；V-44哈明窗-41AVr .V-53不萊克曼窗-5712tt 一 V-74三角形窗-25加V-25凱寃窗(0-57jW-80=7. 865)圖1窗函數(shù)參數(shù)利用窗函數(shù)設(shè)計(jì)低通濾波器，只需利用函數(shù)fir1(N,wn,窗函數(shù)(N+1)，就可求得濾波器的幅頻響應(yīng)和相頻響應(yīng)。漢明窗的應(yīng)用方式為用fir1(N,wn)函數(shù)得到低通濾波器，用w,t=freqz(b,1,512), abs(w)求得幅頻響應(yīng)。

8、矩形窗的應(yīng)用方式為用fir1(N,wn ，boxcar(n)函數(shù)得到低通濾波器，w,t=freqz(b,1,512)，abs(w )求得幅頻響應(yīng)。三角窗的應(yīng)用方式為用fir1(N,wn ，triang(n)函數(shù)得到低通濾波器，w,t=freqz(b,1,512)，abs(w )求得幅頻響應(yīng)。漢寧窗的應(yīng)用為用fir1(N,wn，hanning(n)函數(shù)得到低通濾波器，w,t= freqz (b,1,512)，abs(w )求得幅頻響應(yīng)。布萊克曼窗的方式為fir1(N,wn，blackman(n)函數(shù)得到低通濾波器w,t=freqz(b,1,512)，abs(w )求得幅頻響應(yīng)。凱澤窗應(yīng)用方式為用

9、fir1(N,wnkaiser(n)函數(shù)得到低通濾波器，用w,t=freqz(b,1,512)，abs(w )求得幅頻響應(yīng)。3.3截止頻率調(diào)整模塊本次試驗(yàn)中，規(guī)定截止頻率為0.25 n為了能直觀的顯示不同的截止頻率，低通濾波器的形狀會(huì)有什么不同，我在 GUI上添加了一個(gè)可編輯文本框，把默認(rèn)數(shù)值設(shè)定為0.25,同時(shí)將程序中對(duì)應(yīng)的截止頻率從一個(gè)常數(shù)變?yōu)閺奈谋究蛑腥〉玫臄?shù)。3.4界面的美化為了制作出友好美觀的界面，必須對(duì)波形的顯示區(qū)域、幅度、頻率等進(jìn)行調(diào)整，改變曲線的粗細(xì)程度，對(duì)各個(gè)控件及整體界面進(jìn)行顏色的設(shè)置，并利用其工具欄將各個(gè)控件對(duì)齊。3.5功能模塊圖功能模塊圖如圖2。圖2功能模塊圖4模塊功能

10、分析或源程序代碼4.1單位沖激響應(yīng)及其幅頻響應(yīng)及其代碼單位沖激響應(yīng)的繪制思路比較簡單，就是將一條直線用stem函數(shù)繪制出來，其圖形必須為23個(gè)點(diǎn)。其代碼如下:n=0:22;x=n./n;stem( n,x); title('h( n)');axis(0,25,0,1.3);單位沖激響應(yīng)的幅頻響應(yīng)要用到專門的函數(shù) m文件。該函數(shù)文件可以在主程序中調(diào)用多次，節(jié)省篇幅。函數(shù) m文件代碼如下: fun cti on xk=dft(x n,N) n=0:1:N-1;k=n;WN=ex p(-j*2* pi/N);nk=n'*k;WNn k=WN.A nk;xk=x n* WNnk

11、;4.2不同窗函數(shù)低通幅頻響應(yīng)對(duì)比模塊窗函數(shù)不同，會(huì)對(duì)設(shè)計(jì)出的低通濾波器的性能不同。應(yīng)盡量選取旁瓣小主瓣窄的窗函數(shù)。為了細(xì)致觀察窗函數(shù)不同對(duì)設(shè)計(jì)的低通濾波器的不同影響，本次試驗(yàn)中設(shè)計(jì)了不同窗函數(shù)的低通濾波器的比較環(huán)節(jié)。各個(gè)窗函數(shù)設(shè)計(jì)低通濾波器的思路是相似的，只是其中的窗函數(shù)是不同的。用矩形窗設(shè)計(jì)低通濾波器的代碼:D=str2 nu m(get(ha ndles.edit1,'stn ng');wn=D* pi;N=23;b=fir1(N,w n,boxcar(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'Li neWi

12、dth',2.0);title('矩形窗');用漢明窗設(shè)計(jì)低通濾波器的代碼:D=str2 nu m(get(ha ndles.edit1,'stn ng');wn=D* pi;N=23;b=fir1(N,w n);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) titleC漢明窗');用三角窗設(shè)計(jì)低通濾波器的代碼:D=str2 nu m(get(ha ndles.edit1,'stn ng');wn=D* pi;N=23;b=fir1(N,

13、w n, tria ng(N+1);w,t=freqz(b,1,512);axis(0,3.5,0,1.5);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) axis(0,3.5,0,1.4);titleC三角窗');用漢寧窗設(shè)計(jì)低通濾波器的代碼:D=str2 nu m(get(ha ndles.edit1,'stn ng');wn=D* pi;N=23;b=fir1(N,w n,ha nnin g(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in e

14、Width',2.0) titleC漢寧窗');用布萊克曼窗設(shè)計(jì)低通濾波器的代碼:D=str2 nu m(get(ha ndles.edit1,'stn ng');wn=D* pi;N=23;b=fir1(N,w n,blackma n(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) title（'布萊克曼窗'）;用凱澤窗設(shè)計(jì)低通濾波器的代碼:D=str2 nu m(get(ha ndles.edit1,'stn ng');wn=

15、D* pi;N=23;b=fir1(N,w n,kaiser(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) titleC凱澤窗');4.3截止頻率調(diào)整模塊為了研究低通濾波器的截止頻率對(duì)濾波器的影響，本次實(shí)驗(yàn)設(shè)計(jì)了調(diào)整截止頻率的模D。但塊，用以直觀的顯示不同的截止頻率對(duì)低通濾波器的影響。為此，在控件中添加了一個(gè)可 編輯文本框，用以輸入要調(diào)整的截止頻率。并且將從可編輯文本框獲取的數(shù)值賦予 是要注意的是截止頻率是個(gè)有實(shí)際意義的常數(shù)，不能夠隨意指定其數(shù)值。所以，在文本框 內(nèi)獲取的是n的倍數(shù)。

16、wc=0.25* n也變?yōu)閣c=D*n。另外，從文本框中獲取的的值默認(rèn)為字符型，但計(jì)算時(shí)用到的應(yīng)為數(shù)值型，所以要進(jìn) 行類型的轉(zhuǎn)換。其代碼如下:D=str2 nu m(get(ha ndles.edit1,'stn ng');5調(diào)試過程及結(jié)論5.1調(diào)試過程在繪制h(n)時(shí)，本來打算繪制一條普通的直線，并用stem函數(shù)畫出，但是直線中不含自變量，因此將其從y=1改為y=n./n。在繪制h(n)的幅頻響應(yīng)時(shí)，我誤認(rèn)為dpt為matlab自帶的函數(shù)，所以運(yùn)行時(shí)報(bào)錯(cuò)，dpt 未定義。因此查找資料，發(fā)現(xiàn)該函數(shù)為自己設(shè)定的函數(shù)文件，并且建立了該文件。在繪制低通濾波器時(shí)，應(yīng)用了 freqz(b

17、,1,512)函數(shù)，在m文件中運(yùn)行正確。但在GUI中出現(xiàn)問題，后來發(fā)現(xiàn)，GUI中出現(xiàn)sub plot函數(shù)后就會(huì)不能正常繪圖。而freqz(b,1,512)函數(shù)，本身自帶了同時(shí)繪制相頻和幅頻的命令，因此會(huì)出現(xiàn)錯(cuò)誤。因此，我利用該函數(shù)求出幅頻特性，但是利用plot繪制幅頻特性。在調(diào)整截止頻率的模塊，運(yùn)行時(shí)報(bào)錯(cuò)，原因是忘記轉(zhuǎn)換從文本框獲取的數(shù)據(jù)類型。從文本框獲取了字符型的數(shù)據(jù)，應(yīng)當(dāng)轉(zhuǎn)換成字符型才能參與運(yùn)算。5.2設(shè)計(jì)結(jié)果F面將本次設(shè)計(jì)的結(jié)果按功能展示:chuanghanhuFIR低通濾波器的設(shè)計(jì)0 90 B0 70 60 5040.100 5925OS0 204截止頻率*pi圖1總體設(shè)計(jì)h(nj J

18、 11 IIIh（rO至月h(n)圖2單位沖激響應(yīng)hg）的頻率響應(yīng)圖3單位沖激幅頻響應(yīng)武漢理工大學(xué)專業(yè)課程設(shè)計(jì) 4 （數(shù)字信號(hào)處理）課程設(shè)計(jì)說明書圖5用漢明窗設(shè)計(jì)低通濾波器15矩形窗圖4用矩形窗設(shè)計(jì)低通濾波器漢明窗武漢理工大學(xué)專業(yè)課程設(shè)計(jì) 4 （數(shù)字信號(hào)處理）課程設(shè)計(jì)說明書圖7用漢寧窗設(shè)計(jì)低通濾波器17三角窗0 2lh(n)K璧"密0.25pi00.511.522.633.5截止頻率1.20G0 60 4圖6用二角窗設(shè)計(jì)低通濾波器漢寧®1 41.210 3101.60 4101.21.522 500.5截止頻率呦）伽盤記三瑩寶 II B 1 B 1 B 17 i武漢理工大學(xué)

19、專業(yè)課程設(shè)計(jì) 4 （數(shù)字信號(hào)處理）課程設(shè)計(jì)說明書布菜克曼窗圖9用凱澤窗設(shè)計(jì)低通濾波器18圖8用布萊克曼窗設(shè)計(jì)低通濾波器凱澤窗汶寧空武漢理工大學(xué)專業(yè)課程設(shè)計(jì) 4 (數(shù)字信號(hào)處理)課程設(shè)計(jì)說明書h(n玄迥圖10改變截止頻率6思考題線性相位滿足的條件？線性相位FIR濾波器的特性及其應(yīng)用領(lǐng)域?答：線性相位要求的條件是h(n)必須為實(shí)序列，且滿足h(n)=±h(N-n-1)。例如本設(shè) 計(jì)中使用的單位沖激響應(yīng)就符合實(shí)序列，且N=23，為奇數(shù)，h( n)=h(23-n-1)。在此條件已滿足的前提下共分四種情況。第一種是h(n)=h(N-n-1)，N為奇數(shù)；第二種是 h(n)=h(N-n-1)，N

20、為偶數(shù)；第三種是 h(n)=-h(N-n-1)，N 為奇數(shù)；第四種 h(n)=-h(N-n-1)，N為偶數(shù)；由此易知本次設(shè)計(jì)所要求的為第一種情況。(2)線性相位FIR濾波器的特性為：通過濾波器的信號(hào)中不同頻率的正弦波產(chǎn)生的相位 延遲與正弦波的頻率呈線性關(guān)系。因而，在通過該濾波器后在濾波器通帶內(nèi)的所有信號(hào)頻 率成分，除了由相頻特性決定的延遲外，可以全部保留。正因?yàn)榫€性相位能夠保留所有頻率成分，且產(chǎn)生的相位延遲與頻率呈線性關(guān)系，所以 相位失真的產(chǎn)生是可以控制的。非線性相位濾波器由于產(chǎn)生的相位失真與頻率不呈線性關(guān) 系所以容易造成信號(hào)的色散。其應(yīng)用領(lǐng)域?yàn)椋赫Z音信號(hào)、圖像信號(hào)、視頻信號(hào)要求線性相位，要求

21、信道具有線性相位特性，因此該類信號(hào)在濾波、除噪時(shí)要用到FIR線性相位濾波器。7心得體會(huì)本次課程設(shè)計(jì)要求在 matlab軟件中用窗函數(shù)法設(shè)計(jì)FIR線性相位低通濾波器。線性相位在實(shí)際之中應(yīng)用的十分廣泛。例如圖像、語音信號(hào)要求信道具有線性相位的特性。濾波器分為低通型、高通型、帶通型和帶阻型。而低通濾波器是最基本的濾波器。在本次設(shè)計(jì)中，基本的功能是比較容易實(shí)現(xiàn)的，因此增添了一些附加的功能。指導(dǎo)書中規(guī)定了 h(n)的長度及值,正好符合線性相位的第一類的第一種情況，所以加窗后的低通濾波器理論上具有線性相位。當(dāng)然加窗后處理的低通濾波器由于階段后出現(xiàn)的誤差，其相位不一定為十分嚴(yán)格的直線,但在相當(dāng)寬的一段頻率中

22、是符合嚴(yán)格的線性相位的。在規(guī)定了 h(n)及其長度后，窗函數(shù)的類型就成為影響低通濾波器的性能的主要因素了。因此本次設(shè)計(jì)最重要的一個(gè)附加功能就是加不同的窗函數(shù)后，所對(duì)應(yīng)的低通濾波器的形狀的對(duì)比。通過對(duì)比發(fā)現(xiàn)，加不同的窗函數(shù)后，低通濾波器的衰減到0.9和衰減到0.1對(duì)應(yīng)的頻率的不同，和尾巴衰減的速率的不同。 Matlab中自帶有6種窗函數(shù)，因此制作了六種窗函數(shù)的對(duì)比。由于用腳本 m文件畫圖，運(yùn)行后生成六個(gè)圖形的文件，不方便比較，于是制作了 GUI界面用于顯示繪制出的圖形。GUI界面的制作過程比較繁瑣，并且sub plot函數(shù)在坐標(biāo)軸中繪圖會(huì)出錯(cuò)所以，應(yīng)用freqz函數(shù)繪圖，不能夠同時(shí)繪制相位圖。由

23、于指導(dǎo)書要求的只有幅頻特性，因此將 函數(shù)與abs函數(shù)結(jié)合，繪制出低通濾波器的幅頻響應(yīng)。freqz通過本次設(shè)計(jì)，我了解了學(xué)習(xí) matlab這個(gè)強(qiáng)大的軟件的方法，就是多多實(shí)際,注重自學(xué)。Matlab中的函數(shù)十分便捷，只需要用fir1函數(shù)就可以得到相應(yīng)的低通濾波器。同時(shí),只需要freqz函數(shù)就可以求出低通濾波器的幅頻特性和相頻特性。8參考文獻(xiàn)1劉泉.數(shù)字信號(hào)處理原理與實(shí)現(xiàn).北京：電子工業(yè)出版社，20092張志涌.精通matlab .北京：北京航空航天大學(xué)出版社，201126附錄源程序fun cti on varargout = chua ngha nshu(vararg in) % CHUANGHA

24、NSHU M-file for chua ngha nshu.figCHUANGHANSHU, by itself, creates a new CHUANGHANSHU or raises theexisti ngsin glet on*H = CHUANGHANSHU returns the han die to a new CHUANGHANSHU or thehan die tothe existi ng sin glet on*CHUANGHANSHU('CALLBACK',hObject,eve ntData,ha ndles,.) calls the localf

25、unction named CALLBACK in CHUANGHANSHU.M with the given inputargume nts.CHUANGHANSHUC Prop erty','Value',.) creates a new CHUANGHANSHU or raisestheexist ing sin glet on*.Start ing from the left, property value p airs are app lied to the GUI before chua ngha nshu_O penin gFc n gets called

26、.A n un recog ni zed prop erty n ame or in valid value makes prop erty app licati on stop. All inputs are passed to chuanghanshu_OpeningFcn via varargin.*See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one in sta nee to run (si ngleto n)".% See also: GUIDE, GUIDATA, GUIH

27、ANDLES% Edit the above text to modify the res ponse to help chua ngha nshu% Last Modified by GUIDE v2.5 27-J un-2013 10:14:02% Begi n ini tialization code - DO NOT EDIT gui_S in glet on = 1;mfile name, .gui_State = struct('gui_Name','gui_S in glet on', gui_S in glet on, 'gui_O pe

28、nin gFc n', chua ngha nshu_O penin gFc n,'gui_Out pu tFc n',chua ngha nshu_Out pu tFc n,'gui_LayoutFc n',.'gui_Callback',)；if n argin && ischar(varargi n1)gui_State.gui_Callback = str2fu nc(varargi n1);end if n argoutvarargout1: nargout = gui_mai nfcn( gui_State,

29、vararg in :);elsegui_mai nfcn( gui_State, vararg in :);end % End in itializati on code - DO NOT EDIT% - Executes just before chua ngha nshu is made visible.fun cti on chua ngha nshu_O penin gFc n( hObject, eve ntdata, han dles, varargi n) % This fun cti on has no out put args, see Outpu tFc n.% hObj

30、ect han dle to figure % eve ntdata reserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)% varargincomma nd line argume nts to chua ngha nshu (seeARARGIN)% Choose default comma nd line out put for chua ngha nshu han dles.out put = hObjec

31、t;% Up date han dles structure guidata(hObject, han dles);% UIWAIT makes chua ngha nshu wait for user res ponse (see UIRESUME) % uiwait(ha ndles.figure1);% - Out puts from this fun cti on are returned to the comma nd line.fun cti on varargout = chua ngha nshu_Out pu tFc n( hObject, eve ntdata, han d

32、les)% varargoutcell array for returni ng out put args (see ARARGOUT);% hObjecthan dle to figure% eve ntdatareserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)% Get default comma nd line out put from han dles structurevarargout1 = han

33、dles.out put;% - Executes on butt on p ress in p ushbutt on1.fun cti on p ushbutt on 1_Callback(hObject, eve ntdata, han dles)% hObjecthan dle to pushbutto n1 (see GCBO)% eve ntdatareserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)D=

34、str2 nu m(get(ha ndles.edit1,'stn ng');wn=D* pi;N=23;b=fir1(N,w n,boxcar(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'Li neWidth',2.0);title('矩形窗');% - Executes on butt on p ress in p ushbutt on2.fun cti on p ushbutt on 2_Callback(hObject, eve ntdata, han dles) % hObject

35、han dle to pushbutto n2 (see GCBO) % eve ntdata reserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)D=str2 nu m(get(ha ndles.edit1,'stri ng');wn=D* pi;N=23;b=fir1(N,w n);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in eWid

36、th',2.0) titleC漢明窗');% - Executes on butt on p ress in p ushbutt on3.fun cti on p ushbutt on 3_Callback(hObject, eve ntdata, han dles) % hObject han dle to pushbutto n3 (see GCBO) % eve ntdata reserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user

37、data (see GUIDATA)D=str2 nu m(get(ha ndles.edit1,'stri ng');wn=D* pi;N=23;b=fir1(N,w n, tria ng(N+1);w,t=freqz(b,1,512);axis(0,3.5,0,1.5);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) axis(0,3.5,0,1.4);titleC三角窗');% - Executes on butt on p ress in p ushbutt on4.fun cti on p ushbutt

38、on4 _Callback(hObject, eve ntdata, han dles) % hObject han dle to pushbutto n4 (see GCBO) % eve ntdata reserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)D=str2 nu m(get(ha ndles.edit1,'stri ng');wn=D* pi;N=23;b=fir1(N,w n,ha

39、nnin g(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) titleC漢寧窗');% - Executes on butt on p ress in p ushbutt on5.fun cti on p ushbutt on 5_Callback(hObject, eve ntdata, han dles) % hObject han dle to pushbutto n5 (see GCBO) % eve ntdata reserved - to be defi ned in a

40、 future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)D=str2 nu m(get(ha ndles.edit1,'stri ng');wn=D* pi;N=23;b=fir1(N,w n,blackma n(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) titleC布萊克曼窗'）;% - Executes on butt on p ress in

41、 p ushbutt on6.fun cti on p ushbutt on 6_Callback(hObject, eve ntdata, han dles) % hObject han dle to pushbutto n6 (see GCBO) % eve ntdata reserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA)D=str2 nu m(get(ha ndles.edit1,'stri ng&

42、#39;);wn=D* pi;N=23;b=fir1(N,w n,kaiser(N+1);w,t=freqz(b,1,512);QX=plot(t,abs(w);set(QX,'L in eWidth',2.0) titleC凱澤窗');% - Executes duri ng object creati on, after sett ing all prop erties.fun cti on axes1_CreateFc n( hObject, eve ntdata, han dles)% hObjecthan dle to axes1 (see GCBO) % e

43、ve ntdata reserved - to be defi ned in a future vers ion of MATLAB% han dles empty - han dles not created un til after all CreateFc ns called% Hint: p lace code in Openin gFc n to popu late axesi fun cti on p ushbutt on 13_Callback(hObject, eve ntdata, han dies) % hObject han dle to pushbutto n13 (see GCBO) % eve ntdata reserved - to be defi ned in a future vers ion of MATLAB% han dlesstructure with han dles and user data (see GUIDATA) n=0:22;x=n./n;stem( n,x); ti