实验12_父权指数与父权概率

实验12_父权指数与父权概率
实验12_父权指数与父权概率

实验十二父权指数与父权概率

【实验目的】

掌握父权指数与父权概率的计算方法,了解父权概率的意义。【实验原理】

亲子鉴定中,经过DNA分型后,若争议父亲与孩子之间的基因型不违反孟德尔遗传规律(Mendelian Law),就有两种可能:一种是该争议父亲就是孩子的亲生父亲(生物学父亲,简称生父)。另一种是他是该人群中的随机男子,此人只是偶然具有孩子的生父基因型组合,以下称为随机男子。将这两种可能进行比较就是似然率(likelihood ratio, LR),此数值即为父权指数(paternity index, PI)。显然,父权指数是反映检验结果(基因型组合)不违反孟德尔遗传规律时,对“争议父亲是孩子的生父”这一主的支持强度。依据支持强度,可以评估他们之间是否存在亲生关系。根据Bayes定理,PI可以转换成父权概率(probability of paternity),用以反映争议父亲是孩子生父可能性的高低,使结论容易被理解。

【仪器与方法】

纸、笔、计算器、电脑及计算软件。

【检材】

争议父、母亲、孩子复合STR基因座分型结果图谱。

【实验步骤】

1. 三联体(母子亲生关系已经确定)案件PI和父权概率计算

将已知实验结果的图谱进行分析,分别列出复合STR基因座分型结果,打开亲权鉴定软件,点击亲权鉴定菜单,选择三联体和试剂盒型号,按照方法输入父母子3人的基因型。输入完毕后,在生成结果表中勾选,点击右下角的计算器计算PI。

2. 二联体(父子)案件PI和父权概率计算

将已知实验结果的图谱进行分析,分别列出复合STR基因座分型结果,打开亲权鉴定软件,点击亲权鉴定菜单,选择二联体和试剂盒型号,按照方法输入父子3人的基因型。输入完毕后,在生成结果表中勾选,点击右下角的计算器计算PI。

注意事项

a. 性别基因座(AMEL)中的X和Y在程序中分别用1和2表示,在设置基因座时,基因座序号X请设为1,Y设为2。

b. 等位基因输入时,分隔符可为'/'、'-'、'+',"13/15"或"13-15"或"13+15"都是正确的,纯合子可只输入一个等位基因,如"13/13"可输入为"13",性别基因(AMEL)在输入时可不使用分隔符,如"X/Y"可直接输入为"XY","XX"可简化输入为"X"("YY"或"Y"将不被程序所授),在偶合率计算时,AMEL的Pi值始终为0.5,而在亲权鉴定中,AMEL 不列入计算。

c. 如果出现突变时,需按照司法部颁布的《亲权鉴定技术规》(SF/Z JD0105001-2010)标准方法计算变异PI值。

【实验结果与分析】

三联体(母子亲生关系已经确定)及二联体案件PI和父权概率

计算参照司法部颁布的《亲权鉴定技术规》(SF/Z JD0105001-2010)标准方法,具体如下文所示。

鉴定意见是依据DNA分型结果对是否存在亲权关系进行的判断,一般分为“排除存在亲权关系”和“支持存在亲权关系”两种。

经过累积非父排除率大于99.99%的多个基因座检测,发现有3个以上的基因座不符合遗传规律,可以排除亲权关系。

经过累积非父排除率大于99.99%的多个基因座检测,发现所检测基因座均符合遗传规律,此时必须计算亲权指数PI(即似然率LR),若CPI≥10000,则支持亲权关系。

在p、q、r分别表示等位基因P、Q、R的分布频率时,可根据下图表计算PI值。

图12-1 三联体常染体STR基因座亲权指数计算公式

图12-2 二联体常染体STR基因座亲权指数计算公式

图12-3 三联体常染体不符合遗传规律的亲权指数计算公式【实验注意事项】

1. 在遇到突变时,务必按照常染体不符合遗传规律的亲权指数

计算公式,以免PI计算错误。

【思考题】

1. 使用多少个STR位点的试剂盒才能达到亲权鉴定PI值计算支持和排除的标准?

【附图】

图12-4 父亲样本STR基因型图谱

图12-5 母亲样本STR基因型图谱

图12-6 孩子样本STR基因型图谱

图12-7 DNA相关概率计算软件界面

附:STR基因座的非父排除率参考表

基因座非父排除率基因座非父排除率D3S13580.4806TPOX0.3701 vWA0.6092CSF1PO0.5001 FGA0.7068D7S8200.5576

D8S11790.6882D2S13380.6950

D18S510.7184D19S4330.6554

D21S110.6502D6S10430.7310

D5S8180.5671D12S3910.6790

D13S3170.6000Penta D0.5907

D16S5390.5749Penta E0.7325 TH010.4046

附:19个STR基因座的等位基因分布频率

概率统计实验复习过程

