用MATLAB设计数字FIR低通滤波器(哈明窗)

《用MATLAB设计数字FIR低通滤波器(哈明窗)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、用MATLAB设计数字FIR低通滤波器（哈明窗）设计数字FIR低通滤波器，技术指标如下：Wp=0.2pi.rp=0.25db;Ws=0.3pi,as=50db.求解，如下图：程序：wp=0.2*pi;ws=0.3*pi;trwidth=ws-wp;M=ceil(6.6*pi/trwidth)+1;%朝正无穷方向取整n=[0:1:M-1];wc=(ws+wp)/2%求的截止频率（弧度）hd=ideal_lp(wc,M)%调用计算理想低通滤波器的单位取样响应wham=(hamming(M))';%哈明窗函数h

2、=hd.*wham;[H,w]=freqz(h,1,1000,'whole');%求Z变换频率响应函数H=(H(1:1:501))';w=(w(1:1:501))'mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);grd=grpdelay(h,1,w);deltaw=2*pi/1000;Rp=-(min(db(1:1:wp/deltaw+1)))As=-round(max(db(ws/deltaw+1:1:501)))subplot(1,1,1

3、)subplot(2,2,1),stem(n,hd);title('理想脉冲响应')axis([0M-1-0.10.3]);xlabel('N');ylabel('hd(n)');subplot(2,2,2),stem(n,wham);title('Hamming窗')axis([0M-1-01.1]);xlabel('N');ylabel('w(n)');subplot(2,2,3),stem(n,h);title('实际脉冲响应')axis([0M-1-0.10.3]);xlabel('N');yla

4、bel('h(n)');subplot(2,2,4),plot(w/pi,db);title('相对标尺的幅度响应');gridon;axis([01-10010]);xlabel('frequencyinpiunits');ylabel('decibels');set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])%设置或修改LTI对象的属性值set(gca,'XTickMode','manual','XTick',[-50,0])set(gca,'XTi

5、ckLabelMode','manual','XTickLabels',['50';'0'])子程序：functionhd=ideal_lp(wc,N);a=(N-1)/2;n=[0:1:(N-1)];m=n-a+eps;hd=sin(wc*m)./(pi*m);