华科电信随机过程实验报告.docxVIP

随机过程实验报告院（系）：电子信息与通信学院专业班级：姓名：学号：指导教师：一．实验任务(1)利用线性同余法产生平均分布的随机数。(2)利用平均分布的随机数，根据中心极限定理产生正态分布的随机数。(3)利用平均分布的随机数，根据反函数产生指数分布的随机数。(4)利用平均分布的随机数，产生泊松分布的随机数。(5)计算任意分布的随机过程的均值。(6)计算泊松过程的自相关序列。二．实验环境(1)平台：Windows XP/7.0(2)编程环境：VC6.0(3)编程语言：C++三．实验代码// MyRand.cpp: implementation of the CMyRand class.////////////////////////////////////////////////////////////////////////#include stdafx.h#include random.h#include MyRand.h#include math.h#ifdef _DEBUG#undef THIS_FILEstatic char THIS_FILE[]=__FILE__;#define new DEBUG_NEW#endif//////////////////////////////////////////////////////////////////////// Construction/Destruction//////////////////////////////////////////////////////////////////////CMyRand::CMyRand(){}CMyRand::~CMyRand(){}void CMyRand::MyRandInit(void){N = 0x7FFFFFFF; //2^31-1K = 16807;//7^5seed = 2;}/*函数功能，采用线性同余法，根据输入的种子数产生一个伪随机数，如果种子不变， 则将可以重复调用产生一个伪随机序列 利用CMyRand类中定义的全局变量：S, K, N, Y。 其中K和N为算法参数，S用于保存种子数，Y为产生的随机数*/unsigned int CMyRand::MyRand(unsigned int seed){//添加伪随机数产生代码if(S == seed){Y = (Y * K) % N;}else{S = seed;Y = (seed * K) % N;}return Y;}/*函数功能，产生一个在min~max范围内精度为4位小数的平均分布的随机数*/double CMyRand::AverageRandom(double min,double max){double dResult;dResult = 0;dResult = (double)MyRand(seed)/N;//产生0-1之间的均匀分布dResult = dResult*(max-min)+min;//将0-1之间的均匀分布扩展到min-max之间的均匀分布return dResult;}/*函数功能，在min 到max 范围内产生正态分布的随机数miu,最大概率密度处的随机变量，即产生的随机数中,概率最大的那个sigma*/double CMyRand::NormalRandom(double miu, double sigma, double min, double max){double dResult;dResult = 0;//添加正态分布随机变量产生代码for(int i=0;i12;i++)dResult += AverageRandom(min,max);//产生12个伪随机数，将它们加起来，dResult = (dResult-6)/(max-min); //再减去6，就能近似得到标准正态变量的样本值dResult = dResult*sigma + miu; return dResult;}/*函数功能，产生指数分布的随机数*/double CMyRand::ExpRandom(double lambda, double min, double max){double dResult = 0.0;//添加指数分布随机变量产生代码dResult = -1.0*log(AverageRandom(min,max))/lambda;return dResult;}/*函数功能，产生泊松分布的随机数*/unsigned int CMyRand::PoisonRandom(double lambda, double min, double max){unsigned int dResult;//添加泊松分布随机变量产生代码double