在JavaScript中生成平滑随机趋势(随机漫游)

时间:2014-02-27 03:05:20

标签: javascript algorithm gaussian random-walk

我正在寻找随机游走/随机趋势算法的JavaScript实现。我需要一些能够坚持趋势的东西(所以,只是普通的随机增量),同时仍然保持在一些指定的边界内。 我尝试通过基于加权平均值选择数字来写出一些东西(重量是使用高斯函数计算的),最后是一条稍微平滑的线(不够好)。然后,我采取了一种不那么直接的方法并尝试在互联网上搜索,虽然我发现了一些概述,但没有具体的内容让我理解。

事实证明(我最初并没有意识到这一点),似乎已经有一系列算法使用高斯方程来产生随机趋势。但是,尽管经过了数小时的搜索,我找不到比抽象方程更多的东西,这对我来说毫无用处。我能找到的最好的是this blog,他在那里展示了我正在寻找的随机数据的图片。他列出了方程式,但我不知道那些实际上应该是什么意思(对我来说,它似乎并不像是一个完整的解决方案)。

那里有哪些算法(最好是JavaScript或类似C的实现)来生成这样的数据?

1 个答案:

答案 0 :(得分:9)

以下是我通过阅读您关联的博客提出的建议。据我所知,这是作者为他的第一张图所做的。

CSS

#container {
    min-width: 310px;
    height: 400px;
    margin: 0 auto;
}

HTML

<div id="container"></div>

的Javascript

Box–Muller transform生成Gaussian Random Numbers

var boxMullerRandom = (function () {
    var phase = 0,
        RAND_MAX,
        array,
        random,
        x1, x2, w, z;

    if (crypto && typeof crypto.getRandomValues === 'function') {
        RAND_MAX = Math.pow(2, 32) - 1;
        array = new Uint32Array(1);
        random = function () {
            crypto.getRandomValues(array);

            return array[0] / RAND_MAX;
        };
    } else {
        random = Math.random;
    }

    return function () {
        if (!phase) {
            do {
                x1 = 2.0 * random() - 1.0;
                x2 = 2.0 * random() - 1.0;
                w = x1 * x1 + x2 * x2;
            } while (w >= 1.0);

            w = Math.sqrt((-2.0 * Math.log(w)) / w);
            z = x1 * w;
        } else {
            z = x2 * w;
        }

        phase ^= 1;

        return z;
    }
}());

Random Walk生成器

function randomWalk(steps, randFunc) {
    steps = steps >>> 0 || 100;
    if (typeof randFunc !== 'function') {
        randFunc = boxMullerRandom;
    }

    var points = [],
        value = 0,
        t;

    for (t = 0; t < steps; t += 1) {
        value += randFunc();
        points.push([t, value]);
    }

    return points;
}

帮助函数从随机游走点获取Y值

function getYValues(points) {
    return points.map(function (point) {
        return point[1];
    });
}

为图表生成X图的辅助函数

function generatePlots(howMany) {
    howMany = howMany >>> 0 || 10;
    var plots = [],
        index;

    for (index = 0; index < howMany; index += 1) {
        plots.push({
            name: 'plot' + index,
            data: getYValues(randomWalk())
        });
    }

    return plots;
}

绘制结果图表,使用jQueryhighcharts.js

$('#container').highcharts({
    title: {
        text: 'Random Walk',
        x: -20 //center
    },
    subtitle: {
        text: 'Random Walk',
        x: -20
    },
    xAxis: {
        type: 'linear'
    },
    yAxis: {
        title: {
            text: 'Value'
        },
        plotLines: [{
            value: 0,
            width: 1,
            color: '#808080'
        }]
    },
    tooltip: {
        valueSuffix: ' units'
    },
    legend: {
        layout: 'vertical',
        align: 'right',
        verticalAlign: 'middle',
        borderWidth: 0
    },
    series: generatePlots(10)
});

jsFiddle