%首先，創建一個一維FIR低通濾波器，變換矩陣是可以根據需要定義的，這里定義了一個截
%斷頻
%率為0.6的十階低通濾波器和一個截斷頻率為0.5的十四階低通濾波器
b1=remez(10,[0 0.4 0.6 1],[1 1 0 0]);
[H1,w]=freqz(b1,1,64,'whole');
colormap(jet(64))
plot(w/pi-1,fftshift(abs(H1)))
b2=remez(14,[0 0.3 0.5 1],[1 1 0 0]);
[H2,w]=freqz(b2,1,64,'whole');
colormap(jet(64))
plot(w/pi-1,fftshift(abs(H2)))
%用于頻率變換的兩個一維濾波器的頻率響應見圖7-18和圖7-19
%然后，用默認的McClellan變換生成期望的近似圓對稱的濾波器
h1=ftrans2(b1);
figure,freqz2(h1,[32 32])
h2=ftrans2(b2);
figure,freqz2(h2,[32 32])