产生瑞利分布的随机数 C语言实现

瑞利分布的概率密度函数为

f(x)=\frac{x}{\sigma ^2}e^{-x^2/2\sigma ^2}x>0

瑞丽分布的均值为\sigma \sqrt{\frac{\pi}{2}},方差为(2-\frac{\pi}{2})\sigma^2

首先使用逆变换法产生参数\beta=2的指数分布的随机变量y,其概率密度函数为f(y)=\frac{1}{2}e^{-\frac{y}{2}}

然后通过变换x=\sigma\sqrt{y},产生瑞丽分布的随机变量x,具体的方法如下:

(1)产生均匀分布的随机数u

(2)计算y=-2ln(u)

(3)计算x=\sigma \sqrt{y}

1、头文件 rayleigh.h

#pragma once
#ifndef RAYLEIGH_H_
#define RAYLEIGH_H_


/*
函数功能:	瑞利分布的随机数
输入参数说明:
sigma		瑞利分布的均值
seed	    长整型指针变量, *seed 为伪随机数的种子
*/

double rayleigh_data(double sigma, long int * seed);

#endif // !RAYLEIGH_H

2、rayleigh.c

#include "rayleigh.h"

#include <stdint.h>
#include "string.h"
#include "stdio.h"
#include "uniform.h"
#include <math.h>

double rayleigh_data(double sigma, long int * seed)
{
	double u, x;
	u = uniform_data(0.0, 1.0, seed);
	x = -2.0 * log(u);
	x = sigma + sqrt(x);
	return x;
}

3、主函数 main.c

#include <stdio.h>
#include <string.h>
#include "uniform.h"
#include "rayleigh.h"

int main()
{
	// 产生50个瑞利分布的随机数
	int i, j;
	long int s;

	double x, sigma;
	sigma = 1.0;  s = 13579;

	for (i = 0; i < 10; i++)
	{
		for (j = 0; j < 5; j++)
		{
			x = rayleigh_data(sigma,  &s);
			printf("%13.7f",x);
		}
		printf("\n");
	}	
	getchar();		// 此行代码是为了保持输出窗口,按任意按键关闭弹窗
	return 0;
}

4、调试输出结果

 

  • 0
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

jk99528878769

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值