Linux中高斯分布的参数设置,C++与正态分布

正态分布(Normal distribution)又名高斯分布(Gaussiandistribution)。若随机变量X服从一个数学期望为μ、方差为σ^2的高斯分布,记为N(μ,σ^2)。其概率密度函数为正态分布的期望值μ决定了其位置,其标准差σ决定了分布的幅度。我们通常所说的标准正态分布是μ = 0,σ = 1的正态分布。

ce46c4cbd3f1c5fd62d267f6c9a61380.png

ecf5fe8c60ddd26929f90cb3802ca528.png

f355b21c4b865f1856781aa1813b26aa.png

从上图可以看出,当相差1个方差(σ), 满足要求的面积有68.27%.

当相差2个方差(σ)时,满足要求的面积有95.45.

当相差3个方差(σ)时,满足要求的面积有99.73%.

满足标准正态分的曲线,可以查表来求得正态分布的幅度.(见文后所附表格)

方差(Variance),是各个数据分别与其和的平均数之差的平方的和的平均数,用字母D表示。在概率论和数理统计中,方差用来度量随机变量和其数学期望(即均值)之间的偏离程度。

标准差(StandardDeviation),是离均差平方和平均后的方根,用σ表示。标准差是方差的算术平方根。标准差能反映一个数据集的离散程度。

测试代码:

// NormalDistribution.cpp : Defines the entry point for the console application.

//

#include

#include

#include

#include

#include

#define _USE_MATH_DEFINES

#include

using namespace std;

// 高斯分布随机数系列,默认期望值为0,方差为1

double GaussRand(double dExpect = 0, double dVariance = 1);

double GaussRand(double dExpect, double dVariance)

{

static double V1, V2, S;

static int phase = 0;

double X;

if ( phase == 0 )

{

do

{

double U1 = (double)rand() / RAND_MAX;

double U2 = (double)rand() / RAND_MAX;

V1 = 2 * U1 - 1;

V2 = 2 * U2 - 1;

S = V1 * V1 + V2 * V2;

} while(S >= 1 || S == 0);

X = V1 * sqrt(-2 * log(S) / S);

}

else

{

X = V2 * sqrt(-2 * log(S) / S);

}

phase = 1 - phase;

return (X*dVariance + dExpect);

}

int _tmain(int argc, _TCHAR* argv[])

{

const int DATA_CNT = 100000;

double dArrData[DATA_CNT] = {0};

double dSum = 0;

// 对所有数赋随机数,默认期望值为0,方差为1

srand(GetTickCount());

for (int nIdx = 0; nIdx < DATA_CNT; nIdx++)

{

// 防止计算方差时数值过大

dArrData[nIdx] = GaussRand();

dSum += dArrData[nIdx];

}

// 求平均数

double dAverageData = dSum / DATA_CNT;

// 计算所有的数的方差(各个数据分别与其和的平均数之差的平方的和的平均数)

double dVariance = 0.0;

for (int nIdx = 0; nIdx < DATA_CNT; nIdx++)

{

double dDeviate = dArrData[nIdx] - dAverageData;

dVariance += pow(dDeviate, 2);

}

dVariance /= DATA_CNT;

// 计算标准差(方差的算术平方根,反映一组数据的离散程序)

double dStandardDeviation = sqrt(dVariance);

// 计算0.5个正负标准差之间包含的数字个数

int nDataCnt = 0;

for (int nIdx = 0; nIdx < DATA_CNT; nIdx++)

{

double dDeviate = dArrData[nIdx] - dAverageData;

if (abs(dDeviate) <= 0.5*dStandardDeviation)

{

nDataCnt++;

}

}

cout<

// 计算1个正负标准差之间包含的数字个数

nDataCnt = 0;

for (int nIdx = 0; nIdx < DATA_CNT; nIdx++)

{

double dDeviate = dArrData[nIdx] - dAverageData;

if (abs(dDeviate) <= dStandardDeviation)

{

nDataCnt++;

}

}

cout<

// 计算2个正负标准差之间包含的数字个数

nDataCnt = 0;

for (int nIdx = 0; nIdx < DATA_CNT; nIdx++)

{

double dDeviate = dArrData[nIdx] - dAverageData;

if (abs(dDeviate) <= 2*dStandardDeviation)

{

nDataCnt++;

}

}

cout<

// 计算3个正负标准差之间包含的数字个数

nDataCnt = 0;

for (int nIdx = 0; nIdx < DATA_CNT; nIdx++)

{

double dDeviate = dArrData[nIdx] - dAverageData;

if (abs(dDeviate) <= 3*dStandardDeviation)

{

nDataCnt++;

}

}

cout<

return 0;

}

6a01612ca5ebb2ca641cd9a63c13a372.png

(附)标准正态分布表

φ( -x) = 1–φ(x)

x

0

0.01

0.02

0.03

0.04

0.05

0.06

0.07

0.08

0.09

0

0.500 0

0.504 0

0.508 0

0.512 0

0.516 0

0.519 9

0.523 9

0.527 9

0.531 9

0.535 9

0.1

0.539 8

0.543 8

0.547 8

