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