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