资源描述:
《哈工程数字信号试验六》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数字信号处理实验试验六:离散时间滤波器设计班级:学号:姓名:指导老师:日期:一、IIR数字滤波器设计1.实验原理2.程序试验1clearwp=6*pi*10^3;ws=9*pi*10^3;ap=1;as=15;%Ä£ÄâÔÐ͵ÍͨÂ˲¨Æ÷¼¼ÊõÖ¸±êFs=30*10^3;%²ÉÑùƵÂÊwp1=wp/Fs;ws1=ws/Fs;%Êý×ÖƵÂÊ[N,Wc]=cheb1ord(wp,ws,ap,as,'s');%È·¶¨½×ÊýNºÍ3dB½ØֹƵÂÊWc[b,a]=cheby1(N,ap,Wc,'low','s');%µ÷ÓÃÇбÈÑ©·ò1ÐÍÂ˲
2、¨Æ÷[bz,az]=impinvar(b,a,Fs);%Âö³åÏìÓ¦²»±ä·¨ÊµÏֱ任w0=[wp1,ws1];hx=freqz(bz,az,w0);[H,W]=freqz(bz,az);%ÇóƵÂÊÏìÓ¦dbHx=-20*log10(abs(hx)/max(abs(H)));%Çówp1£¬ws1¶ÔӦ˥¼õplot(W,abs(H));xlabel('Ïà¶ÔƵÂÊ');ylabel('·ùƵ');%»æÖÆͼÏñgrid试验2clearwp=1.4*pi*10^3;ws=pi*10^3;ap=1;as=32;%Ä£ÄâÔÐ͸ßͨÂË
3、²¨Æ÷¼¼ÊõÖ¸±êFs=2*10^3;%²ÉÑùƵÂÊwp1=wp/Fs;ws1=ws/Fs;%Êý×ÖƵÂÊomp1=2*Fs*tan(wp1/2);omp2=2*Fs*tan(ws1/2);%Ä£ÄâÔÐ͵ÍͨƵÂÊÔ¤»û[N,Wc]=cheb1ord(omp1,omp2,ap,as,'s');%È·¶¨½×ÊýNºÍ3dB½ØֹƵÂÊWc[b,a]=cheby1(N,ap,Wc,'high','s');%µ÷ÓÃÇбÈÑ©·ò1ÐÍÂ˲¨Æ÷[bz,az]=bilinear(b,a,Fs);%Ë«ÏßÐÔ·¨ÊµÏֱ任w0=[wp1,
4、ws1];hx=freqz(bz,az,w0);[H,W]=freqz(bz,az);%ÇóƵÂÊÏìÓ¦dbHx=-20*log10(abs(hx)/max(abs(H)));%Çówp1£¬ws1¶ÔӦ˥¼õplot(W,abs(H));xlabel('Ïà¶ÔƵÂÊ');ylabel('·ùƵ');%»æÖÆͼÏñgrid1.实验结果试验1试验2一、窗函数法设计FIR数字滤波器1.实验原理滤波器的理想频率响应函数为Hd(ejω),则其对应的单位脉冲响应为 窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼hd(n)。由于hd(n)往往是无限长序列
5、,且是非因果的,所以用窗函数。w(n)将hd(n)截断,并进行加权处理: h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数H(ejω)为 。用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。2.程序试验1clearN=26;wc=0.4;h1=fir1(25,wc,boxcar(N))[H,W]=freqz(h1,1);subplot(221);plot(W/pi,abs(H));title('¾ØÐÎÕñ·ùÌØÐÔ/dB');xlabel
6、('Ïà¶ÔƵÂÊ');ylabel('H(w)');h2=fir1(25,wc,hanning(N))[H,W]=freqz(h2,1);subplot(222);plot(W/pi,abs(H));title('ººÄþÕñ·ùÌØÐÔ/dB');xlabel('Ïà¶ÔƵÂÊ');ylabel('H(w)');h3=fir1(25,wc,blackman(N))[H,W]=freqz(h3,1);subplot(223);plot(W/pi,abs(H));title('²¼À³¿ËÂüÕñ·ùÌØÐÔ/dB');xlabel('Ïà¶ÔƵÂÊ');ylabe
7、l('H(w)');clearN1=26;wc=0.4;h1=fir1(25,wc,boxcar(N1))[H,W]=freqz(h1,1);subplot(221);plot(W/pi,abs(H));title('Õñ·ùÌØÐÔ/dB');xlabel('Ïà¶ÔƵÂÊ');ylabel('H(w)');N2=51;h2=fir1(50,wc,boxcar(N2))[H,W]=freqz(h2,1);subplot(222);plot(W/pi,abs(H));title('Õñ·ùÌØÐÔ/dB');xlabel('