dol/examples/example5/fft_script.m

   1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   2 %
   3 % Matlab script to compute the reference result for the FFT example.
   4 %
   5 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   6
   7 close all;
   8 clear all;
   9 clc;
  10
  11 wNk = inline('exp(-j * 2 * pi * k / N)', 'N', 'k');
  12
  13 %coefficients produced by c's random number generator in producer_fire()
  14 x = [ -9 + j *  4
  15       -2 + j *  0
  16       -3 + j *  0
  17        8 + j *  3
  18        6 + j *  0
  19       -7 + j * -9
  20       -5 + j *  2
  21        1 + j *  2
  22        8 + j * -4
  23       -6 + j * -4
  24        9 + j * -2
  25       -8 + j * -4
  26       -3 + j *  1
  27       -9 + j *  7
  28       -9 + j * -9
  29       -7 + j *  7
  30        4 + j * -6
  31        5 + j * -5
  32        4 + j *  8
  33        7 + j * -6
  34       -3 + j *  9
  35       -6 + j * -2
  36       -8 + j * -9
  37        2 + j *  7
  38       -5 + j * -7
  39        7 + j * -9
  40        0 + j *  9
  41        5 + j *  1
  42       -8 + j *  1
  43        5 + j *  6
  44       -2 + j * -8
  45        6 + j * -8
  46       -5 + j *  7
  47        5 + j * -9
  48        8 + j * -5
  49       -6 + j *  9
  50        2 + j * -3
  51        3 + j * -6
  52        5 + j *  0
  53       -7 + j *  6
  54        3 + j * -4
  55       -2 + j *  4
  56       -4 + j * -2
  57        9 + j * -8
  58       -3 + j *  3
  59        0 + j * -6
  60        6 + j *  4
  61       -2 + j * -4
  62        9 + j *  0
  63       -8 + j * -9
  64        0 + j * -4
  65        7 + j * -4
  66        8 + j *  6
  67       -2 + j * -1
  68       -8 + j *  6
  69        9 + j * -3
  70        1 + j *  4
  71        3 + j * -3
  72        2 + j * -2
  73        9 + j *  6
  74       -8 + j * -7
  75       -7 + j *  2
  76        8 + j * -1
  77        1 + j *  0
  78       -1 + j * -3
  79        3 + j * -9
  80        1 + j *  9
  81       -2 + j *  8
  82        8 + j * -5
  83        6 + j *  4
  84        8 + j *  3
  85        9 + j *  9
  86        3 + j * -9
  87        4 + j * -5
  88        4 + j * -2
  89        8 + j *  2
  90       -3 + j *  1
  91       -8 + j * -9
  92        4 + j *  6
  93        4 + j *  7
  94        2 + j *  7
  95       -4 + j * -6
  96        6 + j *  9
  97        3 + j * -5
  98        6 + j * -7
  99       -4 + j * -6
 100        3 + j * -1
 101        5 + j *  8
 102        6 + j * -6
 103       -5 + j *  0
 104        2 + j * -2
 105       -2 + j * -5
 106       -4 + j * -8
 107        9 + j * -5
 108        3 + j *  3
 109       -3 + j * -9
 110       -5 + j * -7
 111        6 + j * -6
 112       -9 + j * -3
 113        1 + j * -9
 114        3 + j * -5
 115        3 + j *  9
 116        8 + j * -9
 117       -8 + j *  3
 118        1 + j * -5
 119        9 + j *  4
 120       -6 + j *  7
 121        3 + j *  5
 122        9 + j *  9
 123       -4 + j * -6
 124       -1 + j * -8
 125       -6 + j *  3
 126        4 + j * -3
 127       -6 + j * -7
 128        2 + j * -3
 129       -8 + j *  0
 130        1 + j *  9
 131        9 + j *  0
 132       -3 + j *  2
 133       -2 + j *  7
 134        8 + j * -9
 135        1 + j *  6
 136       -4 + j * -1
 137        5 + j * -1
 138        8 + j *  9
 139        3 + j * -9
 140        8 + j *  3
 141       -3 + j * -2 ];
 142
 143 fft(x(1:8))
 144 fft(x(1:16))
 145 fft(x(1:64))
 146 fft(x)