ssa

1
package movmean
2

3
import (
4
	"fmt"
5
	"testing"
6
)
7

8
func TestMovmean(t *testing.T) {
9
	// https://www.mathworks.com/help/matlab/ref/movmean.html
10
	// A = [4 8 6 -1 -2 -3 -1 3 4 5];
11
	// M = movmean(A,3)
12
	// M = 1×10
13
	// 6.0 6.0 4.3 1.0 -2.0 -2.0 -0.3  2.0 4.0 4.5
14
	A := []float64{4.0, 8.0, 6.0, -1.0, -2.0, -3.0, -1.0, 3.0, 4.0, 5.0}
15
	k := 3
16

17
	result, _ := Movmean(A, k)
18
	fmt.Println(len(result))
19
	for _, res := range result {
20
		fmt.Printf("%.2f ", res)
21
	}
22
}
23