fig742.m

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

%《数字信号处理--原理、实现及应用》第六章 图7.4.3绘图程序:fig743.m
% 电子工业出版社出版  高西全  丁玉美  阔永红合著 2006年8月

clear;close all
fc=1/4; fs=5/16;
%(a)
N=33;wo=[1,10];fo=[0,fc,fs,1];mo=[1,1,0,0];
hn=remez(N,fo,mo,wo);
hw=fft(hn,512);
w=0:511;w=2*w/512;
subplot(2,2,1);
plot(w,20*log10(abs(hw)));grid on;
axis([0,max(w)/2,-60,10]);
%xlabel('ω/π');
ylabel('20lg|Hg(ω)|')
title('(a) N=33,W=[1,10]')
%(b)
N=33;wo=[10,1];fo=[0,fc,fs,1];mo=[1,1,0,0];
hn=remez(N,fo,mo,wo);
hw=fft(hn,512);
w=0:511;w=2*w/512;
subplot(2,2,2);
plot(w,20*log10(abs(hw)));grid on;
axis([0,max(w)/2,-60,10]);
%xlabel('ω/π');
ylabel('20lg|Hg(ω)|')
title('(b) N=33,W=[10,1]')
%(c)
N=33;wo=[1,1];fo=[0,fc,fs,1];mo=[1,1,0,0];
hn=remez(N,fo,mo,wo);
hw=fft(hn,512);
w=0:511;w=2*w/512;
subplot(2,2,3);
plot(w,20*log10(abs(hw)));grid on;
axis([0,max(w)/2,-60,10]);
xlabel('ω/π');ylabel('20lg|Hg(ω)|')
title('(c) N=33,W=[1,1]')
%(d)
N=63;wo=[1,1];fo=[0,fc,fs,1];mo=[1,1,0,0];
hn=remez(N,fo,mo,wo);
hw=fft(hn,512);
w=0:511;w=2*w/512;
subplot(2,2,4);
plot(w,20*log10(abs(hw)));grid on;
axis([0,max(w)/2,-60,10]);
xlabel('ω/π');ylabel('20lg|Hg(ω)|')
title('(d) N=63,W=[1,1]')