理想带阻滤波器

Halcom

1. function H = freqfilter_ideal(M,N,D0,W)
   
2. % 理想带阻滤波器
   
3. % input:
   
4. %     M,N:频域滤波器的尺寸
   
5. %     D0：带阻滤波器的截止频率
   
6. %     W:带宽
   
7. % output:
   
8. %       H：M x N的矩阵，表示频域滤波器矩阵，数据类型为double，
   
9. 
   
10. u = -M/2:M/2-1;
    
11. v = -N/2:N/2-1;
    
12. [U,V] = meshgrid(u,v);
    
13. D = sqrt(U.^2+V.^2);
    
14.         H = double( or( D<(D0-W/2), D>(D0+W/2) ) );
    
15. end
    

复制代码