【改进的多同步挤压变换】基于改进多同步挤压的高分辨率时频分析工具,用于分析非平稳信号(Matlab代码实现)

简介: 【改进的多同步挤压变换】基于改进多同步挤压的高分辨率时频分析工具,用于分析非平稳信号(Matlab代码实现)

💥💥💞💞欢迎来到本博客❤️❤️💥💥


🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。


⛳️座右铭:行百里者,半于九十。


📋📋📋本文目录如下:🎁🎁🎁


目录


💥1 概述


📚2 运行结果


🎉3 参考文献


🌈4 Matlab代码、数据、文章


💥1 概述

文献来源:


844fcc6f23a24595bd2437b6f7fa053a.png


该文提出一种高分辨率时频(TF)分析方法,用于强非平稳信号的分析。通过传统方法生成的TF表示通常过于模糊,无法为此类信号提供精确的特征。最近提出的一种称为多同步挤压变换(MSST)的方法克服了传统方法中存在的大多数问题,这似乎是一种很有前途的工具。但是,MSST仍然存在一个主要问题,即非重分配点问题,这可能导致一些特殊TF点的能量模糊问题。本文主要关注解决这个问题。研究发现,MSST中的此类问题是由重分配步骤的离散过程中的舍入操作引起的。然后采用一种有效的方法来使用简单的策略来解决这个问题。此外,研究中还提供了离散实现。数值分析表明,所提方法能够有效提高与MSST相当的能量浓度。与其他先进方法的比较还表明,所提出的方法在处理强非平稳信号和噪声附加信号方面具有更好的性能。在实验信号分析中,我们进行了三个实验,以验证所提方法在真实世界信号分析中的有效性。


原文摘要:


In this paper, a high-resolution time-frequency (TF) analysis method is presented for the analysis of strongly non-stationary signals. TF representations generated by conventional methods are usually too blurry to provide precise features for such signals. A recently proposed method, called multisynchrosqueezing transform (MSST), overcomes most of the problems that exist in conventional methods, which seems to be a promising tool. However, the MSST still has a major problem, i.e., non-reassigned point problem, which may lead to the blurry energy problem for some special TF points. This paper mainly focuses on resolving this problem. This study finds that such a problem in the MSST is caused by the rounding operation in the discrete procedure of the reassigned step. An effective method is then employed to address this problem using a simple strategy. Additionally, discrete implementation is provided in the study. The numerical analysis shows that our proposed method can effectively improve the energy concentration comparable to the MSST. Comparisons with other advanced methods also show that the proposed method offers better performance in addressing strongly non-stationary signals and noise-added signals. In the experimental signal analysis, we carry out three experiments to validate the effectiveness of the proposed method in the analysis of real-world signals.


有轴承缺陷的旋转机器通常会产生脉冲信号,振动传感器[1],[2],[3],[4],[5],[6],[7]。然而,当机器以可变速度运行时,振动信号处理变得具有挑战性。这是因为在这种情况下测量的信号通常表现出强烈的非平稳特性,即信号的频率随时间而变化很大。非平稳信号分析在旋转机械故障诊断中越来越受到关注。时频(TF)分析(TFA)技术具有应对此类挑战的强大能力。回顾TFA技术在旋转机械故障诊断中的发展的几篇文章可以在[8],[9],[10]中找到。从最近的研究中可以知道,能量浓度是评估TFA方法性能的关键指标[11]。这是因为集中TF表示(TFR)具有更好的表征信号故障特征的能力。然而,受海森堡不确定性原理的限制,传统的TFA方法,例如短时傅里叶变换(STFT),小波变换(WT)和S变换,在处理强时变信号时难以提供集中的结果。近几十年来,旨在克服传统方法缺点的各种新开发的TFA技术引起了广泛关注,例如,重新分配方法(RM)[12],同步挤压变换(SST)[13,14],解调SST(DSST)[15,16],高阶SST [17],[18],[19],同步提取变换(SET)[20,21]和多SST(MSST)[22]。


RM技术旨在根据信号的瞬时频率(IF)和群延迟的局部估计,将数据映射到更接近集中区域的新坐标[12],从而锐化能量涂抹的TFR。这种能量涂抹的TFR通常是由STFT,WT或S变换预先产生的。但是,RM 的映射应用于频谱图或标度图,通常定义为 STFT 或 WT 的平方幅度。这样的映射会丢失信号的相位信息,这也意味着无法从RM结果重建信号。作为一种类似RM的后处理技术,SST仅根据IF的估计值执行TFR映射,从而保留了逆能力。这使得SST在许多领域更具优势,例如机械故障诊断[23,24],地震信号分析[25,26]和呼吸动力学分析[27]。最近的各种研究侧重于在表征强时变信号时进一步增强SST技术的能力[13],[14],[15],[16],[17],[18],[19],[20],[21],[22]。


