资 源 简 介
clear all;
omega=1000;sigma=20;M=3;c=-3;d=8;
Fs=500;T=1.024;
t=0:1/Fs:T-1/Fs;
L=length(t);
signal=sin(omega*t);
noise=zeros(1,L);
t1=T/2:1/Fs:T-1/Fs;
L1=length(t1);
w=0.1*randn(1,L1);
noise(L-L1+1:L)=w;
x=signal+noise;
figure(1);
subplot(411);plot(t,x);
X=fft(x,L);
for s=1:M
omega=c+d*s;
for n=1:L
psi(n)=exp(-(n-s*omega)^2/(4*s*omega)^2);
end
WTx(s,:)=((pi/s)^0.5/sigma)*ifft(X.*psi);
end
subplot(412);plot(t,abs(WTx(1,:))); %输出尺度为s=1
subplot(413);plot(t,abs(WTx(2,:))); %输出尺度为s=2
subplot(414);plot(t,abs(WTx(3,:))); %输出尺度为s=3
