fg_06_18.m

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

% fg_06_18   Comb filter z-plane and dB plots.
N=32; r=.99; s=[2:5];
x=ones(1,201);
[y,b,a]=comb(N,s,r,x);
vr=s/N

% Pole-zero and power gain plots.
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]);
set(gca,'fontsize',16);
for n=1:length(s),
   [dB,v]=power_gain(b(n,:),a(n,:),1000);
   line(v,dB,'color','b'); grid on;
   xlabel('Hz-s'); ylabel('Power gain in dB');
end
set(gca,'fontsize',16);
set(gca,'position',[0.58 0.15 0.327 0.94-.15])