【正文】 T。 fs=1/T。As=100。 %調(diào)用 函數(shù)確定切比雪夫2濾波器階數(shù)[B,A]=cheby2(N,Ap,Wp,39。)。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。Frequency39。)。 Ws1=*pi/T。 %采樣頻率為周期倒數(shù)Ap=1。[N,Wc]= cheb1ord(Wp,Ws,Ap, As,39。stop39。 %調(diào)用 函數(shù)設(shè)計(jì)切比雪夫Ⅰ濾波器W=linspace(0,pi,400*pi)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。)。圖316沖激響應(yīng)比變法設(shè)計(jì)的切比雪夫1型帶阻濾波器運(yùn)行波形用雙線性變換法設(shè)計(jì)切比雪夫1型帶阻濾波器，程序如下：Wp1=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。)。s39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。xlabel(39。Magnitude39。 Ws1=*pi/T。 %采樣頻率為周期倒數(shù)Ap=1。[N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。stop39。 %調(diào)用 函數(shù)設(shè)計(jì)切比雪夫2濾波器W=linspace(0,pi,400*pi)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。)。圖318脈沖響應(yīng)不變法設(shè)計(jì)的切比雪夫2型帶阻濾波器運(yùn)行波形用雙線性變換法設(shè)計(jì)設(shè)計(jì)切比雪夫2型帶阻濾波器，程序如下：Wp1=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。)。s39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。xlabel(39。Magnitude39。它的幅度平方函數(shù)為： ()式中的是雅可比橢圓函數(shù)，是與通帶衰減有關(guān)的參數(shù)，階數(shù)N等于通帶和阻帶內(nèi)最大點(diǎn)和最小點(diǎn)的總和。 Ws=*pi/T。As=100。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。)。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。Frequency/pi39。)。 %設(shè)置采樣周期為1fs=1/T。s39。,39。 %指定一段頻率值 hf=freqs(B,A,W)。 %繪出橢圓數(shù)字高通濾波器的幅頻特性曲線grid on。ylabel(39。用脈沖響應(yīng)不變法設(shè)計(jì)橢圓形帶通濾波器，程序如下：Wp1=*pi/T。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。)。s39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。xlabel(39。Magnitude39。 Ws1=*pi/T。 %采樣頻率為周期倒數(shù)Ap=1。[N,Wn]=ellipord(Wp,Ws,Ap,As,39。bandpass39。 %調(diào)用ellip函數(shù)設(shè)計(jì)橢圓濾波器W=linspace(0,pi,400*pi)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。)。圖322雙線性變換法設(shè)計(jì)的橢圓形帶通濾波器運(yùn)行波形，通帶最大衰減為1dB，阻帶最小衰減為100dB。Ws2=*pi/T。As=100。s39。,39。 %指定一段頻率值 hf=freqs(B,A,W)。 %繪出橢圓數(shù)字帶阻濾波器的幅頻特性曲線grid on。ylabel(39。 Wp2=*pi/T。 %設(shè)置采樣周期為1fs=1/T。Ws=[Ws1,Ws2]。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。)。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。Frequency/pi39。)。在完成課程設(shè)計(jì)的這段時(shí)間里，我收獲頗多。在課程過程中，總是遇到這樣或那樣的問題。畢業(yè)設(shè)計(jì)是不僅對(duì)我專業(yè)知識(shí)和專業(yè)基礎(chǔ)知識(shí)的一次實(shí)際檢驗(yàn)和鞏固，同時(shí)也是走向工作崗位前的一次熱身。為以后的工作積累了經(jīng)驗(yàn)，增強(qiáng)了信心。同時(shí)，濾波器作為數(shù)字信號(hào)處理中必不可少的重要元件仍然處于不斷完善進(jìn)步之中，其性能也在不斷增強(qiáng)。在完成課程設(shè)計(jì)的時(shí)候，我盡量的把課程設(shè)計(jì)和實(shí)際工作有機(jī)的結(jié)合起來，實(shí)踐與理論相結(jié)合。ylabel(39。 %繪出橢圓數(shù)字帶阻濾波器的幅頻特性曲線grid on。 %指定一段頻率值 hf=freqs(B,A,W)。,39。s39。As=100。Ws2=*pi/T。)。Frequency/pi39。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。Ws=[Ws1,Ws2]。 %設(shè)置采樣周期為1fs=1/T。 Wp2=*pi/T。Magnitude39。xlabel(39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。s39。)。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。圖322脈沖響應(yīng)不變法設(shè)計(jì)的橢圓形帶通濾波器運(yùn)行波形用雙線性變換法設(shè)計(jì)的橢圓形帶通濾波器，程序如下：Wp1=*pi/T。)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %調(diào)用ellip函數(shù)設(shè)計(jì)橢圓濾波器W=linspace(0,pi,400*pi)。bandpass39。[N,Wn]=ellipord(Wp,Ws,Ap,As,39。 %采樣頻率為周期倒數(shù)Ap=1。 Ws1=*pi/T。)。Frequency/pi39。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。)。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù)[B,A]=ellip(N,Ap,As,Wn,39。As=100。 Ws=*pi/T。ylabel(39。 %繪出橢圓數(shù)字高通濾波器的幅頻特性曲線grid on。 %指定一段頻率值 hf=freqs(B,A,W)。,39。s39。 %設(shè)置采樣周期為1fs=1/T。雙線性變換法設(shè)計(jì)步驟如下：（1）首先寫出該濾波器的系統(tǒng)函數(shù)（2）利用雙線性變換法轉(zhuǎn)換，數(shù)字濾波器的系統(tǒng)函數(shù)為，通帶最大衰減為1dB，阻帶最小衰減為100dB。圖319雙線性變換法設(shè)計(jì)的切比雪夫2型帶阻濾波器運(yùn)行波形 橢圓形濾波器的設(shè)計(jì)橢圓低通濾波器采用有限零點(diǎn)設(shè)計(jì)來逼近理想的低通濾波器的特性。)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %調(diào)用 函數(shù)設(shè)計(jì)切比雪夫2濾波器W=linspace(0,pi,400*pi)。stop39。[N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。 %采樣頻率為周期倒數(shù)Ap=1。 Ws1=*pi/T。Magnitude39。xlabel(39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。s39。)。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。圖317雙線性變換法設(shè)計(jì)的切比雪夫1型帶阻濾波器運(yùn)行波形用脈沖響應(yīng)不變法設(shè)計(jì)切比雪夫2型帶阻濾波器，程序如下：Wp1=*pi/T。)。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %調(diào)用 函數(shù)設(shè)計(jì)切比雪夫Ⅰ濾波器W=linspace(0,pi,400*pi)。stop39。[N,Wc]= cheb1ord(Wp,Ws,Ap, As,39。 %采樣頻率為周期倒數(shù)Ap=1。 Ws1=*pi/T。Magnitude39。xlabel(39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。s39。)。 %設(shè)置通帶最大和最小衰減Wp=[Wp1,Wp2]。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。用脈沖響應(yīng)不變法設(shè)計(jì)切比雪夫1型帶阻濾波器：程序如下：Wp1=*pi/T。ylabel(39。 grid on。 %指定一段頻率值 hf=freqs(B,A,W)。,39。s39。Ws=[Ws1,Ws2]。Ws2=*pi/T。)。Frequency39。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。)。 %調(diào)用 函數(shù)確定切比雪夫2濾波器階數(shù)[B,A]=cheby2(N,Ap,Wc,39。Ws=[Ws1,Ws2]。 fs=1/T。 Wp2=*pi/T。ylabel(39。 %繪出切比雪夫Ⅰ數(shù)字帶通濾波器的幅頻特性曲線grid on。 %指定一段頻率值 hf=freqs(B,A,W)。,39。s39。Ws=[Ws1,Ws2]。Ws2=*pi/T。)。Frequency39。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。)。 %調(diào)用 函數(shù)確定切比雪夫Ⅰ濾波器階數(shù)[B,A]=cheby1(N,Ap,Wc,39。Ws=[Ws1,Ws2]。 %設(shè)置采樣周期為1fs=1/T。 Wp2=*pi/T。Magnitude39。xlabel(39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。s39。 %設(shè)置通帶最大和最小衰減 [N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。Magnitude39。xlabel(39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。s39。 %設(shè)置通帶最大和最小衰減 [N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。 %設(shè)置歸一化通帶和阻帶截止頻率T=1。Magnitude39。xlabel(39。 %計(jì)算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B