C语言数字信号处理——1.1 均匀分布的随机数

-Aerolite-

已于 2023-04-08 10:13:47 修改

阅读量380

点赞数

分类专栏： C/C++ 数字信号处理文章标签：信号处理 c语言

于 2023-03-15 14:59:52 首次发布

本文链接：https://blog.csdn.net/weixin_43348948/article/details/129554974

版权

C/C++ 数字信号处理专栏收录该内容

1 篇文章 0 订阅

订阅专栏

该文详细介绍了如何在C语言中生成(a,b)区间上均匀分布的随机数，包括数学原理、混合同余法的具体步骤以及程序代码示例。在混合同余法中，使用了特定的参数a=2045,c=1,M=2^20，并通过变换公式将(0,1)区间上的随机数映射到目标区间。

摘要由CSDN通过智能技术生成

一、功能

产生 $(a, b)$ 区间上均匀分布的随机数

二、原理

均匀分布的概率密度函数为：
$f(x)=\begin{cases} \frac{1}{b-a}& , \mathrm{a}\leqslant x\leqslant b\\ 0& , othercase\\ \end{cases}$
通常用 $U (a, b)$ 表示。均匀分布的均值为 $\frac{a+b}{2}$ ，方差为 $\frac{\left( a-b \right) ^2}{12}$ 。
产生均匀分布随机数的方法如下：
首先，由给定的初值 $x_0$ ，用混合同余法，
$\left\{ \begin{aligned} &x_i=\left( ax_{i-1}+c \right) \left( \mathrm{mod} M \right)\\ &y_i=x_i/M\\ \end{aligned} \right.$
产生 $(0, 1)$ 区间上的随机数 $y_i$ 。其中： $a = 2045$ ， $c = 1$ ， $M=2^{20}$ ；然后通过变换 $z_i=a+(b-a)y_i$ 产生 $(a, b)$ 区间上的随机数 $z_i$ 。

三、程序

3.1 dsp.h头文件函数声明

#pragma once
#ifndef __HEAD_H__		
#define __HEAD_H__

#include <math.h>
#include <string.h>

double uniform(double a, double b, long int* seed);  // 均匀分布的随机数

3.2 dsp.c源文件函数定义

#include <math.h>
#include <string.h>
#include "dsp.h"

double uniform(double a, double b, long int* seed) {
	double t;
	*seed = 2045 * (*seed) + 1;
	*seed = *seed - (*seed / 1048576) * 1048576;
	t = (*seed) / 1048576.0;
	t = a + (b - a) * t;
	return t;
}

3.3 main.c中函数调用

#include <stdio.h>
#include "dsp.h"

int main()
{
	double uniform(double a, double b, long int* seed);
	double a = 0.0, b = 1.0, x;
	long int s = 13579;
	for (int i = 0; i < 10; i++){
		for (int j = 0; j < 5; j++) {
			x = uniform(a, b, &s);
			printf("%13.7f", x);
		}
		printf("\n");
	}
}