butter高通滤波

%butterworth高通滤波
p=ones(100);
q=ones(size(p));
d=0.5;  %决定漩涡的大小
n=2;  %决定漩涡的尖点
[f1,f2]=freqspace(size(p),'meshgrid');
r=sqrt(f1.^2+f2.^2);
for i=1:size(p,1)
    for j=1:size(p,2)
        t=r(i,j)/(d*d);
       q(i,j)=t^n/(t^n+1);
    end
end
f=p.*q;

mesh(f1,f2,f);

模拟图像代码：

%butterworth高通滤波
p=ones(100);
q=ones(size(p));
d=0.5;  %决定漩涡的大小
n=2;  %决定漩涡的尖点
[f1,f2]=freqspace(size(p),'meshgrid');
r=sqrt(f1.^2+f2.^2);
for i=1:size(p,1)
    for j=1:size(p,2)
        t=r(i,j)/(d*d);
       q(i,j)=t^n/(t^n+1);
    end
end
f=p.*q;
mesh(f1,f2,f);