§13.6 概率统计实验 [学习目标] 1. 会用Mathematica 求概率、均值与方差; 2. 能进行常用分布的计算; 3. 会用Mathematica 进行期望和方差的区间估计; 4. 会用Mathematica 进行回归分析。 概率统计是最需要使用计算机的领域,过去依靠计算器进行统计计算,由于计算机的普及得以升级换代。本节介绍Mathematica 自带的统计程序包,其中有实现常用统计计算的各种外部函数。 一、 样本的数字特征 1. 一元的情况 Mathematica 的内部没有数理统计方面的功能,但是带有功能强大的数理统计外部程序,由多个程序文件组成。它们在标准扩展程序包集的Statistic 程序包子集中,位于目录 D :\Mathematica\4.0\AddOns\StandardPackages\Statistics 下。通过查看Help ,可以找到包含所需外部函数的程序文件名。 在程序文件DescriptiveStatistics.m 中,含有实现一元数理统计基本计算的函数,常用的有: SampleRange[data] 求表data 中数据的极差(最大数减最小数)。 Median[data] 求中值。 Mean[data] 求平均值∑=n i i x n 1 1。 Variance[data] 求方差(无偏估计)∑=--n i i x x n 12)(11。 StandardDeviation[data] 求标准差(无偏估计)∑=--n i i x x n 1 2)(11。 VarianceMLE[data] 求方差∑=-n i i x x n 1 2)(1。 StandardDeviationMLE[data] 求标准差∑=-n i i x x n 1 2)(1。 实际上程序文件中的函数很多,这里只列出了最常用的函数,其它计算函数可以通过Help 浏览。 例1 给出一组样本值:6.5,3.8,6.6,5.7,6.0,6.4,5.3,计算样本个数、最大值、最小值、均值、方差、标准差等。

概率统计-习地的题目及答案详解(1)

习题一 1.1 写出下列随机试验的样本空间,并把指定的事件表示为样本点的集合: (1)随机试验:考察某个班级的某次数学考试的平均成绩(以百分制记分,只取整数); 设事件A 表示:平均得分在80分以上。 (2)随机试验:同时掷三颗骰子,记录三颗骰子点数之和; 设事件A 表示:第一颗掷得5点; 设事件B 表示:三颗骰子点数之和不超过8点。 (3)随机试验:一个口袋中有5只球,编号分别为1,2,3,4,5,从中取三个球; 设事件A 表示:取出的三个球中最小的号码为1。 (4)随机试验:某篮球运动员投篮练习,直至投中十次,考虑累计投篮的次数; 设事件A 表示:至多只要投50次。 (5)随机试验:将长度为1的线段任意分为三段,依次观察各段的长度。 1.2 在分别标有号码1~8的八张卡片中任抽一张。 (1)写出该随机试验的样本点和样本空间; (2)设事件A 为“抽得一张标号不大于4的卡片”,事件B 为“抽得一张标号为偶数的 卡片”,事件C 为“抽得一张标号能被3整除的卡片”。 试将下列事件表示为样本点的集合,并说明分别表示什么事件? (a )AB ; (b) B A +; (c) B ; (d) B A -; (e) BC ; (f) C B + 。 1.3 设A 、B 、C 是样本空间的事件,把下列事件用A 、B 、C 表示出来: (1)A 发生; (2)A 不发生,但B 、C 至少有一个发生; (3)三个事件恰有一个发生; (4)三个事件中至少有两个发生; (5)三个事件都不发生; (6)三个事件最多有一个发生; (7)三个事件不都发生。 1.4 设}10,,3,2,1{ =Ω,}5,3,2{=A ,}7,5,3{=B ,}7,4,3,1{=C ,求下列事件: (1)B A ; (2))(BC A 。 1.5 设A 、B 是随机事件,试证:B A AB A B B A +=-+-)()(。 1.6 在11张卡片上分别写上Probability 这11个字母,从中任意抽取7张,求其排列结果为ability 的概率。 1.7 电话号码由6位数字组成,每个数字可以是0,1,2,…,9中的任一个数字(但第一位不能为0),求电话号码是由完全不相同的数字组成的概率。 1.8 把10本不同的书任意在书架上放成一排,求其中指定的3本书恰好放在一起的概率。

《数学实验》试题答案

北京交通大学海滨学院考试试题 课程名称:数学实验2010-2011第一学期出题教师:数学组适用专业: 09机械, 物流, 土木, 自动化 班级:学号:姓名: 选做题目序号: 1.一对刚出生的幼兔经过一个月可以长成成兔, 成兔再经过一个月后可以 繁殖出一对幼兔. 如果不计算兔子的死亡数, 请用Matlab程序给出在未来24个月中每个月的兔子对数。 解: 由题意每月的成兔与幼兔的数量如下表所示: 1 2 3 4 5 6 ··· 成兔0 1 1 2 3 5··· 幼兔 1 0 1 1 2 3··· 运用Matlab程序: x=zeros(1,24); x(1)=1;x(2)=1; for i=2:24 x(i+1)=x(i)+x(i-1); end x 结果为x = 1 1 2 3 5 8 13 21 3 4 5 5 89 144 233 377 610 987 1597 2584 4181 6765 1094 6 7711 2865 7 46368 2.定积分的过程可以分为分割、求和、取极限三部分, 以1 x e dx 为例, 利用

已学过的Matlab 命令, 通过作图演示计算积分的过程, 并与使用命令int() 直接积分的结果进行比较. 解:根据求积分的过程,我们先对区间[0,1]进行n 等分, 然后针对函数x e 取和,取和的形式为10 1 i n x i e e dx n ξ=≈ ∑ ? ,其中1[ ,]i i i n n ξ-?。这里取i ξ为区间的右端点,则当10n =时,1 x e dx ?可用10 101 1.805610 i i e ==∑ 来近似计算, 当10n =0时,100 100 1 01 =1.7269100 i x i e e dx =≈ ∑?,当10n =000时,10000 10000 1 1 =1.718410000 i x i e e dx =≈ ∑ ?. 示意图如下图,Matlab 命令如下: x=linspace (0,1,21); y=exp(x); y1=y(1:20); s1=sum(y1)/20 y2=y(2:21); s2=sum(y2)/20 plot(x,y); hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],'b') end syms k;symsum(exp(k/10)/10,k,1,10);%n=10 symsum(exp(k/100)/100,k,1,100);%n=100 symsum(exp(k/10000)/10000,k,1,10000);%n=10000

实验5:概率统计实验

