fg_06_12.m

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

% fg_06_12  First digital design example.
[d,c]=bw_analog_weights(1,6,tan(pi*.3));
[b,a]=bilin(d,c);

% Modified for short course:
[b,a]=bw_weights(1,6,.1);

[dB,nu]=power_gain(b,a);
sp_fig(1,9,4);
subplot(1,2,1);
pz_plot(b,a);
set(gca,'fontsize',16);
set(gca,'position',[0.1 0.1709 0.4-.1 0.7876]);
subplot(1,2,2);
axis([0 .5 -40 0]); grid on;
set(gca,'fontsize',16);
set(gca,'position',[0.5780 0.2 0.3270 0.93-.2])
line(nu,dB,'color','r');
xlabel('Hz-s'); ylabel('dB');