%%%% 二阶频域高滤波器
n=2;
Dcut=10; %设置剪切频率
d0 = 30; %截止频率D0
d1 = 40; %截止频率D1
u = fix(M/2);
v = fix(N/2);
for i = 1:M
for j = 1:N
d = sqrt((i-u)^2+(j-v)^2);
%%%% 计算梯形高通滤波器传递函数
BUTTERH(u,v)=1/(1+(sqrt(2)-1)*(Dcut/d)^2);
EXPOTH(u,v)=exp(log(1/sqrt(2))*(Dcut/d)^2);
if d
h=0;
elseif d<=d1
h=(d-d1)/(d0-d1);
else
h=1;
end
end
仅参考供 望斧正