%用ftrans2函數設計一個通帶在0.1和0.6間（量化頻率，這里1.0對應采樣頻率的一半或π弧度）的近似圓對稱二維通帶濾波器。
%ftrans2函數由一FIR濾波器創建二維濾波器
%首先用信號處理工具箱函數remez設計一個一維FIR通帶濾波器
subplot(121)
colormap(jet(64))
b = remez(10,[0 0.05 0.15 0.55 0.65 1],[0 0 1 1 0 0]);
[H,w] = freqz(b,1,128,'whole');
plot(w/pi-1,fftshift(abs(H)));
%應用函數ftrans2和缺省 McClellan變換矩陣創建近似圓對稱濾波器
subplot(122)
h = ftrans2(b);
freqz2(h)