CN115146221A - 基于多重同步压缩的mstct时频分析方法及系统 - Google Patents
基于多重同步压缩的mstct时频分析方法及系统 Download PDFInfo
- Publication number
- CN115146221A CN115146221A CN202210646074.7A CN202210646074A CN115146221A CN 115146221 A CN115146221 A CN 115146221A CN 202210646074 A CN202210646074 A CN 202210646074A CN 115146221 A CN115146221 A CN 115146221A
- Authority
- CN
- China
- Prior art keywords
- signal
- time
- frequency
- window function
- chirp
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了基于多重同步压缩的MSTCT时频分析方法及系统,该方法包括:通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数;根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信号的chirplet变换和频率重分配算子;根据信号的chirplet变换获取信号的时频表示结果;基于频率重分配算子,对信号的时频表示结果进行多重同步压缩处理,直至压缩结果满足预设条件,输出信号的时频图。通过使用本发明,能够在不提高计算成本的前提条件下实现信号的时频分析结果具有较高的分辨率和能量聚集度。本发明作为基于多重同步压缩的MSTCT时频分析方法及系统,可广泛应用于信号处理技术领域。
Description
技术领域
本发明涉及信号处理技术领域,尤其涉及基于多重同步压缩的MSTCT时频分析方法及 系统。
背景技术
时频分析(TFA)方法为刻画实际应用中常见的非平稳信号的复杂时变特征提供了一种 强有力的方法,因此,研究分析非平稳多分量信号的新方法具有重要的现实意义,传统的TF 分析方法包括线性TF和二次TF,线性TF分析方法的分辨率往往不能达到令人满意的水平, 而二次TF分析方法在分析多分量信号时,其非线性积导致了交叉项,降低了TF表示的可读 性,为了进一步提高同步压缩变换的性能,人们提出了许多集中TF方法,例如多重同步压 缩变换(MSST)方法,但是使用MSST对强调频信号进行时频分析时,容易出现在迭代次 数内无法正确逼近信号真实IF的情况,现有的还有利用信号chirp率(CR)和同步提取多项 式chirplet变换的TFA方法等利用高阶信息来使调制信号的TF结果的能量更集中,然而,这 种高阶的同步压缩变换方法的计算代价较大,且对噪声鲁棒性较差,不能进行精确地重构。
发明内容
为了解决上述技术问题,本发明的目的是提供基于多重同步压缩的chirplet变换(MSTCT) 时频分析方法及系统,能够在不提高计算成本的前提条件下实现信号的时频分析结果具有较 高的分辨率和能量聚集度。
本发明所采用的第一技术方案是:基于多重同步压缩的MSTCT时频分析方法,包括以 下步骤:
通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数;
根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信号的chirplet变 换和频率重分配算子;
根据信号的chirplet变换获取信号的时频表示结果;
基于频率重分配算子,对信号的时频表示结果进行多重同步压缩处理,直至压缩结果满 足预设条件,输出信号的时频图。
进一步,所述通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数这一步 骤,其具体包括:
获取待分析的信号、斜坡函数和窗函数;
根据斜坡函数对窗函数进行预处理,得到预处理后的窗函数;
根据预处理后的窗函数对待分析的信号进行短时傅里叶变换与估计处理,得到信号的估 计瞬时chirp率系数。
进一步,所述根据斜坡函数对窗函数进行预处理,得到预处理后的窗函数这一步骤,其 具体包括:
计算窗函数对时间的一阶偏导数和二阶偏导数;
根据斜坡函数分别与窗函数和窗函数对时间的一阶偏导数进行相乘运算,得到斜坡函数 与窗函数的乘积和斜坡函数与窗函数对时间的一阶偏导数的乘积;
结合窗函数对时间的一阶偏导数、窗函数度对时间的二阶偏导数、斜坡函数与窗函数的 乘积和斜坡函数与窗函数对时间的一阶偏导数的乘积,构建预处理后的窗函数。
进一步,所述根据预处理后的窗函数对待分析的信号进行短时傅里叶变换的表达式如下 所示:
进一步,所述信号的估计瞬时chirp率系数的计算公式如下所示:
上式中,表示信号的估计瞬时chirp率系数,表示使用窗函数得到的短时傅里叶变换结果,表示使用窗函数得 到的短时傅里叶变换结果,表示使用窗函数g(t)得到的短时傅里叶变换结果, 表示使用窗函数得到的短时傅里叶变换结果,表示使 用窗函数Tg(t,ω)=tg(t)得到的短时傅里叶变换结果,表示对(·)取实部,表示对(·)取 虚部。
进一步,所述根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信号 的chirplet变换和频率重分配算子这一步骤,其具体包括:
对短时傅里叶变换后的信号的能量进行计算,得到短时傅里叶变换后信号的能量;
对信号的估计瞬时chirp率系数和短时傅里叶变换后信号的能量进行估计分析,得到信号 的chirp率;
对信号的chirp率进行估计分析,得到信号的chirplet变换;
根据信号的chirplet变换计算信号的频率重分配算子。
进一步,所述信号的chirplet变换的计算公式如下所示:
进一步,所述频率重分配算子的计算表达式如下所示:
本发明所采用的第二技术方案是:基于多重同步压缩的MSTCT时频分析系统,包括:
估计模块,用于通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数;
分析模块,用于根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信 号的chirplet变换和频率重分配算子;
获取模块,用于根据信号的chirplet变换获取信号的时频表示结果;
输出模块,基于频率重分配算子,对信号的时频表示结果进行多重同步压缩处理,直至 压缩结果满足预设条件,输出信号的时频图。
本发明方法及系统的有益效果是:本发明通过窗函数对信号进行预估处理,并根据预估 结果进行信号的chirplet变换,可以提高信号同步压缩的变换性能,通过MSTCT对信号的 chirplet变换进行自调整,能够在同等的迭代计算成本的条件下比现有的变换技术有着更好的 时频分辨率,且输出的信号时频图对噪声的干扰有更好的鲁棒性,具有更高的分辨率。
附图说明
图1是本发明基于多重同步压缩的MSTCT时频分析方法的步骤流程图;
图2是本发明基于多重同步压缩的MSTCT时频分析系统的结构框图;
图3是本发明获取的信号的时域波形图;
图4是通过同步压缩变换方法对信号进行时频分析后得到的时频图;
图5是应用本发明方法对信号进行时频分析后得到的时频图;
图6是通过同步压缩变换方法对含有噪声的信号进行时频分析后得到的时频图;
图7是应用本发明方法对含有噪声的信号进行时频分析后得到的时频图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的详细说明。对于以下实施例中的步骤编 号,其仅为了便于阐述说明而设置,对步骤之间的顺序不做任何限定,实施例中的各步骤的 执行顺序均可根据本领域技术人员的理解来进行适应性调整。
参照图1,本发明提供了基于多重同步压缩的MSTCT时频分析方法,该方法包括以下步 骤:
S1、通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数;
S11、获取待分析的信号、斜坡函数和窗函数;
S111、参照图3,获取待分析的信号,所述待分析的信号包含具有噪声的信号和不含噪 声的信号,从图中可以看出信号f(t)在时域上是一个信号,600秒之前的部分是一个单频信号, 600秒之后的部分是一个非线性的调频信号;
具体地,不含噪声的信号如下所示:
所述具有噪声的信号具体如下所示:
f0(t)=f(t)+n(t)
上式中,f(t)表示不含噪声的信号,f0(t)表示含有噪声的信号,n(t)表示高斯白噪声,t表 示时频分析中的时间变量,其中采样频率为1Hz;
S12、根据斜坡函数对窗函数进行预处理,得到预处理后的窗函数;
S121、预设窗函数与斜坡函数,计算窗函数对时间的一阶偏导数和二阶偏导数;
具体地,输入窗函数,其具体如下所示:
上式中,g(t)表示窗函数,σ2表示方差;
其中σ=0.32,迭代次数M=10;
所述窗函数对时间的一阶偏导数的计算过程如下所示:
上式中,Dg(t)表示窗函数对时间的一阶偏导数;
所述窗函数对时间的二阶偏导数的计算过程如下所示:
上式中,D2g(t)表示窗函数对时间的二阶偏导数;
S122、根据斜坡函数分别与窗函数和窗函数对时间的一阶偏导数进行相乘运算,得到斜 坡函数与窗函数的乘积和斜坡函数与窗函数对时间的一阶偏导数的乘积;
具体地,所述斜坡函数与窗函数的乘积计算过程如下所示:
Tg(t)=tg(t)
上式中,Tg(t)表示斜坡函数与窗函数的乘积;
所述斜坡函数与窗函数对时间的一阶偏导数的乘积的计算过程如下所示:
上式中,DTg(t)表示斜坡函数与窗函数对时间的一阶偏导数的乘积。
S123、结合窗函数对时间的一阶偏导数、窗函数度对时间的二阶偏导数、斜坡函数与窗 函数的乘积和斜坡函数与窗函数对时间的一阶偏导数的乘积,构建预处理后的窗函数。
S13、根据预处理后的窗函数对待分析的信号进行短时傅里叶变换与估计处理,得到信号 的估计瞬时chirp率系数。
具体地,根据预处理后的窗函数对待分析的信号进行短时傅里叶变换,所述预处理后的 窗函数包含原有的窗函数一共具有五种形式,分别对待分析的信号进行短时傅里叶变换,其 具体地变换过程如下所示:
根据傅里叶的变换结果对信号进行估计处理,所述估计处理的具体过程如下所示:
上式中,表示信号的估计瞬时chirp率系数,表示使用窗函数得到的短时傅里叶变换结果,表示使用窗函数得 到的短时傅里叶变换结果,表示使用窗函数g(t)得到的短时傅里叶变换结果, 表示使用窗函数得到的短时傅里叶变换结果,表示使 用窗函数Tg(t,ω)=tg(t)得到的短时傅里叶变换结果,表示对(·)取实部,表示对(·)取 虚部。
S2、根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信号的chirplet 变换和频率重分配算子;
S21、对短时傅里叶变换后的信号的能量进行计算,得到短时傅里叶变换后信号的能量; 具体地,短时傅里叶变换后信号的能量的计算公式如下所示:
S22、对信号的估计瞬时chirp率系数和短时傅里叶变换后信号的能量进行估计分析,得 到信号的chirp率;
具体地,所述chirp率即信号的瞬时频率变化率,描述的是信号瞬时频率的变化趋势与变 化速度,根据信号的估计瞬时chirp率系数并结合短时傅里叶变换后信号的能量对信号的chirp 率进行再次估计,其具体如下所示:
S23、对信号的chirp率进行估计分析,得到信号的chirplet变换;
具体地,基于信号的chirp率估计信号的chirplet变换,其具体如下所示:
S24、根据信号的chirplet变换计算信号的频率重分配算子。
具体地,频率重分配算子的计算过程如下所示:
S3、基于频率重分配算子,对信号的chirplet变换进行多重同步压缩处理,输出信号的时 频图。
S31、根据信号的chirplet变换获取信号的时频表示结果;
S32、基于频率重分配算子,对信号的时频表示结果进行多重同步压缩处理,直至压缩结 果满足预设条件,输出信号的时频图。
具体地,所述多重同步压缩处理即为MSTCT变换,设置迭代次数是人为设置的,迭代 次数越多,信号压缩得越好,但是相应的计算量也更大;反之迭代次数少,则效果没有那么 好,计算量小,当迭代次数达到一定后,信号压缩水平趋于稳定,再增加迭代次数对效果提 升不大,依据频率重分配算子对chirplet变换得到的时频表示进行多重同步压缩,得到具有高 能量聚集度的时频表示结果,其中,MSTCT变换的具体过程如下所示:
上式中,Ts[M](t,ω)表示经过M次MSTCT后的结果,ω、a均表示频率;
其中,上述公式只对a进行积分,而不对ω进行积分。
最后,参照图5和图7,输出信号f(t)经过MSTCT后的时频图,时频图能够让我们同时 观察一个信号的时域信息和频域信息,能够有效刻画非平稳信号的频率成分随时间的变化规 律;时频分析方法可以应用于地震、天文、引力波、雷达、声纳、生物医学、语音和机械工 程等领域,通过时频分析得到的时频图可以对这些领域中的信号处理提供一个非常直观展示, 如图4和图6所示,为通过同步压缩变换方法所得到的信号的时频图,从图中可以看出,用 同步压缩变换方法进行时频分析后得到的时频图虽然在信号的单频部分有着较好的能量聚集 度,但是在信号的非线性调频部分的能量涂抹现象仍然未得到较好解决,且信号的频率变化 地越快,则得到的曲线分辨率越差,这也表明同步压缩变换方法在对被噪声污染的信号进行 时频分析时,噪声鲁棒性差。
参照图2,基于多重同步压缩的MSTCT时频分析系统,包括:
估计模块,用于通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数;
分析模块,用于根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信 号的chirplet变换和频率重分配算子;
获取模块,用于根据信号的chirplet变换获取信号的时频表示结果;
输出模块,基于频率重分配算子,对信号的时频表示结果进行多重同步压缩处理,直至 压缩结果满足预设条件,输出信号的时频图。
上述方法实施例中的内容均适用于本系统实施例中,本系统实施例所具体实现的功能与 上述方法实施例相同,并且达到的有益效果与上述方法实施例所达到的有益效果也相同。
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉 本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换,这些等 同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (9)
1.基于多重同步压缩的MSTCT时频分析方法,其特征在于,包括以下步骤:
通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数;
根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信号的chirplet变换和频率重分配算子;
根据信号的chirplet变换获取信号的时频表示结果;
基于频率重分配算子,对信号的时频表示结果进行多重同步压缩处理,直至压缩结果满足预设条件,输出信号的时频图。
2.根据权利要求1所述基于多重同步压缩的MSTCT时频分析方法,其特征在于,所述通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数这一步骤,其具体包括:
获取待分析的信号、斜坡函数和窗函数;
根据斜坡函数对窗函数进行预处理,得到预处理后的窗函数;
根据预处理后的窗函数对待分析的信号进行短时傅里叶变换与估计处理,得到信号的估计瞬时chirp率系数。
3.根据权利要求2所述基于多重同步压缩的MSTCT时频分析方法,其特征在于,所述根据斜坡函数对窗函数进行预处理,得到预处理后的窗函数这一步骤,其具体包括:
计算窗函数对时间的一阶偏导数和二阶偏导数;
根据斜坡函数分别与窗函数和窗函数对时间的一阶偏导数进行相乘运算,得到斜坡函数与窗函数的乘积和斜坡函数与窗函数对时间的一阶偏导数的乘积;
结合窗函数对时间的一阶偏导数、窗函数对时间的二阶偏导数、斜坡函数与窗函数的乘积和斜坡函数与窗函数对时间的一阶偏导数的乘积,构建预处理后的窗函数。
6.根据权利要求5所述基于多重同步压缩的MSTCT时频分析方法,其特征在于,所述根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信号的chirplet变换和频率重分配算子这一步骤,其具体包括:
对短时傅里叶变换后的信号的能量进行计算,得到短时傅里叶变换后信号的能量;
对信号的估计瞬时chirp率系数和短时傅里叶变换后信号的能量进行估计分析,得到信号的chirp率;
对信号的chirp率进行估计分析,得到信号的chirplet变换;
根据信号的chirplet变换计算信号的频率重分配算子。
9.基于多重同步压缩的MSTCT时频分析系统,其特征在于,包括以下模块:
估计模块,用于通过窗函数对信号进行估计处理,得到信号的估计瞬时chirp率系数;
分析模块,用于根据信号的估计瞬时chirp率系数,对信号的chirp率进行分析,得到信号的chirplet变换和频率重分配算子;
获取模块,用于根据信号的chirplet变换获取信号的时频表示结果;
输出模块,基于频率重分配算子,对信号的时频表示结果进行多重同步压缩处理,直至压缩结果满足预设条件,输出信号的时频图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210646074.7A CN115146221A (zh) | 2022-06-08 | 2022-06-08 | 基于多重同步压缩的mstct时频分析方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210646074.7A CN115146221A (zh) | 2022-06-08 | 2022-06-08 | 基于多重同步压缩的mstct时频分析方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115146221A true CN115146221A (zh) | 2022-10-04 |
Family
ID=83407600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210646074.7A Pending CN115146221A (zh) | 2022-06-08 | 2022-06-08 | 基于多重同步压缩的mstct时频分析方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115146221A (zh) |
-
2022
- 2022-06-08 CN CN202210646074.7A patent/CN115146221A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | A generalized synchrosqueezing transform for enhancing signal time–frequency representation | |
CN106597408B (zh) | 基于时频分析和瞬时频率曲线拟合的高阶pps信号参数估计方法 | |
CN107085140B (zh) | 基于改进的SmartDFT算法的非平衡系统频率估计方法 | |
KR101294681B1 (ko) | 기상 신호 처리장치 및 그 처리방법 | |
Belega et al. | Amplitude and phase estimation of real-valued sine wave via frequency-domain linear least-squares algorithms | |
Damaschke et al. | A fair review of non-parametric bias-free autocorrelation and spectral methods for randomly sampled data in laser Doppler velocimetry | |
Roy et al. | Raised cosine filter-based empirical mode decomposition | |
CN109100687B (zh) | 一种雷达设备lfm脉冲信号pslr参数确定方法 | |
CN112067130A (zh) | 一种多频率分量的幅值提取方法及其系统 | |
Emresoy et al. | Iterative instantaneous frequency estimation and adaptive matched spectrogram | |
Matania et al. | Algorithms for spectrum background estimation of non-stationary signals | |
CN115146221A (zh) | 基于多重同步压缩的mstct时频分析方法及系统 | |
Tan et al. | Identification of a Wiener–Hammerstein system using an incremental nonlinear optimisation technique | |
CN112130138A (zh) | 一种超视距雷达的电离层相位污染校正方法 | |
CN108801296B (zh) | 基于误差模型迭代补偿的传感器频响函数计算方法 | |
CN109460614B (zh) | 基于瞬时带宽的信号时间-频率分解方法 | |
CN112505640B (zh) | 基于参数自适应的扩展b分布脉冲信号时频分析方法 | |
Boudraa | Instantaneous frequency estimation of fm signals by ψb-energy operator | |
CN110579800B (zh) | 一种基于高精度同步挤压变换的地震数据数字处理方法 | |
JP4350336B2 (ja) | スペクトル・アナライザ用の分解能フィルタ | |
Faisal et al. | Suppression of false-terms in wigner-ville distribution using time and frequency windowing | |
Tyagi et al. | Extension of FBSE-EWT for Complex Signal Analysis | |
Wolf et al. | Amplitude and frequency estimator for aperiodic multi-frequency noisy vibration signals of a tram gearbox | |
CN115691537B (zh) | 一种耳机音频信号的分析与处理系统 | |
CN111723699B (zh) | 基于gst和rprg的多分量lfm信号时频线修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |