小波分解synchrosqueezed变换

6视图(30天)
你好,
我想分解信号近似一个细节。像一维小波变换摧毁它。 http://de.mathworks.com/help/wavelet/ug/ex1d2.png
这也可能与小波synchrosqueezed变换?或者是墓场”只有“我可以得到时频图吗? http://de.mathworks.com/help/examples/wavelet_product/synchrosqueezingExample_02.png
谢谢你!

接受的答案

韦恩王
韦恩王 2016年4月15日
你好马克,synchrosqueezing和连续小波变换不能提供完美重建密切相关的输入信号一般像DWT。
你可能想看看MODWPT和/或MODWT。
MODWT将给你一个倍频带分解的输入信号(为逐步精细的八度乐队)与节能和完美重建没有时间分辨率的损失。
例如
负载wecg;
wt = modwt (wecg);
现在,如果你看看wt的能量矩阵的行,和那些你会回到原始信号的能量
总和(var (wt, 1 2))
%比较
var (wecg)
现在,如果你想要逐步建立你的信号从各级“细节”,使用MODWTMRA
mra = modwtmra (wt);
如果你这mra长列总和你最终将重建原始信号。
马克斯(abs (wecg-sum (mra) '))
您可以查看各级的贡献是什么样子:
次要情节(1,1);
情节(wecg);
kk = 6
次要情节(6 1 kk);
情节(mra (kk-1:));
结束
更看到MODWPT MODWPTDETAILS高精度时频分析。

答案(1)

韦恩王
韦恩王 2016年4月15日
你好马克,synchrosqueezed变换不专门提供“细节”像离散小波变换,但它提供的能力做一些非常相似。毕竟,在离散小波变换细节相当于带通项的输入。与逆synchrosqueezed transfom可以重建frequency-localized信号近似。例如
负载multicompsig;
sig = sig1 + sig2;
海温,[F] =墓场(sig, sampfreq);
轮廓(t、F、abs (sst));
包含(“时间”);ylabel (“赫兹”);
网格;
标题(“双组分信号Synchrosqueezed变换”);
现在让我们来重建一个近似,而低频正弦波和比较原始sig1孤立。
xrec = iwsst (sst F[250] 25日);
情节(sig1);持有;
情节(xrec“r”);网格;
同时,看看wsstridge,你可以识别和重建在时频脊。
(冰箱,iridge) = wsstridge (F海温,10日,“NumRidges”2);
持有;
情节(t,冰箱,“k”,“线宽”2);
你可以通过iridge向量或矩阵(矩阵,因为有两个山脊)IWSST简单地转化synchrosqueezed变换沿着山脊。
1评论
马克·休伯
马克·休伯 2016年4月15日
编辑:马克·休伯 2016年4月18日
谢谢你!我不确定我理解是对的。xrec细节和冰箱近似,对吧?但是我想做一个多级一维小波分解(sry基因没有说)小波synchrosqueezed变换。
当我做一个多级一维小波分解我收到以下:
信号s = an + d_1 +……+ d_n
所以,在你的情况中,an =冰箱和d_1 = xrec(不完全相同的(“=”是假的,如果我们是准确的)但这个想法)s =冰箱+ xrec。
但是我怎样才能得到更多更深层次像s = fridge_n + xrec_1 +……+ xrec_n
是一个多级一维小波变换与小波synchrosqueezed变换(这里一样的,不完全是,但只有主意多级一维小波变换相似)可能吗
非常感谢你

登录置评。

类别

找到更多的在人工智能的信号和图像帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!