撰写人姓名:撰写时间:审查人姓名: 实验全过程记录实验 名称概率统计实验 时间2学时 地点数学实验室 姓名学号 同实验者学号 一、实验目的 1、掌握利用MATLAB处理简单的概率问题; 2、掌握利用MATLAB处理简单的数理统计问题。 二、实验内容: 1、熟练掌握几种常用的离散型、连续型随机变量的函数命令; 2、熟练掌握常用的描述样本数据特征的函数命令(如最值、均值、中位数(中值)、方差、标准差、几何平均值、调和平均值、协方差、相关系数等); 3、掌握常用的MATLAB统计作图方法(如直方图、饼图等); 4、能用MATLAB以上相关命令解决简单的数据处理问题; 5、熟练掌握常用的参数估计和假设检验的相关的函数命令; 6、能用参数估计和假设检验等相关命令解决简单的实际问题。 三、实验用仪器设备及材料 软件需求: 操作系统:Windows XP或更新的版本; 实用数学软件:MATLAB 7.0或更新的版本。 硬件需求: Pentium IV 450以上的CPU处理器、512MB以上的内存、5000MB的自由硬盘空间、CD-ROM驱动器、打印机、打印纸等。 四、实验原理: 概率论与数理统计等相关理论 五、实验步骤: 1、对下列问题,请分别用专用函数和通用函数实现。 ⑴X服从[3, 10]上均匀分布,计算P{X≤4},P{X>8};已知P{X>a}=0.4,求a。 p1=unifcdf(4,3,10) p2=1-unifcdf(8,3,10) p11=cdf('unif',4,3,10) p22=1-cdf('unif',8,3,10) unifinv(0.6,3,10) icdf('unif',0.6,3,10) p1 =

数学实验答案-1

