ep721.m

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

%《数字信号处理--原理、实现及应用》第七章 例7.2.1 程序：ep721.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 用窗函数法设计线性相位高通FIRDF
wp=pi/2;ws=pi/4;
DB=wp-ws;           %计算过渡带宽度
N0=ceil(6.2*pi/DB); %根据表7.2.2汉宁窗计算所需h(n)长度N0
N=N0+mod(N0+1,2);   %确保h(n)长度N是奇数
wc=(wp+ws)/2/pi;    %计算理想高通滤波器通带截止频率(关于π归一化）
hn=fir1(N-1,wc,'high',hanning(N));  %调用fir1计算高通FIRDFh(n)
%以下是绘图部分
M=1024;
hk=fft(hn,M);
n=0:N-1;
subplot(2,2,1);stem(n,hn,'.');line([0,30],[0,0])
xlabel('n');ylabel('h(n)');
k=1:M/2;
w=2*(0:M/2-1)/M;
subplot(2,2,2);plot(w,20*log10(abs(hk(k))));
axis([0,1,-80,5]);xlabel('ω/π');ylabel('20lg|Hg(ω)|');
grid on