Digital Signal Processing Using Matlab读书笔记 第三章（三）（P47-P51）

DTFT的属性

在前面的章节中，我们为了展示图像讨论了两个重要的属性，现在我们讨论剩下的有用的属性。

应用

如果序列x(n)是实数且是偶对称，那么X也是实数偶对称。
复数序列的一个类似的性质将会在问题3.7中讨论。



在下面几个例子中，我们将会使用有限区间序列证明其中一些性质
例3.7

x1=rand(1,11);x2=rand(1,11);n=0:10;
alpha=2;beta=3;
k=0:500;w=(pi/500)*k;
X1=x1*(exp(-j*pi/500)).^(n'*k);
X2=x2*(exp(-j*pi/500)).^(n'*k);
x=alpha*x1+beta*x2;
X=x*(exp(-j*pi/500)).^(n'*k);
X_check=alpha*X1+beta*X2;;
error=max(abs(X-X_check));

因为两个傅里叶变换数组的差的模的最大值小于10的-14次方，两个数组在matlab的数字精度的限制下是完全相同的
例3.8

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));

例3.9

n=0:100;
x=cos(pi*n/2);
k=-100:100;
w=(pi/100)*k;
X=x*(exp(-j*pi/100)).^(n'*k);
y=exp(j*pi*n/4).*x;
Y=y*(exp(-j*pi/100)).^(n'*k);
subplot(1,1,1);
subplot(2,2,1);plot(w/pi,abs(X));grid;axis([-1,1,0,60])
subplot(2,2,2);plot(w/pi,angle(X)/pi);grid;axis([-1,1,-1,1])
subplot(2,2,3);plot(w/pi,abs(Y));grid;axis([-1,1,0,60])
subplot(2,2,4);plot(w/pi,angle(Y)/pi);grid;axis([-1,1,-1,1])

由图可知平移了pi/4。