随机过程试验上机报告华科

1、随机模拟与实验上机实验报告 通信0809 潇洒哥 学号 一. 随机模拟的典型步骤： 1、根据问题构建模拟系统 2、仿真系统中各种分布的随机变量3、运行模拟系统，进行统计测量 4、分析数据，输出结果二主要工具基本工具：c、c+等编程模拟、matlab网络模拟：opnet modeler、ns2：大型网络仿真 cassap：数字信号处理；spw：电子系统三编译代码/ myrand.cpp: implementation of the cmyrand class./#include "stdafx.h"#include "random.h"#include &

2、quot;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; /231-1k = 16807; /75seed = 2;/*函数功能，采用线性同余法，根据输入的种子数产生一个伪随机数，如果

3、种子不变， 则将可以重复调用产生一个伪随机序列实现思路：利用cmyrand类中定义的全局变量：s, k, n, y。其中k和n为算法参数，s用于保存种子数，y为产生的随机数，第一次调用检查将seed赋值与s获得y的初值，之后调用选择rand（）函数赋值与y。*/unsigned int cmyrand:myrand(unsigned int seed)if(s != seed)s = seed;y = (seed * k) % n;elsey = (y * k) % n;if(y = 0)y = rand();return y;/*函数功能，产生一个在minmax范围内精度为4位小数的平均分布

4、的随机数实现思路：min和max代表该均匀随机变量的取值范围，精确到小数点后4位。该参数在randomdlg.cpp中onaverage ()和onaveragestat()函数中已经设置，因此只需将该数小数点右移四位强制取整即可获得整数范围内的随机数，采用对（max-min）求余再加上min并将获得的结果缩小10000倍即可。*/double cmyrand:averagerandom(double min,double max)int mininteger = (int)(min*10000);int maxinteger = (int)(max*10000);int randintege

5、r = myrand(seed);int diffinteger = maxinteger - mininteger;int resultinteger = randinteger % diffinteger + mininteger;return resultinteger/10000.0;/*函数功能，在min 到max 范围内产生正态分布的随机数miu,最大概率密度处的随机变量，即产生的随机数中,概率最大的那个sigma实现思路：设由averagerandom函数获得01间隔均匀分布随机数u(0，1)，i=1,2,n，且相互独立，由中心极限定理可知，当n较大时*/double cmyra

6、nd:normalrandom(double miu, double sigma, double min, double max)double dresult;dresult = 0;for(int i=0;i<12;i+)dresult += averagerandom(min,max);dresult = (dresult-6)/(max-min);dresult = dresult*sigma + miu;return dresult;/*函数功能，产生指数分布的随机数实现思路： 1）用averagerandom产生均匀分布随机数ui； 2）计算指数分布随机数：xi=-ln ui

7、/*/double cmyrand:exprandom(double lambda, double min, double max)double dresult = 0.0;while(dresult < 0.01)dresult = averagerandom(min,max);dresult = -1.0*log(dresult) / lambda;return dresult;/*函数功能，产生泊松分布的随机数实现方法：利用递归法*/unsigned int cmyrand:poisonrandom(double lambda, double min, double max)uns

8、igned int dresult;double u = averagerandom(min,max);int i = 0;double p = exp(-1*lambda);double f = p;while(u >= f)p = lambda*p/(i+1);f += p;i +;dresult = i;return dresult;/*函数功能，计算任意分布的随机过程的均值实现过程：用for语句调用任意函数加和循环再求平均*/double cmyrand:ex(void)int m;double ex=0.0;/添加均值计算代码for (m=0; m<1000; m+)ex

9、 += normalrandom(0,1.2, 0,1);return ex/1000;/*函数功能，计算随机过程的自相关序列实现过程：由平稳随机过程数字特征求解的相关原理可得rx(m)=i2e-2|m|; */double* cmyrand:rx(double lambda, int points)int m,i=5;double *rx = (double*)malloc(2*points+1)*sizeof(double);/添加自相关序列产生代码/产生的自相关序列存入rx中，rx可当作数组使用/不要在本函数中释放该数组!for (m=-points; m<=points; m+)rx(m+points) = i * i * exp(-2 * lambda * abs(m);return rx;4 执行截图均匀分布均匀统计正态分布正态统计指数分布指数统计泊松分布泊松统计均值（此时取normalrandom(0,1.2, 0,1)）自相关(此时取int m,i=5;)5 实验心得1. 该课程设计所设计的各种模块和原理与随机过程的理论知识紧密结合。所以这样与课本知识关联紧的课程实际比较容易完成，同时不仅加强了对书本基础知识的理解，而且实践环节得到了很好的锻炼，能够激发实践者对随机的兴趣甚至能够自己动