fg_06_16.m

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

% fg_06_16  Resonator.
r=1-.01;
% zeros.
N=32;
th=6.5*2*pi/N;
vc=6.5/N
beta=0:N-1;
z=exp(2*pi*j*beta/N);
b=poly(z);
% poles.
a=[1 -2*r*cos(th) r^2];
% Adjust max. gain to one.
Hth=gain_f(b,a,th/(2*pi));
b=b/Hth;
[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','b');
xlabel('Hz-s'); ylabel('dB');