ep723.m

高西全2009年最新版<数字信号处理>源码程序

%《数字信号处理--原理、实现及应用》第七章 例7.2.3 程序：ep723.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 调用窗函数法设计一个线性相位FIR带通滤波器
close all;clear;
wlp=0.2*pi;wls=0.35*pi;wus=0.65*pi;wup=0.8*pi;  %设计指标参数赋值
B=wls-wlp;          %过渡带宽度
N=ceil(12*pi/B);    %计算阶数N,ceil(x)为大于等于x的最小整数
wp=[(wls+wlp)/2/pi,(wus+wup)/2/pi];     %设置理想带通截止频率
hn=fir1(N,wp,'stop',blackman(N+1));     %带阻滤波器要求h(n)长度为奇数，所以取N+1
%一下为绘图部分
M=1024;
hk=fft(hn,M);
n=0:N;
subplot(2,2,1);stem(n,hn,'.');
xlabel('n');ylabel('h(n)');
grid on
k=1:M/2+1;
w=2*(0:M/2)/M;
subplot(2,2,2);plot(w,20*log10(abs(hk(k))));
axis([0,1,-90,5]);xlabel('ω/π');ylabel('20lg|Hg(ω)|');
grid on