jamesstein.m

来自「老外写的小波变换的工具箱」· M 代码 · 共 9 行

function xhat = JamesStein(y)
% JamesStein -- James-Stein Shrinker
% Usage:
%    xhat = JamesStein(y)
% Inputs:
%    y      Noisy Data, Standard Deviation of Noise = 1.
% Outputs:
%    xhat   James-Stein Estimate of x
% Notes:
%    dimension(x) > 2 is required.
%
	ss  = sum(y.^2);
	n   = length(y);
	shrnk = max((ss - (n-2))/ss,0);
	xhat = shrnk .* y;

    
    
 
 
%
%  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]