巴特沃斯陷波滤波器

Halcom

1. function H = freqfilter_btw_sink(M,N,u0,v0,D0,n)
   
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. 
    
16. H = 1./(1+ (D0^2./(D1.*D2)).^n );
    
17. 
    
18. end

复制代码