基于matlab的傅里叶频域滤波

clear;clc;figure;  
P=imread（‘zjy.jpg‘）; 
I=rgb2gray（P）;
figureimshow（I）;title（‘a图 原图‘）;  
f=double（I）;[rc]=size（f）;  
F=fft2（f）;G=fftshift（F）;  
d0=15; %半径范围  
n=2;%巴特沃斯阶次  
a=0.5;b=2.0; %高频强调滤波传递函数系数  
mu=floor（r/2）;mv=floor（c/2）;  
for u=1:r  
    for v=1:c  
        d=sqrt（（u-mu）^2+（v-mv）^2）;  
        Hlpbtw=1/（1+0.414*（d/d0）^（2*n））;  
        Hhpbtw=1-Hlpbtw;    Ghpbtw（uv）=Hhpbtw*G（