weixin_45004370的博客

原创 重心频率

求重心频率 因为自己查资料查到头秃才终于知道怎么求重心披绿，网上的资料大多不对，原来这里面涉及到一个功率谱的概念。 以轴承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 %功率谱