uniform.h
#pragma once
class uniform
{
private:
double a, b, generate_num;
int * seed;
int s;
int M, N, i, j;
public:
uniform()
{
M = 1048576;
N = 2045;
}
void generate();
double random_number(double, double, int *);
};
double uniform::random_number(double a, double b, int * seed)
{
(*seed) = N * (*seed) + 1;
(*seed) = (*seed) - ((*seed) / M) * M;
generate_num = static_cast<double>((*seed)) / M;
generate_num = a + (b - a) * generate_num;
return (generate_num);
}
gauss.h
#pragma once
#include "uniform.h"
class gauss
{
private:
double mean, sigma, x, y, generate_num;
int s;
int * seed;
int i, j, m;
public:
gauss() {}
void generate();
double random_number(double, double, int *);
};
double gauss::random_number(double mean, double sigma, int * seed)
{
uniform unif_num;
x = 0;
for (m = 0; m < 12; m++)
{
x += unif_num.random_number(0.0, 1.0, seed);
}
x = x - 6.0;
y = mean + x * sigma;
return (y);
}
lognorm.h
#pragma once
#include <math.h>
#include "gauss.h"
#include "uniform.h"
class lognorm
{
private:
double sigma, u, x, y, generate_num;
int * seed;
int s, i, j;
public:
lognorm() {}
void generate();
double random_number(double, double, int *);
};
double lognorm::random_number(double u, double sigma, int * seed)
{
gauss gau_num;
y = gau_num.random_number(u, sigma, seed);
x = exp(y);
return (x);
}
lognorm.cpp
#include <iostream>
#include <iomanip>
#include "lognorm.h"
using namespace std;
void main()
{
lognorm solution;
solution.generate();
}
void lognorm::generate()
{
cout << "输入对数正态分布的均值:";
cin >> u;
cout << "输入对数正态分布的方差:";
cin >> sigma;
cout << "输入随机数的种子:";
cin >> s;
cout << "生成的随机数结果为:" << endl;
for(i = 0; i < 10; i++)
{
for(j = 0; j < 5; j++)
{
generate_num = random_number(u, sigma, &s);
cout << setw(10) << generate_num;
}
cout << endl;
}
}