《数字信号处理课程设计说明书-数字信号处理DFT对称性验证及应用》.docVIP

武汉理工大学《数字信号处理》课程设计说明书 PAGE 5 - 1.DFT的对称性原理分析 1.1共轭对称序列 长度为的有限长序列,若满足 , （1.1） 称序列为共轭对称序列，一般用来表示。 若满足 , （1.2） 称序列为共轭反对称序列，一般用来表示 把 代入式（1.1）与式（1.2），得 , （1.3） , （1.4） 式（1.3）与式（1.4）说明共轭对称序列与其共轭序列以成偶对称，共轭反对称序列与其共轭序列成奇对称 设一长度为的有限长序列，令 （1.5） （1.6） 则有 （1.7） 说明任一有限长序列，都表示成一个共轭对称序列与共轭反对称序列的和，称为的共轭对称分量，称为的共轭反对称分量。在频域下同样有类似结论 （1.8） 式中 （1.9） （1.10） 1.2有限长序列的对称分量分解及其DFT表示 （1）当x(n)为长度N的复数序列时,有 = ] = （1.11） 同理可得 (1.12) 式（1.11）和（1.12）说明复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量；复序列虚数部分的离散傅立叶变换的共轭反对称分量。 另一方面，由式（1.7）知有限长序列可分解为共轭对称分量与共轭反对称分量，可得其离散傅立叶变换 = （1.13） 同理可得 = （1.14） 上面两式说明复序列共轭对称分量序列的离散傅立叶变换是原来序列离散傅立叶变换的实数部分；复序列共轭对称分量的离散傅立叶变换是原来序列离散傅立叶变换的虚数部分。 离散傅立叶变换的对称性，在求实序列的离散傅立叶变换中有重要作用。例如，有两个实数序列和，为求其离散傅立叶变换，可以分别用和作为虚部和实部构造一个复数序列x(n)，求出x(n)的离散傅立叶变换，然后根据式（1.9）和（1.10）得到的共轭对称分量和，分别对应和，从而实现一次DFT的计算可得到两个序列DFT的高效算法。而DFT可以通过一次快速FFT变换来实现。 （2）当x(n)为长度N的实数序列或纯虚数序列时,有 当x(n)为实序列时，则 又据)的对称性： 当x(n)为纯虚序列时，则 又据)的对称性： 2.对称性分析及流程图 本次课设分两个部分，一个是要验证11点的DFT的对称性，另一个是要用一次快速傅立叶变换FFT实现两个序列的DFT 由于函数ezplot只能画出既存在Symbolic Math Toolbox中又存在于总MATLAB工具箱中的函数，而 circevod（实信号分解为循环偶分量和循环奇分量）和dft(计算离散付利叶变换)仅存在Symbolic Math Toolbox中，因此需要在自己的工作目录work下创建。此后可以直接调用这些函数。11点的DFT的对称性流程图2.1 和一次快速傅立叶变换FFT实现两个序列的DFT流程图2.2 开始开始 开始 开始 图2.2 一次FFT变换实现两序列的DFT结束得到和调用fft函数输入x=+j 图2.2 一次FFT变换实现两序列的DFT 结束 得到和 调用fft函数 输入x=+j n=0:10;x=input(请输入序列x=) 调用circevod，求对称分量 画出对称