一种改进S变换的有限窗长时频分析方法技术

本发明专利技术公开了一种改进S变换的有限窗长时频分析方法，在频率与窗长之间建立基于反正切函数的映射，使得窗长在一定范围内随频率自适应变化，将改进后的窗函数代入S变换，实现高分辨率时频分析，包括以下步骤：信号采样，得到离散序列后进行快速傅里叶变换获得信号频谱，再根据信号频谱特点以及分辨率要求确定控制因子，将得到的控制函数代入窗函数中，窗函数经过快速傅里叶变换后与扩维后的信号频谱相乘，再经傅里叶逆变换得到时频谱。本方法既实现了窗长随信号频率变化而变化，又限定了窗长的变化范围，保证了时频谱各个部分都有很高的分辨率。

【技术实现步骤摘要】
一种改进S变换的有限窗长时频分析方法
本专利技术属于一种雷达信号处理技术，特别是一种改进S变换的有限窗长高精度时频分析方法。
技术介绍
非平稳信号是雷达信号处理中最常见的信号，而时频分析是分析该类信号的重要工具。为了准确分析信号的局部特性，时频分析将一维时域信号映射到二维时频平面，从而获得信号的时频分布。目前，常用的时频分析方法主要包括短时傅里叶变换(STFT)、希尔伯特-黄变换(HHT)、S变换(ST)、小波变换(WT)等。DennisGabor于1946年提出了短时傅里叶变换，其基本思想是通过加窗实现信号的分段傅里叶变换，从而得到信号的时变特性。但STFT所用窗函数固定，与时间和频率无关，是一种单一分辨率的分析方法。1998年，黄锷等人提出的希尔伯特-黄变换，将经验模态分解与Hilbert谱分析相结合，将信号分解为若干固有模态函数(IMF)后映射到时频域中，但其理论与算法还未完善，存在模态混叠、端点效应等问题。而小波变换的思想来源于伸缩与平移方法，是一种窗口面积固定但形状可改变的时频局部化分析方法，能根据高低频信号特点自适应调整时频窗，有着“数学显微镜”美称。但小波基设计难度较大，还有容许性条件的约束，同时存在时频分辨率不足、尺度频率转换复杂等缺陷。为了弥补短时傅里叶变换和小波变换的不足，Stockwell提出了S变换，引入了可变高斯窗函数，且时窗宽度与频率导数成反比。该方法得到的时频谱在低频部分频率分辨率高，在高频部分频率分辨率低，即分辨率可变。但是这种反比关系使得窗函数在局部出现窗长过宽和过窄的问题，导致低频处时间定位失效，高频处频率定位失效。专利申请号为CN201611158226.X，专利技术名称为“一种基于非线性模式分解和自适应最优核的时频分析方法”的中国专利，首先运用非线性模式分解算法，将多分量非平稳信号分解为一组有物理意义的非线性模式分量，再利用自适应最优核的时频分析方法使核函数自适应地随着信号的变化而变化。该方法可以改善时频聚焦能力，局限性在于自适应最优核得到比较困难。专利申请号为CN201610946585.5，专利技术名称为“一种基于反褶积广义S变换的地震频谱成像方法”的中国专利，通过将原始信号与高斯窗各自的魏格纳分布进行二维褶积得到时频谱。该方法可以压制Wigner-Ville分布的交叉项的产生，同时使广义S变换谱获得了较高的时频聚集性，但局限性在于低频和高频处的分辨率不足问题得不到解决。由上可知，现有的时频分析方法还存在不足，需进一步改进来实现高精度的时频分析。
技术实现思路
本专利技术所解决的技术问题在于提供一种改进S变换的有限窗长时频分析方法。实现本专利技术目的的技术解决方案为：一种改进S变换的有限窗长时频分析方法，包括以下步骤：步骤1、对信号进行采样，采样点数为N＝t/T,其中t为信号时长，T为采样周期，得到信号的离散序列h[kT](k＝1,2,…,N)；步骤2、对信号的离散序列h[kT]进行快速傅里叶变换，得到离散信号频谱其中(n＝0,1,2,…,N-1)；步骤3、确定控制因子a、b、c的值，得到窗长控制函数后确定窗函数；步骤4、对窗函数进行快速傅里叶变化得到窗函数频谱；步骤5、将扩维后的信号频谱跟窗函数频谱相乘，再对其进行傅里叶逆变换；步骤6、重复步骤4、步骤5，直至所有频率点全部计算完成，得到高分辨率时频谱。本专利技术与现有技术相比，其显著优点为：1)窗函数可以根据频率自适应做出调整，在低频处有较高的频率分辨率，在高频处有较高的时间分辨率；2)通过反正切函数限制了窗长的变化范围，解决了在低频处由于时窗宽度过宽引起的时间定位不准确问题，以及在高频处由于频窗宽度过宽引起的频率定位不准确问题；3)对于不同类型的信号和不同的分辨率要求，通过调整控制因子取值，都能实现高分辨率时频分析，具有很强的灵活性。下面结合附图对本专利技术作进一步详细描述。附图说明图1是本专利技术一种改进S变换的有限窗长时频分析方法流程图。图2是本专利技术一种改进S变换的有限窗长时频分析方法的算法流程图。图3是本专利技术实施例1的时频分析结果图。图4是本专利技术实施例2的时频分析结果图。图5是本专利技术实施例3的时频分析结果图。具体实施方式结合图1和图2所示，本专利技术的一种改进S变换的有限窗长时频分析方法，包括以下步骤：步骤1、对信号进行采样，采样点数为N＝t/T,其中t为信号时长，T为采样周期，得到信号的离散序列h[kT](k＝1,2,…,N)；步骤2、对信号的离散序列h[kT]进行快速傅里叶变换，得到离散信号频谱其中(n＝0,1,2,…,N-1)；步骤3、确定控制因子a、b、c的值，得到窗长控制函数后确定窗函数，具体步骤为：步骤3-1、确定时窗取值范围[Δtmin,Δtmax]，其中Δtmin为最小时窗长度，Δtmax为最大时窗长度，通过下列不等式确定a和c的取值范围：步骤3-2、在取值范围内确定a和c的值，并取其中fs为采样频率，并将a、b、c的值代入窗长控制函数中：步骤3-3、将窗长控制函数带入到窗函数中，得到改进后的窗函数表达式：步骤4、对窗函数进行快速傅里叶变化得到窗函数频谱，对窗函数进行快速傅里叶变化得到窗函数频谱时，是以为起始频率点，具体公式为：其中，n从0开始取值。步骤5、将扩维后的信号频谱跟窗函数频谱相乘，再对其进行傅里叶逆变换，具体步骤为：步骤5-1、将信号频谱扩维得到扩维后的信号频谱其中m＝0,1,2,…N-1；步骤5-2、将扩维后的信号频谱与窗函数频谱G(m,n)相乘；步骤5-3、对步骤5-2的结果进行傅里叶逆变换，得到步骤6、重复步骤4、步骤5，直至所有频率点全部计算完成，得到高分辨率时频谱。判断所有频率点是否全部计算完成,具体是判断n是否满足大于等于N，若满足，输出时频谱；若不满足，n自加1后重复步骤4、步骤5；全部完成后判断时频谱的分辨率是否满足设定的分辨率要求，若满足，输出时频分析结果；若不满足，调整控制因子的值，重复步骤3、步骤4、步骤5，直至满足设定的分辨率要求，输出高分辨率时频谱。本专利技术一方面实现了窗函数大小随着分析信号的频率进行自适应改变；另一方面，通过反正切函数控制时窗宽度在一定范围内随着频率的增大而减小，实现在时频谱的各个区域都有较高的分辨率。下面以三种信号为例来说明本专利技术的时频分析方法。实施例1仿真信号为两个正弦信号的叠加，信号频率分别为100Hz和400Hz，解析式为：h(t)＝sin(200πt)+sin(800πt)t∈[0,1]信号采样频率fs＝1024Hz,图3为采用反正切函数联合S变换的有限窗长时频方法得到的时频谱。该信号只有两个固定的频率分量，对于频率不随时间变化的信号，只需考虑频率分辨率，因此通过调整控制因子，取a＝1,b＝50,c＝30将频窗宽度控制在较小范围内，保证时频谱有较好的频率分辨率性能。由图3可知，该方法可以实现很好的频率分辨率。实施例2仿真信号为调频斜率为k＝400的线性调频信号，解析式为：信号采样频率fs＝1024Hz，图4为采用反正切函数联合S变换的有限窗长时频方法对线性调频信号的分析结果对于线性调频信号，频率变化较大，所以调整控制因子，取a＝5,b＝50,c＝50可以得到如图的分析结果，解决了原S变换高频处信号发散、能量聚集性差的问题，具有很好的时频性能。实本文档来自技高网...

