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

1、课程设计任务书学生姓名： 专业班级： 电信0801 指导教师： 工作单位： 信息工程学院 题 目：dft对称性的验证及以应用初始条件： 具备数字信号处理的理论知识；具备matlab编程能力；了解dft的对称原理及应用；提供编程所需要的计算机一台要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1、独立编写程序验证dft的对称性2、用dft的对称性用一次fft实现两个序列的fft变换3、完成符合学校要求的设计说明书时间安排：一周，其中3天程序设计，2天程序调试指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日目 录摘要i1 dft基础知识11.1离散傅

2、立叶变换（dft）定义1 1.2复共轭序列的dft 11.3 dft的共轭对称性21.3.1有限长共轭对称序列和共轭反对称序列21.3.2共轭对称性分析32程序设计与分析62.1 n点dft对称性的验证6 2.1.1程序流程图62.1.2程序编写与结果分析72.2用一次fft实现两个序列的dft 132.2.1程序流程图.132.2.2程序编写与结果分析133 课程设计心得体会16参考文献17摘 要有限长序列在数字信号处理是很重要的一种序列，反映它的有限长特点的一种有用工具是离散傅里叶变换(dft)。离散傅里叶变换除了作为有限长序列的一种傅里叶表示法在理论上相当重要之外，而且由于存在着计算离散

3、傅里叶变换的有效快速算法，因而离散傅里叶变换在各种数字信号处理的算法中起着核心的作用。而离散傅立叶变换的对称性，在求实序列的离散傅立叶变换中有重要作用。可以实现一次dft的计算得到两个序列dft的高效算法，而dft可以通过一次快速fft变换来实现。关键词：dft 共轭对称性 matlab 1 dft基础知识1.1离散傅立叶变换（dft）定义有限长序列在数字信号处理是很重要的一种序列，当然可以用z变换和傅里叶变换来研究它，但是，可以导出反映它的有限长特点的一种有用工具是离散傅里叶变换(dft)。离散傅里叶变换除了作为有限长序列的一种傅里叶表示法在理论上相当重要之外，而且由于存在着计算离散傅里叶变

4、换的有效快速算法，因而离散傅里叶变换在各种数字信号处理的算法中起着核心的作用。设x(n)是一个长度为m的有限长序列，则定义x(n)的n点离散傅里叶变换为：正变换：=dft = = 反变换：=idft= 或rn(k)rn(k)x(n)= rn(n) =rn(n)式中，n称为dft变换区间长度，nm。dft隐含有周期性。1.2复共轭序列的dft设是的复共轭序列，长度为n，则（1）已知dft则 dft= 且（2）已知dft则 dft= 1.3 dft的共轭对称性dft有对称性，但由于dft中讨论的序列及其离散傅立叶变换均为有限长序列，且定义区间为0到n-1，所以这里的对称性是指关于n/2点的对称性。

5、下面讨论dft的共轭对称性质。1.3.1 有限长共轭对称序列和共轭反对称序列长度为的有限长序列,若满足 , （1.1） 称序列为共轭对称序列，一般用来表示。若满足 , （1.2）称序列为共轭反对称序列，一般用来表示即=, 0nn-1=, 0nn-1当n 为偶数时，把 代入式（1.1）与式（1.2），得 , （1.3）, （1.4） 式（1.3）与式（1.4）说明共轭对称序列与其共轭序列以成偶对称，共轭反对称序列与其共轭序列以成奇对称。当n为奇数时，把 代入式（1.1）与式（1.2），得 , （1.6） , （1.6） 式（1.5）与式（1.6）说明共轭对称序列与其共轭序列以成偶对称，共轭反对称

6、序列与其共轭序列以成奇对称。设一长度为的有限长序列，令则有 （1.7）这说明任一有限长序列，都表示成一个共轭对称序列与共轭反对称序列的和，在频域下同样有类似结论 （1.8）式中 （1.9） （1.10）1.3.2 共轭对称性分析（1）当x(n)为长度n的复数序列时,有 = = (1.11)同理可得 (1.12)即式（1.11）和（1.12）说明复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量；复书序列虚数部分的离散傅立叶变换是原来序列离散傅里叶变换的共轭反对称分量。另一方面，由式（1.7）知有限长序列可分解为共轭对称分量与共轭反对称分量，即=+ 可得其离散傅立叶变换 =

