实用统计分析实验.docVIP

实验一随机数的生成 1?问题背景 多次重复地抛掷一枚匀质的硬币是一个古老而现实的实验问题，通过分析 “正面向上”出现的频率，我们可以从屮得到许多结论，但要做这个简单而重复 的实验，很多人没有多余的时间或耐心来完成它，现在借助于计算机的帮助，人 人都nJ以在很短的吋间内完成它。因此 借助于计算机进行模拟随机试验，产生 服从各类分布的随机数，通过数据处理和分析，我们可以从屮发现许多有用的规 律，或祥来验证我们理论推导的结论是否正确，木实验的主要H的是产生服从某 种分布的随机数。 实验目的与要求 熟悉常见分布的随机数产生的有关命令； 拿握随机模拟的方法； 提高读者观察实验现象或处理数据方面的能力。 二、实验操作过程 随机数生成的基本原理 生成服从给定分布的随机数，需要首先生成服从均匀分布的随机数，常用的 生成均匀分布随机数的方法是同余法，其递推公式 X] = (ax} + c) mod m 给定初值九0可以迭代出均匀随机数坷，兀2,…，心，将它们进行标准化(此 吋随机数界于0和1之间)或极差标准化(此时随机数界于J和1之间)，可以 得到均匀分布的随机数。 获得均匀分布的随机数之后，可以用多种方法构造基于该随机数的随机变 量，常用的方法是反函数法，即利用随机变量兀的分布函数F(Q的反函数厂「⑴ 来推求随机变量，基本算法是： 产生均匀分布随机数厂 令X,. = F-,(r/)，然后返回。 下而结合正态分布随机变量的生成进行具体介绍: 正态分彳IJ的分彳IJ函数为 2右匚旳[冷(宁)皿 式小“为期望，/为方差，由屮心极限定理，有^^=-(y b n /=1 2 当n=l2吋，可达到较好精度，故x = (J/;.-6)7 + /z, /=1 尢就是基于均匀分布随机数◎的服从正态分布的随机数。 1 ?二项分布的随机数的产生 基本数学原理：设X服从参数为77丿的二项分布。 在MATLAB屮用函数binornd产生参数为n,p的二项分布的随机数，其基本 的调用格式如下： -R=binornd (N, P) %N, P为二项分布的两个参数，生成服从参数为 N, P的二项分布的一个随机数； ? R=binornd (N, P, m,n) %m,n分别表示随机数产生的行数和列数。 例1T产生参数为10,概率为0.5的二项分布的随机数。 产生1个随机数； 产生10个随机数； %产生一个随机数%产生1行10列共10个随机数.%同命令 binornd (10,0.5,2,3).产生 %产生一个随机数 %产生1行10列共10个随机数.%同命令 binornd (10,0.5,2,3). R2=binomd (10,0.5,1,10) R3=binornd(l0,0.5,[2,3]) 2均匀分布的随机数的产生 基本数学原理：设X在区间(a,b)上有从均匀分布。 在MATLAB屮用函数unifmd产生均匀分布的随机数,其基木调用格式如下: ?R=binornd(a,b) %返冋参数为a,b的连续型均匀分布的随机数; R=binornd(a,b,m) %m指定产一 m行m列个随机数； R=binornd(a,b,m,n) %m,n分别表示产生的随机数的行数和列数。 例1-2产生区间(0.1) ±的连续型均匀分布的随机数 产生6X6个随机数； 产生6 (要求2行3列)个随机数。 解只需要命令窗口小依次输入下列命令： random 1 =unifrnd(0,1,6), %产生 6 行 6 列个随机数. randoml=unifmd(0,1,2,3), %产生 2 行 3 列个随机数. 注意 命令iinidmd(N,2,3)产生2行3列个离散型均匀分布的随机数. 3正态分布的随机数的产生 基本数学原理：设X服从参数为“和”2的正态分布。 在MATLAB屮用函数normrnd产屮参数为?的正态分彳|j的随机数，其基 本的调用格式如下： R=normrnd(MU,SIGMA) %返回均值为MU,标准差为SIGMA的正态分 布的随机数，R可以是向量或矩阵； ?R=normrnd(MU,SIGMA,m) %m指定随机数的行数与列数，与R同维数, 产生m行m列个随机数； R=normrnd(MU,SIGMA,m,n) %m,n分别表示R的行数和列数。 例1-3生成满足下列情形的正态分布随机数： 均值和标准差变化； 随机数输出为矩阵； 均值为矩阵。 解(1)在命令窗口小输入： n 1 =normrnd( 1:6,1./(1:6)) %1./(1:6)运算结果是回车后显示: 2 3 4 5 6 n 1 =2.1650 2.3134 3.0250 4.0879 4.8607 6.2827 结果表示：均值〃为1, 2, 3, 4, 5, 6,标准差b对应地为1,的正 2 3 4 5 6