2024年4月11日发(作者：)

四种函数的拟合

1正态分布

#include

//#include "iomanip.h"

#include "iomanip"

#include

#include

using namespace std;

const int A=2045;

const int C=1;

const int M=1048576;//2的20次方

double junyun(double a,double b,long int *seed)//采取的方法是混合同余法其中

A,C,M都是自定义的常数

{

double t;

*seed=A*(*seed)+C;

*seed=(*seed)%M;

t=(*seed)/(double)M;

t=a+(b-a)*t;

return t;

}

double gauss(double mean,double sigma,long int *s)

{

int i;

double x,y;

for(x=0,i=0;i<12;i++)

x+=junyun(0.0,1.0,s);