【正文】 %調(diào)用 ellip 函數(shù)設(shè)計橢圓濾波器 W=linspace(0,pi,400*pi)。s39。bandpass39。)。 [N,Wn]=ellipord(Wp,Ws,Ap,As,39。 %設(shè)置通帶最大和最小衰減 Wp=[Wp1,Wp2]。 %采樣頻率為周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截止頻率 T=1。 Ws1=*pi/T。 用 脈沖響應(yīng)不變法設(shè)計橢圓形帶通濾波器， 程序 如下： Wp1=*pi/T。)。 ylabel(39。Frequency/pi39。 %繪出橢圓數(shù)字高通濾波器的幅頻特性曲線 grid on。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。 %指定一段頻率 值 hf=freqs(B,A,W)。)。,39。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù) [B,A]=ellip(N,Ap,As,Wn,39。s39。 As=100。 %設(shè)置采樣周期為 1 fs=1/T。 Ws=*pi/T。)。 ylabel(39。Frequency/pi39。 %繪出橢圓數(shù)字高通濾波器的幅頻特性曲線 grid on。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。 %指定一段頻率值 hf=freqs(B,A,W)。)。,39。 %調(diào)用函數(shù)確定橢圓濾波器階數(shù) [B,A]=ellip(N,Ap,As,Wn,39。s39。 As=100。 %設(shè)置采樣周期為 1 fs=1/T。 Ws=*pi/T。 雙線性變換法 設(shè)計步驟如下： （ 1）首先寫出該濾波器的系統(tǒng)函數(shù) RCaas asH 1,)( ??? （ 2）利用雙線性變換法轉(zhuǎn)換，數(shù)字濾波器的系統(tǒng)函數(shù) )(1 zH 為 基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計 34 12111121 1 )1(|)()(11 ????? ?????? zazasHzHzzTsa 22,2 21 ????? aTaTaaTaTa 數(shù)字高通濾波器 的設(shè)計 ??? ， ??? ，通帶最大衰減為 1dB，阻帶最小衰減為 100dB。它的幅度平方函數(shù)為： )(11)(222 ????Na JjH ? () 式中的 ()NJ ? 是雅可比橢圓函數(shù)， ? 是與通帶衰減有關(guān)的參數(shù)，階數(shù) N 等于通帶和阻帶內(nèi)最大點和最小點的總和。 圖 319雙線性變換法設(shè)計的切比雪夫 2型帶阻濾波器運行波形 橢圓形 濾波器的設(shè)計 橢圓低通濾波器采用有限零點設(shè)計來逼近理想的低通濾波器的特性。Magnitude39。)。 xlabel(39。 %返回頻率響應(yīng) 基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計 33 plot(W/pi,abs(Hz))。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。 %調(diào)用 函數(shù)設(shè)計切比雪夫 2 濾波器 W=linspace(0,pi,400*pi)。s39。stop39。)。 [N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。 %設(shè)置通帶最大和最小衰減 Wp=[Wp1,Wp2]。 %采樣頻率為周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截止頻率 T=1。 Ws1=*pi/T。 圖 318脈沖響應(yīng)不變法設(shè)計的切比雪夫 2型帶阻濾波器運行波形 用雙線性 變換法設(shè)計設(shè)計切比雪夫 2 型 帶阻濾波器， 程序 如下： Wp1=*pi/T。Magnitude39。)。 基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計 32 xlabel(39。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %計算模擬濾波 器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。 %調(diào)用 函數(shù)設(shè)計切比雪夫 2 濾波器 W=linspace(0,pi,400*pi)。s39。stop39。)。 [N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。 %設(shè)置通帶最大和最小衰減 Wp=[Wp1,Wp2]。 %采樣頻率為 周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截止頻率 T=1。 Ws1=*pi/T。 圖 317雙線性變換法設(shè)計的切比雪夫 1型帶阻濾波器運行波形 用 脈沖響應(yīng)不變法設(shè)計切比雪夫 2 型 帶阻濾波器， 程序 如下： Wp1=*pi/T。Magnitude39。)。 xlabel(39。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。 %調(diào)用 函數(shù)設(shè)計切比雪夫Ⅰ濾波器 W=linspace(0,pi,400*pi)。s39。stop39。)。 [N,Wc]= cheb1ord(Wp,Ws,Ap, As,39。 %設(shè)置通帶最大和最小衰減 Wp=[Wp1,Wp2]。 %采樣頻率為周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截 止頻率 T=1。 Ws1=*pi/T。 基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計 30 圖 316沖激響應(yīng)比變法設(shè)計的切比雪夫 1型帶阻濾波器運行波形 用 雙線性變換法設(shè)計切比雪夫 1 型 帶阻濾波器， 程序 如下： Wp1=*pi/T。Magnitude39。)。 xlabel(39。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。 %調(diào)用 函數(shù)設(shè)計切比雪夫Ⅰ濾波器 W=linspace(0,pi,400*pi)。s39。stop39。)。 [N,Wc]= cheb1ord(Wp,Ws,Ap, As,39。 %設(shè)置通帶最大和最小衰減 Wp=[Wp1,Wp2]。 %采樣頻率為周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截止頻率 T=1。 Ws1=*pi/T。 用 脈沖響應(yīng)不變法設(shè)計 切比雪夫 1 型 帶阻濾波器： 程序如下 ： Wp1=*pi/T。)。 ylabel(39。Frequency39。 grid on。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。 %指定一段頻率值 hf=freqs(B,A,W)。)。,39。 %調(diào)用 函數(shù)確定切比雪夫 2 濾波器階數(shù) [B,A]=cheby2(N,Ap,Wp,39。s39。 As=100。 Ws=[Ws1,Ws2]。 fs=1/T。 Ws2=*pi/T。 Wp2=*pi/T。)。 ylabel(39。Frequency39。 grid on。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。 %指定一段頻率值 hf=freqs(B,A,W)。)。,39。 %調(diào)用 函數(shù)確定切比雪夫 2 濾波器階數(shù) [B,A]=cheby2(N,Ap,Wc,39。s39。 Ws=[Ws1,Ws2]。 As=100。 fs=1/T。 Ws2=*pi/T。 Wp2=*pi/T。)。 ylabel(39。Frequency39。 %繪出切比雪夫Ⅰ數(shù)字帶通濾波器的幅頻特性曲線 grid on。 %調(diào)用雙線性變換法 Hz=freqz(D,C,W)。 %指定一段頻率值 hf=freqs(B,A,W)。)。,39。 %調(diào)用 函數(shù)確定切比雪夫Ⅰ濾波器階數(shù) [B,A]=cheby1(N,Ap,Wc,39。s39。As=100。 Ws=[Ws1,Ws2]。 fs=1/T。 Ws2=*pi/T。 Wp2=*pi/T。)。 ylabel(39。Frequency39。 %繪出切比雪夫Ⅰ數(shù)字帶通濾波器的幅頻特性曲線 grid on。 %調(diào)用脈沖響應(yīng)不變法 Hz=freqz(D,C,W)。 %指定一段頻率值 hf=freqs(B,A,W)。)。,39。 %調(diào)用 函數(shù)確定切比雪夫Ⅰ濾波器階數(shù) [B,A]=cheby1(N,Ap,Wc,39。s39。 Ws=[Ws1,Ws2]。 As=100。 %設(shè)置采樣周期為 1 fs=1/T。 Ws2=*pi/T。 Wp2=*pi/T。 圖 311雙線性變換法設(shè)計的切比雪夫 2型濾波器運行波形 數(shù)字帶通濾波器的 設(shè)計 1 ??? ， 1 ??? ， 2 ??? ， 2 ??? ，通帶最大衰減為 1dB，阻帶最小衰減為 100dB。Magnitude39。)。 xlabel(39。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。 %調(diào)用 函數(shù)設(shè)計切比雪夫 2 濾波器 W=linspace(0,pi,400*pi)。s39。)。 %設(shè)置通帶最大和最小衰減 [N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。 %采樣頻率為周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截止頻率 基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計 24 T=1。 圖 310脈沖響應(yīng)不變法設(shè)計的切比雪夫 2型濾波器運 行波形 用雙線性 變換法設(shè)計切比雪夫 2 型高通 濾波器 ， 程序 如下： Wp=*pi/T。Magnitude39。)。 xlabel(39。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。 %調(diào)用 函數(shù)設(shè)計切比雪夫 2 濾波器 W=linspace(0,pi,400*pi)。s39。)。 %設(shè)置通帶最大和最小衰減 [N,Wc]= cheb2ord(Wp,Ws,Ap, As,39。 %采樣頻率為周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截止頻率 T=1。 圖 39雙線性變換法設(shè)計的切比雪夫 1型濾波器運行波形 用 脈沖響應(yīng)不變法設(shè)計切比雪夫 2 型 高通濾波器 ， 程序 如下： 基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計 23 Wp=*pi/T。Magnitude39。)。 xlabel(39。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=bilinear(B,A,fs)。 %調(diào)用 函數(shù)設(shè)計切比雪夫Ⅰ濾波器 W=linspace(0,pi,400*pi)。s39。high39。)。 %設(shè)置通帶最大和最小衰減 [N,Wc]= cheb1ord(Wp,Ws,Ap, As,39。 %采樣頻率為周期倒數(shù) Ap=1。 %設(shè)置歸一化通帶和阻帶截止頻率 T=1。 圖 38脈沖響應(yīng)不變法 設(shè)計的切比雪夫 1型濾波器運行波形 基于 MATLAB 的 IIR 數(shù)字濾波器的設(shè)計 22 用 雙線性變換法設(shè)計切播學府 1 型高通 濾波器， 程序 如下： Wp=*pi/T。Magnitude39。)。 xlabel(39。 %返回頻率響應(yīng) plot(W/pi,abs(Hz))。 %計算模擬濾波器的幅頻響應(yīng) [D,C]=impinvar(B,A,fs)。 %調(diào)