forked from maddyblue/go-dsp
-
Notifications
You must be signed in to change notification settings - Fork 0
/
spectral.go
47 lines (40 loc) · 1.46 KB
/
spectral.go
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
/*
* Copyright (c) 2012 Matt Jibson <matt.jibson@gmail.com>
*
* Permission to use, copy, modify, and distribute this software for any
* purpose with or without fee is hereby granted, provided that the above
* copyright notice and this permission notice appear in all copies.
*
* THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
* WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
* MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
* ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
* WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
* ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
* OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
*/
// Package spectral provides spectral analysis functions for digital signal processing.
package spectral
// Segment x segmented into segments of length size with specified noverlap.
// Number of segments returned is (len(x) - size) / (size - noverlap) + 1.
func Segment(x []float64, size, noverlap int) [][]float64 {
stride := size - noverlap
lx := len(x)
var segments int
if lx == size {
segments = 1
} else if lx > size {
segments = (len(x)-size)/stride + 1
} else {
segments = 0
}
r := make([][]float64, segments)
for i, offset := 0, 0; i < segments; i++ {
r[i] = make([]float64, size)
for j := 0; j < size; j++ {
r[i][j] = x[offset+j]
}
offset += stride
}
return r
}