主要内容

删除时间局域频率分量

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

T = 0:1/2000:1-1/2000;Dt = 1/2000;X1 = sin(50*pi*t).*exp(-50*pi*(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).²);X4 = 2*sin(140*pi*t).*exp(-80*pi*(t-0.8).^2);X = x1+x2+x3+x4;情节(t, x)网格;标题(叠加信号的

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

获取并显示CWT。

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

图中包含一个axes对象。使用默认莫尔斯小波的解析CWT轴对象包含图像、直线、区域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)

图中包含一个axes对象。标题为“幅值尺度图”的坐标轴对象包含图像、直线、区域3个类型的对象。

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

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

图中包含2个轴对象。标题为Original Signal的Axes对象1包含一个类型为line的对象。axis对象2(标题为Signal,第一个25hz组件被删除)包含一个类型为line的对象。

将重构信号与没有以0.2秒为中心的25赫兹分量的原始信号进行比较。

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

图中包含一个axes对象。坐标轴对象包含两个line类型的对象。这些对象表示逆CWT近似,原始信号没有25赫兹。