求重心频率
因为自己查资料查到头秃才终于知道怎么求重心频率,网上的资料大多不对,原来这里面涉及到一个功率谱的概念。
以轴承phm2012为例
def for_FREmeanf(y):
'''
求取重心频率
'''
pp=list()
for i in y:
N=i.size
Y = fft(i) * 2 / N
Y=np.abs(Y)
f = (Y ** 2) / len(Y)*25600 %功率谱
Y = Y[:Y.size // 2]
f = f[:f.size // 2]
m=sum(Y*f) / sum(Y)
pp.append(m)
return pp