[发明专利]一种短时傅里叶变换窗长的自适应选择方法

权利要求书

1.一种短时傅里叶变换窗长的自适应选择方法，其特征在于，对待分析的振动信号进行分段处理，每段信号采用不同的窗长，信号的分段及每段信号的窗长都根据信号特征自适应得到，方法包括以下步骤：

步骤1：对振动信号进行长度为L的滑动加窗预处理，L的取值为3～16个采样点，计算每一个窗口内信号的短时能量，得到整个信号的短时能量曲线；

步骤2：去除短时能量的平均值和趋势项，并对短时能量数据进行平滑滤波处理，得到滤波后的短时能量曲线；

步骤3：在滤波后的短时能量曲线中寻找各个局部极值点，并分别以能量曲线中最大值绝对值的1/k1和最小值绝对值的1/k2作为极大值和极小值的筛选阈值，其中k1、k2都为大于2.0的一个实数；

步骤4：将所有局部极值点中绝对值小于筛选阈值的极值点去除，仅保留绝对值大于筛选阈值的极大值点和极小值点；

步骤5：计算所有相邻极值点连线的斜率，得到极值点连线的斜率曲线；

步骤6：设定斜率绝对值的阈值th，从左向右依次比较各极值点连线的斜率绝对值，将斜率绝对值小于th的连续范围内的数据作为一段，将斜率绝对值大于th的连续范围内的数据作为另一段，若所有斜率都大于或小于th，则整个信号作为一段，不需再分段；

步骤7：采用最小信息熵法分别估算各段的最优窗长，即对于各段信号，依次计算窗长度从WL1到WL2时的信息熵S(WL)，然后寻找其中S(WL)值最小点对应的窗长，并将此窗长作为该分段内信号的最佳窗长；设短时傅里叶变换的谱图为：

G(t,f)＝|STFT_x(t,f)|² (1)

式(1)中t表示时间序列，f表示频率序列，x表示输入信号；设P(A_u,v)为谱图中第u行、第v列对应的点出现的概率，其计算公式为：

式(2)中m和n表示谱图中共含有m行n列个点；根据公式(1)和(2)，信息熵的计算公式为：

其中，a为大于1的正整数；

步骤8：采用各段的最优窗长分别对各段的数据进行短时傅里叶变换，再将各段结果进行组合，得到全部信号的短时傅立叶变换谱图。