fg_10_23.m

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

% fg_10_23       Decimation & Expansion Spectra
ylab=['|L^1(\omega)|';'|H^1(\omega)|'];
sp_fig(1); p=14;
for i=1:2,
	b=fir_weights(31,i,5,.25);
	H=abs(fft(b,801));
	subplot(3,2,i+2);
	line(0:800,H,'color','k'); grid on;
	axis([0 800 0 1.1]);
	set(gca,'xticklabel',' ');
	set(gca,'yticklabel','0||1.0','fontname','times','fontsize',p);
	dx2=[10 50 40 70 60];
	a2=get(gca,'xtick')-dx2;
	b2=-.13*(ones(1,5));
	xlab2=char('0','\pi/2T','\pi/T','3\pi/2T','2\pi/T');
	text(a2,b2,xlab2,'fontname','times','fontsize',p);
	ylabel(ylab(i,:),'fontname','times','fontsize',p);
	xlabel('\omega (rad/s)');
end