ssa

1
package periodogram_test
2

3
// Входные данные для длины окна = 32
4
func Data32(N int) ([]float64, []float64, []float64) {
5

6
	// Окно Блэкмана-Харриса
7
	w32 := make([]float64, N)
8
	w32[0] = 0.0234200000000000
9
	w32[1] = 0.0200796208567949
10
	w32[2] = 0.0119986516061235
11
	w32[3] = 0.00453856337218171
12
	w32[4] = 0.00521782261016186
13
	w32[5] = 0.0219495235192088
14
	w32[6] = 0.0611985478246830
15
	w32[7] = 0.126474585987881
16
	w32[8] = 0.217470000000000
17
	w32[9] = 0.329974013305730
18
	w32[10] = 0.456501360083246
19
	w32[11] = 0.587419445831711
20
	w32[12] = 0.712282177389838
21
	w32[13] = 0.821092467276898
22
	w32[14] = 0.905301440485947
23
	w32[15] = 0.958471779849594
24
	w32[16] = 0.976640000000000
25
	w32[17] = 0.958471779849594
26
	w32[18] = 0.905301440485947
27
	w32[19] = 0.821092467276898
28
	w32[20] = 0.712282177389838
29
	w32[21] = 0.587419445831712
30
	w32[22] = 0.456501360083247
31
	w32[23] = 0.329974013305730
32
	w32[24] = 0.217470000000000
33
	w32[25] = 0.126474585987881
34
	w32[26] = 0.0611985478246830
35
	w32[27] = 0.0219495235192089
36
	w32[28] = 0.00521782261016190
37
	w32[29] = 0.00453856337218168
38
	w32[30] = 0.0119986516061234
39
	w32[31] = 0.0200796208567949
40

41
	// Входной сигнал
42
	x32 := make([]float64, N)
43
	x32[0] = 1.09754040499941
44
	x32[1] = 0.804647186185957
45
	x32[2] = 0.0975404049994096
46
	x32[3] = -0.609566376187138
47
	x32[4] = -0.902459595000591
48
	x32[5] = -0.609566376187138
49
	x32[6] = 0.0975404049994093
50
	x32[7] = 0.804647186185957
51
	x32[8] = 1.09754040499941
52
	x32[9] = 0.804647186185957
53
	x32[10] = 0.0975404049994098
54
	x32[11] = -0.609566376187137
55
	x32[12] = -0.902459595000591
56
	x32[13] = -0.609566376187138
57
	x32[14] = 0.0975404049994091
58
	x32[15] = 0.804647186185956
59
	x32[16] = 1.09754040499941
60
	x32[17] = 0.804647186185957
61
	x32[18] = 0.0975404049994101
62
	x32[19] = -0.609566376187137
63
	x32[20] = -0.902459595000591
64
	x32[21] = -0.609566376187138
65
	x32[22] = 0.0975404049994071
66
	x32[23] = 0.804647186185956
67
	x32[24] = 1.09754040499941
68
	x32[25] = 0.804647186185957
69
	x32[26] = 0.0975404049994085
70
	x32[27] = -0.609566376187137
71
	x32[28] = -0.902459595000591
72
	x32[29] = -0.609566376187138
73
	x32[30] = 0.0975404049994068
74
	x32[31] = 0.804647186185956
75

76
	// Результат выполнения
77
	rez32 := make([]float64, N/2+1)
78
	rez32[0] = 0.197079184668859
79
	rez32[1] = 0.167207747427803
80
	rez32[2] = 0.119414954305659
81
	rez32[3] = 0.572832117787472
82
	rez32[4] = 1.27047439645810
83
	rez32[5] = 0.588406913455670
84
	rez32[6] = 0.0492587604356405
85
	rez32[7] = 0.000336672614800940
86
	rez32[8] = 1.71105942505116e-32
87
	rez32[9] = 2.50134060820484e-32
88
	rez32[10] = 1.26250044383115e-32
89
	rez32[11] = 1.52094171115659e-32
90
	rez32[12] = 2.20536548117705e-31
91
	rez32[13] = 2.30042433812434e-31
92
	rez32[14] = 7.80670862679593e-32
93
	rez32[15] = 1.50014758619937e-32
94
	rez32[16] = 4.65788399041705e-32
95

96
	return x32, w32, rez32
97
}
98