今天我们主要讲一下多重同步压缩变换multisynchrosqueezing transform(MSST)在信号处理中的应用.
原文:Yu G , Wang Z , Zhao P . Multisynchrosqueezing Transform[J]. Industrial Electronics, IEEE Transactions on, 2018.
多重同步压缩变换 (MSST)通过迭代重分配技术对信号的时频分量能量分布进行多次同步压缩变换SST以提高非稳态信号时频分布中特征信号分量的能量集中度。
先简要看一下SST原理
![](https://img-blog.csdnimg.cn/img_convert/28559a153de47c018d67fe3bbea2a430.jpeg)
MSST是在SST基础上通过多次迭代使得信号时频分量的能量在其瞬时频率附近得到进一步压缩,以获得能量高度集中的时频分布。多重同步压缩变换表达式为:
通过多次迭代,由式(4)表述的瞬时频率将会越来越接近信号的真实瞬时频率,从而获得强时变信号的高分辨率时频谱图。
首先,先导入一个实际的机械振动信号
load('vib_data1.mat')
fs = 12000; %采样频率
N = 1200; % 采样点数
time = (1:N)/fs; % 时间序列
fre = (fs/2)/(N/2):(fs/2)/(N/2):(fs/2); %频率序列
data=data(1:N);
看一下时域图和频谱图
![](https://img-blog.csdnimg.cn/img_convert/816b3a6f0b0a4d292174aaeed19ec642.jpeg)
看一下信号的STFT时频谱图
![](https://img-blog.csdnimg.cn/img_convert/42fb6aac815defa0467447e8ce548567.jpeg)
然后看一下标准的同步压缩变换,即SST时频谱图(基于短时傅里叶变换STFT)
![](https://img-blog.csdnimg.cn/img_convert/204378f296182f92bec9562c20b77296.jpeg)
竖直2阶同步压缩变换时频谱如下(以后会讲,包括水平同步压缩变换)
![](https://img-blog.csdnimg.cn/img_convert/4c7c1139a6eabb44b4d0edc39c5d48a8.jpeg)
竖直4阶同步压缩变换时频谱如下(以后会讲)
![](https://img-blog.csdnimg.cn/img_convert/583d27c25d3d5b677bc7660e5199b221.jpeg)
MSST时频谱图如下
![](https://img-blog.csdnimg.cn/img_convert/846299cc94f679d59dea9c19d98afc8c.jpeg)
顺便看下MSST提取的两条时频脊线
![](https://img-blog.csdnimg.cn/img_convert/bab2ab52c231f3b89178ec0998390d7c.jpeg)
两条时频脊线所对应的两个分解分量
看下连续小波变换CWT时频谱图
![](https://img-blog.csdnimg.cn/img_convert/b3fb61173eb6fe1649557f16af02faf4.jpeg)
同步压缩小波变换的时频谱,下次试试高阶同步压缩小波变换,效果应该会不错
![](https://img-blog.csdnimg.cn/img_convert/83cf762496d7c669c8fcec5727a0c2eb.jpeg)
然后再看一组转子振动数据,其时域波形和频谱如下
![](https://img-blog.csdnimg.cn/img_convert/c9f0f40547a7ed2464d57a5e816a3b7c.jpeg)
其MSST时频谱如下
![](https://img-blog.csdnimg.cn/img_convert/f8420195fcf430b5120e3fcf0a482a42.jpeg)
根据上面的MSST时频谱图,可以得到主要的3个分量
![](https://img-blog.csdnimg.cn/img_convert/ba751f540f0a20c9227445d23beabddb.jpeg)
看一下连续小波变换的CWT时频谱图
![](https://img-blog.csdnimg.cn/img_convert/cde6e9a442b9062a72a1151b2c97cbae.jpeg)
其同步压缩小波变换时频谱图如下
![](https://img-blog.csdnimg.cn/img_convert/eb0a3bb2ce96e1d31d1c42f0f6ac8203.jpeg)
放大一下
![](https://img-blog.csdnimg.cn/img_convert/54d55b5beecead0e38d595fa58276dc9.jpeg)