ex3p8.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 9 行

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