fig725.m

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

%《数字信号处理--原理、实现及应用》第七章 例7.2.2 程序：ep722.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月
% 用kaiser窗函数设计线性相位低通FIRDF
wc=0.5;N=30;
wn=bartlett(N+1);
hn=fir1(N,wc,wn);  %调用kaiser计算低通FIRDF的h(n)
%以下是绘图部分
M=1024;
wk=fft(wn,M);
hk=fft(hn,M);
n=0:N;
subplot(2,2,1);stem(n,wn,'.');
xlabel('n');ylabel('w_B(n)');axis([0,30,0,1.2]);
k=1:M/2;
w=2*(0:M/2-1)/M;
subplot(2,2,2);plot(w,20*log10(abs(wk(k))/max(abs(wk))));
axis([0,1,-80,5]);xlabel('ω/π');ylabel('20lg|W_B_g(ω)|');
grid on
subplot(2,2,3);stem(n,hn,'.');
xlabel('n');ylabel('h(n)');
subplot(2,2,4);plot(w,20*log10(abs(hk(k))));
axis([0,1,-40,5]);xlabel('ω/π');ylabel('20lg|H_g(ω)|');
grid on