C中正态分布的错误答案

Question

我使用此链接中给出的方程应用了正态分布函数： http://en.wikipedia.org/wiki/Normal_distribution

我为正态分布编写了这段代码

float m=2.0;
float s=0.5;
float x[3]={1.0, 2.0, 3.0};
float xs[3];
const float pi = 3.141;

for (int i=0; i<3; i++)
{
    xs[i]=(1/s*sqrt(2*pi))*exp(-(pow(x[i]-m,2)/2*pow(s,2)));
    printf("\n%f",xs[i]);
}

此代码的答案是 4.42、5.01、4.42

我在 Matlab 中有相同的代码

x_s_new=[1 2 3];
x_s=2+0.5.*randn(1,9);
x_s=x_s_new+0.5.*randn(1,3);
plot(x_s_new)

但是 matlab 中的答案是 0.8 , 1.9 , 3.7

谁能告诉我哪里出错了？

我想使用 C 应用正态分布

谢谢:)

Answer 1

您的 Matlab 代码与您的 C 代码不同。 C 代码计算正态分布在某些点的概率密度函数的值。 Matlab 代码根据正态分布生成随机数。

C代码对应的Matlab代码为

m = 2;
s = 0.5;
x = [1 2 3];

for i = 1 : 3
    xs(i) = (1/s*sqrt(2*pi)) * exp(-( (x(i)-m)^2/2*s^2));
    fprintf('%f\n',xs(i));
end

并给出相同的结果

4.424183
5.013257
4.424183

但是，有一个错误，因为 / 在 Matlab 和 C 中只适用于紧邻的下一个操作数。正确的应该是

xs(i) = 1/(s*sqrt(2*pi)) * exp(-( (x(i)-m)^2/(2*s^2)));

相应地在 C 中。

要将使用randn 的代码翻译成C - 据我所知，C 中没有标准函数可以生成正态分布的随机数，您需要找到包含此类函数的库，或者从@ 开始自己构建它987654326@.