随机变量生成算法——产生离散型随机变量样本值的方法

随机变量生成算法——产生离散型随机变量样本值的方法

算法介绍
设离散型随机变量X具有分布律
现在来产生X的随机数：
先产生伪随机数u，令

由于

所以X具有给定的分布律。

下面将引用浙大版概统（第四版）中的例题为例

1、计算过程

2、C++代码实现

#include
int main()
{
     
    float x[10001],y[10001],z[10001],u[10001];
    x[0]=1;
    y[0]=2;
    z[0]=3;
    for(int i=0;i<10000;i++)
    {
     
        x[i+1]=fmod(171*x[i],30269);
        y[i+1]=fmod(170*y[i],30307);
        z[i+1]=fmod(172*z[i],30323);
        u[i+1]=fmod(x[i+1]/30269+y[i+1]/30307+z[i+1]/30323,1);
    }
    FILE *fp;
    fp=fopen("/*输入自定义路径*/","w");
    for(int i=1;i<=10000;i++)
    {
     
        if(u[i]<0.2)
            fprintf(fp,"1 ");
        else if(u[i]<0.35)
            fprintf(fp,"2 ");
        else if(u[i]<0.6)
            fprintf(fp,"3 ");
        else
            fprintf(fp,"4 ");
    }//模拟算法部分
    fclose(fp);
    return 0;
}

3、使用matlab画直方图

u=[…(从C++程序生成文件中导入数据)]
histogram(u)

4、生成直方图预览