《随机模拟与实验上机报告.docxVIP

随机模拟与实验上机报告一随机模拟的典型步骤1、根据问题构建模拟系统2、仿真系统中各种分布的随机变量3、运行模拟系统，进行统计测量4、分析数据，输出结果二随机模拟的主要工具基本工具：C、C++等编程模拟、matlab 网络模拟：OPNET Modeler、NS2：大型网络仿真CASSAP：数字信号处理；SPW：电子系统三实验原理伪随机数产生函数：unsigned int CMyRand::MyRand(unsigned int seed)实现思路：利用线性同余法，设置y0，即设置种子yn=kyn-1(mod N)，un=yn/N即（ran0）N=231-1，k=75，周期≈2×109平均分布随机数产生函数：double CMyRand::AverageRandom (double min,double max)实现思路：C语言中的rand函数可以产生0－1间隔均匀分布随机数，要产生min~max范围内精度为四位小数的随机数，则只需将该数小数点右移四位强制取整即可获得整数范围内的随机数，采用对（max-min）求余再加上min并将获得的结果缩小10000倍即可。正态分布随机数产生函数double CMyRand::NormalRandom(double miu, double sigma, double min, double max)实现思路：利用中心极限定理，通过多个独立的均匀随机变量之和构成正态随机变量,产生均值为miu，方差为sigma的正态分布的随机数x设Ui～U(0，1)，i=1,2,…,n，且相互独立，由中心极限定理可知，当n较大时实验中取n=15,E==0.5,D=再用变换公式x=Z*sigma+miu可以得到正态分布随机数指数分布随机数产生函数double CMyRand::ExpRandom(double lambda, double min, double max)实现思路：利用逆变换法产生参数为λ的指数分布，其指数分布函数为其反函数为因此，X的模拟方法为，先产生均匀分布随机数{ui}；再计算指数分布随机数：泊松分布随机变数产生函数unsigned int CMyRand::PoisonRandom (double lambda, double min, double max)实现思路：由于，所以采用递推法构成泊松分布产生[0,1]均匀分布随机数u若uF,令X=I,停止转向3)计算随机过程的自相关序列double* CMyRand::Rx(double lambda, int points)实现思路：平稳随机过程数字特征求解的相关原理RX(m)=I2e-2λ|m|; KX(m)= RX(m)-m2X(1)E(X(n))= I*P{X(n)=+I}+(-I)*P{X(n)=-I}=0(2)当m0时四代码编译// 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^5srand((int)time(0));seed=3;//rand();//seed = 2;seed++;}/*函数功能，采用线性同余法，根据输入的种子数产生一个伪随机数，如果种子不变，则将可以重复调用产生一个伪随机序列。利用CMyRand类中定义的全局变量：S, K, N, Y。其中K和N为算法参数，S用于保存种子数，Y为产生的随机数。*/unsigned int CMyRand::MyRand(unsigned int seed){if(S=seed){Y=(Y*K)%N;}else{ seed+