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