7、（1.13）同理可得 = （1.14）即上面两式说明复序列共轭对称分量序列的离散傅立叶变换是原来序列离散傅立叶变换的实数部分；复序列共轭对称分量的离散傅立叶变换是原来序列离散傅立叶变换的虚数部分。综上可得到有限长复序列的dft 的共轭对称性质如下将有限长序列x(n)分成实部与虚部，即： 则：将有限长序列x(n)分成共轭对称部分和共轭反对称部分，即=+，则：（2）当x(n)为长度n的实数序列或纯虚数序列时,有当x(n)为实序列时，则 又据)的对称性：有当x(n)为纯虚序列时，则 又据)的对称性：有 离散傅立叶变换的对称性，在求实序列的离散傅立叶变换中有重要作用。例如，有两个实数序列和，为求其离散

8、傅立叶变换，可以分别用和作为虚部和实部构造一个复数序列x(n)，求出x(n)的离散傅立叶变换，然后根据式（1.9）和（1.10）得到的共轭对称分量和，分别对应和，从而实现一次dft的计算可得到两个序列dft的高效算法。而dft可以通过一次快速fft变换来实现。2程序设计与分析本次课设计分两个部分，一个是要验证n点的dft的对称性，另一个是要用一次快速傅立叶变换fft实现两个序列的dft2.1 n点dft对称性的验证2.1.1程序流程图由于函数ezplot只能画出既存在symbolic math toolbox中又存在于总matlab工具箱中的函数，而gedc（实信号分解为循环偶分量和循环奇分量

