【正文】 有關(guān)的參數(shù)，階數(shù)N等于通帶和阻帶內(nèi)最大點和最小點的總和。階數(shù)和系統(tǒng)函數(shù)由通帶的截止頻率，通帶內(nèi)最大衰減，阻帶截止頻率及阻帶內(nèi)的最小衰減等參數(shù)確定。雙線性變換法設(shè)計步驟如下：（1）首先寫出該濾波器的系統(tǒng)函數(shù)（2）利用雙線性變換法轉(zhuǎn)換，數(shù)字濾波器的系統(tǒng)函數(shù)為，通帶最大衰減為1dB，阻帶最小衰減為100dB。用脈沖響應(yīng)不變法設(shè)計橢圓形高通濾波器，程序如下：Wp=*pi/T。 Ws=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置采樣周期為1fs=1/T。 %采樣頻率為周期倒數(shù)Ap=1。As=100。 %設(shè)置通帶最大和最小衰減 [N,Wn]=ellipord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。high39。,39。s39。)。 %調(diào)用ellip函數(shù)設(shè)計橢圓濾波器W=linspace(0,pi,400*pi)。 %指定一段頻率值 hf=freqs(B,A,W)。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %繪出橢圓數(shù)字高通濾波器的幅頻特性曲線grid on。xlabel(39。Frequency/pi39。)。ylabel(39。Magnitude39。)。圖320脈沖響應(yīng)不變法設(shè)計的橢圓形高通濾波器運行波形用雙線性變換法設(shè)計的橢圓形高通濾波器，程序如下：Wp=*pi/T。 Ws=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置采樣周期為1fs=1/T。 %采樣頻率為周期倒數(shù)Ap=1。As=100。 %設(shè)置通帶最大和最小衰減 [N,Wn]=ellipord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。high39。,39。s39。)。 %調(diào)用ellip函數(shù)設(shè)計橢圓濾波器W=linspace(0,pi,400*pi)。 %指定一段頻率值 hf=freqs(B,A,W)。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %繪出橢圓數(shù)字高通濾波器的幅頻特性曲線grid on。xlabel(39。Frequency/pi39。)。ylabel(39。Magnitude39。)。圖321雙線性變換法設(shè)計的橢圓形高通濾波器運行波形，，通帶最大衰減為1dB，阻帶最小衰減為100dB。用脈沖響應(yīng)不變法設(shè)計橢圓形帶通濾波器，程序如下：Wp1=*pi/T。 Wp2=*pi/T。 Ws1=*pi/T。Ws2=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置采樣周期為1fs=1/T。 %采樣頻率為周期倒數(shù)Ap=1。As=100。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。Ws=[Ws1,Ws2]。[N,Wn]=ellipord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。bandpass39。,39。s39。)。 %調(diào)用ellip函數(shù)設(shè)計橢圓濾波器W=linspace(0,pi,400*pi)。 %指定一段頻率值 hf=freqs(B,A,W)。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %繪出橢圓數(shù)字帶通濾波器的幅頻特性曲線grid on。xlabel(39。Frequency/pi39。)。ylabel(39。Magnitude39。)。圖322脈沖響應(yīng)不變法設(shè)計的橢圓形帶通濾波器運行波形用雙線性變換法設(shè)計的橢圓形帶通濾波器，程序如下：Wp1=*pi/T。 Wp2=*pi/T。 Ws1=*pi/T。Ws2=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置采樣周期為1fs=1/T。 %采樣頻率為周期倒數(shù)Ap=1。As=100。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。Ws=[Ws1,Ws2]。[N,Wn]=ellipord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。bandpass39。,39。s39。)。 %調(diào)用ellip函數(shù)設(shè)計橢圓濾波器W=linspace(0,pi,400*pi)。 %指定一段頻率值 hf=freqs(B,A,W)。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %繪出橢圓數(shù)字帶通濾波器的幅頻特性曲線grid on。xlabel(39。Frequency/pi39。)。ylabel(39。Magnitude39。)。圖322雙線性變換法設(shè)計的橢圓形帶通濾波器運行波形，，通帶最大衰減為1dB，阻帶最小衰減為100dB。用沖擊響應(yīng)不變法設(shè)計橢圓形帶阻濾波器，程序如下：Wp1=*pi/T。 Wp2=*pi/T。 Ws1=*pi/T。Ws2=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置采樣周期為1fs=1/T。 %采樣頻率為周期倒數(shù)Ap=1。As=100。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。Ws=[Ws1,Ws2]。[N,Wn]=ellipord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。stop39。,39。s39。)。 %調(diào)用ellip函數(shù)設(shè)計橢圓濾波器W=linspace(0,pi,400*pi)。 %指定一段頻率值 hf=freqs(B,A,W)。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %繪出橢圓數(shù)字帶阻濾波器的幅頻特性曲線grid on。xlabel(39。Frequency/pi39。)。ylabel(39。Magnitude39。)。圖323脈沖響應(yīng)不變法設(shè)計的橢圓形帶阻濾波器運行波形雙線性變換法設(shè)計橢圓型帶阻濾波器，程序如下：Wp1=*pi/T。 Wp2=*pi/T。 Ws1=*pi/T。Ws2=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置采樣周期為1fs=1/T。 %采樣頻率為周期倒數(shù)Ap=1。As=100。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。Ws=[Ws1,Ws2]。[N,Wn]=ellipord(Wp,Ws,Ap,As,39。s39。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。stop39。,39。s39。)。 %調(diào)用ellip函數(shù)設(shè)計橢圓濾波器W=linspace(0,pi,400*pi)。 %指定一段頻率值 hf=freqs(B,A,W)。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %繪出橢圓數(shù)字帶阻濾波器的幅頻特性曲線grid on。xlabel(39。Frequency/pi39。)。ylabel(39。Magnitude39。)。圖324雙線性變換法設(shè)計的橢圓形帶阻濾波器運行波形總 結(jié)課程設(shè)計是我們作為學(xué)生在學(xué)習(xí)階段的最后一個環(huán)節(jié)，是對所學(xué)基礎(chǔ)知識和專業(yè)知識的一種綜合應(yīng)用，是一種綜合的再學(xué)習(xí)、再提高的過程，這一過程對學(xué)生的學(xué)習(xí)能力和獨立思考及工作能力也是一種培養(yǎng)，同時畢業(yè)設(shè)計的水平也反映了大學(xué)教育的綜合水平。在完成課程設(shè)計的時候，我盡量的把課程設(shè)計和實際工作有機的結(jié)合起來，實踐與理論相結(jié)合。這樣更有利于自己能力的提高。在完成課程設(shè)計的這段時間里，我收獲頗多。通過這次課程設(shè)計使我了解到MATLAB功能強大，使用方便，已經(jīng)廣泛地應(yīng)用在濾波器的設(shè)計與仿真等各個方面，并且其仍然處于不斷的發(fā)展之中，功能不斷增強，應(yīng)用更為開放。同時，濾波器作為數(shù)字信號處理中必不可少的重要元件仍然處于不斷完善進步之中，其性能也在不斷增強。課程設(shè)計過程中，不僅鞏固了我的基礎(chǔ)理論知識，而且使我各個方面的能力都有很大的提高。在課程過程中，總是遇到這樣或那樣的問題。有時發(fā)現(xiàn)一個問題的時候，需要做大量的工作，花大量的時間才能解決。為以后的工作積累了經(jīng)驗，增強了信心。既讓我們懂得了怎樣把理論應(yīng)用于實際，又讓我們懂得了在實踐中遇到的問題怎樣用理論去解決，它更是自己綜合運用所學(xué)知識,發(fā)現(xiàn)、提出、分析和解決實際問題,鍛煉實踐能力的重要環(huán)節(jié)。畢業(yè)設(shè)計是不僅對我專業(yè)知識和專業(yè)基礎(chǔ)知識的一次實際檢驗和鞏固，同時也是走向工作崗位前的一次熱身。參考文獻[1] ：上海交通大學(xué)出版社，1998[2] ．北京：電子工業(yè)出 [3] ．西安：西安電子科技大學(xué).[4] ．北京：清華大學(xué)出版社，2003[5] 、理論、算法與實現(xiàn)[M] ．北京：清華大學(xué)出版社，1997[6] 陳希林．一種LabWindows/CVI與MATLAB混合編程的實現(xiàn)方法[J].微計算機信息，2005[7] ．北京：電子工業(yè)出版社，2006[8] ．西安：西北工業(yè)大學(xué)出版社，1999[9] ．西安電子科技大學(xué)出版社 第三版，200843