光谱二阶差分Gabor展开法土壤铜铅污染鉴别

土壤是人类生存环境的重要载体,因此,土壤重金属污染问题一直备受关注。随着遥感技术的发展,高光谱遥感在土壤重金属研究中取得了大量的成果,但是,基本上是根据土壤中有机质、铁、粘土矿物等的光谱吸收特征和反演土壤中重金属含量,而不能够区分土壤重金属污染光谱之间的微弱差异。通过盆栽土壤不同浓度铜(Cu)、铅(Pb)污染实验得到不同浓度Cu和Pb污染下盆栽土壤光谱曲线、土壤含水率和有机质含量,提出了一种光谱二阶差分Gabor展开方法探测不同浓度Cu和Pb污染下土壤光谱曲线之间的微弱差异。以二阶差分为基础,首先将土壤光谱转换为稀疏光谱,然后结合土壤稀疏光谱与Gabor展开理论,在频率域中检测不同浓度土壤重金属污染光谱之间的微弱差异,因此,摆脱了单纯通过土壤光谱反射率信息反演土壤重金属含量的研究,而是对土壤重金属污染光谱信息进行时频分析,最终达到检测土壤重金属污染瞬时光谱存在的目的。结果表明：受Cu和Pb污染的盆栽土壤光谱二阶差分Gabor展开系数尺度及等高线分布有较大的差异,Cu污染的盆栽土壤光谱二阶差分Gabor展开系数尺度分布存在两个较高的峰值,且等高线在第1 800~3 600项之间稀疏分布,Pb污染的盆栽土壤光谱二阶差分Gabor展开系数尺度分布存在一个较高的峰值,且等高线在第3 200~3 600项之间密集分布;二阶差分Gabor展开法检测的土壤Cu和Pb污染结果与土壤Cu和Pb含量、土壤含水率、土壤有机质是密切相关的,由于土壤Cu和Pb含量、有机质含量、含水率的不同,土壤Cu和Pb污染二阶差分Gabor展开光谱尺度分布而不同。根据相关性分析结果,分别将土壤Cu和Pb污染划分为三组：Cu(50)~Cu(300),Cu(400)~Cu(800),Cu(1 000)以上;Pb(50)以下,Pb(100)~Pb(300),Pb(400)~Pb(1 200)。     

HHT时频分析土壤光谱的重金属铜离子污染信息提取模型

土壤中不同浓度Cu2+含量映射到土壤光谱上的信息量十分微弱,并且这些高光谱数据中也存在着难以避免的噪声,因而本研究的关键是如何在土壤光谱复杂的噪声环境中提取微弱Cu2+信息。经验模态分解算法(EMD)能够有效去除高光谱数据中的噪声,且EMD是Hilbert变换对“非线性非稳定”信号时频分析的前提,当引入Huang变换后,可利用Hilbert-Huang变换(HHT)模型时频分析高光谱数据以实现降噪处理与信息提取。通过时频的HHT分析不同浓度Cu2+污染下的土壤光谱,完成从原始光谱经EMD分解出各本征模态函数(IMF)分量的包络线、调制信号和频谱等曲线中挖掘土壤光谱的Cu2+污染信息。研究结果表明,相同浓度Cu2+污染时的土壤光谱HHT时频分析结果相同,不同浓度时则不同,所以也可依据IMF分量反演土壤Cu2+含量。因此,高光谱数据的HHT时频分析能为土壤光谱的信息挖掘、光谱诊断和Cu2+含量反演等提供一种新的方法和思路。     

烃源岩可溶烃、热解烃、总有机碳含量与光谱指数相关性研究

