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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
function [y,c]= stft(x, varargin)
//Compute the short-time Fourier transform of the vector X
//Calling Sequence
//Y = stft (X)
//Y = stft (X, WIN_SIZE)
//Y = stft (X, WIN_SIZE, INC)
//Y = stft (X, WIN_SIZE, INC, NUM_COEF)
//Y = stft (X, WIN_SIZE, INC, NUM_COEF, WIN_TYPE)
//[Y,C] = stft (X)
//[Y,C] = stft (X, WIN_SIZE)
//[Y,C] = stft (X, WIN_SIZE, INC)
//[Y,C] = stft (X, WIN_SIZE, INC, NUM_COEF)
//[Y,C] = stft (X, WIN_SIZE, INC, NUM_COEF, WIN_TYPE)
//Parameters
//X: Real scalar or vector
//WIN_SIZE: Size of the window used
//INC: Increment
//WIN_TYPE: Type of window
//Description
//Compute the short-time Fourier transform of the vector X with NUM_COEF coefficients by applying a window of WIN_SIZE data points and an increment of INC points.
//
//Before computing the Fourier transform, one of the following windows is applied:
//
//"hanning" -> win_type = 1
//
//"hamming" -> win_type = 2
//
//"rectangle" -> win_type = 3
//
//The window names can be passed as strings or by the WIN_TYPE number.
//
//The following defaults are used for unspecified arguments:WIN_SIZE= 80, INC = 24, NUM_COEF = 64, and WIN_TYPE = 1.
//
//Y = stft (X, ...)' returns the absolute values of the Fourier coefficients according to the NUM_COEF positive frequencies.
//
//'[Y, C] = stft (x, ...)' returns the entire STFT-matrix Y and a 3-element vector C containing the window size, increment, and window type, which is needed by the 'synthesis' function.
funcprot(0);
lhs= argn(1);
rhs= argn(2);
if(rhs <1 | rhs>5)
error("Wrong number of input arguments");
end
if(lhs<1 | lhs>2)
error("Wrong number of output arguments");
end
select(rhs)
case 1 then
select(lhs)
case 1 then
y= callOctave("stft", x);
case 2 then
[y,c]= callOctave("stft", x);
end
case 2 then
select(lhs)
case 1 then
y= callOctave("stft", x,varargin(1));
case 2 then
[y,c]= callOctave("stft", x, varargin(1));
end
case 3 then
select(lhs)
case 1 then
y= callOctave("stft", x,varargin(1), varargin(2));
case 2 then
[y,c]= callOctave("stft", x,varargin(1), varargin(2));
end
case 4 then
select(lhs)
case 1 then
y= callOctave("stft", x,varargin(1), varargin(2), varargin(3));
case 2 then
[y,c]= callOctave("stft", x,varargin(1), varargin(2), varargin(3));
end
case 5 then
select(lhs)
case 1 then
y= callOctave("stft", x,varargin(1), varargin(2), varargin(3), varargin(4));
case 2 then
[y,c]= callOctave("stft", x,varargin(1), varargin(2), varargin(3), varargin(4));
end
end
endfunction
|
