ftranfir.m

来自「是一本教科书的各个章节的源程序代码,适合初学数字图像的人」· M 代码 · 共 21 行

    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);
  figure(1)
  plot(f,m,'--',w/pi,abs(h0))%如图4-17所示
  
  %采用频率变换法产生一个20阶的Parks-McClellan二维带通FIR滤波器
  h=ftrans2(b);
  figure(2)
  freqz2(h,[32 32]);
  
  f=imread('lena256.bmp');
  figure(3);
  imshow(f);
  fh=filter2(h,f);
  figure(4)
  imshow(fh);