在线性TFA算法框架下提出了SST技术,例如WT,STFT和S变换。然而,受线性TF原子缺点的限制,线性TFA方法不能很好地处理强时变信号。相应的TF结果经常受到能量模糊问题的困扰。SST也面临同样的问题。解调技术设计非线性TF原子来表征时变信号,可以有效克服线性TFA方法的问题。此外,基于解调技术的新型SST方法在提高能量浓度方面显示出广阔的潜力。这种技术通常称为DSST方法[15,16]。然而,解调技术必须根据信号的先验信息设计非线性TF原子[9]。在实践中,很难甚至不可能提前确定真实世界信号的基本信息。这阻碍了DSST技术的工程应用。


提出了高阶SST方法来处理高度调频(FM)信号,该方法不需要事先提供有关信号的任何信息[17],[18],[19]。理论上,SST假设分析的信号应该是纯谐波信号。这意味着 SST 仅适用于处理微弱的 FM 信号。为了改善这种情况,在更复杂的信号模型上建立了高阶SST的框架,例如线性FM信号和高阶多项式FM信号。在无噪声情况下,高阶SST可以为强时变信号提供高度集中的结果。然而,最近的研究发现,高阶SST方法对噪声非常敏感。使用这种方法通常很难获得具有高噪声的信号的满意结果[22]。


SET方法旨在仅保留与信号时变特征密切相关的TF系数[20,21]。此外,SET消除了大多数弱相关的TF系数。因此,SET结果比SST结果更集中。但是,SET仅提供信号的近似重建。重建的性能随着信号非线性度的增加而降低。


最近发表的一篇论文介绍了MSST方法,该方法采用迭代程序来提高SST的能量浓度。MSST允许对信号进行完美的重建,不需要先验信息。[22]中的研究表明,在处理无噪声和加噪声信号时,MSST可以提供比RM,DSST和高阶SST更集中的结果。MSST方法似乎是实现理想TFR(ITFR)的有前途的工具[28]。然而,[22]的讨论部分指出,MSST方法中存在一个主要问题,即阻碍TF特征集中表征的非重分配点问题。


📚2 运行结果


212b43eb86144b3489629d108a275cda.png

9dd4fd6d0d904965aa0e6fa08150377a.png


部分代码:

function [Ts] = IMSST_Z(x,hlength,num)
% Computes the IMSST (Ts)  of the signal x.
% INPUT
%    x      :  Signal needed to be column vector.
%    hlength:  The length of window function.
%    num    :  iteration number.
% OUTPUT
%    Ts     :  The SST
%    tfr     :  The STFT
[xrow,xcol] = size(x);
if (xcol~=1),
    error('X must be column vector');
end;
if (nargin == 2),
    num=10;
end
if (nargin == 1),
    num=10;
    hlength=xrow/8;
end
hlength=hlength+1-rem(hlength,2);
ht = linspace(-0.5,0.5,hlength);ht=ht';
% Gaussian window
h = exp(-pi/0.32^2*ht.^2);
[hrow,~]=size(h); Lh=(hrow-1)/2;
N=xrow;
t=1:xrow;
tfr= zeros (round(N/2),N) ;
omega = zeros (round(N/2),N-1);
omega2 = zeros (round(N/2),N);
Ts = zeros (round(N/2),N);
%Compute STFT
for icol=1:N,
    ti= t(icol); tau=-min([round(N/2)-1,Lh,ti-1]):min([round(N/2)-1,Lh,xrow-ti]);
    indices= rem(N+tau,N)+1;
    rSig = x(ti+tau,1);
    tfr(indices,icol)=rSig.*conj(h(Lh+1+tau));
end;
tfr=fft(tfr);
tfr=tfr(1:round(N/2),:);
%2D IF of the SST
for i=1:round(N/2)


🎉3 参考文献

部分理论来源于网络,如有侵权请联系删除。



🌈4 Matlab代码、数据、文章

