数字信号处理 Lecture 6: Properties of Discrete Fourier Transformation 杨再跃

Slides:



Advertisements
Similar presentations
2 和 5 的倍数的特征 运动热身 怎样找一个数的倍数? 从小到大写出 2 的倍数( 10 个): 写出 5 的倍数( 6 个) 2 , 4 , 6 , 8 , 10 , 12 , 14 , 16 , 18 , 20 5 , 10 , 15 , 20 , 25 , 30.
Advertisements

第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
练一练: 在数轴上画出表示下列各数的点, 并指出这些点相互间的关系: -6 , 6 , -3 , 3 , -1.5, 1.5.
信号与系统 第三章 傅里叶变换 东北大学 2017/2/27.
§3.4 空间直线的方程.
《解析几何》 -Chapter 3 §7 空间两直线的相关位置.
第3章 离散傅立叶变换 DFS DFS的性质 DFT DFT的性质 圆周卷积 利用DFT计算线性卷积 频率域抽样.
1.2 信号的描述和分类.
數位訊號處理 第4章 離散時間訊號與LTI系統之傅利葉分析
第五章 二次型. 第五章 二次型 知识点1---二次型及其矩阵表示 二次型的基本概念 1. 线性变换与合同矩阵 2.
内容提要 傅立叶级数 傅立叶变换 典型信号的傅立叶变换 周期信号的傅立叶变换 抽样信号的傅立叶变换 抽样定理 第二章 傅立叶变换( FT )
第2章 时域离散信号和系统的频域分析 教学内容包括: 序列的傅立叶变换定义及性质 Z变换的定义与收敛域 利用z变换分析信号和系统的频域特性.
《高等数学》(理学) 常数项级数的概念 袁安锋
第一章 绪论.
第一章 行列式 第五节 Cramer定理 设含有n 个未知量的n个方程构成的线性方程组为 (Ⅰ) 由未知数的系数组成的n阶行列式
第二节 微积分基本公式 1、问题的提出 2、积分上限函数及其导数 3、牛顿—莱布尼茨公式 4、小结.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
2-7、函数的微分 教学要求 教学要点.
第7章 离散信号的频域分析 离散Fourier级数 离散Fourier变换 第3章 连续信号的频域分析 连续Fourier级数
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
实验二 FFT与DFT计算时间的比较及圆周卷积代替线性卷积的有效性实验
Signals and Systems Lecture 28
信号与系统基础 (二) 王烁
第2章 Z变换 Z变换的定义与收敛域 Z反变换 系统的稳定性和H(z) 系统函数.
第3章 离散傅里叶变换(DFT) 及其快速算法(FFT)
Discrete Fourier Transforms
元素替换法 ——行列式按行(列)展开(推论)
熟悉傅里叶变换的性质 熟悉常见信号的傅里叶变换 了解傅里叶变换的MATLAB实现方法. 熟悉傅里叶变换的性质 熟悉常见信号的傅里叶变换 了解傅里叶变换的MATLAB实现方法.
计算机数学基础 主讲老师: 邓辉文.
§2 求导法则 2.1 求导数的四则运算法则 下面分三部分加以证明, 并同时给出相应的推论和例题 .
Biomedical signal processing
Chapter 3 Discrete Fourier-Transform (Part Ⅰ)
数字信号处理 Lecture 4: Analysis of Discrete-time System 杨再跃
绿色圃中小学教育网 比例 比例的意义 绿色圃中小学教育网
实验一: 信号、 系统及系统响应 1、实验目的 1 熟悉连续信号经理想采样前后的频谱变化关系, 加深对时域采样定理的理解。
第一章 函数与极限.
数列.
2.1.
第四章习题.
线段的有关计算.
晶体管及其小信号放大 -单管共射电路的频率特性.
线 性 代 数 厦门大学线性代数教学组 2019年4月24日6时8分 / 45.
晶体管及其小信号放大 -单管共射电路的频率特性.
熟悉傅里叶变换的性质 熟悉常见信号的傅里叶变换 了解傅里叶变换的MATLAB实现方法. 熟悉傅里叶变换的性质 熟悉常见信号的傅里叶变换 了解傅里叶变换的MATLAB实现方法.
§2-1现实生活中的问题与函数的概念 例2.钟表问题
用计算器开方.
2019/5/2 实验一 离散傅立叶变换的性质及应用 实验报告上传到“作业提交”。 08:20:28.
实验一 熟悉MATLAB环境 常用离散时间信号的仿真.
2019/5/4 实验三 离散傅立叶变换的性质及应用 06:11:49.
正切函数的图象和性质 周期函数定义: 一般地,对于函数 (x),如果存在一个非零常数T,使得当x取定义域内的每一个值时,都有
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
函 数 连 续 的 概 念 淮南职业技术学院.
3.1 离散傅里叶变换的定义 3.2 离散傅里叶变换的基本性质 3.3 频率域采样 3.4 DFT的应用举例
2019/5/11 实验四 FIR滤波器的特性及应用 05:31:12.
2019/5/11 实验三 线性相位FIR滤波器的特性 05:31:30.
1.设A和B是集合,证明:A=B当且仅当A∩B=A∪B
2019/5/21 实验一 离散傅立叶变换的性质及应用 实验报告上传到“作业提交”。 11:21:44.
1.4.3正切函数的图象及性质.
1.4.3正切函数的图象及性质.
第二节 函数的极限 一、函数极限的定义 二、函数极限的性质 三、小结 思考题.
§2 方阵的特征值与特征向量.
第6章 离散信号与系统的频域分析 6.1 周期信号的离散时间傅里叶级数 6.2 非周期信号的离散时间傅里叶变换
24.4弧长和扇形面积 圆锥的侧面积和全面积.
定义5 把矩阵 A 的行换成同序数的列得到的矩阵,
第三节 数量积 向量积 混合积 一、向量的数量积 二、向量的向量积 三、向量的混合积 四、小结 思考题.
三角 三角 三角 函数 余弦函数的图象和性质.
§4.5 最大公因式的矩阵求法( Ⅱ ).
第三章 图形的平移与旋转.
Presentation transcript:

