低通滤波器参数.docxVIP

0 0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 频率(kHz) 低通滤波器参数：Fs=8000, fp=2500, fs=3500, Rp=ldB, As = 30dB,其他滤波器可以通过与低通之间的 映射关系实现。 勰模拟滤波器 %巴特沃斯一一滤波器设计 wp二2*p i*2500；ws=2*pi *3500；Rp=1；As=30； [N,wc]=buttord(wp,wstRp,As,1 s')%计算率波器的阶数和3dB截止频率 [B,A]=butter(N,wc, * s') ；%计算滤波器系统函数分子分母多项式 fk=0：800/512:8000；wk=2*pi*fk； Hk=f reqs(B,A?wk)； figure plot(fk/1000,20*logl0(abs(Hk)))； grid on.xlabel (r 频率(kHz) ') .ylabel (‘忙畐度(dB)') title(,巴特沃斯模拟滤波器') axis([0t4,-35,5]) 巴特沃斯模拟滤波器 5 0 -5 -10 -15 -20 -25 -30 -35 %% %切比雪夫I一一滤波器设计 wp=2*pi*2500；ws=2*pi*3500；Rp=1；As=30； [N1, wp 1 ] =cheblord (wp, ws, Rp, Ast d) %计算切比雪夫 旳阶数和通带边界频率 [Bl,Al] =cheby 1 (Nl,Rp,wpl,' s1) ；%计算滤波器系统函数分子分母多项式 fk二0:800/512:8000；wk=2*pi*fk； Hk=freqs(Bl,Alt wk)；figure, plot(fk/1000,20*logl0(abs(Hk)))； grid on.xlabel (f 频率(kHz) ') .ylabel ('幅度(dB)') titleC切比雪夫I模拟滤波器') axis([0t4,-35,5]) %% 切比雪夫橫拟滤液器 10?25■■ 一 10 ? 25 ■ ■ 一 o 3035 ? ? 5 %切比雪夫II一一滤波器设计 wp=2*pi*2500；ws=2*pi*3500；Rp=1；As=30； [N2, wso] =cheb2ord (wp fws,Rp,As, 1 s1)%计算切比雪夫 莠的阶数和通带边界步 [B2, A2] =cheby2 (N1, Rp. wso,' s °) ； %计算滤波器系统函数分子分土 fk二0:800/512:8000；wk=2*pi*fk； Hk=f reqs(B1,A1.wk)；figure, plot(fk/1000,20*logl0(abs(Hk)))； grid on.xlabel (r 频率(kHz) *) Fylabel (r 度(dB)') title(,切比雪夫II模拟滤波器') axis([0,4,-35.5]) 切比雪天II模拟滤波器 5 0.5 1.5 2 2.5 3 3.6 硕率(kHz) ?30 i： 10 1520 ■ ■ 25 QDP) am %% %椭圆一一滤波器设计 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') ;$计算滤波器系统函数分子分母多项式 fk二0:800/512:8000；wk=2*pi*fk； Hk=freqs(Bl,Alt wk);figure, plot(fk/1000,20*logl0(abs(Hk)))； grid on.xlabel 频率(kHz) '),ylabel ('：畐度(dB)') axis([0,4,-35.5]), titleC 椭圆模拟滤波器') 椭回模拟滤液器 %计算频率响应 %计算频率响应 %求幅频特性 %化为分贝值 %检验通带波动 %检验最小阻带衰减 -脉冲响应法°) %% %数字滤波器 %脉冲响应法滤波器设计 fp=2500；fs=3500；Fs=8000； ^=2*fp/Fs；ws=2*fs/Fs；%求归一化数字通带截止频率，求归一化数字阻带起始频率 deltaw^ws-wp；%求过渡带宽 NO=ceil (6. 6/deltaw);%求窗 口长度 N=N0+mod(N0+lt2)； %确保窗口长度N为奇数 n=N-l；%求出滤波器的阶数n wi= (ws+wp)/2； %求：務的截止频 b=firl(n,wn)%利用firl函数求出滤波舉的系数 [Hk,w] = freqz(b,1)； mag = abs(Hk)； db = 20*1og10(mag/max(mag)