1.(1) [1 2 3 4;0 2 -1 1;1 -1 2 5;]+(1/2).*([2 1 4 10;0 -1 2 0;0 2 3 -2]) 2. A=[3 0 1;-1 2 1;3 4 2],B=[1 0 2;-1 1 1;2 1 1] X=(B+2*A)/2 3. A=[-4 -2 0 2 4;-3 -1 1 3 5] abs(A)>3 % 4. A=[-2 3 2 4;1 -2 3 2;3 2 3 4;0 4 -2 5] det(A),eig(A),rank(A),inv(A) 求计算机高手用matlab解决。 >> A=[-2,3,2,4;1,-2,3,2;3,2,3,4;0,4,-2,5] 求|A| >> abs(A) ans = ( 2 3 2 4 1 2 3 2 3 2 3 4 0 4 2 5 求r(A) >> rank(A) ans =

4 求A-1 《 >> A-1 ans = -3 2 1 3 0 -3 2 1 2 1 2 3 -1 3 -3 4 求特征值、特征向量 >> [V,D]=eig(A) %返回矩阵A的特征值矩阵D 与特征向量矩阵V , V = - + + - - + - + - + - + D = { + 0 0 0 0 - 0 0 0 0 + 0 0 0 0 - 将A的第2行与第3列联成一行赋给b >> b=[A(2,:),A(:,3)'] b = 《 1 - 2 3 2 2 3 3 -2

1. a=round(unifrnd(1,100)) i=7; while i>=0 i=i-1; b=input('请输入一个介于0到100的数字:'); if b==a ¥ disp('You won!'); break; else if b>a disp('High'); else if b

高中数学统计与概率知识点

高中数学统计与概率知识点(文) 第一部分:统计 一、什么是众数。 一组数据中出现次数最多的那个数据,叫做这组数据的众数。 众数的特点。 ①众数在一组数据中出现的次数最多;②众数反映了一组数据的集中趋势,当众数出现的次数越多,它就越能代表这组数据的整体状况,并且它能比较直观地了解到一组数据的大致情况。但是,当一组数据大小不同,差异又很大时,就很难判断众数的准确值了。此外,当一组数据的那个众数出现的次数不具明显优势时,用它来反映一组数据的典型水平是不大可靠的。 3.众数与平均数的区别。 众数表示一组数据中出现次数最多的那个数据;平均数是一组数据中表示平均每份的数量。 二、.中位数的概念。 一组数据按大小顺序排列,位于最中间的一个数据(当有偶数个数据时,为最中间两个数据的平均数)叫做这组数据的中位数。 三.众数、中位数及平均数的求法。 ①众数由所给数据可直接求出;②求中位数时,首先要先排序(从小到大或从大到小),然后根据数据的个数,当数据为奇数个时,最中间的一个数就是中位数;当数据为偶数个时,最中间两个数的平均数就是中位数。③求平均数时,就用各数据的总和除以数据的个数,得数就是这组数据的平均数。 四、中位数与众数的特点。 ⑴中位数是一组数据中唯一的,可能是这组数据中的数据,也可能不是这组数据中的数据; ⑵求中位数时,先将数据有小到大顺序排列,若这组数据是奇数个,则中间的数据是中位数;若这组数据是偶数个时,则中间的两个数据的平均数是中位数; ⑶中位数的单位与数据的单位相同; ⑷众数考察的是一组数据中出现的频数; ⑸众数的大小只与这组数的个别数据有关,它一定是一组数据中的某个数据,其单位与数据的单位相同;(6)众数可能是一个或多个甚至没有;(7)平均数、众数和中位数都是描述一组数据集中趋势的量。

概率统计实验报告(三)剖析

线性回归实验报告(三) 实验目的:通过本次实验,了解matlab和spss在非参数检验中的应用,学会用matlab和spss做非参数假设检验,主要包括单样本和多样本非参数假设检验。 实验内容: 1.单样本假设检验; 2.多样本假设检验. 实验结果与分析: 1.单样本K-S儿童身高 操作步骤: ⑴分析-非参数检验-旧对话框-1-样本KS; ⑵将“周岁儿童身高”变换到检验变量列表,由于样本量太少,点击精确按钮,选择精确检验方法; ⑶回到K-S检验对话框,点击选项按钮,设置输出参数,勾选描述性和四分位数; ⑷输出检验结果。 从图形特征上看,儿童身高的分布非常接近正态分布,但是仍需要用K-S来检验

诊断。 结论:K-S检验统计量Z值为0.936,显著性为0.344,大于显著性水平0.05,所以不能拒绝原假设,认为周岁儿童的身高服从正态分布。 2.单样本游程——电缆 操作步骤: ⑴分析-非参数检验-旧对话框-游程; ⑵将“耐电压值”变换到检验变量列表; ⑶回到游程检验对话框,点击选项按钮,设置输出参数,勾选描述性和四分位数; ⑷输出检验结果。

结论:中位数渐进显著性为0.491,平均数和众数为1,大于显著性水平0.05,所以不能拒绝原假设,所以该组电缆耐电压值是随机的。 3.多独立样本——儿童身高 操作步骤: ⑴分析-非参数检验-旧对话框-K个独立样本检验; ⑵将“周岁儿童身高”变换到检验变量列表;将“城市标志”变换到分组变量,设置分组变量范围; ⑶回到多独立样本检验对话框,点击选项按钮,设置输出参数,勾选描述性和四分位数; ⑷输出检验结果。

结论:多个样本的K-W检验,即秩和检验目的是看各总体的位置参数是否一样,渐近显著性值为0.003,小于显著性水平0.05,所以拒绝原假设,因而四个城市儿童身高的分布存在显著性差异。 4.多样本配对——促销方式 操作步骤: ⑴分析-非参数检验-旧对话框-K个相关样本检验; ⑵将“促销形式1”、“促销形式2”、“促销形式3”变换到检验变量列表; ⑶回到多个关联样本检验对话框,点击选项按钮,设置输出参数,勾选描述性和四分位数; ⑷输出检验结果。

概率统计实验报告

概率统计实验报告 班级16030 学号16030 姓名 2018 年1 月3 日

1、 问题概述和分析 (1) 实验内容说明: 题目12、(综合性实验)分析验证中心极限定理的基本结论: “大量独立同分布随机变量的和的分布近似服从正态分布”。 (2) 本门课程与实验的相关内容 大数定理及中心极限定理; 二项分布。 (3) 实验目的 分析验证中心极限定理的基本结论。 2、实验设计总体思路 2.1、引论 在很多实际问题中,我们会常遇到这样的随机变量,它是由大量的相互独立的随机 因素的综合影响而形成的,而其中每一个个别因素在总的影响中所起的作用是微小的,这种随机变量往往近似的服从正态分布。 2.2、 实验主题部分 2.2.1、实验设计思路 1、 理论分析 设随机变量X1,X2,......Xn ,......独立同分布,并且具有有限的数学期望和方差:E(Xi)=μ,D(Xi)=σ2(k=1,2....),则对任意x ,分布函数 满足 该定理说明,当n 很大时,随机变量 近似地服从标准正 态分布N(0,1)。因此,当n 很大时, 近似地服从正 态分布N(n μ,n σ2). 2、实现方法(写清具体实施步骤及其依据) (1) 产生服从二项分布),10(p b 的n 个随机数, 取2.0=p , 50=n , 计算n 个随 机数之和y 以及 ) 1(1010p np np y --; 依据:n 足够大,且该二项分布具有有限的数学期望和方差。 (2) 将(1)重复1000=m 组, 并用这m 组 ) 1(1010p np np y --的数据作频率直方图进 行观察. 依据:通过大量数据验证随机变量的分布,且符合极限中心定理。

matlab数学实验复习题(有标准答案)

复习题 1、写出3 2、i nv(A)表示A的逆矩阵; 3、在命令窗口健入 clc,4、在命令窗口健入clea 5、在命令窗口健入6、x=-1:0.2:17、det(A)表示计算A的行列式的值;8、三种插值方法:拉格朗日多项式插值,分段线性插值,三次样条插值。 9、若A=123456789?? ????????,则fliplr (A)=321654987?????????? A-3=210123456--??????????A .^2=149162536496481?????????? tril(A)=100450789?????????? tri u(A,-1)=123456089??????????diag(A )=100050009?????????? A(:,2),=2 58A(3,:)=369 10、nor mcd f(1,1,2)=0.5%正态分布mu=1,s igm a=2,x =1处的概率 e45(@f,[a,b ],x0),中参数的涵义是@fun 是求解方程的函数M 文 件,[a,b ]是输入向量即自变量的范围a 为初值,x0为函数的初值,t 为输出指定的[a,b],x 为函数值 15、写出下列命令的功能:te xt (1,2,‘y=s in(x)’

hold on 16fun ction 开头; 17 ,4) 3,4) 21、设x 是一向量,则)的功能是作出将X十等分的直方图 22、interp 1([1,2,3],[3,4,5],2.5) Ans=4.5 23、建立一阶微分方程组? ??+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来) 二、写出运行结果: 1、>>ey e(3,4)=1000 01000010 2、>>s ize([1,2,3])=1;3 3、设b=ro und (unifrnd(-5,5,1,4)),则=3 5 2 -5 >>[x,m]=min(b);x =-5;m=4 ,[x,n ]=sort(b ) -5 2 3 5 4 3 1 2 mea n(b)=1.25,m edian(b)=2.5,range(b)=10 4、向量b如上题,则 >>an y(b),all(b<2),all(b<6) Ans =1 0 1 5、>>[5 6;7 8]>[7 8;5 6]=00 11 6、若1234B ??=???? ,则 7、>>diag(d iag (B ))=10 04 8、>>[4:-2:1].*[-1,6]=-4 12 9、>>acos(0.5),a tan(1) ans = 1.6598 ans=

高中数学专题――概率统计专题.

