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