fg_10_30.m

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

% fg_10_30       Time-freq. plane   (DCT)
sp_fig(1); p=14; nx=10; ny=9;
axes('position',[.2 .1 .6 .6]);
axis([0 nx 0 ny]);
% x-axis
set(gca,'xtick',[0:nx],'xticklabel',' ')
d=[.15 .3 .4 .4 .4 .4 .4 .4 .4 .4 .4];
a=get(gca,'xtick')-d;
b=-.4*(ones(1,nx+1));
lab=char('0','NT','2NT',' ','4NT',' ','6NT',' ','8NT',' ',' ');
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 0 0 0 0 0 0 .5];
a=-.3*(ones(1,ny))-d;
lab=char('0','\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]*ones(1,9); x=[1; 1]*[0:8];
line(x,y,'color','k');