专题二概率统计专题 【命题趋向】概率与统计是高中数学的重要学习内容,它是一种处理或然问题的方法,在工农业生产和社会生活中有着广泛的应用,渗透到社会的方方面面,概率与统计的基础知识成为每个公民的必备常识.概率与统计的引入,拓广了应用问题取材的范围,概率的计算、离散型随机变量的分布列和数学期望的计算及应用都是考查应用意识的良好素材.在高考试卷中,概率与统计的内容每年都有所涉及,以解答题形式出现的试题常常设计成包含离散型随机变量的分布列与期望、统计图表的识别等知识为主的综合题,以考生比较熟悉的实际应用问题为载体,以排列组合和概率统计等基础知识为工具,考查对概率事件的识别及概率计算.解答概率统计试题时要注意分类与整合、化归与转化、或然与必然思想的运用.由于中学数学中所学习的概率与统计内容是最基础的,高考对这一部分内容的考查注重考查基础知识和基本方法.该部分在高考试卷中,一般是2—3个小题和一个解答题. 【考点透析】概率统计的考点主要有:概率与统计包括随机事件,等可能性事件的概率,互斥事件有一个发生的概率,古典概型,几何概型,条件概率,独立重复试验与二项分布,超几何分布,离散型随机变量的分布列,离散型随机变量的期望和方差,抽样方法,总体分布的估计,正态分布,线性回归等.【例题解析】 题型1 抽样方法 -)中,在公证部门监督下按照随机抽取的方法确【例1】在1000个有机会中奖的号码(编号为000999 定后两位数为的号码为中奖号码,该抽样运用的抽样方法是() A.简单随机抽样B.系统抽样C.分层抽样D.以上均不对 分析:实际“间隔距离相等”的抽取,属于系统抽样. 解析:题中运用了系统抽样的方法采确定中奖号码,中奖号码依次为:088,188,288,388,488,588,688,788,888,988.答案B. 点评:关于系统抽样要注意如下几个问题:(1)系统抽样是将总体分成均衡几个部分,然按照预先定出的规则从每一部分抽取一个个体,得到所需要的样本的一种抽样方法.(2)系统抽样的步骤:①将总体中的个体随机编号;②将编号分段;③在第一段中用简单随机抽样确定起始的个体编号;④按事先研究的规则抽取样本.(3)适用范围:个体数较多的总体. 例2(2008年高考广东卷理3)某校共有学生2000名,各年级男、女生人数如表.已知在全校学生中随机抽取1名,抽到二年级女生的概率是0.19.现用分层抽样的方法在全校抽取64名学生,则应在三年级抽取的学生人数为() A.24B.18C.16D.12 Array 分析:根据给出的概率先求出x的值,这样就可以知道三年级的学生人数,问题就解决了. x=?=,这样一年级和二年级学生的解析:C 二年级女生占全校学生总数的19%,即20000.19380 +++=,三年级学生有500人,用分层抽样抽取的三年级学生应是总数是3733773803701500 64 50016 ?=.答案C. 2000 点评:本题考查概率统计最基础的知识,还涉及到一点分析问题的能力和运算能力,题目以抽样的等可能性为出发点考查随机抽样和分层抽样的知识. 例3.(2009江苏泰州期末第2题)一个社会调查机构就某地居民的月收入调查了10000人,并根据所得数据画了样本的频率分布直方图(如下图).为了分析居民的收入与年龄、学历、职业等方面的关系, 2500,3500(元)月收入段应抽要从这10000人中再用分层抽样方法抽出100人作进一步调查,则在[) 出人.

概率论与数理统计数学实验

概率论与数理统计数学实验 目录 实验一几个重要的概率分布的MATLAB实现 p2-3实验二数据的统计描述和分析 p4-8实验三参数估计 p9-11实验四假设检验 p12-14实验五方差分析 p15-17实验六回归分析 p18-27

实验一 几个重要的概率分布的MATLAB 实现 实验目的 (1) 学习MATLAB 软件与概率有关的各种计算方法 (2) 会用MATLAB 软件生成几种常见分布的随机数 (3) 通过实验加深对概率密度,分布函数和分位数的理解 Matlab 统计工具箱中提供了约20种概率分布,对每一种分布提供了5种运算功能,下表给出了常见8种分布对应的Matlab 命令字符,表2给出了每一种运算功能所对应的Matlab 命令字符。当需要某一分布的某类运算功能时,将分布字符与功能字符连接起来,就得到所要的命令。 例1 求正态分布()2,1-N ,在x=处的概率密度。 解:在MATLAB 命令窗口中输入: normpdf,-1,2) 结果为: 例2 求泊松分布()3P ,在k=5,6,7处的概率。 解:在MATLAB 命令窗口中输入: poisspdf([5 6 7],3) 结果为: 例3 设X 服从均匀分布()3,1U ,计算{}225P X .-<<。 解:在MATLAB 命令窗口中输入: unifcdf,1,3)-unifcdf(-2,1,3) 结果为:

例4 求概率995.0=α的正态分布()2,1N 的分位数αX 。 解:在MATLAB 命令窗口中输入: norminv,1,2) 结果为: 例5 求t 分布()10t 的期望和方差。 解:在MATLAB 命令窗口中输入: [m,v]=tstat(10) m = 0 v = 例6 生成一个2*3阶正态分布的随机矩阵。其中,第一行3个数分别服从均值为1,2,3;第二行3个数分别服从均值为4,5,6,且标准差均为的正态分布。 解:在MATLAB 命令窗口中输入: A=normrnd([1 2 3;4 5 6],,2,3) A = 例7 生成一个2*3阶服从均匀分布()3,1U 的随机矩阵。 解:在MATLAB 命令窗口中输入: B=unifrnd(1,3,2,3) B = 注:对于标准正态分布,可用命令randn(m,n);对于均匀分布()1,0U ,可用命令rand(m,n)。

实验二极限与连续数学实验课件习题答案

