fourierdenoise.m

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

function recon = FourierDeNoise(y)
% FourierDeNoise -- Hard thresholding in Fourier Domain
%  Usage
%    recon = FourierDeNoise(y)
%  Inputs
%    y        1-d complex valued signal; 
%             noise should be complex-valued white noise
%  Outputs
%    recon    reconstruction after sqrt(2 log(n)) thresholding
%
	n = length(y);
	fy = fft(y);
	af = abs(fy);
	ty = af > (sqrt(2 * log(n)) .* sqrt(n));
	cy = fy .* ty;
	recon = ifft(cy);
    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:42 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]