烃源岩的可溶烃含量(S1)、热解烃含量(S2)和总有机碳含量(TOC)是进行烃源岩研究和评价的重要指标,以往主要通过实验室分析方法获得,时效性得不到保证。为了探讨野外快速测量这些指标的可行性,采集了乌兹别克斯坦克孜勒库姆地块10件烃源岩样品,运用地面ASD光谱仪和碳硫分析仪,分别测定了样品的可见光-近红外-短波红外(Vis-NIR-SWIR)反射光谱和S1, S2和TOC含量参数。运用光谱拟合微分,获得了上述样品反射光谱的一阶导数、二价导数等微分数据;并运用线性回归方程拟合和拟合结果统计分析,研究了这些数据同S1, S2和TOC含量之间的相关性。结果表明:烃源岩S1含量与原始光谱、光谱一阶导数、二阶导数分别有1 973, 511和239个波段相关,相关系数绝对值最大值分别为0.612, 0.823和0.889,相关性程度分别为弱、较强、极强;烃源岩S2含量与原始光谱、光谱一阶导数、二阶导数分别有2, 144和205个波段相关,相关系数绝对值最大值分别为0.561, 0.867和0.926,相关性程度分别为弱、较强、极强;烃源岩TOC含量与原始光谱无相关性,与光谱一阶导数、二阶导数分别有18、 180个波段相关,相关系数绝对值最大值分别为0.882和0.879,相关性都极强。此外, S1含量与光谱二阶导数的2 012 nm波段拟合优度最高,达0.790;S2含量与光谱二阶导数的363 nm波段拟合优度最高,达0.858;TOC含量与光谱一阶导数的2 480 nm波段拟合优度最高,达0.777。结果说明,利用可见光-近红外-短波红外反射光谱反演烃源岩S1, S2和TOC含量具有比较好的效果。为研究野外烃源岩的原位检测奠定了重要的基础。     

类标准化土壤样品Cd含量高光谱定量反演

针对土壤Cd高光谱遥感反演精度低、特征波段难以有效识别的问题,以湖南省某典型重金属污染矿区为例,开展类标准化Cd污染土壤样品和自然污染土壤样品Cd含量高光谱定量反演对比研究。创新之处主要在于以相对清洁的背景土壤作为内控样本,通过实验室定量添加Cd标准溶液制备符合自然污染规律的类标准化Cd污染土壤样品方案的提出与对比实验的开展,研究实验过程还包括自然污染土壤样品野外采集,类标准化样品与自然样品土壤重金属、有机质含量及350～2 500 nm光谱反射率等的测定,光谱全要素主成分逐步回归土壤Cd含量反演建模。研究发现,类标准化土壤样品Cd含量高光谱反演模型精度(adjR2=0.87)明显高于野外自然污染土壤样品Cd含量高光谱反演建模结果(adjR2=0.39)。土壤Cd含量与土壤光谱反射率间确实存在一定的响应关系,但各光谱波段对土壤Cd含量及其变化产生响应的程度存在差异,其中1 000,2 000以及2 300 nm波段光谱响应信号较强。该研究创新引入的类标准化Cd污染土壤样品制备方法有助于深入探索重金属Cd污染土壤光谱特征响应规律并发现重金属Cd含量及其变化的真正指示性特征波段,可为多要素混淆污染模式下模拟反演土壤重金属含量提供先验知识。     

微分算法的艾比湖湿地自然保护区土壤有机质多光谱建模

针对以往利用高光谱数据来来反演土壤有机质（SOM）的可行性与可靠性,结合微分处理对光谱数据信息提取的高效性,提出了直接对多光谱遥感影像进行微分处理就可得出SOM建模研究,旨在为今后SOM速测提供参考。采用Landsat 8_OLI 多光谱遥感影像数据,对多光谱遥感影像进行辐射定标、几何校正、大气校正、镶嵌和裁剪,运用IDL软件对影像进行一阶微分处理和二阶微分处理,发现一阶微分图像能够更好地表达地物的真实情况,更好地区别水体与土壤。原始遥感影像包含大量的信息其中还包括噪声,通过微分处理后的遥感影像剔出了原始影像中反射率值突兀变化的部分。在研究区采用五点法采集土壤样品。室内实验用重铬酸钾氧化-容量法测得SOM数据。多光谱数据结合地面实测SOM数据,分析SOM与多光谱数据反射率的关系,发现一阶微分处理后的遥感数据与SOM含量的相关性存在敏感波段,说明一阶微分处理可以将原始遥感图像数据在多光谱范围内的一些隐含的土壤有机质信息释放出来。选取相关性高的数据建立基于原始遥感数据、一阶微分数据、二阶微分数据的单波段多光谱线性模型和多波段多光谱线性模型,选取最优模型来估算和反演土壤有机质含量。结论如下：（1）通过对原始影像进行微分处理发现,微分处理后的影像变化明显,一阶微分处理的影像噪声降低,更加突出了影像中土壤有机质隐藏的信息。二阶微分处理的影像抑制了土壤有机质信息。（2）原始遥感影像各波段数据对土壤有机质含量的相关性较低,一阶微分处理后的遥感影像数据反映出土壤有机质敏感波段即部分波段数据相关性明显高于原始数据,二阶微分处理后的遥感影像各波段数据对土壤有机质含量的相关性较弱。（3）多波段建模效果要优于单波段建模;一阶微分多波段模型预测精度最优,其模型的决定系数和模型拟合的决定系数分别为0.898和0.854,该模型对估算研究区内的SOM含量效果较好;综合比较了单波段模型和多波段模型的拟合精度,发现无论在单波段模型还是多波段模型一阶微分处理后的模型都具有更好的预测能力。（4）基于一阶微分多波段模型对研究区SOM进行反演,反演结果与实际情况相符合,对干旱区SOM含量制图提供了切实可行的方法和参考。     