天水师范学院数学与统计学院 实验报告 实验项目名称极限与连续 所属课程名称数学实验 实验类型上机操作 实验日期 2013-3-22 班级 10数应2班 学号 291010836 姓名吴保石 成绩

【实验过程】(实验步骤、记录、数据、分析) 1.数列极限的概念 通过计算与作图,加深对极限概念的理解. 例2.1 考虑极限3321 lim 51 x n n →∞++ Print[n ," ",Ai ," ",0.4-Ai]; For[i=1,i 15,i++,Aii=N[(2i^3+1)/(5i^3+1),10]; Bii=0.4-Aii ;Print[i ," ",Aii ," ",Bii]] 输出为数表 输入 fn=Table[(2n^3+1)/(5n^3+1),{n ,15}]; ListPlot[fn ,PlotStyle {PointSize[0.02]}] 观察所得散点图,表示数列的点逐渐接近直线y=0 .4 2.递归数列 例2.2 设n n x x x +==+2,211.从初值21=x 出发,可以将数列一项项地计算出来,这样定义的数列称为 数列,输入 f[1]=N[Sqrt[2],20]; f[n_]:=N[Sqrt[2+f[n-1]],20]; f[9] 则已经定义了该数列,输入 fn=Table[f[n],{n ,20}] 得到这个数列的前20项的近似值.再输入 ListPlot[fn ,PlotStyle {PointSize[0.02]}] 得散点图,观察此图,表示数列的点越来越接近直线2y =

例2.3 考虑函数arctan y x =,输入 Plot[ArcTan[x],{x ,-50,50}] 观察函数值的变化趋势.分别输入 Limit[ArcTan[x],x Infinity ,Direction +1] Limit[ArcTan[x],x Infinity ,Direction -1] 输出分别为2 π 和2π-,分别输入 Limit[sign[x],x 0,Direction +1] Limit[Sign[x],x 0,Direction -1] 输出分别为-1和1 4.两个重要极限 例2.4 考虑第一个重要极限x x x sin lim 0→ ,输入 Plot[Sin[x]/x ,{x ,-Pi ,Pi}] 观察函数值的变化趋势.输入 Limit[Sin[x]/x ,x 0] 输出为1,结论与图形一致. 例2.5 考虑第二个重要极限1 lim(1)x x x →∞+,输入 Limit[(1+1/n)^n ,n Infinity] 输出为e .再输入 Plot[(1+1/n)^n ,{n ,1,100}] 观察函数的单调性 5.无穷大 例2.6 考虑无穷大,分别输人 Plot[(1+2x)/(1-x),{x ,-3,4}] Plot[x^3-x ,{x ,-20,20}] 观察函数值的变化趋势.输入 Limit[(1+2x)/(1-x),x 1] 输出为-∞ 例2.7 考虑单侧无穷大,分别输人 Plot[E^(1/x),{x ,-20,20},PlotRange {-1,4}] Limit[E^(1/x),x 0,Direction +1] Limit[E^(1/x),x 0,Direction -1] 输出为图2.8和左极限0,右极限∞.再输入 Limit[E^(1/x),x 0] 观察函数值的变化趋势. 例2.8 输入 Plot[x+4*Sin[x],{x ,0,20Pi}] 观察函数值的变化趋势. 输出为图2 .9.观察函数值的变化趋势,当x →∞时,这个函数是无穷大,但是,它并不是单调增加.于是,无并不要求函数单调 例2.9 输入

概率论与数理统计实验报告

概率论与数理统计 实验报告 概率论部分实验二 《正态分布综合实验》

实验名称:正态分布综合实验 实验目的:通过本次实验,了解Matlab在概率与数理统计领域的应用,学会用matlab做概率密度曲线,概率分布曲线,直方图,累计百分比曲线等简单应用;同时加深对正态分布的认识,以更好得应用之。 实验内容: 实验分析: 本次实验主要需要运用一些matlab函数,如正态分布随机数发生器normrnd函数、绘制直方图函数hist函数、正态分布密度函数图形绘制函数normpdf函数、正态分布分步函数图形绘制函数normcdf等;同时,考虑到本次实验重复性明显,如,分别生成100,1000,10000个服从正态分布的随机数,进行相同的实验操作,故通过数组和循环可以简化整个实验的操作流程,因此,本次实验程序中要设置数组和循环变量。 实验过程: 1.直方图与累计百分比曲线 1)实验程序 m=[100,1000,10000]; 产生随机数的个数 n=[2,1,0.5]; 组距 for j=1:3 for k=1:3 x=normrnd(6,1,m(j),1); 生成期望为6,方差为1的m(j)个 正态分布随机数

a=min(x); a为生成随机数的最小值 b=max(x); b为生成随机数的最大值 c=(b-a)/n(k); c为按n(k)组距应该分成的组数 subplot(1,2,1); 图形窗口分两份 hist(x,c);xlabel('频数分布图'); 在第一份里绘制频数直方图 yy=hist(x,c)/1000; yy为各个分组的频率 s=[]; s(1)=yy(1); for i=2:length(yy) s(i)=s(i-1)+yy(i); end s[]数组存储累计百分比 x=linspace(a,b,c); subplot(1,2,2); 在第二个图形位置绘制累计百分 比曲线 plot(x,s,x,s);xlabel('累积百分比曲线'); grid on; 加网格 figure; 另行开辟图形窗口,为下一个循 环做准备 end end 2)实验结论及过程截图 实验结果以图像形式展示,以下分别为产生100,1000,10000个正态分布随机数,组距分别为2,1,0.5的频数分布直方图和累积百分比曲线,从实验结果看来,随着产生随机数的数目增多,组距减小,累计直方图逐渐逼近正态分布密度函数图像,累计百分比逐渐逼近正态分布分布函数图像。

