采样数据.txt 是前后2部分数据组成的,每一部分都是稳态的,可分别进行fft/apfft较正
   前段数据的fft/apfft谱如图3所示.
   图3b中兰色为N=1024阶fft对数振幅谱a1,红色为apfft对数振幅谱a2,a1和a2的频谱形状相同的信号为稳态信号,
   图3a中兰色为fft相位谱a1,红色为apfft相位谱,在5个振幅以峰值处apfft相位谱都是水平的,说明都只有一个频率成分,都可以校正.

校正后的5个成分的频率f,振幅a和相位p值为
f1 =18.117  a1=869.53  p1=187.02
f2 =113.28  a2=296.71  p2=84.166
f3 =317.07  a3=1467.4  p3=20.983
f4 =620.88  a4=343.78  p4=30.945
f5 =1026.3  a5=7.7855  p5=27.063

 

在20009-4-11日志<如何读fft/apfft谱(一切尽在图中)>中,强调"要作分析加校正,先看它的时域波形,取一段都有数据的做分析".还遇到也是振动论坛中的一个贴子中的数据,是约5000采样,但前1000个是接上去的,分析取数据就要以1000以后才对.反正fft/apfft谱分析中,只要看fft振幅谱a1和apfft振幅谱a2形状是否重合,即可判定是否稳态信号.

         还有一个拼接数据实例, 2008-7-29在网页

http://www.chinavib.com/forum/redirect.php?fid=47&tid=41727&goto=nextoldset

上的out5.txt数据文件,是4组数据,每组4580个.也早就发现数据out5.txt要从1000后取分析才对,当时放下了.前几天发现有一个10000个数据的文件采样数据.txt 由前后两部分组成的,想起out5.txt的前1000个数据和后面的不一样.今从计算机中找到原处理图表,就写出来供参考.

用N=1024阶fft/apfft谱分析笫4组数据从100开始的笫4组数据如图1.

   图1a中兰色为fft相位谱a1,红色为apfft相位谱,在20和39振幅峰值处apfft相位谱都不是水平的,说明都不只有一个频率成分,都不可以校正.
    图1b中兰色为N=1024阶fft对数振幅谱a1,红色为apfft对数振幅谱a2,除了20频率处a1和a2重合, 在最大29频率处a1和a2不相同,66 和75处a1和a2的频谱形状大不相同的,说明信号为非稳态信号. 不能用fft/apfft谱分析和校正.

图1  笫4组数据从100开始的N=1024阶fft/apfft谱分析图

      用N=1024阶fft/apfft谱分析笫4组数据从1100开始的笫4组数据如图2
      图2a中兰色为fft相位谱a1,红色为apfft相位谱,在20和39等振幅峰值处apfft相位谱都是水平的,说明都只有一个频率成分,都可以校正.所以整个数据是稳态的,
           图2b中兰色为N=1024阶fft对数振幅谱a1,红色为apfft对数振幅谱a2, 全部振幅峰值处a1和a2重合,说明信号为稳态信号.能用fft/apfft谱分析和校正.
       校正后三个主要的频率成分的3个参数如下
      fm =  1.8777e+001    3.7572e+001    7.5149e+001
      Am =  1.3070e-001    6.6257e-001    8.5679e-002
   phase = -1.8790e+000    1.9721e+000    2.4062e+000

 

          图2 笫4组数据从1100开始的N=1024阶fft/apfft谱分析图
    至於重构,把各个频率成分相加起来即可.如图3中,图3(a)中为1000后的原波形图,3(b)中兰色为原波形,红色为重构波形.二者吻合,反映fft/apfft法分析的3个频率成分的9个参数正确.图3(c)是重构误差,很小的.