example2_11.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 8 行

x=rand(1,11);n=0:10;
k=0:500; w=(pi/500)*k;
X=x*(exp(-j*pi/500)).^(n'*k);
[y,m]=sigshift(x,n,2);
Y=y*(exp(-j*pi/500)).^(m'*k);
error=max(abs(Y-exp(-j*2*w).*X));
disp(error);