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

matlab fir低通滤波器 用MATLAB设计数字FIR低通滤波器(哈明窗) 用MATLAB设计数字FIR低通滤波器（哈明窗） 设计数字FIR低通滤波器，技术指标如下： Wp=0.2pi.rp=0.25db; Ws=0.3pi,as=50db. 求解，如下图： 程序： wp=0.2*pi;ws=0.3*pi; trwidth=… 1、滤波器简介 从广义上讲，任何对某些频率（相对于其他频率来说）进行修正的系统称为滤波器。严格地讲，对输入信号通过一定的处理得到输出信号，这个处理通常是提取信号中某频率范围内的信号成分，把这种处理的过程称为滤波。实现滤波处理的运算电路或设备称为滤波器… 中北大学 课程设计说明书 2010/2011 学年第（二）学期 学 专 院： 业： 信息商务学院 信息与通信工程 顾雯琪 、谢雅馨、刘晓荣、刘海亮 学 生 姓 名： 学 号01、02、06、19 信息处理实践： 低通数字 FIR … 用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=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) subplot(2,2,1),stem(n,hd);title(‘理想脉冲响应’) axis([0 M-1 -0.1 0.3]);xlabel(‘N’);ylabel(‘hd(n)’); subplot(2,2,2),stem(n,wham);title(‘Hamming窗’) axis([0 M-1 -0 1.1]);xlabel(‘N’);ylabel(‘w(n)’); subplot(2,2,3),stem(n,h);title(‘实际脉冲响应’) axis([0 M-1 -0.1 0.3]);xlabel(‘N’);ylabel(‘h(n)’); subplot(2,2,4),plot(w/pi,db);title(‘相对标尺的幅度响应’);grid on; axis([0 1 -100 10]);xlabel(‘frequency in pi units’);ylabel(‘decibels’); set(gca,’XTickMode’,’manual’,’XTick’,[0,0.2,0.3,1])%设置或修改LTI对象的属性值 set(gca,’XTickMode’,’manual’,’XTick’,[-50,0]) set(gca,’XTickLabelMode’,’manual’,’XTickLabels’,[‘50’;’0’]) 子程序： function hd=ideal_lp(wc,N); a=(N-1)/2; n=[0:1:(N-1)]; m=n-a+eps; hd=sin(wc*m)./(pi*m); 用MATLAB设计数字FIR低通滤波器（哈明窗） 设计数字FIR低通滤波器，技术指标如下： Wp=0.2pi.rp=0.25db; Ws=0.3pi,as=50db. 求解，如下图： 程序： wp=0.2*pi;ws=0.3*pi; trwidth=… 用MATLAB设计数字FIR低通滤波器（哈明窗） 设计数字FIR低通滤波器，技术指标如下： Wp=0.2pi.rp=0.25db; Ws=0.3pi,as=50db. 求解，如下图： 程序： wp=0.2*pi;ws=0.3*pi; trwidth=… 用MATLAB设计数字FIR低通滤波器（哈明窗） 设计数字FIR低通滤波器，技术指标如下： Wp=0.2pi.rp=0.25db