C# 普通随机数 [重复]

Question

我想创建一个函数，它接受Double mean、Double deviation 并返回一个正态分布的随机数。

示例：如果我传入 5.00 作为平均值，2.00 作为偏差，68% 的时间我会得到一个介于 3.00 和 7.00 之间的数字

我的统计数据有点弱……。任何人都知道我应该如何处理这个问题？我的实现将是 C# 2.0，但只要数学函数是标准的，您可以随意用您选择的语言回答。

我认为this 可能实际上是我正在寻找的。任何帮助将其转换为代码？

提前感谢您的帮助。

Answer 1

请参阅此 CodeProject 文章：Simple Random Number Generation。代码很短，它从均匀分布、正态分布和指数分布中生成样本。

Answer 2

您可能对 Math.NET，特别是 Numerics 包感兴趣。

警告：数字包面向 .NET 3.5。如果您的目标是早期版本，您可能需要使用 Iridium 包...

Answer 3

这是一些返回两个值（rand1 和 rand2）的 C，只是因为算法有效地这样做了。它是Box-Muller transform 的极性形式。

void RandVal (double mean1, double sigma1, double *rand1, double mean2, double sigma2, double *rand2)
{
double u1, u2, v1, v2, s, z1, z2;

do {
    u1 = Random (0., 1.);  // a uniform random number from 0 to 1
    u2 = Random (0., 1.);
    v1 = 2.*u1 - 1.;
    v2 = 2.*u2 - 1.;
    s = v1*v1 + v2*v2;
} while (s > 1. || s==0.); 

z1 = sqrt (-2.*log(s)/s)*v1;
z2 = sqrt (-2.*log(s)/s)*v2;
*rand1 = (z1*sigma1 + mean1);
*rand2 = (z2*sigma2 + mean2);
return;

}

Answer 4

这个库也不错：

.NET random number generators and distributions

Answer 5

抱歉，我没有任何代码给您，但我可以将您指向 some algorithms on Wikipedia。我猜你选择的算法取决于你想要它有多准确以及它需要多快。

Answer 6

对于那些引用这个问题的人，一个简单的解决方案可能是：

Random rand = new Random();
double normRand  = alglib.invnormaldistribution(rand.NextDouble())

根据需要按 mu 和 sigma 缩放。
alglib 库位于www.alglib.net

Answer 7

MetaNumerics 库，也是 .NET，将超级快速地计算正态分布（以及来自统计数据的任何其他内容）。查看功能页面了解更多详细信息。 Codeplex 页面在这里：http://metanumerics.codeplex.com/。

Answer 8

数学网

来自第二个最佳答案

    public static double GenerateRandomVariant(double mean,double deviation,System.Random rand=null, int factor=1)
    {

        rand = rand ?? new Random();
        double randNormal=(MathNet.Numerics.Distributions.Normal.Sample(rand, mean , deviation));
        return factor * randNormal;

    }

Box-Mueller 变换

通过链接获得最佳答案（快两倍？）

你/yoyoyoyosef Random Gaussian Variables

    public static double GenerateRandomVariant(double mean, double deviation, System.Random rand=null, int factor = 1)
    {
        rand = rand ?? new Random();
        double u1 = 1.0 - rand.NextDouble(); //uniform(0,1] random doubles
        double u2 = 1.0 - rand.NextDouble();
        double randStdNormal = Math.Sqrt(-2.0 * Math.Log(u1)) *
                     Math.Sin(2.0 * Math.PI * u2); //random normal(0,1)
        double randNormal=(
                     mean +  deviation * randStdNormal); //random normal(mean,stdDev^2)
        return randNormal * factor;
    }

Answer 9

我知道这篇文章有点老了，但我想分享我昨天创建的一个小项目。 我认为更简单的方法是使用 C++ 11 并在托管 C++ 中创建一个 .dll。 有一个link 源代码和一个包含已编译 dll 的 zip。

还有我制作的代码：

// NormalDistributionRandom.h
#include <random>

#pragma once

using namespace System;

namespace NormalDistribution 
{

    class _NormalDistributionRandom
    {
        std::default_random_engine engine;
        std::normal_distribution<double> distribution;

    public:
        _NormalDistributionRandom(double mean, double deviation) : distribution(mean, deviation)
        {
        }

        double Next()
        {
            return distribution(engine);
        }
    };

    public ref class NormalDistributionRandom
    {
    private:

        void* Distribution;

    public:

        NormalDistributionRandom( double mean, double deviation)
        {
            Distribution = new _NormalDistributionRandom(mean, deviation);
        }

        double Next()
        {
            return ((_NormalDistributionRandom*)Distribution)->Next();
        }
        ~NormalDistributionRandom()
        {
            this->!NormalDistributionRandom();
        }
    protected:
        !NormalDistributionRandom()
        {
            if (Distribution != nullptr)
            {
                delete (_NormalDistributionRandom*)Distribution;
                Distribution = nullptr;
            }
        }
    };

}