码迷,mamicode.com
首页 > 其他好文 > 详细

瑞利分布的随机数

时间:2019-10-14 21:03:17      阅读:280      评论:0      收藏:0      [点我收藏+]

标签:概率   sig   随机   spl   函数   --   使用说明   blog   tps   

一、功能

产生瑞利分布的随机数。

二、方法简介

瑞利分布的概率密度函数为
\[ f(x) = \frac{x}{\sigma ^{2} }e^{-x^{2}/2\sigma ^{2}} \ x > 0 \]
瑞利分布的均值为\(\sigma \sqrt{\frac{\pi }{2}}\),方差为\(\left ( 2 - \frac{\pi }{2} \right )\sigma ^{2}\)

首先用逆变换法产生参数\(\beta = 2\)的指数分布的随机变量\(y\),其概率密度函数为\(f(y) = \frac{1}{2} e^{-\frac{y}{2}}\);然后通过变换\(x = \sigma \sqrt{y}\),产生瑞利分布的随机变量\(x\),具体方法如下:

  1. 产生均匀分布的随机数\(u\),即\(u \sim U(0,1)\)
  2. 计算\(y = - 2 \ ln(u)\)
  3. 计算\(x = \sigma \sqrt{y}\)

三、使用说明

是用C语言实现产生瑞利分布随机数的方法如下:

/************************************
    sigma   ---瑞利分布的参数sigma
    seed    ---随机数种子
************************************/
#include "math.h"
#include "uniform.c"

double rayleigh(double sigma, long int *s)
{
    u = uniform(0.0, 1.0, s);
    x = -2.0 * log(u);
    x = sigma * sqrt(x);
    return(x);
}

uniform.c文件参见均匀分布的随机数

瑞利分布的随机数

标签:概率   sig   随机   spl   函数   --   使用说明   blog   tps   

原文地址:https://www.cnblogs.com/liam-ji/p/11631644.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!