fg_10_31.m

英文书《Digital Signal Processing with Examples in MATLAB》附带的MATLAB实例

% fg_10_31       Time-freq. plane   (Equal bands)
sp_fig(1); p=14; nx=10; ny=9;
axes('position',[.2 .2 .6 .6]);
axis([0 nx 0 ny]);
% x-axis
set(gca,'xtick',[0:nx],'xticklabel',' ')
d=[.1 .15 .2 0 .2 0 .2 0 .2 0 0];
a=get(gca,'xtick')-d;
b=-.4*(ones(1,nx+1));
lab=char('0','T','2T',' ','4T',' ','...',' ','NT',' ',' ');
text(a,b,lab,'fontsize',p);
text(nx/3,-1.2,'Time (t)','fontsize',p);
% y-axis
set(gca,'ytick',[0:ny-1],'yticklabel',' ')
b=get(gca,'ytick');
d=[0 1.0 1.2 0 0 0 0 0 .5];
a=-.3*(ones(1,ny))-d;
lab=char('0','\pi/NT','2\pi/NT',' ',' ',' ',' ',' ','\pi/T');
text(a,b,lab,'fontsize',p);
h=text(-1,ny/3,'Frequency (\omega)');
set(h,'fontsize',p,'rotation',90);
% Horizontal lines
x=[0; 8]*ones(1,9); y=[1; 1]*[0:8];
line(x,y,'color','k');
x=[8; 10]*ones(1,9); y=[1; 1]*[0:8];
line(x,y,'color','k','linestyle','--');
% Vertical lines
y=[0; 8]; x=[8; 8];
line(x,y,'color','k');