mifig51.m

老外写的小波变换的工具箱

% mifig51 -- De-Noising of Gaussian White Noise
%
% Panel (a): Heaviside, MIPT
% Panel (b): Heaviside, AIPT
% Panel (c): Doppler, MIPT
% Panel (d): Doppler, AIPT
% 

scriptname = 'mifig51';

J = 9;
Len = 3^J;
domain = linspace(0,1,Len);

z1 = randn(1,Len); z2 = randn(1,Len);
doppler = 10*MakeSignal('Doppler',Len) + z1; % Doppler
s = 8100;
heav = 4 .* [zeros(1,s) ones(1,Len-s)] + z2; % Heaviside

clear z1 z2

subplot(2,3,1), plot(domain, doppler), title('Doppler+Gaussian White Noise')
subplot(2,3,4), plot(domain, heav), title('Heaviside+Gaussian White Noise')

subplot(2,3,2)
L = 3; D = 2;

%%%% Set Thresholds
ThresType = 'theoretical'
% ThresType = 'simulation'

if strcmp(ThresType,'theoretical')
  GaussT = MIPTThreshold('Gaussian',J,L,'theoretical');
elseif strcmp(ThresType,'simulation')
  % GaussT = MIPTThreshold('Gaussian',J,L,'simulation');
  % Takes too long to do it real time.
  load GaussT
end

wc = FMIPT(doppler,L,D);
for j=(L+1):J,
  wc(triad(j)) =  HardThresh(wc(triad(j)), GaussT(j));
end
doppler_mt = IMIPT(wc,L,D);
plot(domain,doppler_mt)
title('(a) Denoised by MIPT, D=2')

pack

subplot(2,3,5)
wc = FMIPT(heav,L,D); 
for j=(L+1):J,
  wc(triad(j)) =  HardThresh(wc(triad(j)), GaussT(j));
end
heav_mt = IMIPT(wc,L,D);
plot(domain,heav_mt)
title('(c) denoised by MIPT, D=2')

pack

% AI
subplot(2,3,3)
wc = FAIPT(doppler,L,D);
for j=(L+1):J,
  wc(triad(j)) =  HardThresh(wc(triad(j)), GaussT(j));
end
ddoppler_at = IAIPT(wc,L,D);
plot(domain,ddoppler_at)
title('(c) Denoised by AIPT, D=2')

pack

subplot(2,3,6)
wc = FAIPT(heav,L,D);
for j=(L+1):J,
  wc(triad(j)) =  HardThresh(wc(triad(j)), GaussT(j));
end
dheav_at = IAIPT(wc,L,D);
plot(domain,dheav_at)
title('(d) Denoised by AIPT, D=2')

pack

for i=[1,2,3], subplot(2,3,i), axis([0,1,-9,9]); end
for i=[4,5,6], subplot(2,3,i), axis([0,1,-4,9]); end

pack

RegisterPlot(scriptname, '{R}')

%
% Copyright (c) 1996.  David Donoho and Thomas P.Y. Yu
%
    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:41 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]