数据是由long9999  2009-5-25 10:51  在主题  滚动轴承外圈故障信号处理（求助）中提供的data130.mat,采样点数是121991，采样频率fs=12000，故障频率是107HZ，

http://www.chinavib.com/forum/thread-82283-1-1.html

data130.mat的N=12000阶fft/apfft谱分析如图一所示. 图一c的对数振幅谱中,兰色fft和红色apfft谱中显示了在109,216,324,432,539,647有6个重合双峰(其它省略),它们对应的图一b中apfft红色相位谱由於太密,看不清, 在图一a1-6中将6个峰值附近(峰值-10到峰值+10)的相位谱放大了

从图一a1可见在108-109, 从图一a2可见在214-218, 从图一a3可见在324-325, 从图一a4可见在431-433, 从图一a5可见在538-540, 从图一a6可见在646-647的相位谱都是水平的,即6个小图正中总有一段水平,图一表明fft/apfft双峰顶一平顶的特性正确.也说明都是单一成分的平稳信号,可以校正.

            图一 数据data130.mat的N=12000阶fft/apfft谱分析图

利用fft/apfft校正法可算出前6个谐波的正确频率,振幅和相位如下表一,  第1行为峰值序号, 第2行为校正振幅值Am, 第3行为校正相位值Pm, 第4行为校正频率值fm,  即6个谐波的校正频率值,  第5行为基频倍频值fm0, 即107.62*[1 2 3 4 5 6], 从表一可见fm=fm0, 即校正谐波频率成倍频关系, 说明校正各谐波频率值正确. 基频107.62Hz,与楼主提供相符.

                   表一  滚动轴承外圈故障信号校正后的谐波参数

1   峰值序号      109       216        324         432          539        647
2  校正振幅值Am   0.084753  0.064767   0.041433    0.028696     0.034052   0.028387

3  校正相位值pm   300.05    238.81     195.71      166.8        114.82     55.178

4  校正频率值fm   107.62    215.28     322.89      430.51       538.19     645.82
5  基频倍频值fm0  107.62    215.24     322.86      430.48       538.1      645.72

 

    数据文件data130.mat的采样点总数是121991,很大,而N=12000阶fft/apfft谱分析只须23999个数据,所以可以在data130.mat中取不同起点作谱分析,其6个峰值点都相同,频率校正值和振幅校正值差不多,而apfft相位谱虽差别很大,但双峰总对上一段水平相位.

    图二每一行6个小图为不同起点数据段峰值附近的相位谱,峰值分别109,216,324,432,539,647,仔细观察,每一个小图中间总有一个小平台.这是apfft相位谱的一个特点

                   图二 不同数据段峰值附近的apfft相位谱