本帖最后由 wdhd 于 2016-9-18 13:25 编辑
out5.txt数据文件是本主题2008-7-29提供的,是4组数据,每组4580个.早就发现数据out5.txt要从1000后取分析才对,当时放下了.前几天发现有一个10000个数据的文件由前后两部分组成的,想起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)是重构误差,很小的.
图3 原信号和重构信号误差图 是什么原因造成的还要分析, 楼主也不会知道,一般也是从它处找来的.,能提供出来真心致以谢意
[ 本帖最后由 zhwang554 于 2009-4-24 21:49 编辑 ] |