不同铅浓度土壤XRF光谱的频率域去噪及规律研究

近年来,全球工业快速发展和城市化的推进引发了一系列环境问题,其中土壤重金属铅(Pb)污染引起了科研人员的广泛关注。X射线荧光(XRF)光谱技术具有成本低、分析速度快、适合大面积监测等优势,已被广泛应用在土壤污染检测与生态环境保护等多个领域,发展前景广阔,因此充分挖掘土壤XRF光谱信息具有很强的现实意义,可为高效地进行土壤污染检测、土壤环境生态指标参量反演及矿区重金属污染早期预警等提供解决方案。目前关于XRF光谱研究大多集中在土壤重金属浓度测量结果的精度评价和环境质量评估,而对于土壤XRF光谱差异特征变化的深入分析研究较少。时频分析可将时域复杂信号变换到频率域空间,从频率域角度检测光谱、信号中存在的异常信息,是一种有效的谱差异特征变化检测方法,其中,谐波分析(HA)可用于电磁信号去噪;平滑伪魏格纳分布(SPWVD)可筛选合适的基函数,突出信号时频局域性细节特征。首先采用HA方法探究不同Pb浓度土壤XRF光谱的去噪效果,然后利用光谱的SPWVD研究实地采样土壤样品的去噪XRF光谱局部规律。研究结果表明:当谐波分解次数为400时,土壤XRF光谱去燥效果较好并节省了时间,且保留了光谱的特征。土壤样品Pb浓度与XRF光谱的SPWVD在400和600～700波段序列上的频率峰值的分布具有一定的规律性,根据此规律性可识别该地区土壤的Pb浓度超标:全部实地采样土壤中,能够识别75%的Pb浓度不超标的样品,在波段序列为400附近有1个较高的频率峰(频率小于400 Hz)或有2个很强的频率峰值(频率大于400 Hz);全部实地采样土壤中,能够识别79.17%的Pb浓度超标的样品,在波段序列为400附近有1个很强的频率峰(频率大于400 Hz),且在600～700波段序列之间有3个层次明显的频率峰分布,据此推断该地区土壤Pb浓度超标的XRF光谱特征波段区间为6.42和9.42～10.92 keV。因此,研究通过引入时频分析法实现了对土壤XRF光谱频域甄析和可视化显现,为深度挖掘Pb污染谱规律特征及异常信息提供一种新思路。     

关联类标准化样品特征波段的矿区土壤重金属Pb高光谱反演

目前土壤重金属Pb的高光谱反演实验样本多源于研究区采样,忽略了自然样本中复杂组分的光谱混杂问题,所建立的模型精度偏低,重金属光谱响应机理解释仍不明确。基于此,选取湖南省某铅锌矿区采集相对清洁背景土壤,运用控制变量方法,制作Pb类标准化样品40个,同时采集矿区86个自然污染样本。首先利用类标准化样品建立偏最小二乘回归(PLSR)定量反演模型,提取Pb光谱反射特征波段,然后基于特征波段进行矿区自然污染样本土壤Pb定量回归建模,对比全波段建模结果,验证特征波段的有效性。实验结果表明:(1)在背景土壤性质稳定的前提下,制作的类标准化样品具有总体一致的光谱形态,同时Pb含量作为引起反射率变化的唯一驱动因子证实了土壤中重金属Pb光谱响应信号的存在;(2)相较于自然样本光谱反演结果,类标准化样品能显著提升土壤Pb含量模型估算能力,决定系数(R■)为0.85,相对分析误差(RPD)2.30,有效排除其他组分对于光谱建模的干扰,模型性能良好;使用变量重要性投影指标(VIP)与PLSR系数提取的类标准化样品特征波段集中在970～1 000, 1 700～2 080和2 220～2 400 nm区间;(3)与全波段建模对比,基于特征波段的矿区自然污染样本的PLSR建模结果,R■由0.32提升至0.55, RPD由1.20提升至1.44,表明特征波段能够较好表征Pb的响应信号,波段提取有助于过滤噪声,减少数据冗余同时提高模型反演精度。提出的基于类标准化样品特征波段反演自然土壤样本Pb含量,有效解决自然样本成分复杂,重金属光谱信号微弱的问题,丰富了土壤重金属遥感监测理论与应用实例。     

