ex.m
上传用户:speoil
上传日期:2022-06-23
资源大小:224k
文件大小:1k
源码类别:

波变换

开发平台:

Matlab

  1. x1=wavread('late_cafe');
  2. x2=x1(1:65536);
  3. s1=auread('late_cafe_out');
  4. s2=s1(1:65536);
  5. %fs=128;
  6. %t=(1:128)/128;
  7. %c1=sin(2*pi*t*8);
  8. %c2=sin(2*pi*t*10);
  9. %c3=sin(2*pi*t*60);
  10. %c4=c1+c2+c3;
  11. %plot(c4)
  12.    h = daubcqf(6); 
  13.    %[s,N] = makesig('HiSine',65536);
  14.    %n = randn(1,N);
  15.     %x = x2 + n/10;     % (approximately 10dB SNR)
  16.      % figure;plot(x2);hold on;plot(s2,'r');
  17.      %Denoise x with the default method based on the DWT
  18.       [xd,xn,opt1] = denoise(x2,h);
  19.      figure;plot(x2);hold on;plot(xd,'r');
  20.     %Denoise x using the undecimated (LSI) wavelet transform
  21.     [yd,yn,opt2] = denoise(x2,h,1,[]);
  22.     figure;plot(x2);hold on;plot(yd,'r');
  23.    fft1=fft(yd,65536);
  24.    fft2=fft(x2,65536); 
  25.    fft3=fft(s2,65536);
  26.    w=(0:32767)/32768*(20000/2);
  27.     figure;
  28.     plot(w,abs(fft1(1:32768)),'c');
  29.     hold on;
  30.     figure;
  31.     plot(w,abs(fft2(1:32768)),'r');
  32.     hold on;
  33.     plot(w,abs(fft3(1:32768)),'y');
  34.     %legend('denoise','noisy','orignal');
  35.     %wavwrite(yd,'clean')