THESE/codes/wave/allcode/den2.m

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