1、低通濾波器參數：Fs=8000 , fp=2500 , fs=3500 ,Rp=1dB , As=30dB,其他濾波器可以通過與低通之間的 映射關系實現。%嗾擬濾波器%巴特沃斯一觥波器設計wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;N,wc=buttord(wp,ws,Rp,As,s ) %計算率波器的階數和3dB截止頻率B,A=butter(N,wc, s );%計算濾波器系統函數分子分母多項式fk=0:800/512:8000;wk=2*pi*fk;Hk=freqs(B,A,wk);figureplot(fk/1000,20*log10(abs(Hk);gri

2、d on ,xlabel( 頻率(kHz ) ),ylabel(幅度(dB) title(巴特沃斯模擬濾波器) axis(0,4,-35,5)巴特沃斯模擬謔波器500爪511.522533.54頻率(kHz)%切比雪夫I一觥波器設計wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;-1D-15如-30-35N1,wp1=cheb1ord(wp,ws,Rp,As,s ) %計算切比雪夫濾波器的階數和通帶邊界頻B1,A1=cheby1(N1,Rp,wp1,s );%計算濾波器系統函數分子分母多項式fk=0:800/512:8000;wk=2*pi*fk;Hk=freqs(B

3、1,A1,wk);figure,plot(fk/1000,20*log10(abs(Hk);grid on ,xlabel( 頻率(kHz ) ),ylabel(幅度(dB)title(切比雪夫I模擬濾波器)axis(0,4,-35,5)%切比雪關睇揪旎波券切比雪關睇揪旎波券%切比雪夫II一瞬波器設計wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;N2,wso=cheb2ord(wp,ws,Rp,As,s ) %計算切比雪夫濾波器的階數和通帶邊界頻率B2,A2=cheby2(N1,Rp,wso,s );%計算濾波器系統函數分子分母多項式fk=0:800/512:800

4、0;wk=2*pi*fk;Hk=freqs(B1,A1,wk);figure,plot(fk/1000,20*log10(abs(Hk);grid on ,xlabel( 頻率(kHz ) ),ylabel(幅度(dB) title(切比雪夫II模擬濾波器) axis(0,4,-35,5)城幸（kHz）切比雷夫切比雷夫II模狙濾渺幕模狙濾渺幕%圓一觥波器設計wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;N,wpo=ellipord(wp,ws,Rp,As,s )%計算濾波器的階數和通帶邊界頻率B,A=ellip(N,Rp,As,wpo,s );%計算濾波器系統函數分

5、子分母多項式fk=0:800/512:8000;wk=2*pi*fk;Hk=freqs(B1,A1,wk);figure,plot(fk/1000,20*log10(abs(Hk);grid on ,xlabel( 頻率(kHz ) ),ylabel(幅度(dB)axis(0,4,-35,5),title(橢圓模擬濾波器)%數字濾波器%脈沖響應法濾波器設計fp=2500;fs=3500;Fs=8000;wp=2*fp/Fs;ws=2*fs/Fs; deltaw=ws-wp;%求過渡帶寬N0=ceil(6.6/deltaw);%窗 口長度N=N0+mod(N0+1,2);%確保窗口長度N為奇數n

6、=N-1; %求出濾波器的階數nwn=(ws+wp)/2;球濾波器的截止頻率b=fir1(n,wn)%用fir1函數求出濾波器的系數Hk,w = freqz(b,1); mag = abs(Hk); db =20*log10(mag/max(mag); dw =pi/512; %關于pi歸一化Rp = -(min(db(1:wp*pi/dw+1) As =-(max(db(ws*pi/dw+1:512)figure,plot(0:pi/511:pi,db),gridaxis(0,4.0,-80,5),title( %求歸一化數字通帶截止頻率，求歸一化數字阻帶起始頻率%計算頻率響應%求幅頻特性%

7、化為分貝值%檢驗通帶波動%檢驗最小阻帶最減on數字濾波器一驪沖響應法)舞字流濮舞字流濮器器. .一府一府同成涇同成涇%fir1窗函數法fp=2500;fs=3500;Fs=8000;rs=30;wp=2*fp*pi/Fs;ws=2*fs*pi/Fs;始頻率Bt=ws-wp;%求過渡帶寬alpha=0.5842*(rs-21)A0.4+0.07886*(rs-21);M=ceil(rs-8)/2.285/Bt);球出濾波器的階數wc=(ws+wp)/2/pi;%求歸一化數字通帶截止頻率，求歸一化數字阻帶起%計算kaiser窗的控制參數hk=fir1(M,wc,kaiser(M+1,alpha)H

8、k,w = freqz(hk,1); mag = abs(Hk); db= 20*log10(mag/max(mag); db1=db;figure,plot(0:pi/511:pi,db1),gridaxis(0,4.0,-80,5),title( %求濾波器的截止頻率并關于pi咖用fir1%歸一化函數求出濾波器的系數計算頻率響應%求幅頻特性%化為分貝值on數字濾波器一fr1窗函數法)%頻率采樣法fp=2500;fs=3500;Fs=8000;rs=30;wp=2*fp*pi/Fs;ws=2*fs*pi/Fs;始頻率Bt=ws-wp;%熒過渡帶寬m=1;alpha=0.5842*(rs-21

9、)A0.4+0.07886*(rs-21);N=ceil(m+1)*2*pi/Bt;%求出濾波器的階數N=N+mod(N+1,2); Np=fix(wp/(2*pi/N); Ns=N-2*Np-1;Hk=ones(1,Np+1),zeros(1,Ns),ones(1,Np); wc=(ws+wp)/2/pi;%濾波器的截止頻率并關于hk=fir1(M,wc,kaiser(M+1,alpha) Hk,w = freqz(hk,1); mag =abs(Hk); db = 20*log10(mag/max(mag); db1=db;figure,plot(0:pi/511:pi,db1),grid

10、 axis(0,4.0,-80,5),title(%求歸一化數字通帶截止頻率，求歸一化數字阻帶起pi咖用firl%計算kaiser窗的控制參數歸一化函數求出濾波器的系數計算頻率響應%求幅頻特性%化為分貝值on數字濾波器一嗷率采樣法)%用等波紋最佳逼近法設計FIR數字濾波器Fs=8000;f=2500,3500;m=1,0;rp=1;rs=30;delta1=(10A(rp/20)-1)/(10A(rp/20)+1);delta2=10A(-rs/20);rip=delta1,delta2;M,fo,mo,w=remezord(f,m,rip,Fs);姻界頻率為模擬頻率時必須加入采樣頻率M=M+1;%古算的M直達不W要求，家hn=remez(M,fo,mo,w);Hk,w =