add tests
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / fft.tst
1 // =============================================================================\r
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab\r
3 // Copyright (C) 2008 - DIGITEO\r
4 //\r
5 //  This file is distributed under the same license as the Scilab package.\r
6 // =============================================================================\r
7 \r
8 a = [1;2;3];\r
9 n=size(a,'*');\r
10 r1 = norm(1/n*exp(2*%i*%pi*(0:n-1)'.*.(0:n-1)/n)*a -fft(a,1));\r
11 if r1 > 10 *%eps then pause,end\r
12 \r
13 r2 = norm(exp(-2*%i*%pi*(0:n-1)'.*.(0:n-1)/n)*a -fft(a,-1));\r
14 if r2 > 100 *%eps then pause,end\r
15  \r
16 sample_rate = 1000;\r
17 t = 0:1/sample_rate:0.6;\r
18 N = size(t,'*');\r
19 s = sin(2*%pi*50*t)+sin(2*%pi*70*t+%pi/4)+grand(1,N,'nor',0,1);\r
20 y = fft(s);\r
21 if or(size(y) <> [1 601]) then pause,end\r
22 \r
23