土壤导数光谱小波去噪与有机质吸收特征提取

使用高光谱仪ASDFieldSpec在波长范围400～1000nm内,采集有机质含量水平不同的土壤反射光谱数据,之后求取一阶导数(一阶导数光谱)并进行小波去噪;通过数值积分计算、相关分析,从去噪后的一阶导数光谱中提取吸收面积表征有机质含量变化。结果显示:(1)土壤一阶导数光谱中含有大量噪声,致使光谱曲线轮廓和有机质的吸收特征难以识别。(2)当小波分解层数J=3时,小波去噪在平滑土壤一阶导数光谱曲线的同时,较好的保持了曲线上多个敏感波段的光谱响应特征。(3)从去噪后的一阶导数光谱中提取的吸收面积S(538,586)与有机质含量的相关系数为-0.8763,较好的反映了土壤有机质含量的变化情况,可以用于土壤有机质含量测算。     

甘蔗初压汁锤度近红外光谱分析的波段优选

用近红外透射光谱技术和偏最小二乘法建立甘蔗初压汁锤度的定量分析模型。用多元散射校正和Savitzky-Golay平滑化方法对原谱、一阶导数谱和二阶导数谱进行预处理。选取400—1100nm、1100—1400nm、1550—1850nm、2100—2350nm以及全谱（400—2500nm）5个波段,每个波段分别采用原光谱、一阶导数谱、二阶导数谱。同时调整Savitzky-Golay平滑点数和PLS因子数,通过多次PLS数值实验比较,按照预测效果确定每个模型的最优平滑点数、因子数,再从中选优。结果表明,采用二阶导数谱的1550-1850nm波段的定标效果最好,模型的预测相关系数、预测均方根偏差、相对预测均方根偏差分别为0.896,0.406,1.95%。1550—1850nm波段可以代替全谱波段得到好的定标效果,为设计小型专用近红外分析仪器提供依据。     

廿蔗初压汁锤度近红外光谱分析的波段优选

用近红外透射光谱技术和偏最小二乘法建立甘蔗初压汁锤度的定量分析模型.用多元散射校正和Savitzky-Golay平滑化方法对原谱、一阶导数谱和二阶导数谱进行预处理.选取400-1100nm、1100-1400nm、1550-1850nm、2100-2350nm以及全谱(400-2500nm)5个波段,每个波段分别采用原光谱、一阶导数谱、二阶导数谱.同时调整Sayitzky-Golay平滑点数和PLS因子数,通过多次PLS数值实验比较,按照预测效果确定每个模型的最优平滑点数、因子数,再从中选优.结果表明,采用二阶导数谱的1550-1850nm波段的定标效果最好,模型的预测相关系数、预测均方根偏差、相对预测均方根偏差分别为0.896,0.406,1.95%.1550-1850nm波段可以代替全谱波段得到好的定标效果,为设计小型专用近红外分析仪器提供依据.     

Hilbert-Huang变换分析THz脉冲信号的时频特性

提出一种基于Hilbert-Huang变换的THz时域光谱分析方法,将THz时域脉冲信号分解成有限数目的单分量信号之和,利用Hilbert变换求得瞬时频率来获得幅值的时频分布——Hilbert-Huang变换谱,实现了通过水蒸气的THz脉冲信号的时频分析,揭示了THz波与水蒸气相互作用的频谱时域分布特性,并与基于小波变换的时频图进行了对比分析。结果表明,该方法可以同时提高THz脉冲时频分布的时间分辨率和频率分辨率,具有局部化分析和自适应选择的特点,还能直观地表现出各频率成分之间的相对时间延迟。     

