概率统计实验指导书.docx
- 文档编号:27709876
- 上传时间:2023-07-04
- 格式:DOCX
- 页数:151
- 大小:1.27MB
概率统计实验指导书.docx
《概率统计实验指导书.docx》由会员分享,可在线阅读,更多相关《概率统计实验指导书.docx(151页珍藏版)》请在冰豆网上搜索。
概率统计实验指导书
《概率论与数理统计实验指导书》
计算机学院
2010年10月
目录
前言……………………………………………………………………………………3
实验一MATLAB的基本使用方法……………………………………………………5
实验二概率分布(概率密度)、分布函数和上分位点的数值计算………………16
实验三随机数的生成…………………………………………………………….27
实验四统计图及概率密度与分布函数作图的综合性实验……………………..32
实验五逆累积分布函数——求概率表达式{X≤C}中待定参数…………………44
实验六随机变量的数字特征……………………………………………………..48
实验七正态分布的综合性实验……………………………………………………55
实验八参数估计…………………………………………………………………..66
实验九假设检验的综合性实验…………………………………………………..73
附录1MATLAB的固有常数与常用数学函数……………………………………..82
附录2MATLAB数理统计工具箱简介……………………………………………..84
附录3MATLAB常见问题及其解决方法………………………………………….89
前言
概率论与数理统计是研究随机现象并找出其统计规律的一门学科,它广泛应用于社会、经济、科学等各个领域.随着社会生产力与科学技术的发展,这门学科的理论和应用也得到了迅速发展,特别是计算机技术及数学软件的发展使得我们不需要过多担心统计分析,即参数估计、假设检验、方差分析和回归分析等问题中的复杂的计算;也不需要过多担心大量的统计数据带来的计算量等问题.
当今社会是一个信息高度发达、人们的社会经济活动日益频繁的社会,大量的信息、数据需要人们处理.如何从这些海量的信息中提取有用的信息,用来指导人们的社会实践活动,越发显得必要而迫切,从而为数理统计提供了日益广阔的舞台.社会实践对数理统计的日益广泛而迫切的需求,对我们的教学活动提出了这样的要求:
加强数理统计的教学,充实其内容,为社会实践提供更好的服务.但要将这一要求体现到数理统计的教学中颇为困难.这是因为,目前一般工科院校均将概率论与数理统计列为一门课程.这样做的优越性自不待言,它能让学生清楚地体会二者的密切关系,将两者的思想方法融会贯通.但弊端也由此而来,由于将二者列为一门课,分配给它们的课时就相对较少,这使得教与学双方均感到这门课教学困难,学完之后也是感到没有完全理解和掌握,应用起来自然也感到较为困难.如何解决这一问题,增加课时固然是一个选择,但在目前各门学科的课时均在压缩的大趋势下不太现实,剩下的选择只能是向先进的教学方式要效益.
掌握数学软件的一些基础的操作,无疑会给每一个概率统计工作者提供了极大的方便.目前的一些概率统计新编教材也都或多或少地增加了部分数学软件内容.在概率统计课程教学中介绍数学软件的一些相关用法已成为教学改革发展的趋势.考虑到概率统计这门课的学时较紧,学生的数学软件基础不尽相同,如何在较短的时间内让学生能使用某一数学软件处理相关的概率统计问题已成为一个教改研究问题.
概率论与数理统计这门课的课时一般安排46学时,其中概率论与数理统计部分的学时分配大致是:
概率论30学时,数理统计16学时.现在我们要加强数理统计统计的教学,虽然可以适当地压缩一下概率论的学时,但概率论重要而且难学,因而压缩的空间有限.如何在此基础上较大幅度地充实数理统计的教学内容而又不致使教学效果受到影响甚或是提高?
引入数学工具软件MATLAB,将大量繁重的计算任务交由MATLAB处理,应当是一个出路.
将MATLAB引入概率统计的教学后,概率统计中的数据处理数值计算变得轻而易举,使得我们可以将精力集中于讲清处理问题的思想方法,极大地提高教学效率.用MATLAB软件辅助《概率论与数理统计》课程的教学《概率统计》是研究随机现象统计规律的一门数学学科,该课程在处理问题的思想方法上跟学生已学过的其他数学课程有着很大的差异,学生学习时感到难以掌握,根据多年的教学实践,在教学过程中要注意这门课程的特殊性,即把培养学生掌握概率统计的基本思想方法,以及解决实际问题的能力放在首位,而解决实际问题需要进行大量的数值计算.为解决以上问题,我们可以利用MATLAB辅助教学,在MATLAB7.x版本中,仅统计工具箱(StatisticToolbox)中的函数就达200多个,功能已足以赶超任何其他专用的统计软件.在应用上,MATLAB具有其它软件不可比拟的操作简单、接口方便、扩充能力强等优势.
MATLAB特点
目前的数学软件有很多,但常用的有MATLAB、Mathematica、Maple、Mathcad、SAS、SPSS等.它们的功能都很强大,但各有侧重.其中Mathematica的符号计算能力较强,MATLAB的矩阵处理能力突出,而SAS、SPSS主要侧重于数据处理与分析.很重要的一点是它们都具有很强的统计分析能力,并且随着软件版本的提高,各项功能日趋完善,操作使用也更简便.
MATLAB是MathWorks公司开发的一款以数值计算为主要特色的数学工具软件,它所带的统计工具箱几乎囊括了诸如参数估计、假设检验、方差分析、回归分析等数理统计的所有领域,并且统计工具箱中的命令调用格式极为简单方便,对财务会计院系的同学来说,掌握起来无需费多大力气.而对机械工程院、信号与电子技术院系、计算机科学工程院系等理工科系部的同学,通过统计工具箱初步了解MATLAB后还可以进一步挖掘其强大的功能,对学习其它理工类课程也极有帮助.
利用MATLAB软件辅助《概率与数理统计》教学是基于MATLAB有如下特点:
1.操作简单易学
MATLAB的基本数据结构是矩阵,它的表达式与数学计算中使用的形式十分相似,便于学习和使用.一般学生即使没有学过MATLAB,在老师的引导下,也可以在几个小时内就学会操作MATLAB.另外,计算机进入课堂的目的是辅助教学,帮助“教师教好,学生学好”该课程,不应该把大量的课时花费在掌握计算机软件的使用与编程上,应要求计算机软件是配角,决不能让它成为课程中的主角,MATLLAB可以达到该目的.
2.功能强大实用
MATLAB提供了统计工具箱,有大量的概率统计函数可直接进行计算.每当教学中增加《概率论与数理统计》新的方法和公式时,教师和学生无需编程就可以在该软件上实现,简化计算过程的繁杂与查表工作.例如:
各种概率密度函数,分布函数的计算,求数学期望、方差和相关系数等,直接调用这些函数即可方便地得到结果.
3.画图方便迅速
MATLAB可方便迅速地用绘图命令plot、plot3画各种二维、三维图形,也有专用的绘制各种统计图形的函数,可以节省大量的时间与精力.用MATLAB画图比教师在黑板上画图要准确,比事先制作的多媒体课件更灵活、生动,达到化抽象为直观的效果,大大帮助学生理解和学习概率统计的各种方法.
4.加快实际应用
《概率论与数理统计》的产生和发展都与实际紧密相联,离开了实际,这门学科就失去了意义与活力.在教学过程中,教师应尽可能以最新的实际例子来教学,使堂上教学效果更加生动有趣.而堂上教师能如此方便、快速地应用MATLAB软件就得到统计分析的结果,会进一步增加学生学习的兴趣以及用《概率与数理统计》知识和MATLAB软件解决实际问题的信心,使学生达到学以致用的目的.
实验一MATLAB的基本使用方法
一、实验问题
1.问题背景
概率论与数理统计是研究大量随机现象统计规律的一门数学学科.如何对现实中的随机现象进行模拟和处理数据,成为概率论与数理统计实验课程的重要内容.在各种数据处理软件中,MATLAB以其功能强大、操作方便著称,赢得了广大用户的青睐.本实验学习MATLAB的经常使用的操作,掌握这些基本操作将大大提高进行实验的效率.
2.实验目的与要求
(1)熟练掌握MATLAB软件的基本操作;
(2)熟练掌握MATLAB中数据输入的基本方法;
(3)熟练掌握数据加、减、乘和除四则运算的基本方法;
(4)熟练掌握函数求导数、求微分和积分运算的基本方法;
(5)熟悉与排列、组合有关的操作命令.
二、实验操作过程
1.MATLAB的基本操作
(1)启动与退出
通常安装MATLAB的计算机,在其桌面上都有MATLAB的图标.双击MATLAB图标,就可以启动MATLAB.也可以从“开始”菜单中启动.
启动MATLAB后,在MATLAB的主窗口中有几个小窗口.最常用的窗口是:
命令窗口(CommandWindow);命令历史窗口(CommandHistory);工作空间窗口(Workspace).见图1-1.
图1-1MATLAB启动画面
在MATLA中,主要的操作都在命令窗口中进行.在命令窗口中运行过的命令存储在命令历史窗口中.运行命令产生的结果存储在工作空间窗口中.命令窗口是MATLAB中最重要的窗口.命令窗口中有命令提示符“>>”,所有的命令都在命令提示符后面输入.
对命令历史窗口中存储的命令,可以用三种方式重新使用:
①在命令历史窗口中双击该命令;②在命令历史窗口中,选定命令后,再回车,就会重新运行该命令;③把命令从命令历史窗口中拖拉到命令窗口中,经过修改,再回车.
退出MATLAB也有多种方法.方法一,最简单的是单击MATLAB主窗口右上角的关闭按钮退出MATLAB.方法二,也可以在命令窗口中输入:
exit,再回车.方法三,在命令窗口中输入:
quit,回车后就可以退出MATLAB.方法四,从MATLAB主窗口左上角的File菜单中的Exit项,也能退出MATLAB.方法五,单击左上角的MATLAB图标,选中其中的退出项,退出MATLAB.方法六,用快捷键,同时按下Alt和F4两键.
(2)MATLAB的常用命令
在MATLAB中,最常用的命令有:
help,clc,clear.
(a)命令help
这是MATLAB中使用最多的一个命令.用它可以查寻命令或函数的使用方法.比如,要知道正弦函数sin的使用方法,只要在命令窗口中输入:
helpsin,回车即可显示出正弦函数的使用方式.
(b)命令clc
命令clc用来清空命令窗口.在命令窗口中输入:
clc,再回车,即可清空命令窗口.
(c)命令clear
命令clear用来清空工作空间窗口.在命令窗口中输入:
clear,再回车,就可以清空工作空间窗口.
2.数据的输入
在MATLAB中,所有的数据都是按矩阵的形式处理的,即便是一个标量,也看作是一行一列的矩阵.
(1)标量的输入
对于标量数据,只要在命令窗口中直接输入即可.
例1-1在命令窗口中输入:
a=4%将数值4赋给变量a.
回车后显示:
a=
4
在工作空间窗口中,可以看到变量a的图标,在命令历史窗口可以看到已经输入的命令:
a=4.
(2)行向量的输入
(a)直接输入:
数据放在方括号“[]”内,其间加逗号“,”或空格分开.
例1-2在命令窗口中输入:
a1=[1,3,6,8]%将行向量(1368)赋给变量a1.
回车后显示:
a1=
1368
(b)等差数列:
以确定的步长等分区间,得到等差数列.如果向量中的数据构成等差数列,则可以用冒号算符来创建.
例1-3在命令窗口中输入:
a2=1:
0.5:
3%将区间[1,3]以0.5为步长等分,赋给变量a2.
回车后显示:
a2=
1.00001.50002.00002.50003.0000
当步长为1时,还可以省略步长.
(3)列向量的输入
(a)直接输入:
数据放在方括号“[]”内,其间加分号“;”分行.
例1-4在命令窗口中输入:
b1=[1;3;6;8]%将列向量(1368)'赋给变量b1.
回车后显示:
b1=
1
3
6
8
(b)把行向量转置成列向量:
加转置运算符号“'”.
例1-5在命令窗口中输入:
b2=[1,3,6,8]'%将行向量(1368)转置后赋给变量b2.
回车后显示:
b2=
1
3
6
8
(4)矩阵的直接输入
简单的矩阵可以直接输入.其行间数据用逗号“,”或空格分隔,用分号“;”分行.
例1-6在命令窗口中输入:
A=[1,2,3;4,5,7]
回车后显示:
A=
123
457
注意:
在MATLAB中,无论是向量,还是矩阵,直接输入的时候都是用方括号“[]”括了进来.在方括号中的数据,如果是用逗号“,”分隔的,则数据在同一行中;如果是用分号“;”分隔的,则数据在不同行中.
(5)生成矩阵的函数
在MATLAB中,有许多的函数可以生成矩阵,常用的有ones,zeros,eye.这三个函数的调用方式类似.
(a)全1阵函数
用函数ones(n,m),可以生成n行m列元素全是1的矩阵.
例1-7在命令窗口中输入:
A1=ones(3)%生成3行3列的元素都是1的矩阵.
回车后显示:
A1=111
111
111
(b)全0阵函数
用函数zeros(n,m),可以生成n行m列元素全是0的矩阵.
例1-8在命令窗口中输入:
A2=zeros(2,3)%生成2行3列的元素都是0的矩阵.
回车后显示:
A2=000
000
(c)单位阵函数
用函数eye(n,m),可以生成n行m列单位矩阵.它和线性代数中讲的单位矩阵(要求方阵)含义不同.
例1-9在命令窗口中输入:
A3=eye(2,3)%生成2行3列的单位矩阵.
回车后显示:
A3=100
010
注意:
MATLAB中的单位阵的意义更广泛,不一定是方阵.
3.数组加、减、乘、除四则运算及其幂、开方、指数与对数运算
(1)数组运算
①数组与标量的四则运算
数组与标量之间的四则运算是指数组中的每个元素与标量进行加、减、乘、除运算.
例1-10对数组进行乘、除与加、减一个数的运算.
在命令窗口中输入:
x=[134;265;324];
a=2*x-2
c=x/2
回车后显示:
a=
046
2108
426
c=
0.50001.50002.0000
1.00003.00002.5000
1.50001.00002.0000
②数组间的四则运算
在MATLAB中,数组间进行四则运算时,参与运算的数组必须具有相同的维数,加、减、乘、除运算是按元素的方式进行的.其中,数组间的乘、除运算符号为“.*”,“./”或“.\”.注意,运算中的小点号不能少,否则将不会按数组运算规则进行.若没有小点号,将按矩阵的乘和求逆矩阵运算,关于矩阵间的四则运算将在下面讨论.
例1-11进行数组间的加、减法、乘法与除法运算.
在命令窗口中输入:
a=[134;265;324];
b=[231;412;453];
c=a+b
d=a./b%注意比较没有小点号时的d=a/b矩阵运算.
回车后显示:
c=
365
677
777
d=
0.50001.00004.0000
0.50006.00002.5000
0.75000.40001.3333
由于数组的除法运算有点特殊,为了便于读者使用,我们对数组的除法运算规则总结如下:
(a)数组间的除法运算为参与运算的数组中对应元素相除,结果数组与参与运算的数组大小相同.
(b)数组与标量的除法运算为数组中的每个元素与标量相除,结果数组与参与运算的数组大小相同.
(c)数组的除法运算符号有两个,即左除号“./”与右除号“.\”,它们的关系是:
a./b=b.\a.
③数组的幂运算
在MATLAB中,数组的幂运算与矩阵的幂运算完全不同.数组的幂运算符号为“.^”(注意运算符中的小点号),用来表示元素对元素的幂运算.而矩阵的幂运算符号为“^”.
例1-12进行数组与数的幂运算.
在命令窗口中输入:
a=[134;265;324];
c=a.^2
回车后显示:
c=
1916
43625
9416
为了便于比较,下面列出矩阵的幂运算.
例1-13与数组幂运算比较,进行矩阵的幂运算.
a=[134;265;324];
c=a^2
c=
192935
295258
192938
例1-14进行数组与数组的幂运算.
在命令窗口中输入:
a=[134;2,6,5;32,4];
b=[231;412;453];
c=a.^b
回车后显示:
c=
1274
16625
813264
上面两数组的幂运算为数组中各对元素间的运算.
④数组的开方运算、指数运算与对数运算
由于在MATLAB中,数组的运算实质上是数组内部每个元素的运算,因此数组的开方运算、指数运算与对数运算与标量运算完全一样,运算函数分别为“sqrt”,“exp”,“log”等.
例1-15进行数组的开方运算.
在命令窗口中输入:
a=[194;251636];
c=sqrt(a)
回车后显示:
c=
132
546
数组的对数运算、指数运算与数组的开方运算形式完全一样.
4.矩阵的基本运算
矩阵的基本运算包括矩阵的四则运算、矩阵与标量的运算、矩阵的幂运算、指数运算、对数运算、开方运算以及矩阵的逆运算、行列式运算等.下面仅对矩阵的四则运算、矩阵与标量的运算进行说明.
(1)矩阵的四则运算
矩阵的四则运算与前面讲的数组运算基本相同,但也有一些差别.其中,矩阵的加、减运算与数组的加、减运算完全相同,要求进行运算的两个矩阵的大小完全相同,使用的运算符号也是“+”与“-”.
例1-16进行矩阵加减运算.
在命令窗口中输入:
a=[12;35;26];
b=[24;18;90];
c=a+b
回车后显示:
c=
36
413
116
设矩阵A是一个i×j大小的矩阵,则要求与之相乘的矩阵B必须是一个j×k大小的矩阵,此时A与B矩阵才能进行相乘.矩阵的乘法运算使用的运算符号是“*”.
例1-17进行矩阵乘法运算.
在命令窗口中输入:
a=[12;35;26];
b=[241;890];
c=a*b%注意比较d=b*a,可见a*b≠b*a.d=b*a
回车后显示:
c=
18221
46573
52622
d=
1630
3561
当然,矩阵乘法也可以像数组乘法那样,进行矩阵元素的相乘,此时要求进行相乘的两矩阵大小完全相同,用的运算符号为“.*”.
例1-18进行矩阵乘法“*”运算,比较矩阵元素间乘法“.*”运算.
在命令窗口中输入:
a=[120;25-1;410-1];
c=[124;2510;0-1-1];
d=c.*a%注意比较e=a.*c,可见a.*c=c.*a.
e=a.*c
回车后显示:
d=
140
425-10
0-101
e=
140
425-10
0-101
在MATLAB中,矩阵的除法运算有两个运算符号,分别为左除“\”与右除“/”.矩阵的右除运算速度要慢一点,而左除运算可以避免奇异矩阵的影响.对于方程Ax=b,若此方程为超定方程,则使用除法运算符“\”与“/”可以自动找到使误差Ax-b的平方和最小的解.若此方程为不定方程,则使用除法运算符“\”与“/”求得的解至多有Rank(A)(矩阵A的秩)个非零元素,而且求得的解是这种类型的解中范数最小的一个.
例1-19进行矩阵除法运算:
解矩阵方程Ax=b.
在命令窗口中输入:
a=[213420;57820;211417;343138];
b=[10203040]';
x=b\a%方程x=A-1b,A存在逆矩阵.
回车后显示:
x=
0.7667
1.1867
0.8767
上例的方程Ax=b为超定情况.注意,结果矩阵X是列向量形式.
例1-20进行矩阵除法运算:
解矩阵方程Ax=b.
在命令窗口中输入:
a=[2134205;78202114;17343138];%A为3行4列矩阵.
b=[102030]';
x=b\a%对于方程Ax=b,A不存在逆矩阵.
回车后显示:
x=
1.6286
1.2571
1.1071
1.0500
上例的方程Ax=b为不定情况.它有三个方程、四个未知量,理论上有无穷多解.这里的解是使解中范数最小的一个.
(2)矩阵与标量的四则运算
矩阵与标量间的四则运算和数组与标量间的四则运算完全相同,即矩阵中的每个元素与标量进行加、减、乘、除四则运算.需要说明的是,当进行除法运算时,标量只能做除数.
例1-21进行矩阵与标量的四则运算.
在命令窗口中输入:
b=[213420;782021;173431];
c=b+2
d=b/2
回车后显示:
c=
233622
802223
193633
d=
10.500017.000010.0000
39.000010.000010.5000
8.500017.000015.5000
5.函数求导数、微分和积分
(1)数值微分与符号微分
微分是高等数学中
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 概率 统计 实验 指导书