0.551 7

0.555 7

0.559 6

0.563 6

0.567 5

0.571 4

0.575 3

0.2

0.579 3

0.583 2

0.587 1

0.591 0

0.594 8

0.598 7

0.602 6

0.606 4

0.610 3

0.614 1

0.3

0.617 9

0.621 7

0.625 5

0.629 3

0.633 1

0.636 8

0.640 4

0.644 3

0.648 0

0.651 7

0.4

0.655 4

0.659 1

0.662 8

0.666 4

0.670 0

0.673 6

0.677 2

0.680 8

0.684 4

0.687 9

0.5

0.691 5

0.695 0

0.698 5

0.701 9

0.705 4

0.708 8

0.712 3

0.715 7

0.719 0

0.722 4

0.6

0.725 7

0.729 1

0.732 4

0.735 7

0.738 9

0.742 2

0.745 4

0.748 6

0.751 7

0.754 9

0.7

0.758 0

0.761 1

0.764 2

0.767 3

0.770 3

0.773 4

0.776 4

0.779 4

0.782 3

0.785 2

0.8

0.788 1

0.791 0

0.793 9

0.796 7

0.799 5

0.802 3

0.805 1

0.807 8

0.810 6

0.813 3

0.9

0.815 9

0.818 6

0.821 2

0.823 8

0.826 4

0.828 9

0.835 5

0.834 0

0.836 5

0.838 9

1

0.841 3

0.843 8

0.846 1

0.848 5

0.850 8

0.853 1

0.855 4

0.857 7

0.859 9

0.862 1

1.1

0.864 3

0.866 5

0.868 6

0.870 8

0.872 9

0.874 9

0.877 0

0.879 0

0.881 0

0.883 0

1.2

0.884 9

0.886 9

0.888 8

0.890 7

0.892 5

0.894 4

0.896 2

0.898 0

0.899 7

0.901 5

1.3

0.903 2

0.904 9

0.906 6

0.908 2

0.909 9

0.911 5

0.913 1

0.914 7

0.916 2

0.917 7

1.4

0.919 2

0.920 7

0.922 2

0.923 6

0.925 1

0.926 5

0.927 9

0.929 2

0.930 6

0.931 9

1.5

0.933 2

0.934 5

0.935 7

0.937 0

0.938 2

0.939 4

0.940 6

0.941 8

0.943 0

0.944 1

1.6

0.945 2

0.946 3

0.947 4

0.948 4

0.949 5

0.950 5

0.951 5

0.952 5

0.953 5

0.953 5

1.7

0.955 4

0.956 4

0.957 3

0.958 2

0.959 1

0.959 9

0.960 8

0.961 6

0.962 5

0.963 3

1.8

0.964 1

0.964 8

0.965 6

0.966 4

0.967 2

0.967 8

0.968 6

0.969 3

0.970 0

0.970 6

1.9

0.971 3

0.971 9

0.972 6

0.973 2

0.973 8

0.974 4

0.975 0

0.975 6

0.976 2

0.976 7

2

0.977 2

0.977 8

0.978 3

0.978 8

0.979 3

0.979 8

0.980 3

0.980 8

0.981 2

0.981 7

2.1

0.982 1

0.982 6

0.983 0

0.983 4

0.983 8

0.984 2

0.984 6

0.985 0

0.985 4

0.985 7

2.2

0.986 1

0.986 4

0.986 8

0.987 1

0.987 4

0.987 8

0.988 1

0.988 4

0.988 7

0.989 0

2.3

0.989 3

0.989 6

0.989 8

0.990 1

0.990 4

0.990 6

0.990 9

0.991 1

0.991 3

0.991 6

2.4

0.991 8

0.992 0

0.992 2

0.992 5

0.992 7

0.992 9

0.993 1

0.993 2

0.993 4

0.993 6

2.5

0.993 8

0.994 0

0.994 1

0.994 3

0.994 5

0.994 6

0.994 8

0.994 9

0.995 1

0.995 2

2.6

0.995 3

0.995 5

0.995 6

0.995 7

0.995 9

0.996 0

0.996 1

0.996 2

0.996 3

0.996 4

2.7

0.996 5

0.996 6

0.996 7

0.996 8

0.996 9

0.997 0

0.997 1

0.997 2

0.997 3

0.997 4

2.8

0.997 4

0.997 5

0.997 6

0.997 7

0.997 7

0.997 8

0.997 9

0.997 9

0.998 0

0.998 1

2.9

0.998 1

0.998 2

0.998 2

0.998 3

0.998 4

0.998 4

0.998 5

0.998 5

0.998 6

0.998 6

x

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

3

0.998 7

0.999 0

0.999 3

0.999 5

0.999 7

0.999 8

0.999 8

0.999 9

0.999 9

1.000 0

(附)正态分布概率表

dc82d5c9964d184a4e09dd624fd66bbe.png

e2cc749c8e242e6dd2ab9681ed30aa58.png

31adb1d62da54a836e1e80876771517f.png

将C语言梳理一下,分布在以下10个章节中:

0b1331709591d260c1c78e86d0c51c18.png

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值