f=[0 0.1 0.2 0.3 0.4 0.6 0.7 0.8 0.9 1];
 m=[1 1 1 1 0 0 1 1 1 1];

%產生一個20階的具有最小二乘線性相位的一維FIR帶阻濾波器。
b=firls(20, f , m, 'Hilbert');
[h0,w]=freqz(b,1,512,2);
plot(f,m,'--',w, abs(h0));
   %如圖4-21所示, 給出了理想(虛線)和實際(實線)的頻率響應曲線
  
%  在一維濾波器函數取理想帶阻濾波器時，采用頻率抽樣法產生一個20階的具有最小   %  二乘線性相位的FIR帶阻濾波器
  h1=fsamp2(f,f,m,[5 5]);
  figure(2)
  freqz2(h1,[32 32]);  %如圖4-22所示

%  在一維濾波器函數取實際帶阻濾波器時，采用頻率抽樣法產生一個20階的具有最小   %  二乘線性相位的FIR帶阻濾波器
  ff=1/512:1/512:1;
  h2=fsamp2(ff,ff,h0',[5 5]);
  figure(3)
  freqz2(h2,[32 32]);  %如圖4-23所示

  %利用前面設計的濾波器進行圖像分析
  % 讀取原始圖像
  f=imread('lena256.bmp');
  %如圖4-19所示
  
  %利用h1進行濾波
  fh1=filter2(h1,f);
  figure(4)
  imshow(fh1); %如圖4-24所示

  %利用h2進行濾波
  fh2=filter2(h2,f);
  figure(5)
  imshow(fh2); %如圖4-24所示