文档介绍:数字滤波器实验报告数字滤波器的设计实验
导读:就爱阅读网友为您分享以下“数字滤波器的设计实验”的资讯,希望对您有所帮助,的支持!
wm = 2*pi*k./N;
Hd = Hr.*exp(-j*(N-1)*wm*);
hh = real(ifft(Hd));
w1 = linspace(0,pi-,1000);
H1 = freqz(hh,1,w1);
amgH = abs(H1);
dbH = 20*log10(amgH);
angH = angle(H1);
subplot(2,2,2);
plot(w1/pi,dbH),grid on;
title('幅度响应');
xlabel('w/pi'),ylabel('|H( e…w)|/dB')
axis([0 1 -100 20]);
subplot(2,2,3);
plot(w1/pi,angH),grid on;
title('相频特性');
xlabel('w/pi'),ylabel('Φ')
subplot(2,2,4);
plot(w1/pi,amgH),grid on;
title('幅频特性');
xlabel('w/pi'),ylabel('|H( e…w)|')
(5)、用雷米兹(Remez)交替算法设计(3)中的滤波器,并比较(3)、(4)、(5)三种方法的结果。程序如下:
f=[0 1];
A=[0 0 1 1 0 0];
weigh=[1 8 1];
b=remez(40,f,A,weigh);
%
[h,w]=freqz(b,1,512,'whole');
magh=abs(h);
dbh=20*log10(magh);
phah=angle(h);
subplot(2,1,1)
plot(w/(2*pi),dbh);grid on;
xlabel('w/(2*pi)'),ylabel('|H( e…w)|/dB'),title('幅频特性'); axis([0 1 -70 5])
subplot(2,1,2)
plot(w/(2*pi),phah);grid on;
xlabel('w/(2*pi)'),ylabel('Φ(w)'),title('相频特性');
四、实验思考
1. 定性地说明用本实验程序设计的FIR滤波器的3dB截止频率在什么位置?它等于理想频率响应Hd(ejω)的截止频率吗?
答:-,大致等于理想频率响应Hd(ejω)的截止频率。
2. 如果没有给定h(n)的长度N,而是给定了通带边缘截止频率ωc和阻带临界频率ωp,以及相应的衰减,你能根据这些条件用窗函数法设计线性相位FIR低通滤波器吗?
答:能,可以根据通带边缘截止频率ωc和阻带临