fg_06_13.m

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

% fg_06_13  Second digital design example.
[d,c]=bw_analog_weights(3,6,tan(pi*.3),tan(pi*.4));
[b,a]=bilin(d,c);

% Modified for short course:
[b,a]=ch_weights(3,6,-30,.1,.2);


[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.15 0.45-.1 0.94-.15])

subplot(1,2,2);
axis([0 .5 -40 0]); grid on;
set(gca,'fontsize',16);
set(gca,'position',[0.58 0.15 0.327 0.94-.15])
line(nu,dB,'color','r');
xlabel('Hz-s'); ylabel('dB');