【正文】 率采樣法偏離設計指標最明顯，阻帶衰減最小，而且設計比采用窗函數(shù)法復雜。只有適當選取過渡帶樣點值，才會取得較好的衰耗特性； （3）利用等波紋切比雪夫逼近法則的設計可以獲得最佳的頻率特性和衰耗特性，具有通帶和阻帶平坦，過渡帶窄等優(yōu)點。 綜上所述，F(xiàn)IR 濾波器很容易實現(xiàn)具有嚴格線性相位的系統(tǒng), 使信號經(jīng)過處理后不產(chǎn)生相位失真,舍入誤差小,而且穩(wěn)定，因此越來越受到廣泛的重視。 以下是一些對于動態(tài)試驗的總結：(1) 當測試現(xiàn)場環(huán)境惡劣時，除采用前置濾波器外，最好加一級數(shù)字濾波；(2) 對 FIR數(shù)字濾波器而言，選擇合適的窗函數(shù)對濾波的效果比單純增加窗寬N 顯著得多，而且不會增加計算量；(3) FFT可減小濾波過程中卷積的運算量。根據(jù)以上原理編制的濾波程序，實際應用中效果良好。與其他高級語言的程序設計相比，MATLAB環(huán)境下可以更方便、快捷地設計出具有嚴格線性相位的FIR 濾波器，節(jié)省大量的編程時間，提高編程效率，且參數(shù)的修改也十分方便，還可以進一步進行優(yōu)化設計。相信隨著版本的不斷提高，MATLAB 在數(shù)字濾波 器技術中必將發(fā)揮更大的作用。同時，用MATLAB計算有關數(shù)字濾波器的設計參數(shù)，如H(z)、h(n)等，對于數(shù)字濾波器的硬件實現(xiàn)也提供了一條簡單而 準確的途徑和依據(jù)。參考文獻[1] 孫強，運用 MATLAB實現(xiàn)數(shù)字濾波器的設計 [J]，電腦學習，2022，6(1)：2527.[2] 劉松強，數(shù)字信號處理系統(tǒng)及其應用[M]，清華大學出版社，1996.[3] 王世一，數(shù)字信號處理—修訂版[M]，北京理工大學出版社，1997.[4] 竺小松，張頌，數(shù)字FIR濾波器的設計與實現(xiàn)[J] ，設計與應用，2022，5(2) ：2224.[5] 陳桂明，應用 MATLAB語言處理數(shù)字信號與數(shù)字圖像 [M]，科學出版社，2022.[6] 鄭君里等，信號與系統(tǒng)（下冊）[M]，北京高等教育出版社，1981.[7] 胡廣書，數(shù)字信號處理—理論、算法與實現(xiàn)[M]，清華大學出版社，1997，235243.[8] 張宗橙，張玲華，曹雪虹，數(shù)字信號處理與應用[M]，東南大學出版社，1997.[9] 張志涌等，精通 MATLAB[M]，北京航空航天大學出版社，2022.[10] 張葛祥，李娜， MATLAB仿真技術與應用[M]，清華大學出版社，2022.[11] 熊欣，基于 MATLAB的FIR數(shù)字濾波器設計[J]，中國科技論文在線，2022，2(3)：3031.[12] 陳懷琛，數(shù)字信號處理教程—MATLAB釋義與實現(xiàn)[M]，電子工業(yè)出版社，2022.[13] 錢月，王春艷，數(shù)字FIR濾波器的設計與仿真實現(xiàn)[J]，商丘師范學院學報，2022，4(3)： 2325.[14] 黃文梅，熊桂林，信號分析與處理[M]，國防科技大學出版社，2022，241249.[15] 劉敏，魏玲， MATLAB通信仿真與應用[M]，國防工業(yè)出版社，2022，242246.[16] 肖忠祥，數(shù)據(jù)采集原理[M]，西北工業(yè)大學出版社，2022.[17] 馬越，閆清東，郝洪濤，基于窗函數(shù)設計的FIR 濾波器在車輛動態(tài)試驗中的應用[J] ，微計算機信息，2022，3(2)：3536.[18] 李壽柏，胡業(yè)林，MATLAB在FIR濾波器設計中的應用[J]，機電工程技術，2022，6(3)： 3839.[19] 張亞妮，基于 MATLAB 的數(shù)字濾波器設計[J]，遼寧工程技術大學學報，2022，4(3)： 3335.[20] Papadamou，Stephanos A New MATLAB—based Toolbox for Mutual Funds[M]．Proceedings of the IASTED Intemational Conference onComputers and Advanced Technology in Educmion，2022．致 謝四年的讀書生活在這個季節(jié)即將劃上一個句號，而于我的人生卻只是一個逗號，我將面對又一次征程的開始。四年的求學生涯在師長、親友的大力支持下，走得辛苦卻也收獲滿囊，在論文即將答辯之際，實在有太多的需要感謝。首先要感謝給與我?guī)椭恼撐膶?，席燕輝老師。您治學的嚴謹，專業(yè)學識的淵博，為我營造了一種良好的精神氛圍，也使我能在遇到問題時找到解決的方案。對于論文寫作的指導,經(jīng)由您悉心的點撥，再經(jīng)思考后的領悟，常常讓我能及時找到問題，及時修改不足之處。 再來，感謝我的同學們，感謝劉瑛、張競、曹濱在畢業(yè)設計這段時間對我的關心和幫助，我們大家之間的信息交流為我們在論文寫作時能有一個良好的環(huán)境，讓我明白在這互幫互助的學習氛圍中，我們都能彼此得到進步。在論文完成之際，我的心情無法平靜，從開始進入課題到論文的順利完成，有多少可敬的師長、同學、朋友給了我無言的幫助，在這里請接受我誠摯謝意！ 同時也感謝學院為我提供良好的做畢業(yè)設計的環(huán)境。 最后，再一次感謝所有在畢業(yè)設計中曾經(jīng)幫助過我的良師益友和同學，以及在設計中被我引用或參考的論著的作者。附 錄頻率取樣法設計的源程序如下：Fs=50。N=25。tao=(n1)/2。Wc=2*15/50。M=fix(wc/(2/n)+1)。Omega=[0:n1]*pi/n。Abs_h=[ones(1,m),zeros(1,n2*m1),ones(1,m1)]。H=abs_h*exp(j*tao*omega)。H=ifit(h)。Figure(1)。Subplot(2,1,1)。stem(omega,abs(h))。Subplot(2,1,2)。stem(omega,abs(h))。Hh=fit(h,512)。Hh_db=20*log10(abs(hh))。Figure(2)。Omega=[0:255]*pi/256。Set(gca,’tick’,[(120,),3,0])。T=(0。99)/fs。S=sin(2*pi*t*10)+sin(2*pi*t*10)+sin(2*pi*t20)。Sf=fft(s,512)。Figure(3)。Subplot(3,1,1)。plot(*([1:length(y)]1,y)。Subplot(3,1,3)。plot(w1,abs(yf(1:256)))。頻率采樣結構參數(shù)的源程序如下:function[c,b,a]=fitfs(h)%c:包含各并聯(lián)部分增益的行向量%b:包含按行排列的分子系數(shù)矩陣%a:包含按行排列的分母系數(shù)矩陣%h:FIR 濾波器的單位脈沖響應向量N=length(h)。H=fft(h,N)。%求脈沖響應的 DFTMagH=abs(H)。angH=angle(H)。L=(N1)/2。a1=[1,1,0]。c1=[real(H(1))]。k=[1:L]c。b=zeros(L,2)。a=ones(L,3)。%初始化矩陣 b,aa(1:L,2)=2*cos(2*pi*k/N)。a=[a。a1]。%設置 a 矩陣的第二列b(1:L,1)=cos(angH(2:L+1))c。%設置 b 矩陣的第一列b(1:L,2)=(cos((angH(2:L+1))c(2*pi*k/N)))。%設置 b 矩陣的第二列e=[2*magH(2:L+1),cl]c。%設置增益矩陣將上面程序中求得的單位脈沖響應向量 h 代入到 fitfs子函數(shù)中求得各參數(shù)如下：c= a=b=高通濾波器的設計程序主要部分如下：　　As=70。　　ωs=*π?！　ˇ豴=*π　　tr_width=ωpωs?！　?　　　　　　　　　　　　 　　M=ceil(()*2*π/(*tr_width)+1)+1?！? disp([‘濾波器的長度為 ’，num2str(M)])。　　beta=*()。 計算凱澤窗的 β 值　　n=[0:1:M1]?！　isp([’線性相位斜率為’，num2str(beta)])?！　_kai=(kaiser(M,beta))’。 求凱澤窗函數(shù)　　ωc=(ωs+ωp)/2?！　d=ideal_lp(π,M)ideal_lp(ωc,M)。 求理想脈沖響應　　h=hd.*w_kai。 設計的脈沖響應為理想脈沖響應與窗函數(shù)乘積　　[db,mag,pha,grd,ω]=freqz_m(h, [1])。　　delta_ω=2*π/1000?！　p=(min(db(ωp/delta_ω+1:1:501)))?！　isp([’實際通帶波動為’,num2str(Rp)])。 　　　以下為作圖程序　　As=round(max(db(1:1:ωs/delta_ω+1)))?！　isp([’最小阻帶衰減為’,num2str(As)])?！　?subplot(1,1,1)?！　?subplot(2,2,1)?！　?stem(n,hd)?！　?title(’理想脈沖響應 ’)。　　　 axis([0 M1 ])?！　?ylabel(’hd(n)’)?！　?subplot(2,2,2)?！　?stem(n,w_kai)?！　?title(’凱澤窗’)?！　?axis([0 M1 0 ])。　　 ylabel(’wd(n)’)。　 　 subplot(2,2,3)。　 stem(n,h)?！　?title(’實際脈沖響應 ’)?！　?axis([0 M1 ])。　　 xlabel(’n’)。ylabel(’h(n)’)?！　?subplot(2,2,4)?！　?plot(ω/π,db)?！　?title(’幅度響應 /dB’)。　　 axis([0 1 100 10])。　　 grid?！　?xlabel(’以 π 為單位的頻率’)。　　 ylabel(’分貝數(shù)/dB’)。低通濾波器的設計程序主要部分如下：ωp=*π?！　ˇ豷=*π。 　　　tr_width=ωsωp。 計算過渡帶寬　　M=ceil(*/tr_width)+1。 按漢寧窗計算濾波器長度　　disp([’濾波器的長度為’,num2str(M)])。　　n=0:M1?！　ˇ豤=(ωs+ωp)/2。 截止頻率取為兩邊緣頻率的平均值　　hd=ideal_lp(ωc,M)。 求理想脈沖響應　　w_han=(hanning(M))’。 求漢寧窗函數(shù)　　h=hd.*w_han。 設計的脈沖響應為理想脈沖響應與窗函數(shù)乘積　　[db,mag,pha,grd,ω]=freqz_m(h, [1])。 以下為作圖語句　　delta_ω=2*π/1000?！　p=(min(db(1:1: ωp/delta_ω+1)))?！　isp([’實際通帶波動為’,num2str(Rp)])。 以下為作圖程序　　As=round(max(db(ωs/delta_ω+1:1:501)))?！　isp([’最小阻帶衰減為’,num2str(As)])?！　　　ubplot(221)　　　　stem(n,hd)?！　　　itle(’理想沖擊響應 ’),　　　　axis([0 M1 ])?！　　　label(’hd(n)’)?！　　　ubplot(222)　　　　stem(n,w_han)?！　　　itle(’漢寧窗’),　　　　axis([0 M1 0 ])。　　　　ylabel(’wd(n)’)。　　　　subplot(223)　　　　stem(n,h)?！　　　itle(’實際沖擊響應 ’),　　　　axis([0 M1 ])。　　　　xlabel(’n’)?！　　　label(′h(n)′)?！　　　ubplot(224)。　　　　plot(ω/π,db)?！　　　itle(′ 幅度響應（ db)′)?！　　　xis([0 1 100 10]),　　　　grid?！　　　label(′以 π 為單位的頻率′)?！　　　label(′分貝數(shù)′