一种基于能量重心法的扫频卷积变换时频分析方法技术

本发明专利技术提供了一种基于能量重心法的扫频卷积变换时频分析方法，包括以下步骤：步骤一、对实验信号做希尔伯特变换得到信号的解析信号；步骤二、构建一个频率从0到二分之一采样频率范围内变化的复数正弦扫频序列；步骤三、对解析信号进行移动加窗操作，得到一系列的加窗信号段；步骤四、将加窗信号与复数正弦扫频序列进行卷积；步骤五、利用能量重心法对卷积序列能量谱进行校正，可以得到最大能量峰值位置在复数正弦扫频序列中的坐标；步骤六、通过能量峰值坐标计算出信号的瞬时频率；应用本技术方案可实现分析平稳信号和非平稳信号时，计算精度高，计算时间短。计算时间短。

【技术实现步骤摘要】
一种基于能量重心法的扫频卷积变换时频分析方法


[0001]本专利技术涉及扫频卷积变换时频分析方法
，特别是一种基于能量重心法的扫频卷积变换时频分析方法。

技术介绍

[0002]在现代科学研究和实际应用中，对信号的采集、描述和处理是不可缺少的重要组成部分。自然界和实际应用中存在大量非平稳信号，如地震波和导波信号、生物电信号、机械振动信号、雷达波信号和语音信号等。因此，为有效研究这些非平稳信号，时频分析得到了广泛关注。
[0003]目前，现有的技术中存在大量的时频分析方法。例如短时傅里叶变换(STFT,Short
‑
Time Fourier transform)、小波变换(WT,Wavelet transform)、线调频小波变换(LMWT,Linear modulated wavelet transform)、基于Wigner
‑
Ville的时频分析方法和自适应非参数时频分析法(Adaptive non
‑
parametric time
‑
frequency analysis)等方法。
[0004]其中短时傅立叶变换采用固定窗，故它的时频分辨率固定不变，受限于Heisenberg测不准的原理，它的时频分辨率不能同时达到最优；小波变换是典型的线性变换，其本质上是对信号进行加窗的线性变换，同样受限于Heisenberg测不准的原理，它的时频分辨率不能同时达到最优；线调频小波变换由于采用线性变换核，故不适合分析非线性调频信号；基于Wigner
‑
Ville的时频分析方法由于受到交叉项的干扰，在分析非平稳信号时不能真实反映信号的时频特征。

技术实现思路

