st.m

来自「随机信号处理中的经典去噪算法集合」· M 代码 · 共 39 行

function [out,wcoef,wcoefrest]=ST(Noisy,L,qmf,sigma)
%  ST -- Soft Threshold Applied to Wavelet Coefficients.
%  Usage 
%    [out,wcoef,wcoefrest] = ST(Noisy,L,qmf,sigma)
%  Inputs
%    Noisy      1-d signal. length(y)= 2^J
%    L      Low-Frequency cutoff for shrinkage (e.g. L=4)
%               Should have L << J!
%    qmf    Quadrature Mirror Filter for Wavelet Transform
%               Optional, Default = Symmlet 8.
%    sigma  Standard deviation of additive Gaussian White Noise.
%  Outputs 
%    out     	estimate, obtained by applying soft thresholding on
%          	 wavelet coefficients
%    wcoef		Wavelet Transform of input signal	
%    wcoefrest    	Wavelet Transform of estimate
%


  n=length(Noisy) ;
  wcoef = FWT_PO(Noisy,L,qmf) ;
  thresh=0.5*sigma*sqrt(2*log(n));
  wcoefrest = SoftThresh(wcoef,thresh);
  out    = IWT_PO(wcoefrest,L,qmf);

% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]

    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%