sinc_decimate.m

来自「高阶sigma-delta调制器设计matlab工具包, 半波带滤波器设计工具」· M 代码 · 共 9 行

function y = sinc_decimate(x,m,r)
% y = sinc_decimate(x,m,r)	Decimate x by m-th order sinc of length r.
x = x(:)';
for i=1:m
	x = cumsum(x);
	x = [x(1:r) x(r+1:end)-x(1:end-r)]/r;
end
y = x(r:r:end);