数学实验四(概率论)_6

数学实验四(概率论) 一.用MATLAB 计算随机变量的分布 1.用MA TLAB 计算二项分布 当随变量(),X B n p 时,在MATLAB 中用命令函数 (,,)Px binopdf X n p = 计算某事件发生的概率为p 的n 重贝努利试验中,该事件发生的次数为X 的概率。 例1 在一级品率为0.2的大批产品中,随机地抽取20个产品,求其中有2个一级品的概率。 解 在MATLAB 中,输入 >>clear >> Px=binopdf(2,20,0.2) Px = 0.1369 即所求概率为0.1369。 2.用MA TLAB 计算泊松分布 当随变量()X P λ 时,在MATLAB 中用命令函数 (,)P poisspdf x lambda = 计算服从参数为lambda 的泊松分布的随机变量取值x 的概率。用命令函数 (,)P poisscdf x lambda = 计算服从参数为lambda 的泊松分布的随机变量在[]0,x 取值的概率。 例2 用MATLAB 计算:保险公司售出某种寿险保单2500份.已知此项寿险每单需交保费120元,当被保人一年内死亡时,其家属可以从保险公司获得2万元的赔偿(即保额为2万元).若此类被保人一年内死亡的概率0.002,试求: (1)保险公司的此项寿险亏损的概率; (2)保险公司从此项寿险获利不少于10万元的概率; (3)获利不少于20万元的概率. 利用泊松分布计算. 25000.0025np λ==?= (1) P(保险公司亏本)= ()()15 250025000(3020)1(15)10.0020.998k k k k P X P X C -=-<=-≤=- ?∑ =15 5 051! k k e k -=-∑ 在MATLAB 中,输入 >> clear >> P1=poisscdf(15,5) P1 = 0. 9999 即 15 5 05! k k e k -=∑= P1 =0.9999 故 P(保险公司亏本)=1-0.9999=0.0001

数学实验(MATLAB版韩明版)5.1,5.3,5.5,5.6部分答案

