hdplot.m

MATLAB及在电子信息课程中的应用例题索引与程序

%《MATLAB在电子信息课程中的应用》第七章通用绘图程序Hdplot
% 幅频特性绘制
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 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