fg_07_15.m

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

% fg_07_15         Waveform segments.
sp_fig(1);
set(gcf,'defaultaxesfontsize',16);
set(gcf,'defaultaxesvisible','off');

N=2^10;
randn('seed',1234);
xr=randn(1,N);
[b,a]=ch_weights(1,8,-40,.15);
y=filters(b,a,xr);
x=y/max(abs(y));

subplot(11,1,1);
axis([0 N -1 1]); 
line(0:N-1,x,'color','r');
h=text(-60,0,'x_k');
set(h,'fontsize',16);
%p1=get(gca,'position')

for p=0:1,
   subplot(11,1,3+p);
   axis([0 N -1 1]);
%  p34=get(gca,'position')
   n=[0:N/2-1]+p*N/2;
   line(n,x(n+1),'color','b');
   h=text(-60,0,['x_{',int2str(p+1),'k}']);
   set(h,'fontname','times','fontsize',16);
end

for p=0:4,
   subplot(11,1,6+p);
   axis([0 N -1 1]);
%  p610=get(gca,'position')
   n=[0:N/2-1]+p*N/8;
   y=x(n+1);
   line(n,x(n+1),'color','g');
   h=text(-60,0,['x_{',int2str(p+1),'k}']);
   set(h,'fontname','times','fontsize',16);
end

set(gcf,'defaultaxesvisible','on');
p=get(gca,'position');
dy=.0757; y1=.8963; y0=y1-10*dy;
axes('position',[p(1) y0 p(3) y1-y0],'color','none');
axis([0 1024 0 10]); grid on;
set(gca,'xtick',[0:128:1024]);
set(gca,'ytick',[0:1:10]);
set(gca,'yticklabel',' ');
xlabel('k');
%h=text(-50,10,'x(kT)')
%set(h,'fontname','times','fontsize',24);