相关文章
|
8月前
|
算法 数据安全/隐私保护 计算机视觉
基于二维CS-SCHT变换和LABS方法的水印嵌入和提取算法matlab仿真
该内容包括一个算法的运行展示和详细步骤,使用了MATLAB2022a。算法涉及水印嵌入和提取,利用LAB色彩空间可能用于隐藏水印。水印通过二维CS-SCHT变换、低频系数处理和特定解码策略来提取。代码段展示了水印置乱、图像处理(如噪声、旋转、剪切等攻击)以及水印的逆置乱和提取过程。最后,计算并保存了比特率,用于评估水印的稳健性。
|
8月前
|
算法
基于小波多普勒变换的回波信号检测matlab仿真
基于小波多普勒变换的回波信号检测matlab仿真
|
8月前
|
机器学习/深度学习 并行计算 算法
霍夫变换椭圆检测(matlab仿真与图像处理系列第2期)
霍夫变换椭圆检测(matlab仿真与图像处理系列第2期)
|
4月前
|
机器学习/深度学习 算法 数据安全/隐私保护
基于NSCT非采样轮廓波变换和CNN网络人脸识别matlab仿真
本项目展示了一种结合非采样轮廓波变换(NSCT)与卷积神经网络(CNN)的人脸识别系统。通过NSCT提取多尺度、多方向特征,并利用CNN的强大分类能力实现高效识别。项目包括ORL人脸库的训练结果对比,提供Matlab 2022a版本下的完整代码及详细中文注释,另有操作步骤视频指导。
|
7月前
|
算法 计算机视觉 异构计算
基于FPGA的图像一维FFT变换IFFT逆变换verilog实现,包含tb测试文件和MATLAB辅助验证
```markdown ## FPGA 仿真与 MATLAB 显示 - 图像处理的 FFT/IFFT FPGA 实现在 Vivado 2019.2 中仿真,结果通过 MATLAB 2022a 展示 - 核心代码片段:`Ddddddddddddddd` - 理论:FPGA 实现的一维 FFT/IFFT,加速数字信号处理,适用于高计算需求的图像应用,如压缩、滤波和识别 ```
|
7月前
|
机器学习/深度学习 算法
基于BP神经网络和小波变换特征提取的烟草香型分类算法matlab仿真,分为浓香型,清香型和中间香型
```markdown 探索烟草香型分类:使用Matlab2022a中的BP神经网络结合小波变换。小波分析揭示香气成分的局部特征,降低维度,PCA等用于特征选择。BP网络随后处理这些特征,以区分浓香、清香和中间香型。 ```
|
8月前
|
算法 数据安全/隐私保护 C++
基于二维CS-SCHT变换和扩频方法的彩色图像水印嵌入和提取算法matlab仿真
该内容是关于一个图像水印算法的描述。在MATLAB2022a中运行,算法包括水印的嵌入和提取。首先,RGB图像转换为YUV格式,然后水印通过特定规则嵌入到Y分量中,并经过Arnold置乱增强安全性。水印提取时,经过逆过程恢复,使用了二维CS-SCHT变换和噪声对比度(NC)计算来评估水印的鲁棒性。代码中展示了从RGB到YUV的转换、水印嵌入、JPEG压缩攻击模拟以及水印提取的步骤。
|
8月前
|
机器学习/深度学习 算法 数据安全/隐私保护
基于DCT变换和位平面分解的数字水印嵌入提取算法matlab仿真
这是一个关于数字水印算法的摘要:使用MATLAB2022a实现,结合DCT和位平面分解技术。算法先通过DCT变换将图像转至频域,随后利用位平面分解嵌入水印,确保在图像处理后仍能提取。核心程序包括水印嵌入和提取,以及性能分析部分,通过PSNR和NC指标评估水印在不同噪声条件下的鲁棒性。
|
8月前
|
算法 数据安全/隐私保护 计算机视觉
基于DCT变换的彩色图像双重水印嵌入和提取算法matlab仿真
**算法摘要:** - 图形展示:展示灰度与彩色图像水印应用,主辅水印嵌入。 - 软件环境:MATLAB 2022a。 - 算法原理:双重水印,转换至YCbCr/YIQ,仅影响亮度;图像分割为M×N块,DCT变换后嵌入水印。 - 流程概览:两步水印嵌入,每步对应不同图示表示。 - 核心代码未提供。
|
8月前
|
机器学习/深度学习 算法 计算机视觉
霍夫变换车道线识别-车牌字符识别代码(matlab仿真与图像处理系列第5期)
霍夫变换车道线识别-车牌字符识别代码(matlab仿真与图像处理系列第5期)