数字信号处理 Lecture 6: Properties of Discrete Fourier Transformation 杨再跃 Email: yangzy@zju.edu.cn 玉泉校区工控新楼507

复习:DFT 为了进行计算机处理,我们对连续信号通过采样离散化,采样的数学描述为: 为了分析x(n)的频率特性,我们对它进行DTFT,得到一个周期性的、连续的频谱曲线,上述过程可以用下图描述: t x(t) x(nTs) Ts - Ts 2Ts t --- 以时间间隔Ts采样 DTFT

为了用计算机在频域下分析信号,我们对x(n)经过DTFT后得到的连续曲线,在其一个周期内继续离散化,具体做法是:把DTFT的一个周期N等分,采样N个频率点上的值,用公式表示为: 上述过程用图表示为: --- DTFT --- DTFT 0 1 2 N-1 DFT 一个周期内N点离散化 采样 定义:我们把X(k)称为长为N的序列x(n)的离散傅里叶变换(DFT),恰好它也是N个点构成的序列,但是要注意,时域下的n值与频域下的k值不是对应的!

因此,我们可以定义频率分辨率: A1. 观察连续信号x(t),它本身有自己的频率f,也就是信号本身的周期为T=1/f A2. 以采样频率fs对x(t)采样后,得到了离散信号x(n),它是由采样后的离散点组成的,相邻两点之间的时间间隔就是采样周期Ts A3. 我们对x(n)进行DTFT,把时间t=nTs的函数x(n) 转换成数字角频率ω的函数X(ω),在这一转换过程中,存在着如下的频率对应关系: A4. 观察DTFT,注意到横轴是数字化角频率,并且具有周期性,周期为2π,在一个周期[0, 2π)内,它给出了信号从 频率=0 连续变化到 频率=fs 的信号频谱(参见第3条的频率对应公式),即DTFT的一个周期[0, 2π)对应了信号时域下的频率[0, fs] A5. 观察DFT,它是离散的,由于离散是通过对[0, 2π)进行N等分得到的,因此相邻离散点之间的间隔为2π/N. A6. 综合4、5两点可以看出,在DFT的给出的频谱图中,相邻两个点之间的间隔表示为数字角频率是2π/N,表示为时域下的实际频率为fs/N 因此,我们可以定义频率分辨率:

B1. 观察连续信号x(t),它本身有自己的频率f,也就是信号本身的周期为T=1/f B2. 以采样频率fs对x(t)采样后,得到了离散信号x(n),它是由采样后的离散点组成的,相邻两点之间的时间间隔就是采样周期Ts B3. 我们可能只截取了周期为T的连续时间信号x(t) 的一小段构成x(n),这一小段的长度为N,我们以N为周期进行序列周期延拓,得到如下关系式: B4. 注意x(t) 具有周期T,而 的周期为N,N值既是序列的点数,也是 的周期值,由于 是一个周期序列,因此我们可以把它展开成傅里叶级数DFS: 其中 B5. 观察DFS的表达式 ,它是描述了一个时域下、周期性的离散信号序列,序列的主值区间取值为第3步得到的x(n) B6.观察DFS的表达式 ,式中旋转因子项的系数 也是一个周期序列,具有如下表达式:

