|
悬臂梁脉冲响应数据的fft/apfft谱分析
采样数据.txt 是前后2部分数据组成的,每一部分都是稳态的,可分别进行fft/apfft较正
前段数据的fft/apfft谱如图1所示.
图1a中兰色为N=1024阶fft对数振幅谱a1,红色为apfft对数振幅谱a2,a1和a2的频谱形状相同的信号为稳态信号,
图1b中兰色为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
对照p1-p5和图1b的相位值,是一一对应的.对照6楼中图1的频率值和振幅值,也是一一对应的.
图1 前段悬臂梁的脉冲响应数据的N=1024阶fft/apfft谱分析
后段数据的fft/apfft谱如图1所示.
图1a中兰色为N=1024阶fft对数振幅谱a1,红色为apfft对数振幅谱a2,a1和a2的频谱形状相同的信号为稳态信号,
图1b中兰色为fft相位谱a1,红色为apfft相位谱,在5个振幅以峰值处apfft相位谱都是水平的,说明都只有一个频率成分,都可以校正.
校正后的5个成分的频率f,振幅a和相位p值为
f1 =19.432 a1=653.14 p1=2.1554
f2 =121.45 a2=214.52 p2 =57.382
f3 =339.63 a3=1019.2 p3 =151.55
f4 =665.3 a4=208.34 p4 =137.15
f5 =1099.3 a5=3.3859 p5 =106.84
对照p1-p5和图1b的相位值,是一一对应的.对照6楼中图1的频率值和振幅值,也是一一对应的.
对照前段和后段的f1-f5,是不同的,后段比前段频率都高,频率越大差距越大.
图2 后段悬臂梁的脉冲响应数据的N=1024阶fft/apfft谱分析
由於采样数据.txt 共有10000个数据,前后分别为5000个数据,可以分别作N=2048阶fft/apfft谱分析. N=1024和N=2048阶二种都作了,比较后5个频率校正值基本相同,所以用N=1024阶fft/apfft即可, 图1b和图2b的水平相位特性很明显.
如果将全部数据做N=4096阶fft/apfft谱分析,得数据的fft/apfft谱如图3所示.
图3a中兰色为N=1024阶fft对数振幅谱a1,红色为apfft对数振幅谱a2,a1和a2的频谱形状大不相同的,说明信号为非稳态信号,
图3b中兰色为fft相位谱a1,红色为apfft相位谱,在5个振幅以峰值处apfft相位谱都不是水平的,说明都不只有一个频率成分,都不可以校正.正是图3发现整个数据是非稳态的.
图3 全部数据悬臂梁的脉冲响应数据的N=4096阶fft/apfft谱分析
[ 本帖最后由 zhwang554 于 2009-4-21 18:49 编辑 ] |
|