c语言产生各种分布函数

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

float Expon(float fl)// 产生负指数分布的随机数
{
float v,u;
u=Drand();
v=-fl*log(u);
return v;
}

float poissn(float la)//poissn（）产生泊松分布随机数
{
int k=0;
float b,t=1.0f,r;
b=exp(-la);
while((t-b)>=0)
{
r=Drand();
t=t*r;
k=k+1;
}
return k;
}

float pgauss(float k,float j)//pgauss（）产生正态分布随机数
{
float v1,v2,s,w,y,sg[2];
do{v1=Drand();
v2=Drand();
s=v1*v1+v2*v2;
}
while(s>=1);
w=sqrt(-2.0*log(s)/s);
sg[0]=v1*w;
sg[1]=v2*w;
y=k*(sg[0]+sg[1])/2+j;
return y;
}

double Drand()//drand（）产生0～1之间的随机数
{
float x;
int i;
for(i=0;i<20;i++)rand();
x=rand();
x=65539*x+1743251541;
x=fmod(x,2147483638);
return(x/2147483638);
}