blob: f5de115f2966238488430478d67a7762b72e6732 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ???? - INRIA - Scilab
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
function analyze (w,fmin,fmax,rate,points)
// Make a frequency plot of the signal w with sampling rate.
// The data must be at least points long.
// The maximal frequency plotted will be fmax, the minimal fmin.
[lhs,rhs]=argn(0);
if ( rhs <= 4 ) then points=8192 ; end ;
if ( rhs <= 3 ) then rate=0; end ;
if ( rhs <= 2 ) then fmax=1500; end ;
if ( rhs <= 1 ) then fmin=10; end ;
defaultrate=22050;
if rate==0; rate=defaultrate; end;
v=w(1:points);
f=abs(fft(v,1));
i=fmin/rate*points:fmax/rate*points;
fr=i/points*rate;
plot2d(fr',f(i)');
endfunction
|
