新MATLAB辅助现代工程数字信号处理 第二版 教学课件 李益华 第1 5章 第2章.ppt

第2章 离散时间信号与系统的频域分析 2.1 离散时间序列的傅里叶变换(DTFT)　　如果信号在频域上是离散的，则信号在时域上就表现为周期性的时间函数。相反，如果信号在时域上是离散的，则其在频域上就表现为周期性的频率函数。因此，一个离散周期序列的频谱一定既是周期的又是离散的。本章先介绍离散时间傅里叶变换(DiscreteTime Fourier Transform，DTFT)，再介绍离散周期序列及其傅里叶级数(DFS)。 2.1.1 DTFT的定义　　序列x(n)的离散时间傅里叶变换(DTFT)定义为 　　式(2.1)和式(2.3)为序列x(n)的傅里叶变换对，即离散时间信号的傅里叶变换对。式(2.1)为正变换，式(2.3)为反变换或称为逆变换。　　式(2.1）又可视为X(ejω)的傅里叶级数展开式，式(2.3)确定的x(n)是傅里叶级数的系数。 2.1.2 DTFT的性质　　1. 线性性　　定义信号x1(n)和x2(n)的DTFT分别是X1(ejω)和X2(ejω)，并令x(n)=ax1(n)+bx2(n)，其中，a和b为常数，则 　　　　X(ejω)=aX1(ejω)+bX2(ejω)　　　　　(2.4)　　2. 时移与频移性　　令y(n)=x(n－n0)，则　　　　　　　　　　Y(ejω)=e－jωn0X(ejω)　　　　　 (2.5) 　　3. 对称性　　将序列x(n)分成实部与虚部，即　　　　　　　　　　x(n)=xR(n)+jxI(n)　　　　　　(2.6)将式(2.6)进行DTFT，得到　　　　　　　　X(ejω)=XR(ejω)+jXI(ejω)　　　　　(2.7)　　X(ejω)的实部XR(ejω)是偶函数，虚部XI(ejω)是奇函数，即　　　　　　　　　　XR(ejω)=XR(e－jω)　　　　　　(2.8)　　　　　　　　　　XI(ejω)=－XI(e－jω)　　　　　 (2.9) 　　4. 时域卷积定理　　设y(n)=x(n)*h(n)，则　　　　　　　　　　Y(ejω)=X(ejω)·H(ejω)　　　(2.10)该定理说明，两序列卷积的DTFT服从相乘的关系。　　线性时不变系统输出y(n)的DTFT等于输入信号x(n)的DTFT乘以系统单位脉冲响应h(n)的DTFT。因此，求系统的输出信号可以在时域用卷积公式计算，也可以在频域按照式(2.10)，求出输出的DTFT，再作逆DTFT求出输出信号。　　5. 频域卷积定理　　若y(n)=x(n)h(n)，则 　　　　　　　　　　　　　　　　　　　　　　　(2.11) 　 　　6. 时域相关定理　　若y(n)是x(n)和h(n)的相关函数，即 　　8. 维纳-辛钦(Wiener-Khinchin)定理　　若x(n)是功率信号，则其自相关函数的傅里叶变换为 　　无论x(n)是实信号还是复信号，其功率谱Px(ejω)始终是ω的实函数，即功率谱失去了相位信息。相关函数和功率谱是描述随机信号的重要统计量。 　　2.2 离散周期信号的傅里叶级数(DFS)2.2.1 DFS的定义　　一个周期为N的周期序列，即　　　　　　　，k为任意整数，N为周期。该序列可用离散傅里叶级数来表示，即可用周期为N的正弦序列来表示 　 故可得到 记　　　　 ，则DFS变换对可写为 　 (2.25) 　　【例 2.1】 设　　　　　　　　，以N=4为周期进行周期延拓，求周期序列的离散傅里叶级数。　　MATLAB程序如下: 　　%MATLAB PROGRAM 2-1　　%实现离散傅里叶级数(DFS)的计算　　%xn代表离散时间序列x(n)，N为离散时间序列x(n)的长度　　%Xk为x(n)的傅里叶级数，且为其主值序列　　xn= [0，1，2，3]; N= 4; %设定序列和周期　　n= [0: 1: N-1]; k= [0: 1: N-1]; 　　%设定n和k　　Wn= exp(-j*2*pi/N); 　　　　　%设定Wn因子　　nk= n′*k; 　　Wnnk = Wn.^nk; 　　　　%计算W矩阵　　Xk= xn*Wnnk; 　　　　%计算DFS的系数Xk　　disp(′xn = ′); disp(xn); 　　disp(′Xk = ′); disp(Xk); %显示计算结果(系数)　　程序运行结果为 　　xn =　　　 0 　1