%《MATLAB在電子信息課程中的應(yīng)用》第七章通用繪圖程序Hdplot
% 幅頻特性繪制
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clf
wc1=0.1;wc2=0.5;wp=0.2;ws=0.3
w=[0,wc1,wc1,wc2,wc2,pi];
Hd=[0,0,1,1,0,0];
subplot(2,1,1),plot(w,Hd);hold on
axis([0,1,0,1.2])
set(gca,'xTick',[0,1])
w1=[0,wp];w2=[wp+0.1,pi];
Hd11=[1.1,1.1];Hd12=[0.9,0.9];
Hd2=[0.1,0.1];
subplot(2,1,2)
plot(w1,Hd11);hold on
plot(w1,Hd12)
fill([0,wp,wp,0,0],[0.9,0.9,1.1,1.1,0.9],'y'),hold on
plot(w1,[1,1])
plot(w2,Hd2);hold on
plot([wp,wp],[0,1.2],'--')
fill([ws,1,1,ws,ws],[0,0,0.1,0.1,0],'y')
set(gca,'xTick',[0,1])
axis([0,1,0,1.2])
gtext('w(pi)'),gtext('w(pi)')
hold off