【技术保护点】
一种改进S变换的有限窗长时频分析方法，其特征在于，包括以下步骤：步骤1、对信号进行采样，采样点数为N＝t/T,其中t为信号时长，T为采样周期，得到信号的离散序列h[kT](k＝1,2,…,N)；步骤2、对信号的离散序列h[kT]进行快速傅里叶变换，得到离散信号频谱

【技术特征摘要】
1.一种改进S变换的有限窗长时频分析方法，其特征在于，包括以下步骤：步骤1、对信号进行采样，采样点数为N＝t/T,其中t为信号时长，T为采样周期，得到信号的离散序列h[kT](k＝1,2,…,N)；步骤2、对信号的离散序列h[kT]进行快速傅里叶变换，得到离散信号频谱其中(n＝0,1,2,…,N-1)；步骤3、确定控制因子a、b、c的值，得到窗长控制函数后确定窗函数；步骤4、对窗函数进行快速傅里叶变化得到窗函数频谱；步骤5、将扩维后的信号频谱跟窗函数频谱相乘，再对其进行傅里叶逆变换；步骤6、重复步骤4、步骤5，直至所有频率点全部计算完成，得到高分辨率时频谱。2.根据权利要求1所述的改进S变换的有限窗长时频分析方法，其特征在于，步骤3中确定控制因子a、b、c的值，得到窗长控制函数后确定窗函数，具体步骤为：步骤3-1、确定时窗取值范围[Δtmin,Δtmax]，其中Δtmin为最小时窗长度，Δtmax为最大时窗长度，通过下列不等式确定a和c的取值范围：步骤3-2、在取值范围内确定a和c的值，并取其中fs为采样频率，并将a、b、c的值代入窗长控制函数中：步骤...

【专利技术属性】
技术研发人员：芮义斌，严丽萍，谢仁宏，李鹏，郭山红，杜禹，吕云涛，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32