multivisu2.m

上課所教所有matlab 小波去燥的教程！~非常適合各位的參考！非常好喔！

function ws = MultiVisu2(wc,L)
%
%%%%%%%%It's My Own Function!!!
%
% MultiVisu2 -- Apply sqrt(2log(n)) Thresholding to Wavelet Coefficients
%  Usage
%    ws = MultiVisu2(wc,L)
%  Inputs 
%    wc    Wavelet Transform of noisy sequence with N(0,1) noise
%    L     low-frequency cutoff for Wavelet Transform
%  Outputs 
%    ws    result of applying VisuThresh to each High Frequency
%          Dyadic Block
%
	[n,J] = dyadlength(wc);
	ws = wc;
    for j=(J-1):-1:L
        [t1,t2]=dyad2HH(j);
        ws(t1,t2)=VisuThresh2(wc(t1,t2),'Hard');
        [t1,t2]=dyad2HL(j);
        ws(t1,t2)=VisuThresh2(wc(t1,t2),'Hard');
        [t1,t2]=dyad2LH(j);
        ws(t1,t2)=VisuThresh2(wc(t1,t2),'Hard');
    end