Digital Signal Processing Using Matlab读书笔记 第三章（一）（P40-P42）

离散时间傅里叶分析

我们已经看到一个线性时不变系统是如何通过它的单位冲激响应表示的，单位冲激响应使得我们可以通过卷积计算任意的输入x(n)对应的响应

该卷积的表示形式是基于能通过按比例以及时延的单位信号的线性结合，每一个基础信号提供了一个新的信号形式，每一形式在不同的系统下既有优点，也有缺点，但是，当系统是线性以及时不变的，只有一个表示形式是最有用的，这就是离散时间傅里叶变换。

离散时间傅里叶变换（DTFT）

如果x(n)是绝对可和的，即

那么它的DTFT如下

X的反向离散时间傅里叶变换(IDTFT)如下

操作符F[]将一个离散信号x(n)转换成了以w为实变量的复连续函数X，w称为数字频率，以弧度为单位。
例3.1


例3.2


因为X是一个复函数，我们要单独表示它的幅频与相频图像，因为w是负无穷到正无穷的实数，在matlab上我们只能表示一部分，使用两个DTFT的重要属性，我们可以将w的取值范围缩小到[0:pi]，在下一部分我们讨论X的有用属性

两个重要性质

周期性

对称性

matlab应用

如果x是无限序列，matlab不能直接通过x计算X，但是我们能通过它来估计X在[0,pi]的频率并且画出幅频以及相频图像
例3.3

>> w=[0:1:500]*pi/500;
>> X=exp(j*w)./(exp(j*w)-0.5*ones(1,501));
>> magX=abs(X);angX=angle(X);
>> realX=real(X);imagX=imag(X);
>> subplot(2,2,1);plot(w/pi,magX);grid
>> subplot(2,2,3);plot(w/pi,angX);grid
>> subplot(2,2,2);plot(w/pi,realX);grid
>> subplot(2,2,4);plot(w/pi,imagX);grid