A1-A6的分析,从对DTFT进行N点抽样的角度,得到DFT定义; B1-B8的述分析,从离散傅里叶级数的角度得到DFT定义; B7.观察 的表达式,可以发现它与前面从DTFT采样得到的DFT有着一样的表达形式,唯一的区别是, 是一个无限长的周期序列,而X(k)是一个N点序列 B8. 如果画出 与X(k)的图,可以发现,后者恰好就是前者主值序列,因此DFT可以被认为是周期序列傅里叶级数的复指数项系数的主值部分,事实上,原始的DFT定义就是根据这一物理含义推导出来的,而不是前面所述对DTFT抽样获得 A1-A6的分析,从对DTFT进行N点抽样的角度,得到DFT定义; B1-B8的述分析,从离散傅里叶级数的角度得到DFT定义; 不论如何得到DFT,最重要的都是要知道为什么要费力得到DFT,即把时域离散信号转化为频域信号,通过计算机分析离散频率点上的信号幅值,试图从此了解原有时域信号的本质特征!

C1. 只要有一个N点序列,就可以用下面的定义找到一个N点DFT分析它的频率响应特性,如果序列的长为L<N,可以补零使长度变为N C2. 反过来,如果用N点DFT还原时域信号,则会得到一个无限长的周期信号,对该信号的主值区间取值就可以得到原有N点离散信号,即频域抽样造成时域周期延拓:

如何应用DFT 时域序列 频域DFT序列 周期与频率: 时域信号周期,信号频率; 时域信号的延拓周期,信号的点数; 采样周期,采样频率; Friday, January 18, 2019 如何应用DFT 时域序列 频域DFT序列 x(n)表示第n个分量值; X(k)表示第k个频谱分量的幅度 周期与频率: 时域信号周期,信号频率; 时域信号的延拓周期,信号的点数; 采样周期,采样频率; 周期与频率: 与模拟频率对应的数字频率; DTFT的周期; DFS的周期; DFT的点数;

Example: 宽为L的矩形序列及其DFT表达 Friday, January 18, 2019 Example: 宽为L的矩形序列及其DFT表达 5. 频率分辨率(抽样间隔): 6. 数据记录长度: 7. 分析频率: 1. 信号频率: 2. 数字角频率: 3. 采样(抽样)频率: 4. 信号最大频率:

如何应用DFT 实际测试一个信号时,为了避免混叠失真,应该尽可能提高采样频率使之必须满足抽样定理: Friday, January 18, 2019 如何应用DFT 实际测试一个信号时,为了避免混叠失真,应该尽可能提高采样频率使之必须满足抽样定理: 在抽样点数N给定的条件下,为了提高频率分辨率,必须降低采样频率,如果不能满足采样定理,则只能降低信号可分析的带宽 在采样频率和信号最高频率给定的情况下,为了提高分辨率,则必须提高采样点数N,从而也会增加数据记录的长度

计算机最小应该记录多长时间的数据(T0=?) 为了计算方便,要求抽样点数是2的次幂,在一次记录中该最少记录多少个点的数据(N=?) Friday, January 18, 2019 例6.1 现对一个信号进行频谱分析,估计该信号最高频率不大于2000Hz,要求频率分辨率不大于5Hz,由于计算机缓存容量限制,不能存储过长的序列。请确定如下指标: 计算机最小应该记录多长时间的数据(T0=?) 为了计算方便,要求抽样点数是2的次幂,在一次记录中该最少记录多少个点的数据(N=?) 抽样后,相邻两个数据点之间的最大时间间隔为多少(T=?) 1. Solution: 2. 3.

序列的对称性 则称xe(n)为偶对称序列 1. 序列xe(n)满足 若两个序列x(n)和y(n)满足 n 则称x(n)和y(n)互为偶对称 Friday, January 18, 2019 序列的对称性 1. 序列xe(n)满足 则称xe(n)为偶对称序列 n 偶对称序列(偶序列) 若两个序列x(n)和y(n)满足 则称x(n)和y(n)互为偶对称 n VS

