离散余弦变换 小波变换 去噪

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