- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
傅立叶变换在丫能谱数据的光滑与去噪中的应用
3.1 丫能谱的时频特征
由能谱仪器测量得到的、被复杂化的核辐射原始线谱,称为仪器谱。线谱变成了具
有高斯分布特征的形状,具有单一能量的线谱,变成连续谱。将原始谱线做一个关于谱频
的傅里叶变换 [8]得到它的谱频图。
6
图 原始谱线
7
图 频频图
3.2 傅立叶变换去噪中低通滤波器的选择
傅立叶变换采用FFT 算法。傅立叶变换法中,关键是滤波函数的选择是否适 当。
在实际应用中,为了抑制滤波器过大的边叶作用,常常采用平滑变化的函数,且使滤波函
数的两端逐渐变小而趋于零,F ()=0 时的 频率称为切断频率
。
MF
( )特别应适当的选择切断频率,切断频率过高,光滑的效果比较差,切断频率过
低,光滑过度 [9 ]。可选的滤波函数有:
1
() 高斯形滤波器(匹配滤波器):
12 2
F () e?
3-1
( )
2 k
N 2
其中 上( 为总道数,一般取 的整数幕)
2
N
是最佳的滤波器,与其信号峰有相同形状的函数。一般: A=1,c :高斯宽度,
(7= H/2.355 H
, 为半宽度。效果:无附加的虚峰,信号的有用信息均集中于峰高的
数值中。
2
() 其它函数:
MFC
0.54 MC)
0.46cos( F() (3-2 )
MFC
(3) 切断频率的确定方法及影响:
1/4 R R=5
从后 谱中找出最大幅度对应的①,即为噪声最大幅度。设定 初始值,例如设
1/4 Rf i
从后 处向前找出大于 的点,此点即为信号起主要作用的点。若找到,从该点向后找第
一个低于①的点,此点即为切断频率 MFC 若找不
到,贝U R=R-0.5,继续寻找。
3.3 傅立叶变换去噪的Matlab 语言实现何
将谱数据放于k.dat 文件中,首先打开文件,将数据其调入存里:
fid=fope n( \spe.dat,r);
date=fsca nf(fid,%f);
status=fclose(fid);
将谱数据做傅里叶变换,使信号从能域变为频域:
y= t(date,8192);
mag=abs(y);
本次做快速傅里叶变换道数已知,若在未知情况下用length 函数,可求出道数 根据实际找出
截断频率:
yj=y;
yj( : )=o ;
文档评论(0)