数字滤波器实验报告数字滤波器的设计实验

《数字滤波器实验报告数字滤波器的设计实验》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、数字滤波器实验报告数字滤波器的设计实验导读:就爱阅读网友为您分享以下“数字滤波器的设计实验”的资讯，希望对您有所帮助，感谢您对92to.com的支持!wm=2*pi*k・/N;Hd=Hi*exp(・j*(N・l)*wm*0・5);hh=real(ifft(Hd));wl=linspace(O，pi・O丄1000);Hl=freqz(hh，l,wl);amgH=abs(Hl);dbH=20*logl0(amgH);angH=angle(Hl);subplot(2,2,2);plot(wl/pi,dbH),gridon;title(

2、9;幅度响应');xlabel('w/pi'),ylabel('

3、H(e—w)

4、/dB')axis([01-10020]);subplot(2,2,3);plot(wl/pi,angH),gridon;title('相频特性');xlabel('w/pi'),ylabel('O')subplot(2,2,4);plot(wl/pi,amgH),gridon;title('幅频特性');xlabel('w/pi')，yl

5、abel('

6、H(e・・・w)

7、')(5)、用雷米兹(Remez)交替算法设计⑶中的滤波器，并比f=[00.380.420.780.821];较(3)、(4)、(5)三种方法的结果。程序如下:A=[001100];weigh=[l81];b=remez(40,f,A,weigh);[h,w]=freqz(b,112,'whole');maghuabs(h)；dbh==20goglo(magh)；phahuang_e(h)；subpoA29191)pof(w、(2^i)a-bh)；gridon;x_abe_

8、(m#39;w、(2^i)/#39;Ly_abe_(^#39;_H(ew)l/dB^#39;)ttHe(0#39;li酋苹甬办#39;)；axisao1.7051)subplox29192)plot(w、(2茨pi)9phah)；gridon;x_abel(^#39JW、(2诜pi)<'Ly_abe_(<'G(w)<')ttHe(<'^^^ff0#39;)；1.定性地说明用本实验程序设计的FIR滤波器的3dB截止频率在什么位置？它等于理想频率响应Hd(ejco)的截止频率吗？答：在0.3-0.5之间，大致

9、等于理想频率响应Hd(ej(o)的截止频率。2.如果没有给定h(n)的长度N,而是给定了通带边缘截止频率(0C和阻带临界频率cop,以及相应的衰减，你能根据这些条件用窗函数法设计线性相位FIR低通滤波器吗？答:能,可以根据通带边缘截止频率oc和阻带临界频率(op,以及相应的衰减与所选用的窗函数来大致判断其长度N,从而设计线性相位FIR低通滤波器。实验四、用MATLAB信号处理工具箱进行滤波器设计训练针对一个含有5Hz.50Hz和100Hz的混和正弦波信号，设计一个FIR带通滤波器并对其进行滤波。利用MATLAB实现的三种方法：程序设计法

10、.FDATool设计法和SPTool设计法。参数要求：釆样频率fs=100Hz,通带下限截止频率fcl=10Hz,通带上限截止频率fc2=20Hz,过渡带宽6Hz,通阻带波动0.01,采用凯塞窗设计。方法1:程序设计法fsamp=100;fcuts=[4102026];mags=[010];devs=[0.010.010.01];[n9wn9beta,ftype]=kaiserord(fcuts9mags^devs9fsamp);n=n+rem(n,2);hh=firl(n,wn?ftype,kaiser(n+l9beta),'

11、noscale');[H，f]=freqz(hh，l，1024，fsamp);figure(l)plot(f,abs(H)),gridon;xlabel('f/Hz')，ylabel('

12、H(ew)

13、'),title('kaiser数字带通滤波器');用一个序列测试滤波器性能figure(2)N=256;n=O:1:N-1;t=n/fsamp;xn=sin(2*pi*5*t)+sin(2*pi*15*t)+sin(2*pi*35*t);subplot(2,2，l)stem(n,xn

14、，'.');gridon;xlabel('n');ylabel('x(n)');title('序列x(n)');axis([0256-22])