4、Q
5、<
6、
7、77(Q)
8、QS<
9、Q
10、<^(1)通带边缘频率Cp；(2)阻带边缘频率(3)通带起伏耳；(4)通带峰值起伏aP=-201og10(l-Jp)[d3],(5)阻带起伏6,最小阻带衰减^=-201og10(<5J[dS]
11、o三.预习要求1、在MATLAB屮，熟悉函数firl>kaiscrord、rcmczord>rcmcz的使用；B二firl(n,Wn/high,,,noscale，)设计滤波器;[n,Wn,beta,ftype]=kaiserord(f,a,dev)估计滤波器阶数；[n,fo,ao,w]=remezord(f,a,dev,fs)计算等波纹滤波器阶数n和加权函数w(3)；B=rcmcz(n,f,a)进行等波纹滤波器的设计。2、阅读扩展练习中的实例，学习FIR滤波器的设计方法及其在MATLAB中的实现;3、给出FIR数字滤波器的冲激响应，绘出它们
12、的幅度和相位频响曲线，讨论它们各自的实现形式和特点。数字滤波器有IIR和FIR两种类型，它们的特点和设计方法不同。四、实验内容:1、用凯塞窗设计一FIR低通滤波器，通带边界频率Op=0.3^,阻带边界频率2=0.5龙，阻带衰减加r不小于50dB。解:首先由过渡带宽和阻带衰减皿來决定凯塞窗的N和0上图给出了以上设计的频率特性，(a)为N=30直接截取的频率特性(b)为凯塞窗设计的频率特性。凯塞窗设计对应的MATLAB程序为：wn=kaiser(30,4.55);nn=[0:1:29];alfa=(30-l)/2;hd=sin(0.4*pi*(n
13、n-alfa))./(pi*(nn-alfa));h=hd.*wn;[hl,wl]=freqz(h,l);或者:b=firl(29,0.4,kaiser(30,4.55));[hl,wl]=freqz(b,l);plot(wl/pi,20*log10(abs(h1)));axis(
14、0,1,-80,10]);grid；xlabelf归一化频率伉);ylabel('幅度/dBJ;还可以使用[n,Wn,beta,ftypeJ=kaiserord(f,a,dev)函数来估计滤波器阶数等，得到凯塞窗滤波器：fcuts=f0.30.51;%归一化频率o
15、mega/pimags=[10];devs=[0.0510A(-2.5)];fn,Wn,beta,ftype]=kaiserord(fcuts,mags,devs);%计算til凯塞窗N,beta的值hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale,);freqz(hh);2、利用MATLAB编程，分别用窗函数法和等波纹滤波器法设计两种FIR数字滤波器，指标要求如下：通带边缘频率：Qp1=0.45,t?QP2=0.65^,通带峰值起伏：ap16、带衰减：as>4C[dB]o1.用窗函数法实现：调用函数[n,wn,bta,ftype]=kaiserord(f,a,dev,fs)参数：f=[0.30.450.650.8]为对应数字频率=0.3^,=0.45^,Qp2=0・65龙，，Q52=0.75龙4二[010]为由f指定的各个频带上的幅值向量，一-般只有0和1表示；和f长度关系为(2*3的长度)一2二(f的长度)devs二[0・010.10870.01]用于指定各个频带输出滤波器的频率响应与其期里幅值之间的最大输岀误差或偏差，长度与a相等，计算公式：阻带衰减误差二為“0亠叫，通带波动衰
17、减误差丛=1-10"°fs缺省值为211Z>>[n,wn,bta,ftype]=kaiserord([0.30.450.650.8],[010],[0.010.10870.01]);%用kaiserord函数估计出滤波器阶数n和beta参数>>hl=firl(n,wn,ftype,kaiser(n+1,bta),'noscale');>>[hhl,wl]=freqz(hl,1,256);>>figurc(l)>>subplot(2,1,1)>>plot(wl/pi,20*logl0(abs(hhl)))>>grid»xlabel('归一化频率
18、w');ylabel('幅度/db');>>subplot(2,1,2)>>plot(wl/pi,angle(hhl))>>grid>>xlabelC归一化频率w'