9、）和dft(计算离散付利叶变换)仅存在symbolic math toolbox中，因此需要在自己的工作目录work下创建。此后可以直接调用这些函数。n点的dft的对称性验证流程图如图2-1所示开始求x序列的共轭对称与反对称分量画出共轭对称与反对称分量图形求出x(k)，xep，xop画出real(x(k) ),imag(x(k) )，xep，xop的图形xep结束图2-1 验证对称性流程图输入x序列n=0:n-12.1.2 程序编写与结果分析首先在目录work下创建gedc的m文件，gedc的m文件是用来生成共轭对称分量与共轭反对称分量的，程序如下：function xec,xoc=gedc(

10、x);n=length(x); n=0:(n-1); xec=0.5*(x + x(mod(-n,n)+1); xoc=0.5*(x - x(mod(-n,n)+1);再是在目录work下创建dft的m文件，dft为离散傅立叶变换，程序如下：function xk=dft(xn,n);n=0:1:n-1;k=0:1:n-1;wn=exp(-j*2*pi/n);nk=n*k;wnnk=wn.nk;xk=xn*wnnk;主程序：（1） n=12，序列为x=2.5 0 1.6 -3 -2 2 1.6 -3 -1 4 4.5 -2 的程序设计与结果分析程序：figure(1)n=0:11;x=inpu

11、t(请输入序列x=);xep,xop=gedc(x);subplot(2,1,1);stem(n,xep);title(共轭对称分量)xlabel(n);ylabel(xep);axis(-0.5,12.5,-3,4);subplot(2,1,2);stem(n,xop);title(共轭反对称分量);xlabel(n);ylabel(xop);axis(-0.5,12.5, -4,4);figure(2)x=dft(x,12) ;xep=dft(xep,12);xop=dft(xop,12);subplot(2,2,1);stem(n,real(x);axis(-0.5,12.5,-10,1

12、0);title( real(x);xlabel(k);subplot(2,2,2);stem(n,imag(x);axis(-0.5,12.5,-17,17);title( imag(x);xlabel(k);subplot(2,2,3);stem(n,xep);axis(-0.5,12.5, -10,10);title(dftxep(n);xlabel(k);subplot(2,2,4);stem(n,imag(xop);axis(-0.5,12.5,-17,17);title(dftxop(n);xlabel(k);结果：图2-2 共轭对称分量与共轭反对称分量图2-3 对称性的验证图形分

13、析：从图2-3可以看出复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量；复数序列虚数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭反对称分量。复序列共轭对称分量序列的离散傅立叶变换是原来序列离散傅立叶变换的实数部分；复序列共轭反对称分量的离散傅立叶变换是原来序列离散傅立叶变换的虚数部分。从而验证了dft的对称性。（2）n=14，序列为 x=(1.2).n 的程序设计与结果分析程序：figure(1)n=0:13;x=input(请输入序列x=);xep,xop=circevod(x);subplot(2,1,1);stem(n,xep);title(共轭对称分量)xl

14、abel(n);ylabel(xep);axis(-0.5,14.5,0,7);subplot(2,1,2);stem(n,xop);title(共轭反对称分量);xlabel(n);ylabel(xop);axis(-0.5,14.5,-6,6);figure(2)x=dft(x,14) ;xep=dft(xep,14);xop=dft(xop,14);subplot(2,2,1);stem(n,real(x);axis(-0.5,14.5,-13,6);title( real(x);xlabel(k);subplot(2,2,2);stem(n,imag(x);axis(-0.5,14.5

15、,-25,25);title(imag(x);xlabel(k);subplot(2,2,3);stem(n,xep);axis(-0.5,14.5,-13,6);title(dftxep(n);xlabel(k);subplot(2,2,4);stem(n,imag(xop);axis(-0.5,14.5,-25,25);title(dftxop(n);xlabel(k);结果：图2-4 共轭对称分量与共轭反对称分量图2-5 对称性的验证图形分析：从图2-5可以看出复数序列实数部分的离散傅立叶变换是原来序列离散傅立叶变换的共轭对称分量；复数序列虚数部分的离散傅立叶变换是原来序列离散傅立叶变换

16、的共轭反对称分量。复序列共轭对称分量序列的离散傅立叶变换是原来序列离散傅立叶变换的实数部分；复序列共轭反对称分量的离散傅立叶变换是原来序列离散傅立叶变换的虚数部分。从而验证了dft的对称性。2.2 用一次fft实现两个序列的dft2.2.1 程序流程图一次快速傅立叶变换fft实现两个序列的dft流程图如图2-5所示。开始输入x=+j调用fft函数得到和结束图2-6 一次fft变换实现两序列的dft2.2.2 程序编写与结果分析程序： x1=input(请输入序列x1=);x2=input(请输入序列x2=);n=input(请输入n=);x=x1+j*x2;x=fft(x,n);k=0:n-1

17、;c=conj(x);xep=0.5*(x+ c(mod(-k,n)+1);xop=-j*0.5*(x- c(mod(-k,n)+1);x1=xepx2=xopsubplot(2,1,1);stem(k,x1);xlabel(k);ylabel(x1);axis(-0.5,7.5,-10,40);subplot(2,1,2);stem(k,x2);xlabel(k);ylabel(x2);axis(-0.5,7.5,-10,40);结果：当运行程序时，会出现提示，按提示输入x1=1 3 5 2 4 6 3 5 2 6，x2=1 2 3 4 5 -5 -4 -3 -2 -1，n=10,程序运行结

18、果如下：x1和x2分别为x1,x2的离散傅立叶变换，x1和x2的图形如图2-7所示图2-7 x1和x2的离散傅立叶变换当直接调用dft时，程序运行结果和上面的是相同的，从而实现了用一次fft实现了两个序列的dft。 3 课程设计心得体会本次课程设计主要是运用本学期所学到的数字信号处理的基础知识来设计一个符合要求的matlab程序来进行dft对称性的验证以及应用，本次设计不仅要求我们要掌握数字信号处理课程的基础知识，还要求我们对matlab编程有深刻的理解和掌握。用新的语言去解决工程问题根本不需要先掌握某一门语言，有效的方法是先了解那门语言的一些基本函数，然后熟悉界面，就可以开始编了。拿到一个课题，不要急于坐在电脑前开始编程，因为当你坐在电脑前都不知道该干什么时，你就是对课题了解得不够。你首先需要的是透彻分析课题，把你要解决的问题写下来和列出各种可能情况。接下来，就考虑看用什么样的算法去解决，等到这一切都定下来后就可以开始着手编程了，如果你