创建一个由指数加权正弦波组成的信号。该信号有两个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)网格在;标题(叠加信号的)
获取并显示CWT。
类(x, 2000);标题(“使用默认莫尔斯小波分析CWT”);
通过消零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)
绘制原始信号和重建信号。
次要情节(2,1,1);情节(t, x);网格在;标题(原始信号的);次要情节(2,1,2);情节(t, xrec)网格在;标题(“去掉第一个25赫兹分量的信号”);
将重构信号与没有以0.2秒为中心的25赫兹分量的原始信号进行比较。
Y = x2+x3+x4;图;xrec情节(t)在情节(t y“r——”网格)在;传奇(“逆CWT近似”,“没有25赫兹的原始信号”);持有从