练习 B的分布规律和分布函数的图形,通过观1、仿照本节的例子,分别画出二项分布()7.0,20 察图形,进一步理解二项分布的性质。 解:分布规律编程作图:>> x=0:1:20;y=binopdf(x,20,; >> plot(x,y,'*') 图像: y x 分布函数编程作图:>> x=0::20; >>y=binocdf(x,20, >> plot(x,y) 图像: 《

1 x 观察图像可知二项分布规律图像像一条抛物线,其分布函数图像呈阶梯状。 2、仿照本节的例子,分别画出正态分布()25,2N的概率密度函数和分布函数的图形,通过观察图形,进一步理解正态分布的性质。 解:概率密度函数编程作图:>> x=-10::10; >> y=normpdf(x,2,5); >> plot(x,y) 图像:

00.010.020.030.040.050.060.070.08x y 分布函数编程作图:>> x=-10::10; >> y=normcdf(x,2,5); ~ >> plot(x,y) 图像:

01x y 观察图像可知正态分布概率密度函数图像像抛物线,起分布函数图像呈递增趋势。 3、设()1,0~N X ,通过分布函数的调用计算{}11<<-X P ,{}22<<-X P , {}33<<-X P . 解:编程求解: >> x1=normcdf(1)-normcdf(-1),x2=normcdf(2)-normcdf(-2),x3=normcdf(3)-normcdf(-3) x1 = x2 = ) x3 = 即:{}6827.011=<<-X P ,{}9545.022=<<-X P ,{}9973.033=<<-X P . 4、设()7.0,20~B X ,通过分布函数的调用计算{}10=X P 与{}10> x1=binopdf(10,20,,x2=binocdf(10,20,-binopdf(10,20, x1 = x2 =

高考数学概率与统计知识点汇编

高中数学之概率与统计 求等可能性事件、互斥事件和相互独立事件的概率 解此类题目常应用以下知识: (1)等可能性事件(古典概型)的概率:P(A)=)()(I card A card =n m ; 等可能事件概率的计算步骤: 计算一次试验的基本事件总数n ; 设所求事件A ,并计算事件A 包含的基本事件的个数m ; 依公式 ()m P A n = 求值; 答,即给问题一个明确的答复. (2)互斥事件有一个发生的概率:P(A +B)=P(A)+P(B); 特例:对立事件的概率:P(A)+P(A )=P(A +A )=1. (3)相互独立事件同时发生的概率:P(A ·B)=P(A)·P(B); 特例:独立重复试验的概率:Pn(k)=k n k k n p p C --)1(.其中P 为事件A 在一次试验中发生的 概率,此式为二项式[(1-P)+P]n 展开的第k+1项. (4)解决概率问题要注意“四个步骤,一个结合”: 求概率的步骤是: 第一步,确定事件性质?? ?? ???等可能事件 互斥事件 独立事件 n 次独立重复试验 即所给的问题归结为四类事件中的某一种. 第二步,判断事件的运算 ?? ?和事件积事件 即是至少有一个发生,还是同时发生,分别运用相加或相乘事件. 第三步,运用公式()()()()()()()()(1) k k n k n n m P A n P A B P A P B P A B P A P B P k C p p -? =???+=+? ??=??=-??等可能事件: 互斥事件: 独立事件: n 次独立重复试验:求解 第四步,答,即给提出的问题有一个明确的答复. 例1. 在五个数字12345,,,,中,。 例2. 若随机取出三个数字,则剩下两个数字都是奇数的概率是 (结果用数值表示). [解答过程]0.3提示:13 35C 33. 54C 10 2P ===?

概率论与数理统计实验报告

概率论与数理统计实验报告 一、实验目的 1.学会用matlab求密度函数与分布函数 2.熟悉matlab中用于描述性统计的基本操作与命令 3.学会matlab进行参数估计与假设检验的基本命令与操作 二、实验步骤与结果 概率论部分: 实验名称:各种分布的密度函数与分布函数 实验内容: 1.选择三种常见随机变量的分布,计算它们的方差与期望<参数自己设 定)。 2.向空中抛硬币100次,落下为正面的概率为0.5,。记正面向上的次数 为x, (1)计算x=45和x<45的概率, (2)给出随机数x的概率累积分布图像和概率密度图像。 3.比较t(10>分布和标准正态分布的图像<要求写出程序并作图)。 程序: 1.计算三种随机变量分布的方差与期望 [m0,v0]=binostat(10,0.3> %二项分布,取n=10,p=0.3 [m1,v1]=poisstat(5> %泊松分布,取lambda=5 [m2,v2]=normstat(1,0.12> %正态分布,取u=1,sigma=0.12 计算结果: m0 =3 v0 =2.1000 m1 =5 v1 =5 m2 =1 v2 =0.0144 2.计算x=45和x<45的概率,并绘图 Px=binopdf(45,100,0.5> %x=45的概率 Fx=binocdf(45,100,0.5> %x<45的概率 x=1:100。 p1=binopdf(x,100,0.5>。 p2=binocdf(x,100,0.5>。 subplot(2,1,1>

plot(x,p1> title('概率密度图像'> subplot(2,1,2> plot(x,p2> title('概率累积分布图像'> 结果: Px =0.0485 Fx =0.1841 3.t(10>分布与标准正态分布的图像 subplot(2,1,1> ezplot('1/sqrt(2*pi>*exp(-1/2*x^2>',[-6,6]> title('标准正态分布概率密度曲线图'> subplot(2,1,2> ezplot('gamma((10+1>/2>/(sqrt(10*pi>*gamma(10/2>>*(1+x^2/10>^(-(10+1>/2>',[-6,6]>。b5E2RGbCAP title('t(10>分布概率密度曲线图'> 结果:

数学实验答案

实验一 %sy1ljq20111668 %第一大题 %1 x=[3,2*pi]; y1=sin(x)+exp(x) %y1= 20.2267 535.4917 %2 x=2:2:10 y2=x.^2+sqrt(2*x) %y2= 6.0000 18.8284 39.4641 68.0000 104.4721 %3 a=2*pi,b=35/180*pi,c=exp(2); y31=sin(a/5)+cos(b)*c y32=tan(b)*cot(a/3) %y31 =7.0038 %y32 =-0.4043 %6 a1=-6.28,a2=7.46,a3=5.37; a11=fix(a1) a21=fix(a2) a31=fix(a3) %a11=-6 %a21=7 %a31=5 %7

y71=abs(a1*a2+a3) y72=a1^2*sqrt(a2*a3/2) %y71 =41.4788 %y72 =176.5066 %8 save sy1 clear %9 load sy1 %10 A=[2 -5 6;8 3 1;-4 6 9]; A1=A' A2=det(A) A3=5*A save sy1 A1 A2 A3 %A1 = 2 8 -4 -5 3 6 6 1 9 %A2 =782 %A3 = 10 -25 30 40 15 5 -20 30 45 %第二大题 %1 X=0:pi/10:2*pi; Y=cos(X);S=[X',Y']

%S = 0 1.0000 0.3142 0.9511 0.6283 0.8090 0.9425 0.5878 1.2566 0.3090 1.5708 0.0000 1.8850 -0.3090 2.1991 -0.5878 2.5133 -0.8090 2.8274 -0.9511 3.1416 -1.0000 3.4558 -0.9511 3.7699 -0.8090 4.0841 -0.5878 4.3982 -0.3090 4.7124 -0.0000 5.0265 0.3090 5.3407 0.5878 5.6549 0.8090 5.9690 0.9511 6.2832 1.0000 %2 a22=input('a22='); b22=input('b22=');

高中数学必修三 概率与统计

高中数学必修三:概率与统计 1.要从已编号(1-50)的50枚最新研制的某型号导弹中随机抽取5枚来进行发射试验,用每部分选取的号码间隔一样的系统抽样方法确定所选取的5枚导弹的编号可能是( ). A.5,10,15,20,25B.3,13,23,33,43C.1,2,3,4,5D.2,4,8,16,32 2.从鱼塘捕得同一时间放养的草鱼240尾,从中任选9尾,称得每尾鱼的质量分别是1.5,1.6,1.4,1.6,1.3,1.4,1.2,1.7,1.8(单位:千克).依此估计这240尾鱼的总质量大约是( ).A.300克B.360千克C.36千克D.30千克 3.以下茎叶图记录了甲.乙两组各五名学生在一次英语听力测试中的成绩(单位:分) 已知甲组数据的中位数为15,乙组数据的平均数为16.8,则,x y的值分别为()A.2,5B.5,5C.5,8D.8,8 4.为了考查两个变量x和y之间的线性关系,甲、乙两位同学各自独立作了10次和15次试验,并且利用线性回归方法,求得回归直线分别为l1,l2,已知两人得的试验数据中,变量x和y的数据的平均值都分别相等,且值分别为s与t,那么下列说法正确的是( ). A.直线l1和l2一定有公共点(s,t)B.直线l1和l2相交,但交点不一定是(s,t) C.必有直线l1∥l2 D.直线l1和l2必定重合 5..设某大学的女生体重y(单位:kg)与身高x(单位:cm)具有线性相关关系,根据一组样本数据(x i,y i)(i=1,2,…,n),用最小二乘法建立的回归方程为$y=0.85x-85.71,则下列结论中不正确的是( ).A.y与x具有正的线性相关关系B.回归直线过样本点的中心(x,y)C.若该大学某女生身高增加1cm,则其体重约增加0.85kgD.若该大学某女生身高为170cm,则可断定其体重比为58.79kg

相关文档
最新文档