96e03f19a79e5afc6f87eeac1bcd1a1c8e8492e7
[scilab.git] / scilab / modules / signal_processing / tests / unit_tests / fft.tst
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) 2008 - Digiteo - unknown author
3 // Copyright (C) 2009 - Digiteo - Vincent LIARD
4 //
5 // This file must be used under the terms of the CeCILL.
6 // This source file is licensed as described in the file COPYING, which
7 // you should have received as part of this distribution.  The terms
8 // are also available at    
9 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
10
11 a = [1;2;3];
12 n=size(a,'*');
13 r1 = norm(1/n*exp(2*%i*%pi*(0:n-1)'.*.(0:n-1)/n)*a -fft([1;2;3],1));
14 if r1 > 10 *%eps then pause,end
15 r2 = norm(exp(-2*%i*%pi*(0:n-1)'.*.(0:n-1)/n)*a -fft(a,-1));
16 if r2 > 100 *%eps then pause,end
17 sample_rate = 1000;
18 t = 0:1/sample_rate:0.6;
19 N = size(t,'*');
20 s = sin(2*%pi*50*t)+sin(2*%pi*70*t+%pi/4)+grand(1,N,'nor',0,1);
21 y = fft(s);
22 if size(y) <> [1 601] then pause,end
23