序列的对称性 则称xo(n)为奇对称序列 2. 序列xo(n)满足 若两个序列x(n)和y(n)满足 n 则称x(n)和y(n)互为奇对称 Friday, January 18, 2019 序列的对称性 2. 序列xo(n)满足 则称xo(n)为奇对称序列 n 奇对称序列(奇序列) 若两个序列x(n)和y(n)满足 则称x(n)和y(n)互为奇对称 n

序列的对称性 任意定义在[0, N-1]上的长度为N的实序列可以被表示称为一个偶序列和一个奇序列的和: 则称xe(n)为偶序列 Friday, January 18, 2019 序列的对称性 1. 序列xe(n)满足 则称xe(n)为偶序列 2. 序列xo(n)满足 则称xo(n)为奇序列 任意定义在[0, N-1]上的长度为N的实序列可以被表示称为一个偶序列和一个奇序列的和:

共轭对称(Conjugate Symmetry) Friday, January 18, 2019 共轭对称(Conjugate Symmetry) 若序列为实序列,则为偶对称(偶序列) 1. 复序列取共轭、折叠后与原序列相等,则称序列为共轭对称序列 若序列为实序列,则为奇对称(奇序列) 2. 复序列取共轭、折叠后与原序列符号相反,则称序列为共轭反对称序列

共轭对称(Conjugate Symmetry) Friday, January 18, 2019 共轭对称(Conjugate Symmetry) 1. 共轭对称序列 2. 共轭反对称序列 任意定义在[0, N-1]上的长度为N的序列可以被表示称为一个共轭对称序列和一个共轭反对称序列的和:

周期共轭对称(Periodic Conjugate Symmetry) Friday, January 18, 2019 周期共轭对称(Periodic Conjugate Symmetry) 序列依次均匀排列在圆周上,序列关于标号为零的点循环对称 1. 定义在[0, N-1]上的长度为N的复序列满足下面条件时,被称为周期(循环)共轭对称序列 2. 定义在[0, N-1]上的长度为N的复序列满足下面条件时,被称为周期(循环)共轭反对称序列

周期共轭对称(Periodic Conjugate Symmetry) Friday, January 18, 2019 周期共轭对称(Periodic Conjugate Symmetry) 1. 周期共轭对称 2. 周期共轭反对称 任意定义在[0, N-1]上的长度为N的序列可以被表示称为一个周期共轭对称序列和一个周期共轭反对称序列的和: 证明?

周期共轭对称(Periodic Conjugate Symmetry) Friday, January 18, 2019 周期共轭对称(Periodic Conjugate Symmetry) 1. 周期共轭实序列xs(n)是关于N/2的对称序列 2. 周期共轭实序列xa(n)是关于N/2反对称序列 如果N为偶数,则

周期共轭对称(Periodic Conjugate Symmetry) Friday, January 18, 2019 周期共轭对称(Periodic Conjugate Symmetry) 1. 周期共轭实序列xs(n)是关于N/2的对称序列 2. 周期共轭实序列xa(n)是关于N/2反对称序列 任意定义在[0, N-1]上的长度为N的实序列可以被表示称为一个对称序列和反对称序列的和:

序列的循环折叠(圆周折叠) 一个有限长序列的循环折叠定义为 操作:将有限长序列以序列中间点(圆周上零标号点)进行反转即可得到; 几何意义:如果原序列顺时针在圆周均布,则循环折叠后将逆时针均布在圆周上;

序列的循环折叠(圆周折叠) x(n) = { 0, 3, 2, 4, 1, 5} N=6 一个有限长序列的循环折叠定义为 例6.2 求下列序列的循环折叠序列h(n) x(n) = { 0, 3, 2, 4, 1, 5} Solution: x'(n) = {0, 5, 1, 4, 2, 3} N=6  

序列的循环移位(圆周移位) 1. 序列的线性移位: 2. 序列的循环移位: 序列左移 序列右移 Friday, January 18, 2019 序列的循环移位(圆周移位) 1. 序列的线性移位: 序列右移   序列左移     2. 序列的循环移位: 1 2 3 4 5 6 7 把有限长序列的元素首尾相连,构成一个圆环,序列的元素在圆环上转动若干位置后构成新的序列 移位2

序列的循环移位(圆周移位) 一个有限长序列的圆周移位定义为 1. 将有限长序列进行周期延拓: 2. 进行移位: 3. 取主值序列

序列的循环移位 n N-1 n 周期延拓 n 左移2 n 取主值 N-1

Friday, January 18, 2019 例6.3 求序列向右循环移位1后的序列: Solution:

Friday, January 18, 2019