DSP实验报告
实验四：
IIR数字滤波器的设计 数字滤波器的设计
实验五：FIR
实验四、IIR数字滤波器的设计
【1】 fc = 0.3kHz，delta = 0.8dB，fr = 0.2kHz，At = 20dB,T = 1ms:设计
一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。 实验结果：
实验代码：
>> Wc = 2*pi*300/1000; Wr = 2*pi*200/1000; Rp = 0.8; Rs = 20;
[N,Wc] = cheb1ord(Wc,Wr,Rp,Rs,'s'); [B,A] = cheby1(N,Rp,Wc,'high','s'); omega = [0:pi/1000:pi];
h = freqs(B,A,omega); gain = 20*log10(abs(h));
plot(omega/(2*pi/1000),gain); 结果分析：
由实验所得关于设计的滤波器的增益曲线来看，当f<200Hz时，衰减大于20dB，当f>300Hz时，衰减趋近于零，满足设计参数要求。
【5】利用双线性变换法设计满足下列指标的切比雪夫型数字带阻滤波 器，并作图验证设计结果：当1kHz<=f<=2kHz时，At>=18dB;当f<=500Hz以及f>=3kHz时，delta<=3dB;采样频率fs = 10kHz。 实验结果
程序代码：
>> W1 =2*10000*tan(2*pi*500/(2*10000));
W2 =2*10000*tan(2*pi*1000/(2*10000)); W3 =2*10000*tan(2*pi*2000/(2*10000)); W4 =2*10000*tan(2*pi*3000/(2*10000)); Wp = [W2,W3]; Ws = [W1,W4];
[N,Wn] = cheb1ord(Wp,Ws,3,18,'s'); [B,A] = cheby1(N,3,Wn,'stop','s'); [num,den] = bilinear(B,A,10000); [h,w] = freqz(num,den); f = w/pi*5000;
plot(f,20*log10(abs(h))); axis([0,3500,-100,10]); 结果分析：
根据设计要求，取要求中的参数值为极限值，所得滤波器增益曲线如上图。由图可知当频率在1kHz到2kHz之间时，增益衰减大于18dB，当频率小于500或大于3000Hz时，增益略小于1，通带波动delta小于3dB,满足设计需求。
实验五、FIR数字滤波器的设计
【1】N = 45，计算并画出矩形窗、汉明窗、布莱克窗的归一化幅度谱，并比较各自的主要特点。
