《用matlab设计数字fir低通滤波器(哈明窗)》由会员分享,可在线阅读,更多相关《用matlab设计数字fir低通滤波器(哈明窗)(3页珍藏版)》请在金锄头文库上搜索。
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=hd.*wham; H,w=freqz(h,1,1000,whol
2、e); %求 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
3、(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);