zinc.m

来自「高阶sigma-delta调制器设计matlab工具包, 半波带滤波器设计工具」· 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;