中文说明：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
English Description：