高斯带阻滤波器

Halcom

1. function H = freqfilter_gaussian(M,N,D0,W)
   
2. % 高斯带阻滤波器
   
3. % input:
   
4. %     指定滤波器的类型’gaussian’
   
5. %     M,N:频域滤波器的尺寸
   
6. %     D0：带阻滤波器的截止频率
   
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 = 1-exp(-(1/2).*((D.^2-D0^2)./(D.*W)).^2);
    
15. end

复制代码