【正文】 byshev I型和Chebyshev II型數(shù)字低通濾波器，要求通帶邊界頻率fp=，通帶最大衰減Rp=；阻帶邊界頻率fs=8kHZ，阻帶最小衰減Rs=30dB,采樣頻率為Fs=20kHZ。Wp=2*Flp/Fs。 %阻帶最小衰減（單位：dB）[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs)。plot(w/pi,20*log10(abs(hw)))。)。切比雪夫I型幅頻響應(yīng)39。grid on。ylabel(39。 %抽樣頻率20KHzFlp=2100。 %歸一化的阻帶截止頻率Rp=。 %返回H(z)的分子分母系數(shù)[hw,w]=freqz(b,a)。xlabel(39。幅度（dB）39。subplot(2,1,2)。ω/π39。)。橢圓濾波器的選擇性相對(duì)前三種是最好的，下降斜度比較大，通帶和阻帶內(nèi)均為等波紋，同樣的性能指標(biāo)，橢圓濾波器可以最低的階數(shù)來(lái)實(shí)現(xiàn)。FDATool可以設(shè)計(jì)幾乎所有的常規(guī)濾波器，包括FIR和IIR的各種設(shè)計(jì)方法。在界面的下半部，用來(lái)設(shè)置濾波器的設(shè)計(jì)參數(shù)；另一部分則是特性區(qū)，在界面的上半部分，用來(lái)顯示濾波器的各種特性。Filter Order(濾波器階數(shù))選項(xiàng)，定義濾波器的階數(shù)，包括Specify Order(指定階數(shù))和Minimum Order(最小階數(shù))。它的具體選項(xiàng)由Filter Type選項(xiàng)和Design Method選項(xiàng)決定。例如設(shè)計(jì)帶通濾波器時(shí)，可以定義Wstop1(頻率Fstop1處的幅值衰減)、Wpass (通帶范圍內(nèi)的幅值衰減)、Wstop2(頻率Fstop2處的幅值衰減)。要得到幅頻、相頻響應(yīng)特性，運(yùn)算量也是很大的。 然后在Frequency Specifications 中選擇Unit 為Hz，給出采樣頻率Fs=1000，；最后在Magnitude Specifications 中選擇Unit 為dB，Apass=1。而且我們?cè)谧詈蟮姆抡婺Ｐ驮O(shè)計(jì)中將載入這一模塊進(jìn)行信號(hào)濾波。首先建立信號(hào)傳輸仿真模塊，先打開(kāi) SIMULINK 工具箱，并且建立一個(gè) Model，在這個(gè)空白 Model 中進(jìn)行環(huán)節(jié)庫(kù)及框圖的建立，在 SIMULINK 菜單下找到 Source，雙擊 Source 圖標(biāo)，將正弦信號(hào)源和常量信號(hào)源拉到 Model 中，然后分別找到 DSP 離散正弦信號(hào)源,乘法 器和加法器以及示波器，然后把需要用來(lái)設(shè)計(jì) IIR 數(shù)字濾波器的模塊 FDATOOL 都拉到 Model 中，把環(huán)節(jié)都布好后，把各環(huán)節(jié)的端口按框圖連接起來(lái)。 總 結(jié)通過(guò)這幾個(gè)實(shí)驗(yàn)，對(duì)設(shè)計(jì)數(shù)字濾波器的整個(gè)過(guò)程有了很好的掌握。致 謝本論文是在雷老師的親切關(guān)懷和悉心指導(dǎo)下完成的。在此，我還要感謝在一起愉快的度過(guò)畢業(yè)設(shè)計(jì)生活的同學(xué)，正是由于你們的幫助和支持，我才能克服一個(gè)一個(gè)的困難和疑惑，直至本文的順利完成。設(shè)計(jì)階數(shù)較高的IIR 濾波器時(shí)，計(jì)算量更大，設(shè)計(jì)過(guò)程中要改變參數(shù)或?yàn)V波器類型時(shí)都要重新計(jì)算。本文利用Matlab的信號(hào)處理工具及其SIMULINK仿真工具包箱能夠方便快捷地設(shè)計(jì)和實(shí)現(xiàn)各種濾波器，使信號(hào)波形更加直觀，并且用FDATOOL工具箱以及專門(mén)用于濾波器設(shè)計(jì)來(lái)設(shè)計(jì)IIR數(shù)字濾波器，充分利用了MATLAB 的交互性好的特點(diǎn)，而且最后用到了SIMULINK 仿真工具箱來(lái)實(shí)現(xiàn)信號(hào)傳輸和濾波器的設(shè)計(jì)，將濾波器的設(shè)計(jì)置于一個(gè)新的平臺(tái)，這對(duì)于研究信號(hào)的傳輸和處理有著極其重要的作用。在這種情況下，濾波器設(shè)計(jì)就要進(jìn)行大量復(fù)雜的運(yùn)算，單純的靠公式計(jì)算和編制簡(jiǎn)單的程序很難在短時(shí)間內(nèi)完成。在論文即將完成之際，我的心情無(wú)法平靜，從開(kāi)始進(jìn)入課題到論文的順利完成，有多少可敬的師長(zhǎng)、同學(xué)、朋友給了我無(wú)言的幫助，在這里請(qǐng)接受我誠(chéng)摯的謝意!最后我還要感謝培，謝謝你們!參考文獻(xiàn)[1] 倪養(yǎng)華．?dāng)?shù)字信號(hào)處理與實(shí)現(xiàn)．上海：上海交通大學(xué)出版社，1998[2] 飛思科技產(chǎn)品研發(fā)中心．MATLAB7輔助信號(hào)處理技術(shù)與應(yīng)用．北京：電子工業(yè)出版社， 2005[3] 樓順天，李博苗．基于MATLAB的系統(tǒng)分析與設(shè)計(jì)一信號(hào)處理．西安：西安電子科技大學(xué)出版社，1998.[4] 張葛祥，李娜．MATLAB仿真技術(shù)與應(yīng)用．北京：清華大學(xué)出版社，2003[5] 胡廣書(shū) 數(shù)字信號(hào)處理、理論、算法與實(shí)現(xiàn)[M] ．北京：清華大學(xué)出版社，1997[6] 陳希林，肖明清．一種LabWindows/CVI與MATLAB混合編程的實(shí)現(xiàn)方法[J].微計(jì)算機(jī)信息，2005[7]劉波．MATLAB信號(hào)處理．北京：電子工業(yè)出版社，2006[8]施陽(yáng)等．MATLAB語(yǔ)言工具箱．西安：西北工業(yè)大學(xué)出版社，1999[9]丁玉美．?dāng)?shù)字信號(hào)處理．西安電子科技大學(xué)出版社 第二版結(jié)束語(yǔ)本文在對(duì)濾波器實(shí)際設(shè)計(jì)中用到了多種設(shè)計(jì)方法，與其他高級(jí)語(yǔ)言的程序設(shè)計(jì)相比，MATLAB 環(huán)境下可以更方便、快捷地設(shè)計(jì)出具有IIR 數(shù)字濾波器，節(jié)省大量的編程時(shí)間，提高編程效率，且參數(shù)的修改也十分方便，還可以進(jìn)一步進(jìn)行優(yōu)化設(shè)計(jì)。從課題的選擇到項(xiàng)目的最終完成，雷老師都始終給予我細(xì)心的指導(dǎo)和不懈的支持。通過(guò)對(duì)數(shù)字帶通濾波器的設(shè)計(jì)，熟悉了MATLAB的運(yùn)行環(huán)境，初步掌握了MATLAB語(yǔ)言在數(shù)字信號(hào)處理中一些基本庫(kù)函數(shù)的調(diào)用和編寫(xiě)基本程序等應(yīng)用；熟悉了濾波器設(shè)計(jì)的一般原理，對(duì)濾波器有了一個(gè)感性的認(rèn)識(shí)；學(xué)會(huì)了數(shù)字高通濾波器設(shè)計(jì)的一般步驟；加深了對(duì)濾波器設(shè)計(jì)中產(chǎn)生誤差的原因以及雙線性變換法優(yōu)缺點(diǎn)的理解和認(rèn)識(shí)。圖38 仿真模型圖39 濾波效果 結(jié)果顯示，仿真模型中顯示 40HZ,75HZ,100HZ 三個(gè)不同頻率的正弦信相加后的波形，干擾后的混合信號(hào)是一個(gè)夾雜多頻率的信號(hào)，我們所需要的有用信號(hào)（40HZ）在傳輸過(guò)程中不可避免的受到其他頻率信號(hào)干擾，我們的目的就是得到所需要的 40HZ 頻段信號(hào)， 盡可能濾除干擾信號(hào)，濾波效果的圖中就是濾波后信號(hào)與原始抽樣信號(hào)的的波形圖，經(jīng)濾波后，信號(hào)周期約為 ,即頻率為 40HZ，說(shuō)明濾波效果很好。它可以處理的系統(tǒng)包括：線性、非線性系統(tǒng)；離散、連續(xù)及混合系統(tǒng)；單任務(wù)、多任務(wù)離散事件系統(tǒng)。通過(guò)菜單選項(xiàng)Analysis 可以在特性區(qū)看到所設(shè)計(jì)的幅頻響應(yīng)、相頻響應(yīng)、沖擊響應(yīng)和零極點(diǎn)配置等特性。第四節(jié) FDATOOL設(shè)計(jì)IIR數(shù)字濾波器我們將以一個(gè)IIR 濾波器的設(shè)計(jì)實(shí)例來(lái)具體說(shuō)明使用MATLAB 工具箱的方便。 Window Specifications選項(xiàng)，當(dāng)選取采用窗函數(shù)設(shè)計(jì)，該選項(xiàng)可定義，它包含了各種窗函數(shù)。采用窗函數(shù)設(shè)計(jì)濾波器時(shí)，由于過(guò)渡帶是由窗函數(shù)的類型和階數(shù)所決定，所以只需定義通帶截止頻率，而不必定義阻帶參數(shù)。如果選擇Minimum Order，則MATLAB根據(jù)所選擇的濾波器類型自動(dòng)使用最小階數(shù)。Design Method(設(shè)計(jì)方法)選項(xiàng)，包括IIR濾波器的Butterwotth(巴特沃思)法、ChebyshevFDATool。第三節(jié) FDATool介紹和界面設(shè)計(jì)FDATool(Filter Designamp。切比雪I型濾波器的幅頻特性在通帶內(nèi)有波動(dòng)，阻帶內(nèi)單調(diào)。ylabel(39。grid on。切比雪夫II型幅頻響應(yīng)39。)。plot(w/pi,20*log10(abs(hw)))。 %阻帶最小衰減（單位：dB）[N,Wn]=cheb2ord(Wp,Ws,Rp,Rs)。Wp=2*Flp/Fs。)。ω/π39。subplot(2,1,2)。幅度（dB）39。xlabel(39。 %返回H(z)的分子分母系數(shù)[hw,w]=freqz(b,a)。 %歸一化的阻帶截止頻率Rp=。 %抽樣頻率20KHzFlp=2100。Fs缺省時(shí)默認(rèn)為2Hz。如果要將歸一化頻率轉(zhuǎn)換為Hz,則將歸一化頻率乘以采樣頻率的一半。(2)、有關(guān)濾波器設(shè)計(jì)當(dāng)中的頻率歸一化問(wèn)題:信號(hào)處理工具箱中經(jīng)常使用的頻率是Nyqu