fsamfir.m

在一维信号中

 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所示