㈠ 如何用c语言生成一个正态分布的样本
调试程序时，随机数种子可以设常数，例如srand(54321);
用 rand() 产生均匀分布随机数 x1,x2
利用瑞利分布得正态分布随机数 y1,y2
再按要求线性缩放一下到[0.01,2] 区间。
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
main(){
#define N 100
double rd[N];
double x1,x2,y1,y2;
double pi2=6.28318530728,mx,mi,ave=0;
int i;
//srand(54321);
srand(time(NULL));
for (i=0;i<=N-2;i=i+2){
x1=1.0*rand()/RAND_MAX;
x2=1.0*rand()/RAND_MAX;
y1= sqrt((-2.0*log(x1))) * cos(pi2*x2);
y2= sqrt((-2.0*log(x1))) * sin(pi2*x2);
rd[i]=y1;
rd[i+1]=y2;
}
mx=rd[0];mi=rd[0];
for (i=0;i<N;i++){
if (rd[i]>mx)mx=rd[i];
if (rd[i]<mi)mi=rd[i];
}
//printf("mi=%lf mx=%lf\n",mi,mx);
for (i=0;i<N;i++) rd[i] = (rd[i]-mi)/(mx-mi+0.001) * (2.0-0.01) + 0.01;
for (i=0;i<N-2;i=i+2) printf("%lf %lf\n",rd[i],rd[i+1]);
return 0;
}
㈡ C++math库里有生成高斯分布随机数的函数么
VC2008 FeturePack1 以后有,参见http://growupsoft.blog.163.com/blog/static/960729200903104720643/
对源态于可以解析表达成C++的特殊概率分布,根据概率论的原理,产生均匀分布的随机分布,而后代入分布函数就可以产生高斯分布的随机数了阿.下面原版转载:
At the core of any pseudorandom number generation software is a routine for generating uniformly distributed random integers.
In C++ TR1 you have your choice of several core generators that it calls “engines.” The following four engine classes are supported in the Visual Studio 2008 feature pack
(微软已经发布了Visual Studio 2008的Services Pack 1，它包含了此前发布的feature pack，以及完整的TR1支持，后来又发布了一个修正：VC 2008 SP1: Problems with STL/TR1 after installing VS2008 SP1，关于：VC9 SP1 Hotfix For The vector<function<FT>> Crash，关于文档：微软TR1文档).
linear_congruential uses a recurrence of the form x(i) = (A * x(i-1) + C) mod M
mersenne_twister implements the famous Mersenne Twister algorithm
subtract_with_carry uses a recurrence of the form x(i) = (x(i - R) - x(i - S) - cy(i - 1)) mod M in integer arithmetic
subract_with_carry_01 uses a recurrence of the form x(i) = (x(i - R) - x(i - S) - cy(i - 1)) mod 1 in floating point arithmetic
Each engine has a seed() method that accepts an unsigned long argument to specify the random number generation seed. It is also possible to set the seed in more detail using template parameters unique to each engine.
微软的C++ TR1可以生成以下分布的随机数：
Generates a Bernoulli distribution.
Generates a binomial distribution.
Generates an exponential distribution.
Generates a gamma distribution.
Generates a geometric distribution.
Generates a normal distribution.
Generates a Poisson distribution.
Generates a uniform integer distribution.
Generates a uniform floating-point distribution.
试验一：在VS2008中档或先建一空的VC++项目文件，然后添加新建项cpp文件行裂伍如下：
#include <random>
#include <iostream>
void main(){
std::tr1::mt19937 eng; // a core engine class：Mersenne Twister generator
std::tr1::normal_distribution<double> dist;
std::tr1::uniform_int<int> unif(1, 52);
for (int i = 0; i < 10; ++i) //产生正态分布的10个随机数
std::cout << dist(eng)<<std::endl;
for(int i = 0; i < 5; ++i) //产生均匀分布的在1到52之间的五个整数随机数
std::cout << unif(eng) << std::endl;
}
在循环中，每循环一次，就调用mt19937 eng一次，产生一个随机数输出。
试验二：关于种子seed
#include <random>
#include <iostream>
#include <time.h>
void main(){
std::tr1::mt19937 eng; // a core engine class：Mersenne Twister generator
std::tr1::normal_distribution<double> dist;
std::tr1::uniform_int<int> unif(1, 52);
eng.seed((unsigned int)time(NULL)); // reseed base engine 设置种子用#include <time.h>, 不能用#include <time>
for (int i = 0; i < 10; ++i) //产生正态分布的10个随机数
std::cout << dist(eng)<<std::endl;
//eng.seed(); // reseed base engine
for(int i = 0; i < 5; ++i) //产生均匀分布的在1到52之间的五个整数随机数
std::cout << unif(eng) << std::endl;
}
试验三：随机数写入文件
#include <random>
#include <iostream>
#include <fstream>
#include <time.h>
using namespace std;
using namespace std::tr1;
void main()
{
mt19937 eng; // a core engine class：Mersenne Twister generator
normal_distribution<double> dist;
uniform_int<int> unif(1, 52);
eng.seed((unsigned int)time(NULL)); // 设置种子用#include <time.h>, 不能用#include <time>
for (int i = 0; i < 10; ++i) //产生正态分布的10个随机数
cout << dist(eng)<<endl;
ofstream fileout("fileout.dat");
for(int i = 0; i < 5; ++i) //产生均匀分布的在1到52之间的五个整数随机数
fileout << unif(eng)<< endl;
fileout.close();
}
试验四：第三方"Mersenne Twister"随机数生成程序使用试验(程序来源：Agner Fog http://www.agner.org/random/)
// 使用说明:从网站下载压缩包，http://www.agner.org/random/randomc.zip
// 展开后,将其中的randomc.h头文件及mersenne.cpp文件Copy到项目文件夹,
// 并将它们加入到项目中，其中包括"Mersenne Twister"的实现
#include <iostream>
#include <time.h>
#include "randomc.h" // define classes for random number generators
using namespace std;
void main()
{
int seed = (int)time(0); // random seed
// choose one of the random number generators:
CRandomMersenne RanGen(seed); // make instance of random number generator
cout<<"\n\nRandom integers in interval from 0 to 99:\n";
for (int i = 0; i < 40; i++) {
int ir = RanGen.IRandom(0,99);
cout<<ir<<" ";
}
cout <<endl;
cout<<"\n\n\n\nRandom floating point numbers in interval from 0 to 1:\n";
for (int i = 0; i < 40; i++) {
float fr = RanGen.Random();
cout<<fr<<" ";
}
cout <<endl;
}
试验五：第三方"Mother-Of-All"随机数生成程序使用试验(程序来源：Agner Fog http://www.agner.org/random/)
// 使用说明:从网站下载压缩包，http://www.agner.org/random/randomc.zip
// 展开后,将其中的randomc.h头文件及mother.cpp文件Copy到项目文件夹,
// 并将它们加入到项目中，其中包括"Mother-Of-All" generator invented by George Marsaglia 的实现
#include <iostream>
#include <time.h>
#include "randomc.h" // define classes for random number generators
using namespace std;
void main()
{
int seed = (int)time(0); // random seed
// choose one of the random number generators:
CRandomMother RanGen(seed); // make instance of random number generator
cout<<"\n\nRandom integers in interval from 0 to 99:\n";
for (int i = 0; i < 40; i++) {
int ir = RanGen.IRandom(0,99);
cout<<ir<<" ";
}
cout <<endl;
cout<<"\n\n\n\nRandom floating point numbers in interval from 0 to 1:\n";
for (int i = 0; i < 40; i++) {
float fr = RanGen.Random();
cout<<fr<<" ";
}
cout <<endl;
}
试验六：第三方"SFMT"随机数生成程序使用试验(程序来源：Agner Fog http://www.agner.org/random/)
重要提示：在编译前，可以修改头文件sfmt.h中的#define MEXP以及下面相应的宏代码：Choose one of the possible Mersenne exponents. Higher values give longer cycle length and use more memory。SFMT利用了SSE2指令，速度最快，但只适合intel系列的部分芯片。
// 使用说明:从网站下载压缩包，http://www.agner.org/random/randomc.zip
// 展开后,将其中的头文件及sfmt.cpp文件Copy到项目文件夹,
// 并将它们加入到项目中，其中包括"SFMT" 的实现
#include <iostream>
#include <time.h>
#include "sfmt.h" // define classes for random number generators
using namespace std;
void main()
{
int seed = (int)time(0); // random seed
// choose one of the random number generators:
CRandomSFMT1 RanGen(seed); //注意可以是CRandomSFMT，是CRandomSFMT0，或CRandomSFMT1
cout<<"\n\nRandom integers in interval from 0 to 99:\n";
for (int i = 0; i < 40; i++) {
int ir = RanGen.IRandomX(0,99);
cout<<ir<<" ";
}
cout <<endl;
cout<<"\n\n\n\nRandom floating point numbers in interval from 0 to 1:\n";
for (int i = 0; i < 40; i++) {
float fr = RanGen.Random();
cout<<fr<<" ";
}
cout <<endl;
}
㈢ 高手进，c语言中如何得到服从正态分布的随机数
最好寻求专门算法
如果不行的话，也可以用大量随机数来模拟，譬如生成1000次0或1，然后求其平均数,可以得到很接近正态分布的.
如果有连续的随机函数，也可以直接求正态分部的积分函数，根据积分函数的反函数来确定位置
㈣ c语言正态分布随机数生成
int n,i,j,k,l,w,g;改为double n,i,j,k,l,w,g;
㈤ 如何用C语言生成[0.01,2]之间符合正态分布的随机数。。。注意是正态分布！！答案采用后再追加50分
# include <stdio.h>
# include <math.h>
# include <stdlib.h>
# include <time.h>
# define MAX_N 3000 /*这个值为N可以定义的最大长度*/
# define N 100 /*产生随机序列的点数，注意不要大于MAX_N*/
# define PI 3.141592653
void randn(double *x,int num)
{
double x1[MAX_N],x2[MAX_N];
int i;
srand((unsigned)time(NULL));
for(i=0;i<N;i++)
{
x1[i]=rand();
x2[i]=rand();
x1[i]=x1[i]/(RAND_MAX+1);
x2[i]=x2[i]/(RAND_MAX+1);
x[i]=sqrt(-2*log(x1[i]))*cos(x2[i]*2*PI);
}
}
void main()
{
double x[N],x_min,x_max;
int i;
FILE *fp;
if((fp=fopen("test.txt","w+"))==NULL)
{
fprintf(stderr,"Can't open the file ");
exit(1);
}
randn(x,N);
x_min=x[0];
x_max=x[0];
for(i=0;i<N;i++)
{
if(x[i]>x_max)
{
x_max=x[i];
}
if(x[i]<x_min)
{
x_min=x[i];
}
}
for(i=0;i<N;i++)
{
x[i]=(x[i]-x_min)/(x_max-x_min)*(2-0.01)+0.01;
}
for(i=0;i<N;i++)
{
printf("%f ",x[i]);
fprintf(fp,"%lf ",x[i]);
if(i%5==4)
{
printf(" ");
}
}
if(fclose(fp)==EOF)
{
printf("Closing error ");
}
}
把生成的数据放入txt文件中，再导入matlab中，查看是否符合正态分布。
得到H1=0，说明确实是正态分布。。。。
㈥ 如何在C语言中生成正态分布的随机数，要源代码~谢谢
随机生成一百个1至100的随机数。
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define N 100
int main(int argc, char *argv[])
{
int i;
int a[N];
srand(time(NULL));
for(i=0;i<N;i++)
a[i]=rand()%100+1;
printf("生成的随机数为:\n");
for(i=0;i<N;i++)
{
printf("%5d",a[i]);
if((i+1)%10==0)
printf("\n");
}
system("PAUSE");
return 0;
}
输出结果如下：
生成的随机数为:
41 15 82 1 23 51 16 96 92 17
86 71 87 69 74 5 50 18 42 52
46 34 52 18 40 74 79 35 22 36
65 94 80 91 18 72 61 79 4 11
61 30 95 55 11 19 38 87 78 52
95 30 99 53 99 99 10 79 70 33
91 85 10 99 47 58 93 41 19 71
56 60 10 24 73 87 18 38 13 73
57 22 91 4 37 60 67 58 85 48
46 7 57 100 73 96 60 44 24 23
请按任意键继续. . .
㈦ 求一个VC++程序，可以产生满足高斯分布的随机数！
float CMycomDlg::Gaussrand(float miu, float sigma)
{
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 (float)(miu+sigma*X);
}