uniform.h文件
#ifndef UNIFORM_H_
#define UNIFORM_H_
/*
函数功能: 产生(a,b)区间上均匀分布的随机数组
输入参数说明:
a 给定区间的下限
b 给定区间的上线
seed 长整型指针变量, *seed 为伪随机数的种子
*/
double uniform_data(double a, double b, long int * seed);
#endif
?
子函数程序uniform.c? ?函数功能:产生(a,b)区间上均匀分布的随机数
#include "uniform.h"
#include <stdint.h>
#include "string.h"
#include "stdio.h"
/*
函数功能: 产生(a,b)区间上均匀分布的随机数
输入参数说明:
a 给定区间的下限
b 给定区间的上线
seed 长整型指针变量, *seed 为伪随机数的种子
*/
double uniform_data(double a, double b,long int * seed)
{
double t;
*seed = 2045.0 * (*seed) + 1;
*seed = *seed - (*seed / 1048576) * 1048576;
t = (*seed) / 1048576.0;
t = a + (b - a) * t;
return t;
}
主函数代码入下
#include <stdio.h>
#include <string.h>
#include "uniform.h"
int main()
{
// 产生50个0-1之间的随机数
double a, b ,x;
int i, j;
long int s;
a = 0.0;
b = 1.0;
s = 13579;
for (i = 0; i < 10; i++)
{
for (j = 0; j < 5; j++)
{
x = uniform_data(a, b, &s);
printf("%13.7f",x);
}
printf("\n");
}
getchar(); // 此行代码是为了保持输出窗口,按任意按键关闭弹窗
return 0;
}
调试输出结果
|