小波分解synchrosqueezed变换
6视图(30天)
显示旧的评论
你好,
这也可能与小波synchrosqueezed变换?或者是墓场”只有“我可以得到时频图吗?
http://de.mathworks.com/help/examples/wavelet_product/synchrosqueezingExample_02.png
谢谢你!
0评论
接受的答案
韦恩王
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高精度时频分析。
0评论
答案(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变换沿着山脊。