ep722.m

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

%《数字信号处理--原理、实现及应用》第七章 例7.2.2 程序：ep722.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 用kaiser窗函数设计线性相位低通FIRDF
fp=1500;fs=2500;rs=40;
wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;
Bt=ws-wp;           %计算过渡带宽度
alph=0.5842*(rs-21)^0.4+0.07886*(rs-21);  %根据(7.2.16)式计算kaiser窗的控制参数α
N=ceil((rs-8)/2.285/Bt); %根据(7.2.17)式计算kaiser窗所需阶数N
wc=(wp+ws)/2/pi;    %计算理想高通滤波器通带截止频率(关于π归一化）
hn=fir1(N,wc,kaiser(N+1,alph));  %调用kaiser计算低通FIRDF的h(n)
%以下是绘图部分
M=1024;
hk=fft(hn,M);
n=0:N;
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