THESE/codes/wave/allcode/bishrink.m

   1 function [w1] = bishrink(y1,y2,T)
   2 % Bivariate Shrinkage Function
   3 % Usage :
   4 %      [w1] = bishrink(y1,y2,T)
   5 % INPUT :
   6 %      y1 - a noisy coefficient value
   7 %      y2 - the corresponding parent value
   8 %      T  - threshold value
   9 % OUTPUT :
  10 %      w1 - the denoised coefficient
  11
  12 R  = sqrt(abs(y1).^2 + abs(y2).^2);
  13 R = R - T;
  14 R  = R .* (R > 0);
  15 w1 = y1 .* R./(R+T);