THESE/codes/wave/allcode/rden2.m

   1 function re = rden2(s,x,t)
   2
   3 % % Example
   4 % s1 = double(imread('st.tif'));
   5 % s = s1(:,:,3);
   6 % x = s + 20*randn(size(s));
   7 % t = 0:5:45;
   8 % e = rden2(s,x,t);
   9 % plot(t,e);
  10
  11 N = length(t);
  12 for k = 1:N
  13     y = denR2D(x,t(k));
  14     re(k) = sqrt(mean(mean((y-s).^2)));
  15 end
  16