【正文】 （42） 在設(shè)計(jì)過(guò)程中，將無(wú)限長(zhǎng)序列變?yōu)橛邢揲L(zhǎng)是通過(guò)時(shí)域加矩形窗乘積實(shí)現(xiàn)數(shù)據(jù)的截?cái)嗟?。時(shí)域乘積對(duì)應(yīng)了頻域卷積，從而對(duì)頻響特征發(fā)生的改變。常見(jiàn)的窗函數(shù)有：矩形窗、三角形(Bartlertt)窗、漢寧(Hanning)窗。海明(Hamming)窗、布拉克曼(Blackman)窗、凱澤(kaiser)窗等，下面介紹幾種常用的窗函數(shù)：　 )()(nwhnd??25 （42/sin)/si()( W: ()()( 2/)1(10wNweeenRR NjRjnjN????幅 度 函 數(shù)窗 譜3） 設(shè)計(jì)一個(gè)窗函數(shù)線性相位 FIR 濾波器，程序及運(yùn)行結(jié)果如下：omegac=。N=81。m=(N1)/2。n=0:2*m+10。h=omegac/pi*sinc(omegac*(nm)/pi)。w=[ones(1,N) zeros(1,length(n)N)]。hd=h.*w。omega=pi:2*pi/300:pi。Hd=freqz(hd,1,omega)。plot(omega,abs(Hd))。圖 42 輸出信號(hào)26 線性相位 FIR 濾波器的雷米茲交換法的仿真設(shè)計(jì)用函數(shù) firs 和 remez 設(shè)計(jì)一個(gè) 50 階 FIR 多通濾波器，濾波器理想頻率響應(yīng)為：f=[0 1] , a=[1 1 0 0 1 1 0 0 1 1 0 0 1 1]。將設(shè)計(jì)的濾波器的幅頻響應(yīng)和理想濾波器幅頻響應(yīng)進(jìn)行比較。%Samp clf。n=20。 f=[0 1]。 a=[1 1 0 0]。 b=firls(n,f,a)。 [h,w1]=freqz(b)。 bb=remez(n,f,a)。 [hh,w2]=freqz(bb)。 figure(1)plot(w1/pi,abs(h),w2/pi,abs(hh),39。r:39。)。 xlabel(39。歸一化頻率39。)。ylabel(39。振幅39。)。legend(39。firls39。,39。remez39。)。 grid on。figure(2)fs=50。t=0:1/fs:2。 f1=5。f2=15。 x1=sin(2*pi*f1*t)+8.*cos(2*pi*f2*t)。subplot(2,1,1),plot(t,x1)title(39。原始信號(hào)39。)y1=filter(b,1,x1)。y2=filter(bb,1,x1)。subplot(2,1,2),plot(t,y1,t,y2,39。r:39。)27legend(39。firls39。,39。remez39。)。 title(39。輸出信號(hào)39。)xlabel(39。時(shí)間/s39。) 0 101一一一一一一一firlsremez圖 43 幅頻特性0 1 21050510 一一一一0 1 221012一一/s一一一一firlsremez圖 44 原始與輸出信號(hào)%Samp clf。n=50。 f=[0 1]。28a=[1 1 0 0 1 1 0 0 1 1 0 0 1 1 ] b=firls(n, f,a)。 [h,w1]=freqz(b)。 bb=remez(n, f,a)。 [hh,w2]=freqz(bb)。 figure(1)plot (w1/pi,abs(h),39。b.39。,w2/pi, abs(hh),39。r:39。,f,a,39。m39。)。xlabel(39。歸一化頻率39。)。ylabel(39。振幅39。)。legend(39。firls39。,39。remez39。,39。理想特性39。)。 figure(2)impz(bb,1),title (39。脈沖響應(yīng)39。) xlabel(39。樣本數(shù)39。)。ylabel(39。幅度39。)0 101一一一一一一一firlsremez一一一一圖 45 幅頻特性 290 5 10 15 20 25 30 35 40 45 0一一一一一一一一一圖 46 脈沖響應(yīng) 結(jié)果分析首先由圖 41 可見(jiàn)，這是一個(gè)第一類線性相位濾波器。濾波器系數(shù) N 為 9（奇數(shù)） ，且 h(n)=h(N1n),幅度特性關(guān)于 對(duì)稱，在 和 處可以取任何???0???值；可以用于實(shí)現(xiàn)低通、高通、帶通和帶阻等各種濾波特性。程序運(yùn)行結(jié)果見(jiàn)圖 41 和圖 43。比較兩種設(shè)計(jì)方法所設(shè)計(jì)的濾波器的幅頻響應(yīng)可見(jiàn)，用 firs 設(shè)計(jì)的濾波器在整個(gè)頻率范圍內(nèi)（包括通帶和阻帶）均具有較好的響應(yīng)，但理想頻率響應(yīng)和實(shí)際頻率響應(yīng)的誤差較大；而函數(shù) remez 設(shè)計(jì)的濾波器在通帶內(nèi)具有等波紋特性，在邊界 和 處及過(guò)渡帶內(nèi)更接近理想頻響。由圖44 可見(jiàn)，輸入 5Hz 和 15Hz 的合成振動(dòng)信號(hào)，5Hz 的信號(hào)對(duì)應(yīng)的采樣頻率為 50Hz的歸一化頻率為 ,15Hz 的信號(hào)對(duì)應(yīng)的歸一化頻率為 ,對(duì)應(yīng)于該濾波器的幅頻響應(yīng)，歸一化頻率為 的振動(dòng)能通過(guò)濾波器，但歸一化頻率為 不能通過(guò)該濾波器，輸出結(jié)果也能看到著一點(diǎn)，只有 5Hz 的振動(dòng)通過(guò)了濾波器。但通過(guò)兩種類型的濾波器的濾波，輸出結(jié)果略有不同，這是由于兩種的濾波器的幅頻特性的略微不同造成的。用 firs 設(shè)計(jì)幅頻響應(yīng)在歸一化頻率為 的頻率略高于 remez 函數(shù)設(shè)計(jì)的輸出。運(yùn)用 firs 函數(shù)設(shè)計(jì)幅頻響應(yīng)在歸一化頻率為 的頻率比 remez 函數(shù)設(shè)計(jì)的振幅響應(yīng)衰減更大，因此運(yùn)用 firs 函數(shù)設(shè)計(jì)輸出更為平滑，即高頻成分更小，而 remez 函數(shù)設(shè)計(jì)的輸出較多的高頻成分，使的輸出結(jié)果不對(duì)稱。程序輸出結(jié)果見(jiàn)圖 45 和圖 46 所示。可見(jiàn) firs 設(shè)計(jì)的濾波器的通帶和阻帶具有較小的波紋，但在整個(gè)頻帶內(nèi)不一致。而 remez 函數(shù)設(shè)計(jì)的濾波器具有較大的30通帶和阻帶波紋，但在正個(gè)頻帶內(nèi)較為一致。圖 43 給出了 remez 函數(shù)設(shè)計(jì)的脈沖響應(yīng)，且具有偶對(duì)稱的關(guān)系，即 b（k）=b（n+2k） 。結(jié) 論我在這次的畢業(yè)設(shè)計(jì)中，設(shè)計(jì)初期遇到了很大困難，主要是軟件的安裝和使用問(wèn)題，但隨著設(shè)計(jì)的深入，這些問(wèn)題已得到初步解決，到目前為止，我已經(jīng)能運(yùn)用MATLAB 進(jìn)行濾波器的仿真。本文由四個(gè)部分組成，由數(shù)字濾波器，數(shù)字濾波器準(zhǔn)則，線性相位 FIR 數(shù)字濾波器，軟件設(shè)計(jì)組成。這四個(gè)部分都是緊密相連的，每一個(gè)部分都是為下一個(gè)部分做鋪墊，正是這樣一環(huán)接一環(huán)的理論加實(shí)踐的設(shè)計(jì)理念，才使得我順利完成設(shè)計(jì)。從這次運(yùn)用 MATLAB 軟件設(shè)計(jì)濾波器，讓我掌握了很多東西，使我懂得 函數(shù)Remez 實(shí)現(xiàn)了 ParkMcClellan 算法，這種算法利用了 Remez 交換算法，使的實(shí)際31頻率響應(yīng)擬合期望頻率響應(yīng)達(dá)到最優(yōu)。從實(shí)際和理想頻率響應(yīng)之間的最大誤差最小化來(lái)看，函數(shù) Remez 設(shè)計(jì)的濾波器是最優(yōu)的，因此成為最優(yōu)濾波器。我對(duì)最優(yōu)化從新有了更深的認(rèn)識(shí)，讓我受益非淺。雖然自己完成了本次畢業(yè)設(shè)計(jì)，我發(fā)現(xiàn)我還有很多東西沒(méi)有掌握，特別是MATLAB 軟件的使用，它以矩陣和向量的運(yùn)算以及運(yùn)算結(jié)果的可視化為基礎(chǔ) ,廣泛應(yīng)用于各個(gè)學(xué)科領(lǐng)域的數(shù)值分析，矩陣計(jì)算、函數(shù)生成、信號(hào)處理、圖形及圖象處理。其用途極其廣泛，對(duì)數(shù)字領(lǐng)域的幫助很大，只有熟練運(yùn)用我才能為自己的設(shè)計(jì)發(fā)揮出最大的潛力，我希望以后自己好好抓緊練習(xí)。還有雖然自己運(yùn)用最大誤差最小化準(zhǔn)則設(shè)計(jì)了最優(yōu)的濾波器，但是其優(yōu)化的參數(shù)還需大量實(shí)踐去證明，所以我要在以后的日子里付出更大的努力，期待自己能有更大更好的發(fā)現(xiàn)。致 謝短暫而漫長(zhǎng)的大學(xué)生活即將結(jié)束，多彩的人生也將拉開新的序幕。在這大學(xué)四年里的最后一次實(shí)驗(yàn)設(shè)計(jì)中，我頗有感觸：首先我要感謝輔導(dǎo)我的吳莉老師，她雖然年輕，但是對(duì)于通信方面的研究也是很深的，也正是因?yàn)槟贻p，才讓吳莉老師在帶課的同時(shí)抽空來(lái)輔導(dǎo)我做畢業(yè)設(shè)計(jì)。我本身的基礎(chǔ)很差，如果沒(méi)有吳莉老師的悉心指導(dǎo)，我想我在設(shè)計(jì)中會(huì)遇到許許多多的困難。其次我想感謝我的學(xué)長(zhǎng)，他們雖然已經(jīng)先踏入了社會(huì)，但是對(duì)我還是經(jīng)常點(diǎn)撥，32不能說(shuō)他們還是像以前那樣精通，但是對(duì)我的課題還是略知一二的。然后我要感謝我的父母，因?yàn)槲以谧霎厴I(yè)設(shè)計(jì)時(shí)遇到了許多困難，有時(shí)甚至堅(jiān)持不下去，如果沒(méi)有他們的鼓勵(lì)和鞭策，我想我也是很難完成此次設(shè)計(jì)的。最后我要感謝和我一組的賀睿同學(xué)，我們不僅是生活中的好兄弟，在做畢業(yè)設(shè)計(jì)時(shí)也是相互幫助，相互學(xué)習(xí)。人生路漫漫，但是我從這最后一次設(shè)計(jì)中學(xué)到了許多。愿所有的老師、同學(xué)們、朋友們合家歡樂(lè)，一生平安！參考文獻(xiàn)[1] MATLAB ：電子工業(yè)出版社，2022.[2] ：人民郵電出版社， 2022.[3] (第二版).北京：清華大學(xué)出版社，2022.[4] 樓順天. MTALAB程序設(shè)計(jì)語(yǔ)言. 西安：西安電子科技大學(xué)出版社， 1998.[5] 劉波，文忠，：電子工業(yè)出版社， 2022.33[6] ：理工大學(xué)出版社，2022.[7] 常華，袁鋼，：清華大學(xué)出版社，2022.[8] ：國(guó)防工業(yè)出版社，1979.[9] ：國(guó)防工業(yè)出版社，1980.[10] ：電子工業(yè)出版社,1996.[11] ：西安電子科技大學(xué)出版社，2022.[12] ：國(guó)防工業(yè)出版社， 2022.[13] 丁玉美，：西安電子科技大學(xué)出版社，2022.[14] Higherorder spectrum factorization in one and two dimensions with applications in sigal modeling and nonminimun phase system Symposium on Circuits and Systems,2022. [15] Antoniou Filters: Analysis and Design. McGrawHill Book Cmpany,1979.[16] Hamming R. Digital Inc,197