dft在信号频谱分析中的应用课设

dft在信号频谱分析中的应用课设

ID:13793454

大小:362.50 KB

页数:7页

时间:2018-07-24

dft在信号频谱分析中的应用课设_第1页
dft在信号频谱分析中的应用课设_第2页
dft在信号频谱分析中的应用课设_第3页
dft在信号频谱分析中的应用课设_第4页
dft在信号频谱分析中的应用课设_第5页
资源描述:

《dft在信号频谱分析中的应用课设》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、设计一DFT在信号频谱分析中的应用一、设计题目DFT在信号频谱分析中的应用二、设计目的1.熟悉DFT的性质。2.加深理解信号频谱的概念及性质。3.了解高密度谱与高分辨率频谱的区别。三、设计原理信号的频谱分析就是计算信号的傅里叶变换。连续信号与系统的傅里叶分析显然不便于直接用计算机进行计算,使其应用受到限制,而DFT是一种时域和频域均离散化的变换,适合数值运算,成为分析离散信号和系统的有力工具。工程实际中,经常遇到的连续信号Xa(t),其频谱函数Xa(jW)也是连续函数。数字计算机难于处理,因而我们采用DF

2、T来对连续时间信号的傅里叶变换进行逼近,进而分析连续时间信号的频谱。四、实现方法离散傅里叶变换是有限长序列的傅里叶变换,它相当于把信号的傅里叶变换进行等频率间隔采样,并且有限长序列的离散傅里叶变换和周期序列的离散傅里叶级数本质是一样的。快速傅里叶变换(FFT)并不是一种新的变换,它是离散傅里叶变换的一种快速算法,并且主要是基于这样的思路而发展起来的:(1)把长度为N的序列的DFT逐次分解成长度较短的序列的DFT来计算。(2)利用WN(nk)的周期性和对称性,在DFT运算中适当的分类,以提高运算速度。(对称

3、性,;周期性,r为任意整数)离散傅里叶变换的推导:离散傅里叶级数定义为(1-1)将上式两端乘以并对n在0~N-1求和可得因为所以这样用k代替m得(1-2)令则(1-2)成为DFS(1-3)(1-1)成为IDFS(1-4)式(1-3)、(1-4)式构成周期序列傅里叶级数变换关系。其中都是周期为N的周期序列,DFS[·]表示离散傅里叶级数正变换,IDFS[·]表示离散傅里叶级数反变换。习惯上,对于长为N的周期序列,把0nN-1区间称为主值区,把称为的主值序列,同样也称为的主值序列。由于,对于周期序列仅有N个独

4、立样值,对于任何一个周期进行研究就可以得到它的全部信息。在主值区研究与是等价的,因此在主值区计算DFS和DFT是相等的,所以DFT计算公式形式与DFS基本相同。其关系为所以离散傅里叶正变换0kN-1离散傅里叶变换(DFT)定义:设有限长序列x(n)长为N(0nN-1),其离散傅里叶变换是一个长为N的频率有限长序列(0kN-1),其正变换为0kN-1()离散傅里叶变换的实质是:把有限长序列当做周期序列的主值序列进行DFS变换,x(n)、X(k)的长度均为N,都是N个独立值,因此二者具有的信息量是相等的。已知

5、x(n)可以唯一确定X(k),已知X(k)可以唯一确定x(n)。虽然离散傅里叶变换是两个有限长序列之间的变化,但它们是利用DFS关系推导出来的,因而隐含着周期性。构造离散傅里叶变换的Matlab实现程序如下:function[Xk]=dft(xn,N)n=[0:1:N-1];k=n;WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^nk;Xk=xn*WNnk快速傅里叶变换(FFT)并不是与DFT不同的另外一种变换,而是为了减少DFT计算次数的一种快速有效的算法共轭对称性:设有限长序列的

6、长度为N,以N为周期的周期延拓列为周期序列的共轭对称分量和共轭反对称分量分别为(1-5)(1-6)同样可以证明,它们满足(1-7)(1-8)则有限长序列的圆周共轭对称分量和圆周共轭反对称分量分别定义为:(1-9)(1-10)由于满足故(1-11)显然,长度为N的有限长序列可以分解为圆周共轭对称分量和圆周共轭反对称分量之和,和的长度皆为N。利用有限长序列与周期序列的共轭对称分量和反对称分量的关系式(1-9)和式(1-10),以及式(1-11)可以推导出DFT的一系列的对称性质(1)DFT式中表示的共轭复序列

7、。证明:DFT又因为所以DFT(2)复序列实部的DFT等于DFT的圆周共轭对称部分,即DFT证明:DFTDFT={DFT+DFT}=利用DFT的对称性可求得的DFT:设则DFT因为所以DFTDFT=五、设计内容1.用MATLAB语言编写计算序列x(n)的N点DFT的m函数文件dft.m。并与MATLAB中的内部函数文件fft.m作比较。程序如下:functionXk=dft(xn,N)iflength(xn)

8、:N-1Xk(1,k+1)=sum(xn.*exp((-1)*j*n*k*(2*pi/N)));End运算量估计:对于N=点序列进行时间抽选奇偶分解FFT计算,需分M级,每级计算N/2个蝶。每一级需N/2次复乘、N次复加,因此总共需要进行:复乘:复加:直接计算N点的DFT,需要次复乘、N(N-1)次复加。N值越大,时间抽选奇偶分解FFT算法越优越。例如当N=2048点时,时间抽选奇偶分解FFT算法比直接计算DFT速度快300多

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。