频谱峰值搜索.doc
《频谱峰值搜索.doc》由会员分享,可在线阅读,更多相关《频谱峰值搜索.doc(6页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、【精品文档】如有侵权,请联系网站删除,仅供学习与交流频谱峰值搜索.精品文档.实验名称:频谱峰值搜索实验目的:1、 熟练掌握离散信号的DFT实现方法;2、 熟练掌握Matlab实现DFT的方法,提高编程实践能力;3、 增强自我学习能力,查阅文献搜索能力;4、 掌握离散信号的时域与频域的对应关系。实验原理:1. 离散复正弦信号的DFT (1)2、Matlab主要函数fft(signal,N); signal:输入信号,N:fft的点数 函数的作用是对输入信号做N点的DFTfftshift(fft(signal,N); 将零频点移到频谱的中间3谱峰搜索算法采用一维黄金分割精搜算法5取代分级搜索过程中
2、的递归精搜。对于一维黄金分割精搜算法, 若函数有且仅有一个极大值位于区间上, 则有当时, 且。此时若在区间上选取一点, 当时, 则且 , 即极大值点将位于三元点组对应的区间上; 否则, 当 时, 则且, 即极大值点将位于三元点组 对应的区间上。在这些三元点组中, 其中间点对应的函数值都是每一轮求解过程中的最大值.这一过程下图所示, 继续对三元点组所对应的区间进行划分, 直到区间足够小, 小到以至于该区间上任何一点都可以表示函数的极大值点。下面是一维黄金分割搜索算法的基本原理。给定三元点组, 假设是,之间的一个分割。即 (2)再假设一个试探点位于区间之间, 且有 (3) 通过选取试探点后, 可以
3、将极大值点压缩到相对长度为的区间。或者是相对长度为的区间上。考虑到搜索极大值的最坏情况,应该使得下式成立: (4)可见试探点应该选为点关于区间的对称点。同样应该保证是之间的一个分割。 (5)把式( 4) 代入到式(5) 中得到如下的二次方程 (6)解得其根为( 考虑到, 舍弃另一根) (7)所以对于给定的三元点组所对应的区间上, 每次选取试探点都是位于较大的一段子区间上, 并且距离原来中间点0.3189660 的位置上。也就是说通过一次这样的试探点的选取可以使得极大值点将位于原区间0.6180340 的更小的区间上。继续这样的计算, 那么区间将变得越来越小, 小到以至于该区间上任何一点都可以表
4、示函数的极大值点。实验步骤:1、 设置输入信号的参数以及DFT变换的点数根据要求,输入信号的模拟频率为,那么采样频率满足即可,为方便观察频率最大值位置,取=2Hz。给定DFT点数为1024点,而为了使的被观察的频谱峰值在频谱图的中央,将抽样时间取在的区间,采样间隔为。其中N=512,满足采样点数为点。这样得到输入信号的表达式为 (8)2、 对信号进行DFT,并画出频谱图。(1)在MATLAB中应用fft(signal,N)对信号signal做N点的FFT;(2)分别应用函数fftshift、abs对DFT结果调整和取绝对值;(3)设置横坐标。根据和的对应关系,得到很坐标的取值范围是-1,1。在
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 频谱 峰值 搜索
限制150内