% ch3example6A.m
Fs=8000; Ts=1/8000;  % 采樣率
f0=500;              % 梳狀濾波器開槽基頻率
bw=60/(Fs/2);        % 歸一化開槽帶寬
ab=-3 ;              % 計算開槽帶寬位置處的衰減分貝值
n=Fs/f0;             % 計算濾波器階數
[num,den] = iircomb(n,bw,ab, 'notch'); % 計算H(z)
freqz(num,den, 4000, 8000); % 作出H(z)的幅頻相頻圖
axis([0 4000 -30 5]);