移除时间局部频率分量

打开实时脚本

创建一个由指数加权正弦波组成的信号。信号有两个25赫兹的分量——一个以0.2秒为中心，一个以0.5秒为中心。它还有两个70hz的组件——一个以0.2秒为中心，一个以0.8秒为中心。前25-Hz和70-Hz分量同时出现。

T = 0:1/2000:1-1/2000;Dt = 1/2000;X1 = sin(50* *t).*exp(-50* *(t-0.2).^2);X2 = sin(50*pi*t).*exp(-100*pi*(t-0.5).^2);X3 = 2*cos(140*pi*t).*exp(-50*pi*(t-0.2).^2);X4 = 2*sin(140*pi*t).*exp(-80*pi*(t-0.8).^2);X = x1+x2+x3+x4;情节(t, x)网格在；标题(叠加信号的）

图中包含一个轴对象。标题为“叠加信号”的axes对象包含一个类型为line的对象。

获取并显示CWT。

类(x, 2000);标题(用默认莫尔斯小波分析CWT）;

图中包含一个轴对象。标题为analytical CWT using Default Morse Wavelet的轴对象包含图像、线、区域类型的3个对象。

通过将CWT系数归零，去除发生在大约0.07秒到0.3秒之间的25hz分量。使用逆CWT (icwt)重建信号的近似值。

[cfs,f] = cwt(x,2000);T1 = .07;T2 = .33;F1 = 19;F2 = 34;cfs(f > F1 & f < F2, t> T1 & t < T2) = 0;Xrec = icwt(cfs);

显示重构信号的CWT。初始的25hz组件被移除。

类(xrec, 2000)

图中包含一个轴对象。标题为Magnitude scalalogram的轴对象包含图像、线、区域类型的3个对象。

绘制原始信号和重建信号。

次要情节(2,1,1);情节(t, x);网格在；标题(原始信号的）;次要情节(2,1,2);情节(t, xrec)网格在；标题(“去除第一个25赫兹分量的信号”）;

图中包含2个轴对象。标题为Original Signal的Axes对象1包含一个类型为line的对象。去掉第一个25-Hz组件的axis对象2包含一个line类型的对象。

将重建信号与不以0.2秒为中心的25-Hz分量的原始信号进行比较。

Y = x2+x3+x4;图;xrec情节(t)在情节(t y“r——”网格)在；传奇(逆CWT近似，“无25赫兹原始信号”）;持有从

图中包含一个轴对象。axis对象包含2个line类型的对象。这些对象表示逆CWT近似，原始信号没有25赫兹。