带阻FIR滤波器

《带阻FIR滤波器》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验总成绩：——————————————————————————装订线—————————————————————————————————报告份数：西安邮电大学数字信号处理课内上机报告专业班级：通工1002班学生姓名：王铁祥学号(班内序号):03101044（07号）2013年6月21日-4-带阻FIR滤波器Matlab代码：functionhd=ideal(w,N);alpha=(N-1)/2;n=[0:N-1];m=n-alpha+eps;hd=sin(w*m)./(pi*m);clc,clear;Wpl=0.2*pi;%根据阻带：

2、0.35~0.65pi，通带：0~0.2pi和0.8pi~piWph=0.8*pi;%确定两个通带截止频率和两个阻带截止频率Wsl=0.35*pi;Wsh=0.65*pi;dert_w=min((Wsl-Wpl),(Wph-Wsh));N=ceil(6.6*pi/dert_w);%根据过度带宽确定Nn=0:1:N-1;Wcl=(Wsl+Wpl)/2;%低通中心频率Wch=(Wsh+Wph)/2;%高通中心频率hd=ideal(pi,N)-ideal(Wch,N)+ideal(Wcl,N);%带通滤波器的原型高通+低通B=hamming

3、(N)';%海明窗-4-h=hd.*B;%加窗后[H,m]=freqz(h,[1],1024,'whole');%获取频率响应mag=abs(H);%幅值db=20*log10((mag+eps)/max(mag));%分贝数pha=angle(H);w=m/pi;%绘图subplot(2,2,1);stem(n,hd);title('理想时域波形');xlabel('n');ylabel('sa(n)');subplot(2,2,2);stem(n,B);title('海明窗');xlabel('n');ylabel('B');su

4、bplot(2,2,3);plot(w,mag);title('幅频响应');xlabel('f/Hz');ylabel('幅度');gridon;-4-axis([0101]);subplot(2,2,4);plot(w,db);title('衰竭特性');xlabel('f/Hz');ylabel('分贝数/db');axis([01-1000]);gridon实验结果：-4-——————————————————————————装订线—————————————————————————————————指导教师评语：实验成绩:指导(辅

5、导)教师:-4-