虫类函数

12个视图(30天)
洗巴Ebra
洗巴Ebra 2018年6月11日
回答: 洗巴Ebra2018年7月7日
我用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”容易与普通语法。请建议。
非常感谢。
最好的

答案(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)
如果我们把一个小得多的规模,我们预计,小波中心频率更高。
1评论
洗巴Ebra
洗巴Ebra 2018年7月4日
亲爱的韦恩
非常感谢你的回答。但我想阴谋在不同尺度的小波函数,以确保他们是正确的。当你的阴谋(psi)(在你的例子),它没有一个正确的形状扩展morlet ! !,因为它似乎是调制与另一个函数,我相信应该是错误的。请帮助。非常感谢。
最好洗巴

登录置评。


洗巴Ebra
洗巴Ebra 2018年7月7日
亲爱的韦恩,
我再次检查了cwt功能和我的数据。“类”Matlab函数绝对错误,当应用于数据,许多问题(主要是证明波动)发生造成致命影响输出! ....说我将使用morlet小波,6 * 32尺度(没有八度= 6)。你可以自己测试它。我花了大量的时间(2 - 3个月)前我终于意识到,太多的问题,我的输出不是从仪表或信号本身,但CWT MATLAB函数! ! ....我不得不最后自己写一个函数做连续小波分析! ! ....
希望mathwork团队纠正这种错误! !…
最好洗巴

标签

s manbetx 845

社区寻宝

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

开始狩猎!