THESE/codes/wave/allcode/dualfilt1.m

   1 function [af, sf] = dualfilt1
   2
   3 % Kingsbury Q-filters for the dual-tree complex DWT
   4 %
   5 % USAGE:
   6 %    [af, sf] = dualfilt1
   7 % OUTPUT:
   8 %    af{i}, i = 1,2 - analysis filters for tree i
   9 %    sf{i}, i = 1,2 - synthesis filters for tree i
  10 %    note: af{2} is the reverse of af{1}
  11 % REFERENCE:
  12 %    N. G. Kingsbury,  "A dual-tree complex wavelet
  13 %    transform with improved orthogonality and symmetry
  14 %    properties", Proceedings of the IEEE Int. Conf. on
  15 %    Image Proc. (ICIP), 2000
  16 % See dualtree
  17 %
  18 % WAVELET SOFTWARE AT POLYTECHNIC UNIVERSITY, BROOKLYN, NY
  19 % http://taco.poly.edu/WaveletSoftware/
  20
  21 % These cofficients are rounded to 8 decimal places.
  22
  23 af{1} = [
  24    0.03516384000000                  0
  25                   0                  0
  26   -0.08832942000000  -0.11430184000000
  27    0.23389032000000                  0
  28    0.76027237000000   0.58751830000000
  29    0.58751830000000  -0.76027237000000
  30                   0   0.23389032000000
  31   -0.11430184000000   0.08832942000000
  32                   0                  0
  33                   0  -0.03516384000000
  34  ];
  35
  36 af{2} = [
  37                   0  -0.03516384000000
  38                   0                  0
  39   -0.11430184000000   0.08832942000000
  40                   0   0.23389032000000
  41    0.58751830000000  -0.76027237000000
  42    0.76027237000000   0.58751830000000
  43    0.23389032000000                  0
  44   -0.08832942000000  -0.11430184000000
  45                   0                  0
  46    0.03516384000000                  0
  47 ];
  48
  49 sf{1} = af{1}(end:-1:1, :);
  50
  51 sf{2} = af{2}(end:-1:1, :);