理想陷波滤波器

Halcom

1. function H = freqfilter_ideal_sink(M,N,u0,v0,D0)
   
2. % 理想陷波滤波器
   
3. % input:
   
4. %     M,N:频域滤波器的尺寸
   
5. %     u0,v0：频率阻止点
   
6. %     D0：带阻滤波器的截止频率
   
7. % output:
   
8. %       H：M x N的矩阵，表示频域滤波器矩阵，数据类型为double，
   
9. u = -M/2:M/2-1;
   
10. v = -N/2:N/2-1;
    
11. [U,V] = meshgrid(u,v);
    
12. D = sqrt(U.^2+V.^2);
    
13. D1 = sqrt( (U-u0).^2 + (V-v0).^2 );
    
14. D2 = sqrt( (U+u0).^2 + (V+v0).^2 );
    
15.         Mask1 = D1<=D0;
    
16.         Mask2 = D2<=D0;
    
17.     Mask = or(Mask1,Mask2);
    
18. H = ones(M,N);
    
19. H(Mask) = 0;
    
20. 
    
21. end

复制代码