产生瑞利分布的随机数 C语言实现_对偶变量mc生rayleigh分布随机数-程序员宅基地

技术标签： c语言 C和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、调试输出结果

本文链接：https://blog.csdn.net/jk99528878769/article/details/124524666

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

三、购物车案例_17(笔记)_购物车案例请求数据地址-程序员宅基地

文章浏览阅读332次。购物车_购物车案例请求数据地址

在unity中DOTween 使用方法_tween number-程序员宅基地

文章浏览阅读603次。DOTween 使用方法_tween number

DSP之TMS320F28335学习总结与笔记（二）————ADC模块_dsp2833x_adc.c-程序员宅基地

文章浏览阅读1.2w次，点赞13次，收藏138次。F28335 ADC模块ADC转换模块A/D转换器（ADC）将模拟量转换为数字量通常要经过四个步骤：采样、保持、量化和编码。采样：将一个时间上连续变化的模拟量转化为时间上离散变化的模拟量。保持：将采样结果存储起来，直到下次采样，这个过程称作保持。一般，采样器和保持电路一起总称为采样保持电路。量化：将采样电平归化为与之接近的离散数字电平，这个过程称作量化。ADC关键指标分辨率：指数字量变化一个最小量时模拟信号的变化量，定义为满刻度与2^n的比值。分辨率又称为精度，通常以数字信号的位_dsp2833x_adc.c

V-rep学习笔记：机器人路径规划2-程序员宅基地

文章浏览阅读1.8k次。　　路径规划问题是机器人学研究的一个重要领域，它是指给定操作环境以及起始和目标的位置姿态，要求选择一条从起始点到目标点的路径，使运动物体(移动机器人或机械臂)能安全、无碰撞地通过所有的障碍物而达到目标位置。路径规划从研究对象上可分为关节式机械臂和移动机器人。一般来讲前者具有更多的自由度，而后者的作业范围要更大一些，这两类对象具有不同的特点，因此在研究方法上略有不同。在V-rep学习笔记：机器人路..._collision pairs

问题：U8的存货同步到本地时，报“数据库可能存在相同的编码”_数据库可能已经存在相同的编码等-程序员宅基地

文章浏览阅读1.5k次。1.环境：U8与本地某程序。跨服务器、跨数据库。2.问题：U8的存货同步到本地时，报“数据库可能存在相同的编码”。3.原因：在确认了U8程序正常、本地程序正常、U8数据库正常、本地数据库正常、链接服务器正常、触发器本身正常后，回到了问题的本身：为何原来的同步执行正常，但是现在的存货保存却异常，怀疑是最新的存货有问题。所以先停了触发器，然后用U8做了一个存货，正常保存后。启用触发器，直接用刚_数据库可能已经存在相同的编码等

EOS与ESD的区别_eos和esd的区别-程序员宅基地

文章浏览阅读7.9k次。摘要：什么是EOS？EOS与ESD的区别是什么？EOS为ELectrical Over Stress的缩写，指所有的过度电性应力。当外界电流或电压超过器件的最大规范条件时，器件性能会减弱甚至损坏。一、什么是EOS？EOS为ELectrical Over Stress的缩写，指所有的过度电性应力。当外界电流或电压超过器件的最大规范条件时，器件性能会减弱甚至损坏。EOS通常产生于：1.电源(AC/DC) 干扰、电源杂讯和过电压。2.由于测试程序切换（热切换）导致的瞬变电流..._eos和esd的区别