THESE/codes/wave/allcode/AntonB.m

   1 function [af, sf] = AntonB
   2
   3
   4 a0 = [
   5    0
   6    0.02674875741081
   7   -0.01686411844287
   8   -0.07822326652899
   9    0.26686411844288
  10    0.60294901823636
  11    0.26686411844287
  12   -0.07822326652899
  13   -0.01686411844287
  14    0.02674875741081
  15    0
  16    0
  17 ];
  18
  19 a1 = [
  20    0
  21    0
  22    0
  23    0.04563588155712
  24   -0.02877176311425
  25   -0.29563588155712
  26    0.55754352622850
  27   -0.29563588155713
  28   -0.02877176311425
  29    0.04563588155712
  30    0
  31    0
  32 ];
  33
  34 s0 = [
  35    0
  36    0
  37    0
  38   -0.04563588155712
  39   -0.02877176311425
  40    0.29563588155712
  41    0.55754352622850
  42    0.29563588155713
  43   -0.02877176311425
  44   -0.04563588155712
  45    0
  46    0
  47 ];
  48
  49 s1 = [
  50    0
  51    0.02674875741081
  52    0.01686411844287
  53   -0.07822326652899
  54   -0.26686411844288
  55    0.60294901823636
  56   -0.26686411844287
  57   -0.07822326652899
  58    0.01686411844287
  59    0.02674875741081
  60    0
  61    0
  62 ];
  63
  64 s0 = 2*s0;
  65 s1 = 2*s1;
  66
  67 aa0 = [
  68    0
  69    0
  70    0.02674875741081
  71   -0.01686411844287
  72   -0.07822326652899
  73    0.26686411844288
  74    0.60294901823636
  75    0.26686411844287
  76   -0.07822326652899
  77   -0.01686411844287
  78    0.02674875741081
  79    0
  80 ];
  81
  82 aa1 = [
  83    0
  84    0
  85    0
  86    0
  87    0.04563588155712
  88   -0.02877176311425
  89   -0.29563588155712
  90    0.55754352622850
  91   -0.29563588155713
  92   -0.02877176311425
  93    0.04563588155712
  94    0
  95
  96 ];
  97
  98 ss0 = [
  99    0
 100    0
 101   -0.04563588155712
 102   -0.02877176311425
 103    0.29563588155712
 104    0.55754352622850
 105    0.29563588155713
 106   -0.02877176311425
 107   -0.04563588155712
 108    0
 109    0
 110    0
 111 ];
 112
 113 ss1 = [
 114    0.02674875741081
 115    0.01686411844287
 116   -0.07822326652899
 117   -0.26686411844288
 118    0.60294901823636
 119   -0.26686411844287
 120   -0.07822326652899
 121    0.01686411844287
 122    0.02674875741081
 123    0
 124 0
 125 0
 126 ];
 127
 128 ss0 = 2*ss0;
 129 ss1 = 2*ss1;
 130
 131
 132 af{1} = [a0 a1];
 133
 134 sf{1} = [s0 s1];
 135
 136 af{2} = [aa0 aa1];
 137
 138 sf{2} = [ss0 ss1];