《数字信号处理： MATLABdft对称性验证以及应用》.docVIP

目 录 1 MATLAB基本操作及常用命令介绍 1 1．1 MATLAB的启动 1 1．2桌面平台 1 1.3 基本平面图形绘制命令 plot 3 2 理论分析 3 2.1实验内容 3 2.2序列对称性的理论验证 3 3 程序验证 5 4 结果分析 8 5 对称性的应用 10 5.1 FFT算法的基本思想 10 5.2 对称性应用的程序实现 12 6 心得体会 15 参考文献 161 MATLAB基本操作及常用命令介绍 1．1 MATLAB的启动 启动MATLAB有多种方式，最常用的方法就是双击系统桌面的MATLAB图标， 也可以在开始菜单的程序选项中选择MATLAB快捷方式。 初次启动MATLAB后，将进入MATLAB默认设置的桌面平台。 1．2桌面平台 默认设置情况下的桌面平台包括6个窗口，分别是MATLAB主窗口、命令窗口、历史窗口、当前目录窗口、发行说明书窗口和工作间管理窗口。下面分别对各窗口做简单介绍。 （1）MATLAB主窗口，与MATLAB的早期版本不同的是，MATLAB 6.1增加了一个主窗口，其他的几个窗口都包含在这个大的主窗口中。 （2）命令窗口， 在主窗口的“View”菜单下选择“Command Window”命令，可以打开或关闭MATLAB的命令窗口。其中“”为运算提示符，表示MATLAB正处在准备状态。当在提示符后输入一段运算式并按Enter键后，MATLAB将给出计算结果，然后再次进入准备状态。 （3）工作间窗口 在主窗口的“View”菜单下选择“Workspace”命令，可以打开或关闭MATLAB的工作间窗口，工作间窗口是MATLAB的重要组成部分。 MATLAB的程序类型有三种 1）脚本M文件： 在命令窗口中输入并执行，它所用的变量都要在工作空间中获取，不需要输入输出参数的调用，退出MATLAB后就释放了。 2）程序M文件： 以.m格式进行存取，包含一连串的MATLAB指令和必要的注解。需要在工作空间中创建并获取变量，也就是说处理的数据为命令窗口中的数据，没有输入参数，也不会返回参数。程序运行时只需在命令窗口中键入其名称即可（不需要后缀名）。 函数M文件（常称为M函数）：与在命令窗口中输入命令一样，函数接受输入参数，然后执行并输出结果。用help命令可以显示它的注释说明。 4）M函数具有标准的基本结构： 1）函数定义行（关键字function）function[out1,out2,..]=filename(in1,in2,..) 输入和输出（返回）的参数个数分别由nargin和nargout两个MATLAB保的变量来给出。2）第一行帮助行以（%）开头，作为lookfor指令搜索的行3）函数体说明及有关注解以（%）开头，用以说明函数的作用及有关内容4）函数体语句：函数体内使用的除返回和输入变量这些在function语句中直接引用的变量以外的所有变量都是局部变量，即在该函数返回之后，这些变量会自动在MATLAB的工作空间中清除掉。如果希望这些中间变量成为在整个程序中都起作用的变量，则可以将它们设置为全局变量。 1.3 基本平面图形绘制命令 plot 1）plot(X,Y) %当X,Y均为实数向量，且为同维向量，则以x的元素为横坐标，y的元素为纵坐标，绘制曲线。若X，Y均为同维同型实数矩阵，X = [X(i)]，Y = [Y(i)]，其中X(i),Y(i)为列向量，则plot(X,Y)依次画出plot(X(i),Y(i))，矩阵有几列就有几条线； 2）plot(X1,Y1,X2,Y2,…)，其中Xi与Yi成对出现，plot(X1,Y1,X2,Y2,…)将分别按顺序取两数据Xi与Yi进行画图。 3）plot(X1,Y1,LineSpec1,X2,Y2,LineSpec2…) 将按顺序分别画出由三参数定义Xi,Yi,LineSpeci的线条。其中参数LineSpeci指明了线条的类型，标记符 2 理论分析 2.1实验内容 1）编写程序验证DFT的正确性 2）利用DFT的对称性用一次FFT实现两个序列的FFT变换 2.2序列对称性的理论验证 (1) 有限长共轭对称序列和共轭反对称序列定义： ?共轭对称序列：x e p ( n ) = x *e p ( N－n ) , 0 ≤n≤ N －1 ?共轭反对称序列：x o p (n) =－x*o p ( N－n ) ,0≤n≤N－1 ?如下图所示：图2-1 共轭对称与共轭反对称序列示意图 ? (2) DFT的共轭对称性1)设x(n)=xr(n)+jxi(n)其中xr(n)=Re[x(n)]=[x(n)+x*(n)]/2 ,xi(n)=Im[x(n)]=[x(n)－x*(n)]/2则?? DFT[xr(n)]=Xep(k)DF