matlab数字信号处理实验报告,数字信号处理matlab实验报告4.doc

数字信号处理matlab实验报告4.doc

实验报告课程名称数字信号处理实验专业班级姓名学号实验名称数字信号处理实验地点实验时间实验成绩实验目的及任务理解信号变换的基本概念理解离散傅立叶变换的基本概念掌握快速傅立叶变换的应用方法掌握离散余弦变换的应用方法掌握Z变换的应用方法了解CHIPZ变换的基本概念掌握HILBEIT变换的初步应用了解倒谱变换的基本概念实验内容与步骤信号变换概述信号是数字信号处理领域中最基本、最重要的概念。而数字信号变换技术，又是对信号进行处理操作的最基本的有效途径之一。因此，数字信号变换技术，便成为数字信号处理领域中专业人员所必须要张我的一项最基本的技能。简单地说，数字信号变换技术就是为了处理操作上的方便和可能，通过数学变换，将一个域内的信号变换映射倒另一个域内的信号的方法。常用的数字信号变换主要有傅立叶变换、离散余弦变换(DCT)、Z变换、CHIRPZ变换、HILBERT变换等。这些变换，都有着各自的理论和其应用背景。MATLAB中的工具箱对这几种典型的变换，都提供了相对应的、具体的应用函数。这可以使得工程人员大大节省无谓的工作量，从而将主要精力放到新技术的创新和研发上面。下面将对这几种变换的含义和应用进行具体的介绍。离散傅立叶变换傅立叶变换是信号分析和处理的重要工具。有限长序列作为离散信号的一种，在数字信号处理种占有着极其重要的位置。对于有限长序列，离散傅立叶变换不仅在理论上有着重要的意义，而且有快速计算的方法－快速傅立叶变换。所以在各种数字信号处理的运算方法中，越来越起到核心的作用。下面，就对离散傅立叶变换及其MATLAB函数应用，结合实际工程实例做说明531傅立叶变换的几种形式1、非周期连续时间信号的傅立叶变换非周期连续时间信号的傅立叶变换可以表示为TXJX＝JXDTETJ逆变换为DJXTXT21在这里，是模拟角频率。可以看到，时域的连续函数造成频域的非周期谱，时域的非周期性造成频域的连续谱。结论非周期连续时间函数对应于一非周期连续频域变换函数。2、周期连续时间信号的傅立叶变换周期为的周期性连续时间信号傅立叶变换是离散频域函数，可表示为TTX21TTJMDEXJX逆变换为DEJTXTJMM这就是经常称之为傅立叶级数的变换形式。在这里，也是模拟角频率。可以看到，时域的连续函数造成频率域的非周期谱，频域函数的离散造成时域函数的周期性。结论周期连续时间函数对应于一非周期离散频域变换函数。3、非周期离散时间信号的傅立叶变换可以表示为NXJEXNJJX逆变换为DEXNJ21在这里，是数字频率，它和模拟角频率的关系为。可以看到，时域的取样T对应于频域的周期延拓，而时域函数的非周期性造成频域的离散谱。结论非周期离散时间函数对应于一周期连续频域变换函数。4、周期离散时间信号的傅立叶变换周期离散时间信号的傅立叶变换－离散傅立叶变换，可以表示为NX102NNKJEKX逆变换为102NKNKNJEX可以看到，时域的取样对应于频域的周期延拓，而时域函数的周期性造成频域的离散谱。结论周期离散时间函数对应于一周期离散频域变换函数。532离散傅立叶变换离散傅立叶级数变换是周期序列，仍不便于计算机计算。但离散傅立叶级数虽是周期序列，却只有个独立的数值，所以它的许多特性可以通过有限长序列延拓来得到。对于N一个长度为的有限长序列，也即只在个点上有非零值，其余NX10NN皆为零，即其他,01NX把序列以为周期进行周期延拓得到周期序列，则有NNX其他,01NXN所以，有限长序列的离散傅立叶变换(DFT)为10,10NNWNXDFTKXNK逆变换为10,10NKXNKINXNNN若将DFT变换的定义写成矩阵形式，则得到XA﹒X，其中DFT变换矩阵A为21111NNWADFTMTX函数用来计算DFT变换矩阵A的函数调用方式A＝DFTMTA(N)返回NN的DFT变换矩阵A。若X为给定长度的行向量，则Y＝XA，返回X的DFT变换Y。AI＝CONJ(DFTMTX(N))/N；返回NN的IDFT变换矩阵AI。DFT的性质两个序列和都是N点有限长序列，设1X2,1DFTKX22X线性式中A，B为任意常数。,2121KBXABNAFT圆周移位一个有限长序列的圆周移位定义XNRMXN式中，表示的周期延拓序列的移位XNXNNXN有限长序列圆周移位后的DFT为}{KXWNRMXDFTKXNNNM3、圆周卷积假设21KXY则有1012102NRMNXNRMNXKIDFTNNNMNMN用表示圆周卷积，则上式可化简为122121NXNXKXIDFTNY4共轭对称性令的共轭复数序列为，则XNKXNXFT用和分别表示序列的实部和虚部，即NXRIXJIR21NXNXIR用和分别表示实部和虚部序列的，即KXRIDFTNXDFTRRKII而且可以证明得到NRKXIIK通常称为的共轭偶部，为的共轭奇部。所以说，对于时KXRKI域、频域的对应关系来说，序列实部对应于的共轭偶部，序列的虚DFTNXXNX部对应于的共轭奇部。K5序列乘积121KXNNXDFT6DFT形式下的帕塞瓦尔定理1010NKNNYXYX由于MATLAB软件本身的特点，序列或向量元素下标从1开始记录，而不是从0开始。因此，上述两式在MATLAB中相应的表达式为NKNNWXKX102,0NK1KNKNX1,N而下面所讨论使用的快速傅立叶变换并不是与不同的另外一种变换，FTDFT而是为减少计算次数的一种快速有效的算法。这种快速算法，主要是利用了下DFTNKNW面两个特性使长序列的分解为更小点数的所实现的。54快速傅立叶变换(FFT)在信号处理中，DFT的计算具有举足轻重的地位，，信号的相关、滤波、谱估计等都要通过DFT来实现。然而，当很大的时候，求一个点的DFT要完成次复数NNN乘法和次复数加法，其计算量相当大。1965年JWCOOLEY和JWTUKEY巧妙地1N利用因子的周期性和对称性，构造了一个DFT快速算法，即快速傅立叶变换FFT。W通过前面的知识，已经知道有限列长为的序列的变换为NNXDFTNKNNWXKX1012,0K其逆变换为10NKNKX1,0NN上机练习1．试用MABLAB求其有限长序列与10801NNX的圆周卷积，(N20)，并画出其结果图。8062NXNCIRCONVT函数FUNCTIONYCIRCONVTX1,X2,NIFLENGTHX1NERROR NSHOULDBIGGERTHANOREQUALTOTHELENGTHOFX1 ENDIFLENGTHX2NERROR NSHOULDBIGGERTHANOREQUALTOTHELENGTHOFX2 X1X1,ZEROS1,NLENGTHX1X2X2,ZEROS1,NLENGTHX2M01N1FORN11NENDYX1H COMMANDWINDOWN10110N20118N20N01N1X108N1X206N2YCIRCONVX1,X2,NSTEMN,Y,G