本发明属于信号处理领域,尤其涉及一种窗长自适应的短时傅里叶变换方法。
背景技术:
1、在信号处理领域,短时傅里叶变换可以同时反映信号在时域和频域的特性,是对非平稳信号进行时频分析的经典方法。由于测不准原理的约束,短时傅里叶变换结果的时间分辨率和频率分辨率之间存在矛盾。短时傅里叶变换窗长的选择决定了时频图的时间分辨率和频率分辨率:当窗长较大时,频率分辨率高而时间分辨率低;当窗长较小时,时间分辨率高而频率分辨率低。因此,短时傅里叶变换窗长的确定在信号处理领域有着重要研究价值与意义。
2、在传统的短时傅里叶变换中,往往根据经验选择窗长,很难得到高时频分辨率的结果。并且一般以固定窗长对整个信号进行时频分析,对于时变信号,固定的时频分辨率难以适应信号在不同时间段的特性,分析效果较差。
技术实现思路
1、发明目的:针对以上现有技术存在的问题,本发明提出一种窗长自适应的短时傅里叶变换方法。该方法基于最小谱质量评价指标准则,根据信号在不同时间段的特性自适应确定傅里叶变换的窗长。本发明突破了窗长固定的短时傅里叶变换的时频分辨率之间的约束,具有良好的时频分析性能。
2、技术方案:为实现上述目的,本发明提出一种窗长自适应的短时傅里叶变换方法,该方法包括以下步骤:
3、步骤1,获取待分析信号,对信号进行镜像延拓;
4、步骤2,改变傅里叶变换的窗长,以各时间点为中心对信号做傅里叶变换得到频谱;
5、步骤3,根据谱质量评价指标选取各时间点的初始窗长;
6、步骤4,在初始窗长的基础上进一步细化窗长间隔,改变傅里叶变换的窗长,以各时间点为中心对信号做傅里叶变换得到频谱;
7、步骤5,根据谱质量评价指标选取各时间点的最终窗长;
8、步骤6,将各时间点的最终窗长对应的傅里叶变换结果进行插值并拼接,以此得到时频图。
9、进一步的,在步骤1中,获取待分析信号x(n),n=1,2,…,n,n为信号采样时刻序号,n为信号采样点数,按照如下公式对信号进行镜像延拓:
10、
11、其中,s(i)为镜像延拓后的信号,i=1,2,…,3n为镜像延拓后的信号采样时刻序号。
12、进一步的,在步骤2中,以2的整数次幂为傅里叶变换的窗长,即l=2m+2,其中,m=1,2,…,m,为窗长序号,m为满足的最大整数,以k为时间间隔,分别以信号s(i)的第n+jk个时间点为中心做傅里叶变换得到频谱,其中,k为小于n的正整数,j=1,2,…,j,为时间点序号,j为不大于n/k的最大整数。
13、进一步的,在步骤3中,根据谱质量评价指标选取各时间点的初始窗长,包括以下步骤:
14、(3.1)谱质量评价指标的计算方法如下:
15、
16、其中,p(j,m)、as(j,m)、am(j,m)、v(j,m)和f(j,m)分别为信号s(i)以第n+jk个时间点为中心、窗长为2m+2时的傅里叶变换的频谱的谱质量评价指标、旁瓣幅度、主瓣幅度、主瓣瞬时带宽和瞬时频率,并规定v(j,m)为频谱归一化幅度的处的带宽;
17、(3.2)按照最小谱质量评价指标的准则选取各时间点的初始窗长为:
18、linitial(j)=2^[2+minitial(j)]
19、其中,
20、
21、上式表示谱质量评价指标p(j,m)取最小值时的窗长序号m。
22、进一步的,在步骤4中,以初始窗长为基础进一步细化窗长间隔,改变窗长进行傅里叶变换,包括以下步骤:
23、(4.1)将各时间点的初始窗长向下向上各延伸一个档作为窗长区间,即
24、
25、(4.2)将窗长区间等间隔取r个点,得到第j个时间点的窗长序列lj={lj(1),lj(2),…,lj(r)},其中,
26、
27、其中,r=1,2,…,r为窗长序号,r为大于2的正整数,
28、(4.3)以信号s(i)的第n+jk个时间点为中心做傅里叶变换,窗长遍历序列lj,将傅里叶变换结果存储为fj(r,f),其表示窗长为lj(r)时的频谱,其中,f为频率。
29、进一步的,在步骤5中,根据谱质量评价指标选取各时间点的最终窗长,包括以下步骤:
30、(5.1)谱质量评价指标的计算方法如下:
31、
32、其中,p(j,r)、as(j,r)、am(j,r)、v(j,r)和f(j,r)分别为信号s(i)以第n+jk个时间点为中心、窗长为lj(r)时的傅里叶变换的频谱的谱质量评价指标、旁瓣幅度、主瓣幅度、主瓣瞬时带宽和瞬时频率,并规定v(j,r)为频谱归一化幅度的处的带宽;
33、(5.2)按照最小谱质量评价指标的准则选取各时间点的最终窗长为:
34、
35、其中,
36、
37、上式表示谱质量评价指标p(j,r)取最小值时的窗长序号r。
38、进一步的,在步骤6中,将各时间点的最终窗长对应的傅里叶变换结果进行插值并拼接得到时频图,包括以下步骤:
39、(6.1)提取各时间点的最终窗长对应的傅里叶变换结果为:
40、sj(f)=fj(rfinal(j),f)
41、(6.2)将上述结果进行线性插值,插值为2m+1个点,得到sj(k),k=1,2,…,2m+1;
42、(6.3)将插值后的结果拼接为:
43、
44、以此得到时频图。
45、有益效果:与现有技术相比,本发明的技术方案具有以下有益技术效果:
46、在传统的短时傅里叶变换中,一般以固定窗长对信号进行时频处理,固定的时频分辨率难以适应信号在不同时间段的特性。本方法以各时间点为中心,通过粗略和细化窗长间隔的结合,两次基于最小谱质量评价指标准则选取合适的傅里叶变换窗长,减小计算量并提升算法性能,能够适应信号的时变特性,具有良好的时频分析效果。
1.一种窗长自适应的短时傅里叶变换方法,其特征在于,该方法包括以下步骤:
2.根据权利要求1所述的一种窗长自适应的短时傅里叶变换方法,其特征在于,在步骤1中,获取待分析信号x(n),n=1,2,…,n,n为信号采样时刻序号,n为信号采样点数,按照如下公式对信号进行镜像延拓:
3.根据权利要求2所述的一种窗长自适应的短时傅里叶变换方法,其特征在于,在步骤2中,以2的整数次幂为傅里叶变换的窗长,即l=2m+2,其中,m=1,2,…,m,为窗长序号,m为满足的最大整数,以k为时间间隔,分别以信号s(i)的第n+jk个时间点为中心做傅里叶变换得到频谱,其中,k为小于n的正整数,j=1,2,…,j,为时间点序号,j为不大于n/k的最大整数。
4.根据权利要求3所述的一种窗长自适应的短时傅里叶变换方法,其特征在于,在步骤3中,根据谱质量评价指标选取各时间点的初始窗长,包括以下步骤:
5.根据权利要求4所述的一种窗长自适应的短时傅里叶变换方法,其特征在于,在步骤4中,以初始窗长为基础进一步细化窗长间隔,改变窗长进行傅里叶变换,包括以下步骤:
6.根据权利要求5所述的一种窗长自适应的短时傅里叶变换方法,其特征在于,在步骤5中,根据谱质量评价指标选取各时间点的最终窗长,包括以下步骤:
7.根据权利要求6所述的一种窗长自适应的短时傅里叶变换方法,其特征在于,在步骤6中,将各时间点的最终窗长对应的傅里叶变换结果进行插值并拼接得到时频图,包括以下步骤: