波变换

开发平台：

Matlab

wbmpen.m：源码内容

							function suggthr = wbmpen(c,l,sigma,alpha,arg)
%WBMPEN Penalized threshold for wavelet 1-D or 2-D de-noising.
%   THR = WBMPEN(C,L,SIGMA,ALPHA) returns global 
%   threshold THR for denoising. THR is obtained by a wavelet 
%   coefficients selection rule using a penalization 
%   method provided by Birge-Massart.
%
%   [C,L] is the wavelet decomposition structure of the signal
%   or image to be de-noised.
%
%   SIGMA is the standard deviation of the zero mean Gaussian 
%   white noise in the de-noising model (see WNOISEST for more
%   information).
%
%   ALPHA is a tuning parameter for the penalty term and it 
%   must be a real number greater than 1. The sparsity of the
%   wavelet representation of the de-noised signal or image 
%   grows with ALPHA. Typically ALPHA = 2.
%
%   THR minimizes the penalized criterion given by:
%   let t* be the minimizer of
%   crit(t) = -sum(c(k)^2,k<=t) + 2*SIGMA^2*t*(ALPHA + log(n/t))
%   where c(k) are the wavelet coefficients sorted in  
%   decreasing order of their absolute value and n is the number   
%   of coefficients; then THR = |c(t*)|.
%
%   WBMPEN(C,L,SIGMA,ALPHA,ARG) computes the global threshold
%   and, in addition, plots three curves:  
%   2*SIGMA^2*t*(ALPHA + log(n/t)), sum(c(k)^2,k<=t) and crit(t).
%   
%   See also WDEN, WDENCMP, WPBMPEN, WPDENCMP.
%   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 26-Oct-98.
%   Last Revision: 14-May-2003.
%   Copyright 1995-2004 The MathWorks, Inc.
%   $Revision: 1.6.4.2 $  $Date: 2004/03/15 22:42:26 $
% Check arguments and set problem dimension.
nbIn = nargin;
if nbIn < 4
  error('Not enough input arguments.');
elseif nbIn > 5
  error('Too many input arguments.');
end
dim = 1; if min(size(l))~=1, dim = 2; end
if dim==1 , last = l(1); else , last = prod(l(1,:)); end
nbcfs  = prod(size(c));
c      = c(last+1:end);
thresh = sort(abs(c));
thresh = thresh(end:-1:1);
rl2scr = cumsum(thresh.^2);
xpen   = [1:length(thresh)];
pen    = 2*sigma^2*xpen.*(alpha + log(nbcfs./xpen));
[dummy,indmin] = min(pen-rl2scr);
suggthr = thresh(indmin);
if nargin==5
    h = figure;
    subplot(311),plot(xpen,pen),xlabel('pen')
    subplot(312),plot(xpen,rl2scr),xlabel('rl2scr')
    subplot(313),plot(xpen,pen-rl2scr),xlabel('crit')
end

						