随机变量生成算法——通过Rayleigh变量(或均匀分布变量)转换为高斯随机变量
算法介绍
由随机过程的理论知两个相互独立的高斯随机变量的均方根是Rayleigh分布的随机变量。
首先,由高斯随机变量的定义可得其概率密度为:
此时,令
即进行如下所示的的坐标变换(即极坐标变换):
由该坐标变换可将原概率密度变换为:
若设R是Rayleigh随机变量,Θ是(0,2pi)均匀分布随机变量。
那么由如下关系式:
最后可得两个相互独立的高斯随机变量的表达式为:
此时称X为同相分量,Y为正交分量。
~~~~ ~~~~ ~~~~ ~~~~ ~~~~ ~~~~ ~~~~ ~~~~ ~~~~ ~~~~ ~~~~
下面将以生成方差为1的高斯变量(即标准正态分布)为例
1、计算过程
由上述公式可得高斯随机变量由两个不相关的均匀分布随机变量的表达式:
或者转化为Rayleigh分布表示形式:
其中, u 1 , u 2 u_1,u_2 u1,u2要求是不相关的均匀分布随机变量, r = − 2 σ 2 l n ( u 1 ) r=\sqrt{-2\sigma^2ln(u_1)} r=−2σ2ln(u1)。
关于如何产生均匀分布随机变量以及如何使用均匀分布随机变量产生Rayleigh分布随机变量,诸位可以参照我之前的几篇文章,这里不做赘述。
随机变量生成算法入门——Wichmann-Hill算法
随机变量生成算法——逆变换法
2、C++代码实现
#include<bits/stdc++.h>
#define PI 3.1415926535897932
{
float re[10001],g[10001];
for(int i=0;i<10000;i++)
{
/*这里代码打不动了,均匀分布直接用随机种子发生器了*/
re[i+1]=sqrt(-2*log(rand()/(RAND_MAX+1.0)));
g[i+1]=re[i+1]*sin(2*PI*(rand()/(RAND_MAX+1.0)));//正交分量
/*g[i+1]=re[i+1]*cos(2*PI*(rand()/(RAND_MAX+1.0)));//同相分量*/
}
FILE *fp;
fp=fopen("/*输入自定义路径*/","w");
for(int i=1;i<=10000;i++)
{
fprintf(fp,"%f ",g[i]);
}
fclose(fp);
return 0;
}
3、使用matlab画直方图
u=[…(从C++程序生成文件中导入数据)]
histogram(u)
4、生成直方图预览