系統仿真与虚拟现实实验报告_3010218158_杨静.docVIP

系统仿真与虚拟现实 实 验 报 告 学院名称 软件学院 专 业 软件工程 学生姓名 杨静 学 号 3010218158 年 班 级 2010级四班 2013年 4月 22日 目　　录 一、随机数生成算法 1、实验内容 4 2、实验原理 4 3、实现方法 5 二、VRML建模试验 1、实验内容 13 2、实验原理 13 3、实现方法 13 三、Matlab基本实验：MATLAB基本操作 1、实验内容 17 2、实验原理 17 3、实现方法 18 四、Matlab进阶实验：单窗口服务模型模拟 26 1、实验内容 26 2、实验原理 26 4、实现方法 26 5、实验结果与结果分析 28 三、总结与展望 28 一、随机数生成算法 1、实验内容 参考教材《系统仿真导论》DEDS01.c程序 利用该程序对如下PMMLCG进行2均匀性检验及独立性检验（=0.05） i-1(mod231-1) 样本点 n = 10000 K = 10 (均匀性) N = 1000 j = 1 （独立性） 修改该程序中的随机数发生器程序，实现对如下组合发生器进行x2均匀性检验及独立性检验（α = 0.05） Zi = 25173Zi-1 + 13849 (mod 65536) Z0 = 1.0 Zi2 = Zi-12 + Zi-22 (mod 128) Z02 = 1.0 Z-12 = 1.0 样本点数：n = 5000 K = 10 （均匀性） N = 1000 j = 1 （独立性） 2、实验原理 均匀性检验：要求由随机数发生器产生的伪随机数Ui能均匀的分布在[0.1]区间内。 最常用的方法是频率检验：将随机数发生器的取值范围[0,1]分成K个互不重叠的等长子区间，由该随机数发生器产生N个随机数Ui（i=1,2，…,N）按照均匀性的要求，随机数在每一个子区间上的概率P=1/K，即在每一个子区间上的随机数的个数理论值为n=N/K。实际上Ui在每一个子区间上的个数不会恰好等于n，而为nj（j=1,2,…,K）这样就会有偏差。 频率检验就是检验实际频率与理论频率之间的偏差大小，一般用χ2检验 ： 显然，若nj = n则χ2等于0，实际频率与理论频率一致，χ2大小反映了偏差程度，也就是随机数的均匀程度。 独立性检验χ2检验 1：原假设H0：给定随机数发生器产生的Ui是独立同分布U（0,1）的随机变量； 2：将[0,1]分成K个等长子区间； 3：由该随机数发生器产生N个随机数Ui； 4：统计计算在每个子区间上随机数的个数nj； 5：计算χ2： 独立性检验： 1：原假设H0：给定随机数发生器产生的Ui是独立同分布U（0,1）的随机变量，即ρ=0； 2：由该随机数发生器产生N个随机数Ui，并计算ρj的均值； 3：若N—j充分大（一般要50），取统计量： 渐进服从标准正态分布N（0,1）； 4：给定检验水平α，记Z1-α为N（0,1）的上1-α的临界点，则当 组合发生器： 1：由Z1产生128个随机数U1，形成数组V=（U11，…，U1128） 2：由Z2产生在[1,128]上分布的随机整数I，返回UI作为组合发生器的第一个随机数； 3：由Z1产生一个随机数取代UI； 4：由Z2生成新的I，依次下去； 4、实现方法 第一部分： /*====================================== The program is used to test the uniform and the dependence of a random generator. */ #include stdio.h #include stdlib.h #include math.h int SUB1(); double SUB2(); int nvalue; double probd[25]; long int a1,c1,m1,x1; int main(void) { int i,m; float mn; int ntypes,nx[20],number,number1; long int z0,a0,c0,m0,x0;