_Matlab在概率統计中的应用.docVIP

PAGE PAGE 96 第8章 Matlab在概率统计中的应用 概率论与数理统计是研究和应用随机现象统计规律性的一门数学科学。其应用十分广泛，几乎遍及所有科学领域、工农业生产和国民经济各部门。本章将利用Matlab来解决概率统计学中的概率分布、数字特征、参数估计以及假设检验等问题。 8.1 数据分析 8.1.1 几种均值 在给定的一组数据中，要进行各种均值的计算，在Matlab中可由以下函数实现。 mean 算术平均值函数。对于向量X，mean (X) 得到它的元素的算术平均值；对于矩阵，mean (X)得到X各列元素的算术平均值，返回一个行向量。 nanmean 求忽略NaN的随机变量的算术平均值。 geomean 求随机变量的几何平均值。 harmmean 求随机变量的和谐平均值。 trimmean 求随机变量的调和平均值。 8.1.2 数据比较 在给定的一组数据中，还常要对它们进行最大、最小、中值的查找或对它们排序等操作。Mtalab中也有这样的功能函数。 max 求随机变量的最大值元素。 nanmax 求随机变量的忽略NaN的最大值元素。 min 求随机变量的最小值元素。 nanmin 求随机变量的忽略NaN的最小值元素。 median 求随机变量的中值。 nanmedian 求随机变量的忽略NaN的中值。 mad 求随机变量的绝对差分平均值。 sort 对随机变量由小到大排序。 sortrows 对随机矩阵按首行进行排序。 range 求随机变量的值的范围，即最大值与最小值的差（极差）。 8.1.3 累和与累积 求向量或矩阵的元素累和或累积运算是比较常用的两类运算，在Matlab中可由以下函数实现。 sum 若X为向量，sum (X)为X中各元素之和，返回一个数值；若X为矩阵，sum (X)为X中各列元素之和，返回一个行向量。 nansum 忽略NaN求向量或矩阵元素的累和。 cumsum 求当前元素与所有前面位置的元素和。返回与X同维的向量或矩阵。 cumtrapz 梯形累和函数。 prod 若X为向量，prod (X)为X中各元素之积，返回一个数值；若X为矩阵，prod (X)为X中各列元素之积，返回一个行向量。 cumprod 求当前元素与所有前面位置的元素之积。返回与X同维的向量或矩阵。 8.2 离散型随机变量的概率及概率分布 8.2.1 几个常见分布 8.2.1.1 二项分布 设随机变量X的分布律为： P {X = k } = Cnk p k (1-p) n - k k = 0, 1, 2, …, n 其中：0 p 1，n为独立重复试验的总次数，k为n次重复试验中事件A发生的次数，p为每次试验事件A发生的概率。则称X服从二项分布，记为X~B (n, p)。 8.2.1.2 Poisson分布 设随机变量X的分布律为： 则称X服从参数为的Poisson分布，记为X~P ()。Poisson是二项分布的极限分布，当二项分布中的n较大而p又较小时，常用Poisson，= np。 8.2.1.3 超几何分布 设一批同类产品共N件，其中有M件次品，从中任取n (n≤N)件，其次品数X恰为k件的概率分布为： 则称次品数X服从参数为 (N, M, n)的超几何分布。超几何分布用于无放回抽样，当N很大而n较小时，次品率在抽取前后差异很小，就用二项分布近似代替超几何分布，其中二项分布的。而且在一定条件下，也可用Poisson分布近似代替超几何分布。 8.2.2 概率密度函数值 无论是离散分布还是连续分布，在Matlab中，都用通用函数pdf或专用函数来求概率密度函数值。而对于离散型随机变量，取值是有限个或可数个，因此，其概率密度函数值就是某个特定值的概率，即利用函数pdf求输入分布的概率。 1. 通用概率密度函数pdf计算特定值的概率 命令：pdf 格式为：Y = pdf (‘name’, k, A) Y = pdf (‘name’, k, A, B) Y = pdf (‘name’, k, A, B, C) 说明：返回以name为分布，在随机变量X = k处，参数为A、B、C的概率密度值；对离散型随机变量X，返回X = k处的概率值，name为分布函数名。 常见的分布有： name = bino（二项分布），hyge（超几何分布），ge