THESE/codes/wave/allcode/dualtree_eg1.m

   1 J = 5;  % J: number of stages\r
   2 \r
   3 % get filters\r
   4 [Faf, Fsf] = FSfarras;\r
   5 [af, sf] = dualfilt1;\r
   6 \r
   7 x = zeros(1,256);  % zero signal\r
   8 \r
   9 % Compute dual-tree complex DWT of zero signal\r
  10 w = dualtree(x, J, Faf, af); \r
  11 % Set a single (real) coefficient to 1\r
  12 w{5}{1}(4) = 1;\r
  13 % Compute the inverse transform \r
  14 y1 = idualtree(w, J, Fsf, sf);\r
  15 \r
  16 % Compute dual-tree complex DWT of zero signal\r
  17 w = dualtree(x, J, Faf, af); \r
  18 % Set a single (imaginary) coefficient to 1\r
  19 w{5}{2}(4) = 1;\r
  20 % Compute the inverse transform \r
  21 y2 = idualtree(w, J, Fsf, sf);\r
  22 \r
  23 % Display real and imaginary parts and magnitude\r
  24 n = [1:256]/256;\r
  25 plot(n,y1,n,y2,n,sqrt(y1.^2+y2.^2))\r
  26 title('COMPLEX 1D WAVELET') \r
  27 xlabel('t');\r
  28 ylabel('\psi(t)');\r