X-Git-Url: http://sraa.de/git/?a=blobdiff_plain;f=dol%2Fexamples%2Fexample5%2Ffft_script.m;fp=dol%2Fexamples%2Fexample5%2Ffft_script.m;h=4a418ea2bfa1fffec3e20636b98809fee5470116;hb=8c411cf24ed0eb889191aaeafd8fa1e69081df42;hp=0000000000000000000000000000000000000000;hpb=dea7a4fb1ed110d3ce6e6d9255103d724bd66c0e;p=jump.git

diff --git a/dol/examples/example5/fft_script.m b/dol/examples/example5/fft_script.m
new file mode 100644
index 0000000..4a418ea
--- /dev/null
+++ b/dol/examples/example5/fft_script.m
@@ -0,0 +1,146 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% Matlab script to compute the reference result for the FFT example.
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+close all;
+clear all;
+clc;
+
+wNk = inline('exp(-j * 2 * pi * k / N)', 'N', 'k');
+
+%coefficients produced by c's random number generator in producer_fire()
+x = [ -9 + j *  4
+      -2 + j *  0
+      -3 + j *  0
+       8 + j *  3
+       6 + j *  0
+      -7 + j * -9
+      -5 + j *  2
+       1 + j *  2
+       8 + j * -4
+      -6 + j * -4
+       9 + j * -2
+      -8 + j * -4
+      -3 + j *  1
+      -9 + j *  7
+      -9 + j * -9
+      -7 + j *  7
+       4 + j * -6
+       5 + j * -5
+       4 + j *  8
+       7 + j * -6
+      -3 + j *  9
+      -6 + j * -2
+      -8 + j * -9
+       2 + j *  7
+      -5 + j * -7
+       7 + j * -9
+       0 + j *  9
+       5 + j *  1
+      -8 + j *  1
+       5 + j *  6
+      -2 + j * -8
+       6 + j * -8
+      -5 + j *  7
+       5 + j * -9
+       8 + j * -5
+      -6 + j *  9
+       2 + j * -3
+       3 + j * -6
+       5 + j *  0
+      -7 + j *  6
+       3 + j * -4
+      -2 + j *  4
+      -4 + j * -2
+       9 + j * -8
+      -3 + j *  3
+       0 + j * -6
+       6 + j *  4
+      -2 + j * -4
+       9 + j *  0
+      -8 + j * -9
+       0 + j * -4
+       7 + j * -4
+       8 + j *  6
+      -2 + j * -1
+      -8 + j *  6
+       9 + j * -3
+       1 + j *  4
+       3 + j * -3
+       2 + j * -2
+       9 + j *  6
+      -8 + j * -7
+      -7 + j *  2
+       8 + j * -1
+       1 + j *  0
+      -1 + j * -3
+       3 + j * -9
+       1 + j *  9
+      -2 + j *  8
+       8 + j * -5
+       6 + j *  4
+       8 + j *  3
+       9 + j *  9
+       3 + j * -9
+       4 + j * -5
+       4 + j * -2
+       8 + j *  2
+      -3 + j *  1
+      -8 + j * -9
+       4 + j *  6
+       4 + j *  7
+       2 + j *  7
+      -4 + j * -6
+       6 + j *  9
+       3 + j * -5
+       6 + j * -7
+      -4 + j * -6
+       3 + j * -1
+       5 + j *  8
+       6 + j * -6
+      -5 + j *  0
+       2 + j * -2
+      -2 + j * -5
+      -4 + j * -8
+       9 + j * -5
+       3 + j *  3
+      -3 + j * -9
+      -5 + j * -7
+       6 + j * -6
+      -9 + j * -3
+       1 + j * -9
+       3 + j * -5
+       3 + j *  9
+       8 + j * -9
+      -8 + j *  3
+       1 + j * -5
+       9 + j *  4
+      -6 + j *  7
+       3 + j *  5
+       9 + j *  9
+      -4 + j * -6
+      -1 + j * -8
+      -6 + j *  3
+       4 + j * -3
+      -6 + j * -7
+       2 + j * -3
+      -8 + j *  0
+       1 + j *  9
+       9 + j *  0
+      -3 + j *  2
+      -2 + j *  7
+       8 + j * -9
+       1 + j *  6
+      -4 + j * -1
+       5 + j * -1
+       8 + j *  9
+       3 + j * -9
+       8 + j *  3
+      -3 + j * -2 ];
+
+fft(x(1:8))
+fft(x(1:16))
+fft(x(1:64))
+fft(x)