- 1、本文档共22页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
概率统计实验指导书2,3
概率统计
实验指导书
理学院实验中心
数学专业实验室编写
2009.12
实验二 统计分析
1 引
1. 问题:湖中有鱼,其数不知。现在请你想一个办法,能将湖中的鱼数大致估计出来。
2. 分析:有两种方法。
[方法一] 设湖中有条鱼。先捕出r条鱼,做上记号后放回湖中(设记号不会消失)。让湖中的鱼充分混合后,再从湖中捕出s条鱼,设其中有T条鱼标有记号,则T是随机变量,且服从超几何分布。应用极大似然估计思想,寻找N,使达到最大,得。于是取作为湖中鱼数的一种估计,其中表示不超过x的最大整数。
[方法二] 用矩估计法.因为T服从超几何分布,其数学期望是,此即捕s条鱼得到有标记的鱼的总体平均数。而现在只捕一次,出现t条有标记的鱼。由矩估计法,令总体一阶原点矩等于样本一阶原点矩,即,也得。
3. 问题的解决:由上面的分析,要想估计出湖中的鱼数,首先需要取到样本数据,然后利用样本数据,采用统计中的点估计法对总体进行估计,其属于统计分析中的一部分。本节重点进行与统计分析相关的实验。
2 实验目的
1、利用常用的统计量描述样本数据的集中和分散程度,并对总体特征进行归纳和分析。
2. 学习用MATLAB对总体均值、方差进行估计。
3. 学习用MATLAB处理假设检验的相关问题。
4. 解决“引”中的实际问题。
3 实验内容
1.使用MATLAB对样本数据进行处理
MATLAB提供了若干对数据进行统计分析的命令,这些命令作用到一个矩阵上会对各列分别作用,得到一个行向量,现将这些命令列举如下:
max 最大分量;
mean 平均值;
std 标准差;
sum 分量和;
product 分量积;
cumsum 元素累和;
min 最小分量;
median 中位数;
sort 按不增次序排序;
hist 直方图;
diff 差分函数;
cumprod 元素累计积
概率论与数理统计实验指导书
第 PAGE 22 页 共 NUMPAGES 22 页
此外,命令corrcoef计算相关系数矩阵,格式为R=corrcoef(X),X为输入矩阵,它的行元素为观测值,列元素为变量,返回相关系数矩阵R,矩阵R的元素为R(i,j);
命令cov计算协方差矩阵,格式为C=cov(X),X若为单个向量,cov(X)返回包含方差的标量;X若为矩阵,X的每一列表示一个变量而行元素为观测值。cov(X)计算结果为协方差矩阵C,矩阵C的元素为C(i,j),它与R(i,j)的关系如下:
。
注 方差命令var(X)与diag(cov(X))所返回的结果相同;标准差命令std(X)等价于sqrt(diag(cov(X)));cov(x,y)(其中x,y为等长度的列向量)与cov([x,y])的计算结果相同。
例1 有一大批糖果,现从中随机地取8袋,称得重量X(单位:g)如下:
505.0 507.0 489.0 502.0 504.0 511.0 488.0 528.0
求这组数据的均值、方差、标准差等。
解 输入下列语句:
data1=[505,507,489,502,504,511,488,528]
mean(data1) %求data1的均值
运行得结果
ans=504.2500
再运行
std(data1) %求data1的标准差S
得结果
ans=12.6463
再运行
hist(data1) %画出data1的直方图
例2 随机地抽查某班10名学生的数学、政治、外语三门课程的考试成绩如表1所示,试求各科的平均分数、标准差、协方差及相关系数矩阵。
表1 学生成绩表
序号
科目122345678910数学
政治
外语78
82
6791
85
7667
81
7257
63
5280
80
6563
72
7669
83
6882
93
8875
82
7885
88
75
解 在MATLAB中,首先建立数据的M文件,从File菜单中选择M-File,即可打开一个编辑器,输入如下语句并以data.m命名。
function data=x
data=[78 82 67 ;91 85 76;67 81 72 ;57 63 52;80 80 65;
63 72 76;69 83 68;82 93 88;75 82 78;85 88 75];
在MATLAB命
文档评论(0)