给定参数N,用频率采样法设计FIR带阻滤波器。.docVIP

西安邮电大学 数字信号处理课内上机报告 专业班级： 学生姓名： 学号(班内序号): 2013 年 06 月 19 日 四、给定参数N，用频率采样法设计FIR带阻滤波器。 clear all; N = 65; Hg = [ones(1,12),zeros(1,8),ones(1,25),zeros(1,8),ones(1,12)];%用于逼近的理想滤波器 Hg(1,11) = 0.5886;Hg(1,12) = 0.1065;Hg(1,20) = 0.1065;Hg(1,21) = 0.5886;%增加过度带优化点 % Hg(1,32) = 0.5886;Hg(1,33) = 0.1065;Hg(1,44) = 0.1065;Hg(1,45) = 0.5886; k = 0:N-1; A = [exp(-j*pi*k*(N-1)/N)];%滤波器相位特性 Hk = Hg.*A;%对理想滤波器采样 hn = ifft(Hk);%计算滤波器的单位取样响应 [H,W] = freqz(hn);%计算滤波器的频率响应 n = 0:N-1; figure(1); plot(n,Hg); axis([0 32 0 1]);grid on; figure(2);