f=[0 0.3 0.4 0.6 0.7 1];
m=[0 0 1 1 0 0];
%產生一個20階的Parks-McClellan一維帶通FIR濾波器
b=remez(20,f,m);
[h0,w]=freqz(b,1,512);
%freqz計算由矢量a(這里a=1)和b構成的數字濾波器：
figure(1),plot(f,m,'--',w/pi,abs(h0));
%理想（虛線）和實際（實線）的頻響曲線
%頻率變換法產生20階的Parks-McClellan二維帶通FIR濾波器
h=ftrans2(b);
figure(2),freqz2(h,[32 32]);    
%利用前面設計的濾波器進行圖像分析
%讀取原始圖像
f=imread('3.jpg');
f=rgb2gray(f);
figure(3),imshow(f);           
fh=filter2(h,f);
figure(4),imshow(fh)