如何实现时域重构？
现有经过趋势项消除及平滑处理后的加速度振动时域信号，在各位大虾们的帮助下，使用频谱比值校正法后的数据如下：XfCorrect1 =73.2019
3.3798
-85.3080
69.1629
1.1844 -105.9674
71.5453
1.0721
139.4704
XfCorrectW1 =
73.2055
4.4180 -109.5894
68.5150
1.4475
-60.9038
71.4856
1.3020
151.7981
XfCorrect2 =
73.1984
3.4438
-87.1826
36.4348
2.3570
-6.9087
69.1557
1.2512 -105.9231
XfCorrectW2 =
73.2037
4.4914 -112.4902
36.4814
3.0521
-7.9999
68.8192
1.4956
129.4980
XfCorrect3 =
73.2056
3.9331
-87.9453
69.1672
1.3885
161.1317
71.5434
1.2511
137.9774
XfCorrectW3 =
73.2076
5.1428 -111.4920
68.5170
1.6659
-65.5076
71.4919
1.5203
146.4510
XfCorrect4 =
73.1910
4.3912
-84.4465
69.1307
1.6531
-93.5139
75.0307
1.3036
81.6861
XfCorrectW4 =
73.2008
5.7094 -112.6594
68.8174
1.9731
129.8562
71.4816
1.7889
148.7992
请问如何在此基础上进行加速度振动信号的时域重构？求各位路过的高手指点一二，在此谢啦~~~~！！ 把各个频率成分相加起来
数据out5.txt前1000和后面部分组成不同
本帖最后由 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 编辑 ]
页:
[1]