带阻FIR滤波器.docVIP

西安邮电大学 数字信号处理课内上机报告 专业班级： 通工1002班 学生姓名： 王铁祥 学号(班内序号): 07号） 2013 年 6 月 21 日 带阻FIR滤波器 Matlab代码： function hd = 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;%根据阻带：0.35~0.65pi，通带：0~0.2pi和0.8pi~pi Wph = 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);%根据过度带宽确定N n = 0:1:N-1; Wcl = (Wsl+Wpl)/2;%低通中心频率 Wch = (Wsh+Wph)/2;%高通中心频率 hd = ideal(pi,N)-ideal(Wch,N)+ideal(Wcl,N);%带通滤波器的原型 高通+低通 B = hamming(N);%海明窗 h = hd.*