虫类函数
4视图(30天)
显示旧的评论
我用cwt用morl的小波函数。这是我用的语法:
[cwtwavex, freqsx] = cwt (x,尺度,‘morl’,‘absglb’, dt)
但似乎有一个错误的函数,因为当我选择一个脉冲函数作为输入(x)和情节中的输出不同的尺度,这应该是基本小波的形状,但是最后的输出主要是morlet函数调制与另一个函数!…这是我写的代码:清除所有;
xlen = 1000;Fs = 2 e4;%采样频率
dt = 1 / f;a0 = 2 ^ (1/32);尺度= 2 * a0。^ (0:6 * 32);
x = 0 (1, xlen + 1);
x (int16 (xlen / 2)) = 1;
wavex = x;[cwtwavex, freqsx] = cwt (x,尺度,‘morl’,‘absglb’, dt);
abs_cwt = abs (cwtwavex);
图(1)% %;情节(x);
图(2);情节(cwtwavex (1:),“r”) ylabel (“cwtx1”);
图(3);情节(cwtwavex (193:), b) ylabel (“cwtx193”);
图(4);网格(abs_cwt);
但如你所见,小波的形状输出最后规模(193)和网状图也很奇怪(许多波动)将源于错误的形状morlet小波的年代(193)(规模)! !
这个问题不存在使用“sym2”!请帮助。
还,似乎有很多限制使用morlet时小波变换函数:不需要“morl”容易与普通语法。请建议。
非常感谢。
最好的
0评论
答案(2)
韦恩王
2018年6月14日
嗨洗巴,首先,我将鼓励你使用新的类接口而不是遗产接口使用。新类接口,您可以使用分析Morlet“爱”。
属于你的当前的问题,我不知道你在期望什么。让我们用δ函数当我们获得类,我们基本上得到各尺度小波(逆转)。
Fs = 2 e4;%采样频率
dt = 1 / f;
a0 = 2 ^ (1/32);
尺度= 2 * a0。^ (0:6 * 32);
x = 0 (e3, 1);
x (500) = 1;
[cwtwavex, freqsx] = cwt (x,尺度,“morl”,dt);
让我们知道检查的频率在193元素的向量
freqsx (193)
这是126赫兹。现在得到相对应的小波尺度。
ψ= cwtwavex(193年);
让我们得到功率谱,我们希望峰值频率是127左右。
[Pxx F] = periodgram (psi, [], 2048, 2 e4);
[~,idx] = max (Pxx);
F (idx)
如果我们把一个小得多的规模,我们预计,小波中心频率更高。