[0005]有鉴于此，本专利技术的目的在于提供一种基于能量重心法的扫频卷积变换时频分析方法，利用扫频信号覆盖了所考虑信号的所有频率分量的特性，该方法是将信号与扫频信号进行卷积，并利用能量重心法对其进瞬时频率的计算。该方法适用于分析平稳信号和非平稳信号，计算精度高，计算时间短。
[0006]为实现上述目的，本专利技术采用如下技术方案：一种基于能量重心法的扫频卷积变换时频分析方法，包括以下步骤：
[0007]步骤一、对实验信号做希尔伯特变换得到信号的解析信号；
[0008]步骤二、根据采样信号的采样频率，构建一个频率从0到二分之一采样频率范围内变化的复数正弦扫频序列；
[0009]步骤三、对解析信号进行移动加窗操作，得到一系列的加窗信号段；
[0010]步骤四、将加窗信号与复数正弦扫频序列进行卷积，当加窗信号移动到与复数正弦扫频序列具有相同频率序列范围时，当加窗信号靠近同频序列时共振峰值逐渐增大至最高点，达到峰值后再逐渐衰减；
[0011]步骤五、利用能量重心法对卷积序列能量谱进行校正，可以得到最大能量峰值位置在复数正弦扫频序列中的坐标；
[0012]步骤六、通过能量峰值坐标计算出信号的瞬时频率；
[0013]根据能量重心法计算信号的瞬时频率，所述计算信号瞬时频率是根据以下方程组得到：
[0014][0015]其中，f为信号的瞬时频率，Y
i
＝y(x)2，y(x)为汉宁窗的主瓣函数，n＝0,1,2,
···
,N
‑
1；N为采样点数。
[0016]在一较佳的实施例中，所述对信号做希尔伯特变换包括步骤：
[0017]对实验信号s(t)做希尔伯特变换可得到其计算公式为：
[0018][0019]故实验信号s(t)的解析信号z(t)的计算公式为：
[0020][0021]其中j为一个虚数单位。
[0022]在一较佳的实施例中，所述复数正弦扫频序列包括步骤：
[0023]构建一个频率从0到二分之一采样频率范围内变化的复数正弦扫频序列其定义式为：
[0024][0025]其中t为时间，f0为初始频率，u0为调频率。
[0026]在一较佳的实施例中，所述的汉宁窗W(n)可定义为：
[0027][0028]其中，n＝0,1,2,
···
,N
‑
1；N为采样点数。
[0029]因此汉宁窗的主瓣函数为：
[0030][0031]在一较佳的实施例中，所述将加窗信号段与复数正弦扫频序列进行卷积的计算公式CCTFT(t,ω)为：
[0032][0033]其中，w(σ)为归一化对称汉宁窗函数，z(τ)为实验信号，s(t)的解析信号。
[0034]在一较佳的实施例中，所述卷积序列的能量通过对卷积序列绝对值的计算得到，能量计算公式E
CCTFT(t,ω)
为：
[0035][0036]在一较佳的实施例中，所述能量重心法的计算公式为：
[0037][0038]与现有技术相比，本专利技术具有以下有益效果：
[0039](1)本专利技术不需要进行傅立叶变换，不受傅里叶变换只能反映信号时域和频域的全局特征的限制，可以对非平稳信号的局部特征进行分析。
[0040](2)本专利技术实现了瞬时频率的计算，结合了能量重心法，适用于平稳信号和非平稳信号，计算精度高，计算时间短。
附图说明
[0041]图1为本专利技术优选实施例的一种基于能量重心法的扫频卷积变换时频分析方法的方法原理图。
具体实施方式
[0042]下面结合附图及实施例对本专利技术做进一步说明。
[0043]应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属
的普通技术人员通常理解的相同含义。
[0044]需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式；如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0045]一种基于能量重心法的扫频卷积变换时频分析方法，包括以下步骤：
[0046]步骤一、对实验信号做希尔伯特变换得到信号的解析信号；
[0047]步骤二、根据采样信号的采样频率，构建一个频率从0到二分之一采样频率范围内变化的复数正弦扫频序列；
[0048]步骤三、对解析信号进行移动加窗操作，得到一系列的加窗信号段；
[0049]步骤四、将加窗信号与复数正弦扫频序列进行卷积，当加窗信号移动到与复数正弦扫频序列具有相同频率序列范围时，当加窗信号靠近同频序列时共振峰值逐渐增大至最高点，达到峰值后再逐渐衰减；
[0050]步骤五、利用能量重心法对卷积序列能量谱进行校正，可以得到最大能量峰值位置在复数正弦扫频序列中的坐标；
[0051]步骤六、通过能量峰值坐标计算出信号的瞬时频率；
[0052]根据能量重心法计算信号的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于能量重心法的扫频卷积变换时频分析方法，其特征在于，包括以下步骤：步骤一、对实验信号做希尔伯特变换得到信号的解析信号；步骤二、根据采样信号的采样频率，构建一个频率从0到二分之一采样频率范围内变化的复数正弦扫频序列；步骤三、对解析信号进行移动加窗操作，得到一系列的加窗信号段；步骤四、将加窗信号与复数正弦扫频序列进行卷积，当加窗信号移动到与复数正弦扫频序列具有相同频率序列范围时，当加窗信号靠近同频序列时共振峰值逐渐增大至最高点，达到峰值后再逐渐衰减；步骤五、利用能量重心法对卷积序列能量谱进行校正，可以得到最大能量峰值位置在复数正弦扫频序列中的坐标；步骤六、通过能量峰值坐标计算出信号的瞬时频率；根据能量重心法计算信号的瞬时频率，所述计算信号瞬时频率是根据以下方程组得到：其中，f为信号的瞬时频率，Y
i
＝y(x)2，y(x)为汉宁窗的主瓣函数，n＝0,1,2,
···
,N
‑
1；N为采样点数。2.根据权利要求1所述的一种基于能量重心法的扫频卷积变换时频分析方法，其特征在于，所述对信号做希尔伯特变换包括步骤：对实验信号s(t)做希尔伯特变换可得到其计算公式为：故实验信号s(t)的解析信号z(...

【专利技术属性】
技术研发人员：钟剑锋，钟嘉杰，钟舜聪，刘东明，
申请(专利权)人：福州大学，
类型：发明
国别省市：