matlab三大低通滤波器的简单使用

下边为三个滤波器的使用，截至频率半径值分别设置30，90，160，，先对灰度图像进行加高斯噪声，然后进行滤波处理

F=fftshift(F);%将变换的原点移到频率矩形的中心  

[M,N]=size(f);  
  
%截至半径30理想低通滤波    
h1=zeros(M,N);  
for i=1:M  
    for j=i:N  
        if(sqrt(((i-M/2)^2+(j-N/2)^2))<30)  
            h1(i,j)=1;  
          
        end  
    end  
end  
G1=F.*h1;  
G1=ifftshift(G1);  
g1=real(ifft2(G1)); %截至半径90理想低通滤波    
h2=zeros(M,N);  
for i=1:M  
    for j=i:N  
        if(sqrt(((i-M/2)^2+(j-N/2)^2))<90)  
            h2(i,j)=1;  
          
        end  
    end  
end  
G2=F.*h2;  
G2=ifftshift(G2);  
g2=real(ifft2(G2)); %截至半径160理想低通滤波