铜污染土壤光谱特征及弱差信息的LH-PSD甄测模型

土壤受重金属污染后，会影响农作物及食品安全，危及人体健康，因此寻找快速、高效甄测土壤重金属污染信息的方法尤为关键。传统化学分析方法存在过程繁杂、费时耗力等缺点，而高光谱遥感因光谱分辨率高、信息量大、快速无损等特点在环境监测等应用方面优势明显。由于电磁遥感信号反射、辐射过程复杂，通过仪器获取的土壤高光谱数据难以直接解析出重金属污染信息，因而，研究并寻求一种能够有效挖掘土壤重金属污染信息的方法对高光谱遥感监测污染意义重大。不同浓度铜(Cu)污染会使土壤理化性质改变，引发土壤光谱产生微弱变化，该研究目的是对Cu污染土壤光谱的特征及弱差信息进行识别、提取与分析，进而挖掘光谱中的重金属污染信息。采用包络线去除(CR)对光谱进行预处理，通过定义局部极大值均值(LMM)与半波高(HWH)，结合时频分析的短时傅里叶变换(STFT)及能量谱密度(PSD)，构建LH-PSD甄测模型。通过模型对极相似土壤光谱进行处理，所获PSD分布图使光谱间的微弱差异可视化显现，并显著区分了相似光谱，验证了模型对光谱特征及弱差信息的甄别能力和有效性。同时应用该模型，对不同Cu污染梯度的土壤实验光谱进行重金属污染信息的提取与分析，研究结果表明，LH-PSD甄测模型中，LMM与HWH可有效提取光谱间差异特征并以阶梯状显现。经模型处理后得到的可视化PSD分布图能直观定性判别土壤是否受重金属Cu污染，即当土壤受重金属Cu污染后，相同采样频数下，在频率为100与600 Hz附近PSD分布会出现明显空缺分离，随着Cu污染浓度的增加，在100~600 Hz之间PSD的分布呈逐渐稀疏态势。能量值E可定量化监测土壤Cu污染程度，即随着土壤中Cu污染浓度的增加，E值呈下降趋势，且与Cu含量的相关系数达到-0.910 5，显著相关。为检验模型的可靠性，研究结合栽种玉米作物的土壤光谱，经LH-PSD甄测模型对其进行分析，所得可视化的PSD分布图结果与实验分析中基本一致，且能量值E的监测结果与土壤中Cu含量相关系数达到-0.973 9，相关性显著，验证了模型的可靠性。因此，LH-PSD甄测模型实现了对土壤光谱从光谱域到时频域的甄析，为深度挖掘重金属污染的光谱特征及弱差信息提供一种新思路。     

Resonance identification in loudspeaker driver units: A comparison of techniques

A number of techniques with the scope of identifying loudspeaker cone resonances have been examined. Namely, the waterfall plot, the wavelet transform and the empirical mode decomposition scheme were compared on the basis of time-frequency resolution and damping estimation. The commonly used waterfall plot is only acceptable at the upper range of the acoustic spectrum. The wavelet transform is especially well suited for the analysis of transient signals from loudspeakers and is a significant improvement over the waterfall method. The newly developed EMD scheme has the highest potential in separating the modal components. By application of the EMD both instantaneous amplitude and frequency can be accurately determined.     

EMD时频分析拉曼光谱和近红外光谱

用时频方法分析拉曼光谱和近红外光谱。经验模态分解光谱成为特征模态分量,模态分量比重计算显示拉曼光谱能量均匀分布于各个分量,而近红外光谱的低阶特征模态分量只承载了较少的原光谱有效信息。真实光谱和数值实验均显示,经验模态分解视拉曼光谱为调幅信号,具有高频能量吸附特性; 视近红外光谱为调频信号,在一阶特征模态分量中可以较好实现高频窄带解调。一阶特征模态分量希尔伯特变换显示,经验模态分解拉曼光谱时易出现模态混叠现象。进一步在时频域分析玉米叶片近红外光谱,经验模态分解后截掉低能量的一、二阶分量,用剩余特征模态分量重构光谱信号,均方根误差为1.001 1,相关系数为0.981 3,两个指标反映出重构精度较高; 分解趋势项表明在近红外光波段,吸光度随着波长的减小呈现递增趋势; 特征模态分量的希尔伯特变换显示,657 cm-1是碱胁迫光谱特有频率,可作为碱胁迫光谱特征频率来辨识。     

