zhwang554 发表于 2009-3-30 08:27

脉动压力数据的谱校正

    楼主提供的脉动压力数据及其频谱可用离散频谱校正法求出其校正后的频率值和振蝠值. 画出校正后的脉动压力数据的功率谱. 这是二年前的老数据,但是真正的实际数据,值得分析与校正.    脉动压力数据的fft/apfft谱如图1a所示, 兰色为fft对数振幅谱,红色为apfft对数振幅谱, 在f>23后fft和apfft振幅谱形状类同,表示它们是稳态正弦信号,可以校正.   图1b为相位谱, 兰色为fft相位谱,红色为apfft相位谱,从红色相位谱可见,在23<f<50频段内,apfft相位谱是水平阶梯形状, 每一个振幅峰值对应的相位谱都是水平的,表明每一峰值只有一个频率成分,可以校正. 只有最后一个f=48处峰值对应的红色相位线是斜的,表明此处不是一个频率成分,不能用单峰校正法校正.
    图1b红色的apfft的水平相位谱除了不须校正外,还有一个用途,就是判断峰值处是否是单一频率成分,每个振幅谱峰值对应的相位谙谱若是水平的,则这个峰值只有一个频率成分,若一个振幅谱峰值对应的相位谱不是水平的,是斜线,则这个峰值有二个以上频率成分,这个判断准则在图1中可看得很清楚.这是apfft水平相位谱的一个有实用价值的性能.
    从图1a及b还可看出,f<23的fft振幅谱和apfft振幅谱形状不重合,apfft相位谱也不是水平阶梯形状,表明这一频段不是稳态信号,无法校正.
                               图1 脉动压力数据的fft和apfft的振幅谱和相位谱
用fft/apfft校正法在23<f<50频段内得到的校正频率,振幅,相位值如下表及图2

   校正频率值          校正振幅值          校正相位值            频差
f1=2.3399e+001   a1=2.7823e-003      p1=1.7839e+002
f2=2.4507e+001   a2=3.1431e-003      p2=1.0102e+002         f2-f1=1.1080
f3=2.5665e+001   a3=3.3657e-003      p3=1.6957e+002         f3-f2=1.1583
f4 =2.6878e+001    a4=3.8814e-003      p4=3.4439e+002         f4-f3= 1.2125
f5 =2.8150e+001    a5=3.3750e-003      p5=5.4287e+001         f5-f4=1.2720
f6=2.9481e+001    a6=2.3415e-003       p6=3.0470e+001         f6-f5=1.3309
f7=3.0876e+001   a7=2.3896e-003      p7=3.1526e+002         f7-f6=1.3951
f8=3.2339e+001   a8=1.9805e-003      p8=2.5739e+001         f8-f7=1.4634
f9 = 3.3865e+001   a9=1.8092e-003      p9=1.3274e+002         f9-f8=1.5257
f10=3.5466e+001   a10=1.8853e-003    p10=7.4403e+001   f10-f9=1.6014
f11 =3.7143e+001a11=1.5073e-003   p11=2.6500e+002      f11-f10=1.6767
f12 =3.8900e+001a12=1.4167e-003   p12=2.3868e+002      f12-f11=1.7573
f13 =4.0741e+001a13=1.0840e-003   p13=1.6831e+002      f13-f12=1.8404
f14 =4.2667e+001a14=9.8551e-004   p14=1.7911e+002      f14-f13=1.9263
f15 =4.4685e+001a15=8.6638e-004   p15=1.5750e+002      f15-f14=2.0178
f16=4.6799e+001   a16=7.0902e-004   p16=2.3486e+002      f16-f15=2.1139
f17=4.9202e+001   a17=3.9827e-004   p17=3.5487e+002      f17-f16=2.038


    图2a为23<f<50之间二个振蝠峰的频差曲线(上表第4列),它是一个均匀上升线,只有最后一点f=44处例外(上面己指出此处不是一个频率成份,不能用单峰校正法校正).频差曲线是一个均匀上升线在生理上如何理解,其增长率有物理意义吗?虽然从未校正的振蝠谱图1a以及沙发贴子的3个功率谱图中也可以看出二峰间频率间隔在增大,但校正后的图2a将增大细节给出了,相当均匀有规则(但不是线性)的增大,还可有具体数据.
    图2b为校正后的各频率成分的振幅值,在10(-3)数量级, 曲线起伏和 图1a的振幅谱一致

                        图2 脉动压力谱的频差和校正振蝠值

[ 本帖最后由 zhwang554 于 2009-3-30 09:09 编辑 ]

zhwang554 发表于 2009-3-30 10:35

补充16楼的贴子--脉动压力功率谱

    16楼的贴子中图2b的振幅值图中,横座标是23<f<50频段内17个频率成份的序号.
    如果把横座标改为0<f<50Hz频率轴,则得如图3校正后脉动压力数据的功率谱估计.这个功率谱的横座标频率是校正后的频率值,纵座标是校正后振幅值的平方的对数.
    如果校正正确的话,这是脉动压力数据在23<f<50频段内去噪后的真实功率谱,它和沙发贴子中的3种功率谱比较一下到顶有意思的.除了注意纵座标的分贝数,还要注意横座标的频率值,在图3中横座标的频率值都是校正后的频率值,每个频率成份只有一条谱线.而不是有的方法由於泄漏形成的一堆谱线.从图3还可先见,二条谱线间的间隔随频率增加而增加.




                        图3 校正后的脉动压力数据的功率谱


[ 本帖最后由 zhwang554 于 2009-3-30 10:43 编辑 ]

lonsvm 发表于 2012-2-28 09:11

学习了,调试看看
页: 1 [2]
查看完整版本: 大家帮忙看看我这个谱分析的结果!谢谢!