正式版

wzp

2021-05-13 7d694a9113118daec5be7ac224dab46a3b20f106

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
'use strict';
 
var assert = require('assert-plus');
 
/**
 * Compute the exponential weighted moving average of a series of values.
 * The time at which you insert the value into `Ewma` is used to compute a
 * weight (recent points are weighted higher).
 * The parameter for defining the convergence speed (like most decay process) is
 * the half-life.
 *
 * e.g. with a half-life of 10 unit, if you insert 100 at t=0 and 200 at t=10
 *      the ewma will be equal to (200 - 100)/2 = 150 (half of the distance
 *      between the new and the old value)
 *
 * @param {Number} halfLifeMs parameter representing the speed of convergence
 * @param {Number} initialValue initial value
 * @param {Object} clock clock object used to read time
 *
 * @returns {Ewma} the object computing the ewma average
 */
function Ewma(halfLifeMs, initialValue, clock) {
    assert.number(halfLifeMs, 'halfLifeMs');
    assert.ok(!Number.isNaN(halfLifeMs), 'halfLifeMs can not be NaN');
    assert.optionalNumber(initialValue, 'initialValue');
    assert.ok(!Number.isNaN(initialValue), 'initialValue can not be NaN');
 
    if (clock !== undefined) {
        assert.func(clock.now, 'clock.now');
    }
 
    this._decay = halfLifeMs;
    this._ewma = initialValue || 0;
    this._clock = clock || Date;
    this._stamp = (typeof initialValue === 'number') ? this._clock.now() : 0;
}
 
module.exports = Ewma;
 
Ewma.prototype.insert = function insert(x) {
    assert.number(x, 'x');
    assert.ok(!Number.isNaN(x), 'x can not be NaN');
 
    var self = this;
    var now = self._clock.now();
    var elapsed = now - self._stamp;
    self._stamp = now;
 
    // This seemingly magic equation is derived from the fact that we are
    // defining a half life for each value. A half life is the amount of time
    // that it takes for a value V to decay to .5V or V/2. Elapsed is the time
    // delta between this value being reported and the previous value being
    // reported. Given the half life, and the amount of time since the last
    // reported value, this equation determines how much the new value should
    // be represented in the ewma.
    // For a detailed proof read:
    // A Framework for the Analysis of Unevenly Spaced Time Series Data
    // Eckner, 2014
    var w = Math.pow(2, -elapsed / self._decay);
    self._ewma = w * self._ewma + (1.0 - w) * x;
};
 
Ewma.prototype.reset = function reset(x) {
    assert.number(x, 'x');
    assert.ok(!Number.isNaN(x), 'x can not be NaN');
 
    var self = this;
    self._stamp = self._clock.now();
    self._ewma = x;
};
 
Ewma.prototype.value = function value() {
    var self = this;
    return self._ewma;
};