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

文档介绍：数字滤波器实验报告数字滤波器的设计实验导读:就爱阅读网友为您分享以下“数字滤波器的设计实验”的资讯,希望对您有所帮助,的支持!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('幅度响应');xlabel('w/pi'),ylabel('|H(e—w)|/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'),ylabel('|H(e・・・w)|')(5)、用雷米兹(Remez)交替算法设计⑶中的滤波器,并比f=[];较(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_(m#39;w、(2^i)/#39;Ly_abe_(^#39;_H(ew)l/dB^#39;)ttHe(0#39;li酋苹甬办#39;);)subplox29192)plot(w、(2茨pi)9phah);gridon;x_abel(^#39JW、(2诜pi)<'Ly_abe_(<'G(w)<')ttHe(<'^^^ff0#39;);?它等于理想频率响应Hd(ejco)的截止频率吗?答:-,大致等于理想频率响应Hd(ej(o)的截止频率。(n)的长度N,而是给定了通带边缘截止频率(0C和阻带临界频率cop,以及相应的衰减,你能根据这些条件用窗函数法设计线性相位FIR低通滤波器吗?答:能,可以根据通带边缘截止频率oc和阻带临界频率(op,以及相应的衰减与所选用的窗函数来大致判断其长度N,从而设计线性相位FIR低通滤波器。实验四、,设计一个FIR带通滤波器并对其进行滤波。利用MATLAB实现的三种方法:。参数要求:釆样频率fs=100Hz,通带下限截止频率fcl=10Hz,通带上限截止频率fc2=20Hz,过渡带宽6Hz,,采用凯塞窗设计。方法1:程序设计法fsamp=100;fcuts=[4102026];mags=[010];devs=[];[n9wn9beta,ftype]=kaiserord(fcuts9mags^devs9fsamp);n=n+rem(n,2);hh=firl(n,wn?ftype,kaiser(n+l9beta),'noscale');[H,f]=freqz(hh,l,1024,fsamp);figure(l)plot(f,abs(H)),gridon;xlabel('f/Hz'),ylabel('|H( ew)|'),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,'.');gridon;xlabel('n');ylabel('x(n)');title('序列x(n)');axis([0256-22]);X=fft(xn^N