1、数字信号处理,Lecture 4: Analysis of Discrete-time System杨再跃Email: 玉泉校区工控新楼 弧度转到了弧度,即覆盖了整个z平面,描述系统特性的工具,Sunday, September 23, 2018,5,离散时间傅里叶变换,离散时间序列x(n)的傅里叶变换,也称为离散时间傅里叶变换(Discrete-time Fourier Transformation, DTFT),定义为:,Sunday, September 23, 2018,6,上述变换的逆变换,(Inverse Discrete-time Fourier Transformation,
2、IDTFT),定义为:,研究DTFT的目的是:计算离散时间信号的频谱,计算离散时间系统的频率响应;序列x(n)收敛(绝对可和)是该序列DTFT存在的充分条件,离散时间傅里叶变换,DTFT是是数字频率的连续函数;DTFT的定义表达式也可以被认为是频谱的傅里叶级数;DTFT是一个以2为周期的周期函数;IDTFT可以被认为是上述傅里叶级数的系数;,Sunday, September 23, 2018,7,8,离散时间傅里叶变换,离散信号DTFT,连续信号FT,连续函数;是数字频率的周期函数;,连续函数;不具有周期性;,考虑序列的z变换,函数变量在复平面内取值,离散时间傅里叶变换,当r=1时,z变换的
3、定义与DTFT一致,即序列在单位圆上的z变换就是它的离散时间傅里叶变换;序列的z变换收敛域包含单位圆时,它的DTFT存在;,离散时间信号的频谱,Sunday, September 23, 2018,10,DTFT是分析离散时间序列的重要工具,通常X(e-j)是一个以为自变量的复函数,具有幅值和相位特性;函数的模(幅值) |X(e-j)|称为信号的幅度谱,辐角()=argX(e-j)称为相位谱,由于X(e-j)是周期为2的周期函数,因此一个周期就足以代表序列的频谱,通常取-, )或者0, 2);为了扩大幅度谱的动态显示范围,幅度谱常采用对数形式表示:20log |X(e-j)|,单位是dB,相位
4、的单位是度或弧度,通常取主值相位表示,即-() 频域共轭+折叠,时域折叠 - 频域折叠,常用DTFT性质,时域卷积定理,两边取傅里叶变换,LTI系统的零状态响应,18,已知复指数序列是LTI系统的特征函数,因此,LTI系统冲击响应的DTFT可由描述系统的差分方程直接得到,用差分方程求解DTFT,例4.4 已知LTI系统由下面的差分方程定义,求:(1) ;(2) 求输入为下面的余弦信号时稳态输出;,19,Solution:,离散系统的频率响应,Sunday, September 23, 2018,20,频率响应是分析系统动态特性的有效工具,可以由如下两种方法来定义:,1. 单位冲击响应的离散时间
5、傅里叶变换;2. 离散系统对全频段正弦输入(或复指数序列)的稳态响应;,频率响应通过幅频响应(曲线图)和相频响应(曲线图)来直观的分析系统特性。,上述二者的操作方法不同,但实际上是等价的!,频率响应(Frequency Response),定义1:离散时间LTI系统的单位冲击响应的傅里叶变换也称为系统的频率响应,称为系统幅频特性(幅度特性、幅度响应),,称为系统的相频特性,22,例 4.5 求下面序列的DTFT:,系统的激励是单位冲击序列时,它的频谱覆盖了于是系统的单位冲击响应可以看成对各次谐波的滤波的总的效果,因此可以用冲击响应的频谱来反映系统特征;,例4.6 已知系统的单位冲击响应为: ,
6、画出幅值与相位的频率响应;如果对上述系统施加输入信号为0.1u(n),求系统的稳态输出,23,Solution:,24,横坐标归一化处理,单位为,频率响应(Frequency Response),定义2:离散时间LTI系统的对全频段正弦信号响应的傅里叶变换也称为系统的频率响应,任何复杂信号都可被看做是不同幅度和相位的正弦波线性叠加而成,这些正弦波包含了丰富的频率信息,因此上述定义也反映系统对各次谐波的传输能力;,频率响应(Frequency Response),定义2:离散时间LTI系统的对全频段正弦信号响应的傅里叶变换也称为系统的频率响应,把正弦信号输入系统,然后逐渐增加输入信号的频率,观察
7、系统的输出在幅值和相位上的变化,特别是与输入信号相比有何改变,由此可以判断系统对输入的传输性能。,Tip:这就是为什么系统的传递函数可以被定义为频域内输出和输入的比值。,频率响应(Frequency Response),定义2:离散时间LTI系统的对全频段正弦信号响应的傅里叶变换也称为系统的频率响应,Sunday, September 23, 2018,28,Solution: (1),Sunday, September 23, 2018,29,Solution: (2) 以 为例说明:,Sunday, September 23, 2018,30,Solution: (2),信号的采样,模拟信
8、号通过采样(抽样)离散化,采样器相当于一个电子开关,它每隔T秒短暂地闭合一次,将模拟信号接通,读入一个采样值,开关的闭合时间为;,Sunday, September 23, 2018,31,信号的采样,理想采样条件下,无穷小,通常开关的闭合时间T,可以近似认为满足理想抽样;采样不可避免的会对原始信号包含的信息造成丢失,同时,由于采样器本身的特点,会在信号上增加一些不必要的信息,需要分析清楚这些影响;数学方法上看,理想采样相当于对信号乘以单位抽样序列;,Sunday, September 23, 2018,32,信号的采样,Sunday, September 23, 2018,33,信号的采样,
9、抽样除了把信号离散化外,还具有其他特性,可以通过研究抽样后信号的频谱进行分析,Sunday, September 23, 2018,34,信号的采样,Sunday, September 23, 2018,35,根据频域卷积定理,两信号在时域相乘的傅立叶变换等于两信号分别的傅立叶变换的卷积,Sunday, September 23, 2018,36,对比模拟信号的频谱 与抽样信号的频谱 可以发现,抽样后在频域内造成了周期延拓,周期为,Sunday, September 23, 2018,37,做如下替换:,Sunday, September 23, 2018,38,如果信号本身所包含的最大频率0
10、不大于采样频率的一半,则在一个延拓周期内,采样后的信号仍然可以保留原来模拟信号的频谱,二者仅在幅值上相差一个常数倍,采样定理,如果 ,在区间 内有:,Sunday, September 23, 2018,39,如果使得采样频率变小,使得信号本身所包含的最大频率0大于采样频率的一半,则在s/2处发生频率分量的重叠,造成信息不能分辨,这一现象称为频域混叠,采样定理,把s/2称为折叠频率( 奈奎斯特频率 ),当信号频谱超过这一频率时会对称折叠回相应的部分,使得信息混合;,Sunday, September 23, 2018,40,如果使得采样频率变小,使得信号本身所包含的最大频率0大于采样频率的一半
11、,则在s/2处发生频率分量的重叠,造成信息不能分辨,这一现象称为频域混叠,把s/2称为折叠频率,当信号频谱超过这一频率时会对称折叠回相应的部分,使得信息混合;实际抽样前,为了防止混叠发生,可加入截止频率为s/2的前置预滤波器;,41,信号的还原,如果采样频率满足奈奎斯特定理,可以把抽样后的信号还原为模拟信号,Sunday, September 23, 2018,42,对上述信号进行傅里叶逆变换:,Sunday, September 23, 2018,43,内插公式,44,信号的还原,在物理实现上,抽样后的信号经过一个理想的低通滤波器,去除周期延拓的部分,就可以还原信号;,Sunday, September 23, 2018,45,实际低通滤波存在过渡带,可以提高采样频率来避免过渡带的影响;,例4.8 已知 ,求下面函数的采样信号及其DTFT,46,Solution:,Sunday, September 23, 2018,47,Assignment 1:构造一个连续时间信号,该信号包含3-5个小于100Hz的整数频率分量,写出信号的函数表达式,画出函数图;设定合适的采样频率,对信号进行采样,得到一个离散时间序列,截取合适长度的序列,画出离散点图;用DTFT画出信号的频谱图; 交作业时间: 3月25日晚发送word文件至,