手动生成伪随机数C++源代码

最新推荐文章于 2024-07-25 00:15:00 发布

Chauncy__xu

最新推荐文章于 2024-07-25 00:15:00 发布

阅读量713

点赞数

分类专栏：智能算法文章标签：算法 c++

本文链接：https://blog.csdn.net/qq_44384577/article/details/104440796

版权

智能算法专栏收录该内容

16 篇文章 10 订阅

订阅专栏

乘同余法产生(0,1)均匀分布的伪随机数

乘同余法的计算公式为：
$x_{i}=A x_{j-1}(\operatorname{Mod} M)$
或者
$x_{i}=A x_{i-1}-\left[A x_{i-1} / M\right] M$

 /*date:2020年2月22日
   Author:Chauncy_xu*/
 #include<iostream>
 #include<math.h>
  using namespace std;
  int main()
  {
      int n;
      cout<<"伪随机数个数："<<endl;
      cin>>n;
      int number[n];
      double pse_number[n];
      number[0]=1;
      int A=3,M=16;
      for(int i=1;i<n;i++)
      {
          number[i]=(A*number[i-1])%M;
      }
      for(int i=0;i<n;i++)
      {
           pse_number[i]=float(number[i])/float(M);
      }
      for(int i=0;i<n;i++)
      {
          cout<<pse_number[i]<<endl;
      }
      return 0;
  }

逆变法生成负指数分布的伪随机数

逆变法

密度函数： $f(x)=\left\{\begin{array}{l} {1,0 \leq x \leq 1} \\ {0, 其它} \end{array}\right.$
分布函数：
$F(x)=P_{x}\{X \leq x\}=\int_{-\infty}^{x} f(x) d x=x$

逆变法的步骤：

已知F(x)，或由f(x)求F(x)
即 $F(x)=\int_{-\infty}^{x} f(x) d x$ 令 $y = F (x)$
推导 $x=F^{-1}(y)$
产生 $y\in U(0,1)$
用 $x=F^{-1}(y)$ 得到 $x_{0},x_{1},x_{2},\cdots,x_{i-1},x_{i}$

逆变法产生负指数分布的伪随机数

首先产生 $\in U(0,1)$ ，然后用 $X=-\frac{1}{\lambda} \mathrm{ln} U$

 /*date:2020年2月22日
   Author:Chauncy_xu*/
 #include<iostream>
 #include<math.h>
  using namespace std;
  int main()
  {
      int n;
      cout<<"伪随机数个数："<<endl;
      cin>>n;
      int number[n];
      double pse_number[n];
      double com_index[n];
      number[0]=1;
      int A=3,M=16,lambda;
      cout<<"输入负指数分布的参数："<<endl;
      cin>>lambda;
      for(int i=1;i<n;i++)
      {
          number[i]=(A*number[i-1])%M;
      }
      for(int i=0;i<n;i++)
      {
           pse_number[i]=float(number[i])/float(M);
           com_index[i]=(-1)*log(pse_number[i])/lambda;
      }
      for(int i=0;i<n;i++)
      {
          cout<<com_index[i]<<endl;
      }
      return 0;
  }