应用经验模态分解和独立成分分析解译岩矿光谱

岩矿光谱由多种矿物光谱混合而成,解译岩矿光谱能够得到岩矿的组分信息,且该方法具有快速、方便、不损坏样品的特点。经验模态分解(empirical mode decomposition, EMD)不能直接分离出混合信号中的源信号,独立成分分析(independent component analysis, ICA)要求混合信号数目不小于其所包括的源信号数目。将EMD和ICA两种方法相融合,首先用EMD分解混合信号得到本征模态函数(intrinsic mode function, IMF),再选择一定数目的IMF与混合信号一起组成ICA的输入数据矩阵,经过ICA运算可以获取单一混合信号中的源信号信息,克服了EMD和ICA两种方法各自的缺陷。研究表明,综合应用EMD和ICA方法可以获取单一混合信号中的源信号信息,混合信号中源信号含量越大,得到的源信号近似值越理想。参与ICA分离的IMF数目决定了分离得到的源信号近似值的数目,并且选择的IMF与混合信号相关系数越大,得到的源信号近似值越理想。运用该方法定量分析岩矿光谱,可以获取组成岩矿的矿物信息,比较适用于野外作业岩矿的快速分析鉴定及成分初步分析。     

闪电回击场多尺度分析

利用聚合经验模式分解(EEMD)方法对闪电负地闪回击电场波形进行了分析并得到一系列具有物理意义的内禀模态函数 (IMF)。 结果表明,IMF1是高频噪声项,趋势项R主要是静电场作用的结果。 IMF2~IMF11表现出闪电通道的不同尺度放电特征。通过计算IMF分形维数,发现越高阶的IMF对应着越大的分形维数,表明越高阶的IMF对应着越曲折的闪电放电通道。相比于基于经验模式分解（EMD）的希尔伯特-黄变换（HHT）谱和小波谱,基于EEMD的HHT谱能改善模态混叠现象,EEMD更适合于观察像闪电这种非平稳非线性信号。     

Time-frequency analysis of time-varying modulated signals based on improved energy separation by iterative generalized demodulation

Energy separation algorithm is good at tracking instantaneous changes in frequency and amplitude of modulated signals, but it is subject to the constraints of mono-component and narrow band. In most cases, time-varying modulated vibration signals of machinery consist of multiple components, and have so complicated instantaneous frequency trajectories on time-frequency plane that they overlap in frequency domain. For such signals, conventional filters fail to obtain mono-components of narrow band, and their rectangular decomposition of time-frequency plane may split instantaneous frequency trajectories thus resulting in information loss. Regarding the advantage of generalized demodulation method in decomposing multi-component signals into mono-components, an iterative generalized demodulation method is used as a preprocessing tool to separate signals into mono-components, so as to satisfy the requirements by energy separation algorithm. By this improvement, energy separation algorithm can be generalized to a broad range of signals, as long as the instantaneous frequency trajectories of signal components do not intersect on time-frequency plane. Due to the good adaptability of energy separation algorithm to instantaneous changes in signals and the mono-component decomposition nature of generalized demodulation, the derived time-frequency energy distribution has fine resolution and is free from cross term interferences. The good performance of the proposed time-frequency analysis is illustrated by analyses of a simulated signal and the on-site recorded nonstationary vibration signal of a hydroturbine rotor during a shut-down transient process, showing that it has potential to analyze time-varying modulated signals of multi-components.     

基于EMD的谱峭度方法在滚动轴承故障检测中的应用

传统谱峭度方法通常采用基于短时傅里叶变换(Short Time Fourier Transform,STFT)的峭度图方法来实现。针对STFT不能保证对瞬态脉冲这种高度非平稳信号最优的分解效果的缺点,提出一种基于经验模式分解(Empirical Mode Decomposition ,EMD)的谱峭度方法。该方法首先利用EMD和Hilbert变换得到信号的时频分布,然后将信号的时频分布按照不同层数分成若干频段,通过计算各频段的峭度值得到相应的峭度图,再根据峭度最大原则选择滤波频段进行带通滤波,最后对滤波信号采用包络分析确定故障信息。实验结果表明：相比传统基于STFT的谱峭度方法,本文方法更能准确的获得轴承加速度信号的故障特征频率信息。