zinc.m

来自「这是一个新的用于图像处理的工具箱」· M 代码 · 共 15 行

function mag = zinc(f,m,n)
% mag = zinc(f,m=64,n=1)	Calculate the magnitude response
% of a cascade of n mth-order comb filters at frequencies f.

if nargin<3
	n = 1;
    if nargin<2
	    m = 64;
    end
end

mag = ones(size(f));
nonzero = mod(f,1) ~=0;
mag = abs( sinc(m*f) ./ sinc(f) ).^n;